场地响应及路径衰减联合反演方法软件实现及应用.docxVIP

场地响应及路径衰减联合反演方法软件实现及应用 朱新运 摘要：地丧波哀减及场地响应研究是地震学研究领域重要的前沿课题C基干谱比法.联合反演求解场地响 应及地震波衰减参数，这种算法理论上是可行的?据此，基干MATLAB开发平台开发研制了地震波衰减参数 与台站场地响应联合反演软件。木文介绍软件设汁的思路.理论原理.方法实现及资料处理过程，简述软 件模块功能及使用方法。该软件计算过程稳定、可靠?界【加设计友好.可用于日常地丧资料的分析，也可 作为地丧科学研究的工具。 关键诃场地响应 衰减参数联合反演软件 1引言 在对耦合了路径、场地、仪器、地震輻射的方向性等多种响应成分的地丧记录(陈运泰 等，2000)进行必要限定的情况下il?算波谱比，可以扣除耦合的几何扩散、仪器响应、地震 辐射方向性，但是场地响应与非均匀性、非弹性路径衰减则难以简单地扣除的。为了获得较 为准确的衰减系数，Chun等(1987)提岀的“逆向双台谱比法”。该方法是计算Lg波衰减参数 比较可靠的方法(Chung, 2002)，它通过两台站连线两端延长线位置上发生的成对地震的记 录谱比计算衰减参数。由于这种方法在资料选择上要求比较苛刻,影响了该方法的推广应用。 Chung等(2005)对逆向双台谱比法”进行了改进，改进后的“逆向双台谱比法”不再严 格要求地震对以及台站对处在一条直线上，只要求地震对及台站对之间的距离满足一立的条 件即可。这一改进扩大了地虎资料的选择范用，但对资料的应用仍然受到一立的限制。 场地响应描述的是近地表地层介质对地震地而运动的放大作用，研究台站场地响应 (Bonilla et al, 1997; Riepl et al, 1998;Langston, 1979; Nath et al, 2000,2001) 的方法很多，但基本上是通过路径衰减项简化求得场地响应的「根^Shearer (1999)的研 究结果，场地响应值与近地表地层介质的介质密度及地震波速度有关，与地震波传播距离无 关，而地震衰减参数是地震波路径衰减的描述，仅与距离有关，由于地震Lg波的衰减对地震 辐射模式不敏感(Mitra et al, 2006),因而使用Lg波作为基本资料进行相关研究时，台站 方位可以忽略。据此，朱新运等(朱新运等，2007, ZHUXin_Yun, et al, 2007;2012;2013;2014) 以地丧Lg波台站记录谱为观测数据，在资料选择上，使用Ms5级以上地箴，为保证通过谱比 消去震源因子时不致引入较大误差，限左台站对地丧方位角变化不大，同时为保证衰减参数 稳泄，限宦地震波穿过构造相对单一的区域，利用路径衰减与台站震中距有关、而场地响应 与震中距无关这一特点，通过联合反演将场地响应与路径衰减互相分开，得到了较好的结果。 本文中作者对场地响应与路径衰减参数联合反演方法论证，同时介绍了在MATLAB平台下 研制的路径衰减与场地响应联合反演软件，及使用该软件研究了华北、华东等区域的地喪波 衰减及场地响应参数成果。软件设汁构件合理、界面美观、计算功能完善、使用方便，可以 作为地震波相关研究的实用工具。 2理论与方法 台站记录的地震波的频谱A(/,d)可以表示为(Hasegawa, 1985)： A(f、d) = So(f)R(JQG(dW(f,d)I(f)Si(门 (1) 式中，/表示频率：So(/)表示震源因子项；/?(/,)表示宸源辐射图型因子：d表示箴 中距：/(/)表示仪器效应：Sf(/)表示场地响应：G(d)表示几何扩散因子， G(〃)= d (2) “是表示几何扩散特性的参数，对于而波波，“ = -0.5,体波= 与(/,〃)表示介质 衰减项，它可以表示为： (3)以/,〃)=宀 (3) 式中，了是衰减系数, (11)如鮎(/) = (11) 如鮎(/) = bij (/) -/(/)? (13) QU (4) 在式（4）中，（/表示波速度:。 假定参与运算的地丧数为加,￡ = 1,2….，〃八需要评估的台站数为/I,两个台站（第 i台与第丿台，jj = 123,???‘）对同一个地喪（第￡个地箴）的记录谱的比可表示为： 仏O,仏）So心）心￡k.）G（d小讥f （门Si心） “、 4.J （几如 J）Sq （f）Rkj （/, q J ）G（%?4 J 兀（/）旳（/） 在上式中，对于同一地震，丧源因子项之比可以消去：震源辐射图型因子R（f、e）与震源机 制解有关，对较大地震选取资料时可考虑丧中距离与台站间距离相比很大，以至于震源对台 站之间方位角变化很小，贝J： 」（/?」）心.丿（/禺丿，对较小地震则可看作点源，震源 辐射因子影响不大，从而，在此情况下震源辐射图型因子项之比可以消去；仪器效应人（/） 与///）是已知的；A（fd）是台站