基于MATLAB的场地效应及路径衰减联合反演软件研制.DOC

场地响应及路径衰减联合反演方法软件实现及应用 朱新运 摘要：地震波衰减及场地响应研究是地震学研究领域重要的前沿课题。基于谱比法，联合反演求解场地响应及地震波衰减参数,这种算法理论上是可行的，据此，基于MATLAB开发平台开发研制了地震波衰减参数与台站场地响应联合反演软件。本文介绍软件设计的思路、理论原理，方法实现及资料处理过程，简述软件模块功能及使用方法。该软件计算过程稳定、可靠，界面设计友好，可用于日常地震资料的分析，也可作为地震科学研究的工具。 关键词 场地响应 衰减参数 联合反演 软件 1 引言 在对耦合了路径、场地、仪器、地震辐射的方向性等多种响应成分的地震记录（陈运泰等，2000）进行必要限定的情况下计算波谱比，可以扣除耦合的几何扩散、仪器响应、地震辐射方向性，但是场地响应与非均匀性、非弹性路径衰减难以简单扣除的。为了获得较为准确的衰减系数，Chun等(1987)提出的“逆向双台谱比法”该方法是计算Lg波衰减参数较可靠的方法（Chung,2002），通过两台站连线两端延长线位置上发生的成对地震的记录谱比计算衰减参数。由于这种方法在资料选择上要求，影响了该方法的推广应用。Chung等（2005）对“逆向双台谱比法”进行了改进，改进“逆向双台谱比法”不再严格要求地震对及台站对在一直线上，只要求地震对及台站对之间的距离满足一定的条件即可这一改进扩大了地震资料的选择范围，但对资料的应用仍然。场地响应描述的是近地表地层介质对地震地面运动的放大作用，研究台站场地响应（Bonilla et al, 1997; Riepl et al, 1998;Langston, 1979; Nath et al, 2000,2001）的方法很多，但基本上是通过路径衰减项简化得场地响应的。根据Shearer（）研究结果，场地响应值与近地表地层介质的介质密度及地震波速度有关，与地震波传播距离无关，而地震衰减参数是地震波路径衰减的描述,仅与距离有关，由于地震Lg波衰减对地震射模式不敏感(Mitra et al2006)，因而使用Lg波作为基本资料进行相关研究时，台站方位可以忽略据此，朱新运等（）以地震Lg波台站记录谱为观测数据，利用路径衰减与台站震中距有关、而场地响应与震中距无关这一特点，通过联合反演将场地响应与路径衰减互相分开，得到了较好的果。 2 理论与方法 台站记录的地震波的频谱可以表示为（Hasegawa,1985）： （1） 式中，表示频率；表示震源因子项；表示震源辐射图型因子；表示震中距；表示仪器效应；表示场地响应；表示几何扩散因子， （2） 是表示几何扩散特性的参数，对于面波波，，体波；表示介质衰减项，它可以表示为： （3） 式中，是衰减系数， （4） 在式（4）中，表示波速度；。 假定参与运算的地震数为，，需要评估的台站数为，两个台站（第台与第台，）对同一个地震（第个地震）的记录谱的比可表示为： （5） 在上式中，对于同一地震，震源因子项之比可以消去；震源辐射图型因子与震源机制解有关，对较大地震选取资料时可考虑震中距离与台站间距离相比很大，以至于震源对台站之间方位角变化很小，则：，对较小地震则可看作点源，震源辐射因子影响不大，从而，在此情况下震源辐射图型因子项之比可以消去；仪器效应与是已知的；是台站记录谱；由式（2）给出。将式（3）代入式（5）整理后得： （6） 对式（6）两边取自然对数可以得： （7） 式中， （8） (9) （10） （11） （12） 则： （13） 选取资料时考虑研究时间段内，在研究区即地震波射线穿过的区域内，没有发生大的构造活动，衰减参数保持恒定；台站场地响应仅与频率有关，两台站对地震效应的幅度的比值仅与频率有关；在参与评估的台站中，对应