[发明专利]基于天然地震震源辐射强度修正的观测信号幅值测量方法

权利要求书

1.一种基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于，包括以下步骤：

1)将一次天然地震在地震台站j上记录的以count为强度单位的原始波形数据去除仪器响应，得到以所测量物理量单位为强度单位的波形数据；

2)利用幅值测量所需要的数字滤波器，对去除仪器响应后的波形数据进行滤波；

3)指定测量时间，对滤波后的波形数据测量相应时间窗内的初始波形信号幅值所述波形信号幅值为峰峰信号幅值或者均方根信号幅值；

4)从地震监测机构查询该地震的震源矩张量M及相应的震源持续时间；

5)计算所述地震台站j与震源之间的距离，以及从震源看向地震台站j的方位角；

6)根据震源所在区域的介质结构模型，利用地震图计算程序计算所述地震台站j所对应震中距的格林函数；

7)在以震源为中心、以所述地震台站j与震源之间的距离为半径的圆周上，从0°方位角开始，均匀间隔设置N个模拟地震台站；利用格林函数和震源矩张量M合成N个模拟地震台站的理论地震图；

8)利用步骤2)和步骤3)中的数字滤波器、时间窗及测量方法，测量N个模拟地震台站的理论地震图波形信号幅值A_i，并根据计算N个模拟地震台站的理论地震图平均幅值即震源平均辐射强度；

9)根据从震源看向地震台站j的方位角，合成地震台站j的理论地震图，并利用步骤2)和步骤3)中的数字滤波器、时间窗及测量方法，测量相应的理论地震图波形信号幅值即地震台站j的震源辐射强度；

10)利用得到修正后的波形信号幅值

其中，为步骤3)得到的初始波形信号幅值；

为步骤8)得到的震源平均辐射强度；

为步骤9)得到的地震台站j的震源辐射强度。

2.根据权利要求1所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

所述地震台站j为区域台站，步骤6)中所述格林函数的计算步骤为：

i.从地震公报中查询得到地震震源深度，查询台站信息得到地震台站j的深度；

ii.计算格林函数。

3.根据权利要求2所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

所述区域台站为震中距小于1500km的地震台站。

4.根据权利要求1所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

所述地震台站j为远震台站，步骤6)中所述格林函数的计算步骤为：

i.根据地震台站j的震中距，查询iaspei91走时表得到相应的P波或S波慢度值p；

ii.利用公式i_s＝arcsin(p·α)计算地震台站j的P波或S波射线在震源处的离源角i_s；

其中，p为P波或S波慢度值；α为震源处的P波或S波的波速；

iii.从地震公报中查询得到地震震源深度；

iv.利用公式得到地震台站j的等效深度H；

其中，H_s为震源深度；△为格林函数计算的等效震中距，对于深度小于30km的地震，△取为400km；i_s为P波或S波射线在震源处的离源角；

v.计算格林函数。

5.根据权利要求4所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

所述远震台站为震中距大于2000km的地震台站。

6.根据权利要求1至5任一所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

所述N＝72。

7.根据权利要求6所述的基于天然地震震源辐射强度修正的观测信号幅值测量方法，其特征在于：

步骤6)中所述地震图计算程序为CPS或SEM。