张丽娜
福建省地震局 福建 福州 350003
由于福建地震台网日常工作中需CAP方法快速并产出可靠的震源机制解,如本省发生5.0级以上地震时,或若地震震级在5.0级以下而中国地震局启动应急响应时,2小时内完成主震震源机制解的计算。然而,CAP方法受台站布局、台站拾取数量等因素的影响,结果也存在着差异性。本研究探索台站布局等因素,为快速产出较为可靠的震源机制解提供参考,缩短震源机制解产出用时。
CAP方法基本思想是将地震记录波形分割为Pn1和面波波段截开,分别赋予不同的权重,利用格点搜索的方法进行地震震源机制的反演。设是台站记录到的去除仪器响应后的地震波形,是相应的理论波形。一个双力偶震源的理论合成位移可以表示为,式中分别对应垂直走向滑动断层、垂直倾向动断层和45°倾向滑动断层3种最基本的断层类型。标量地震矩,为格林函数,为震源的辐射花样信息,为台站方位角,分别为待求震源机制解的断层走向、倾角和滑动角。对于中的未知量,若已知格林函数,就可以通过求解式子:。将初步得出的震源机制解与源函数、格林函数按照上式计算之后合成出理论图,将此理论图和实测图对比拟合。
CAP反演运算中,理论波形的得出依赖于各个台站位置的格林函数计算,目前计算格林函数的方法有多种,使用F-K方法( 频率-波数法)在水平层状模型中计算高频理论地震图时能够得到多种震相,可以计算出完整的理论地震图,并且对于不同震源类型、震中距都有普适性[1]。……