[发明专利]一种次声信号与地震事件的自动关联方法

权利要求书

1.一种次声信号与地震事件的自动关联方法，其特征在于：包括以下步骤：

S1、确定待关联地震事件对应的次声信号可监测震中距mag_th：首先提取待关联地震事件的体波震级mb、近震震级ML，选取其中的最大值m作为台站筛选震级，确定待关联地震事件对应的次声信号可监测震中距mag_th与m的关系如下：

根据上述关系，确定次声监测台网中符合震级-震中距的具有监测数据的可关联次声台站列表1；

S2、从历史次声事件数据库中搜索与待关联地震事件距离小于特定震中距阈值的参考事件，确定该地区历史次声事件曾关联的次声台站，形成可关联次声台站列表2，综合可关联次声台站列表1和可关联次声台站列表2，形成待关联次声台站清单；

S3、对于每一个待关联次声台站，分别计算其与待关联地震事件间的震中距dist_i,j，台站至地震事件的方位角az_i,j；其中i为待关联地震事件序号，j为待关联次声台站序号；

利用震中距dist_i,j计算信号频率上限f max_i,j、不同等效声速条件下的信号理论走时，选取走时误差delt_time_i,j；

fmax_i,j＝4-0.055*dist_i,j

t1_i,j＝111.1949*dist_i,j/v_h

t2_i,j＝111.1949*dist_i,j/v_m

t3_i,j＝111.1949*dist_i,j/v_l

其中v_h、v_m、v_l分别表示等效声速的上限、平均值及下限；dealt_time1、dealt_time2分别为次声信号走时误差识别低误差系数及走时误差识别高误差系数；dealt_time1、dealt_time2根据需求预先设定；

S4、关联次声台站信号方位角误差；

S5、根据地震事件发生时间确定次声信号的检索起止时间t_start_i,j和t_end_i,j；

对于次声信号检索起始时间，若t1_i,j-0.5*dealt_time_i,j>0：

t_start_i,j＝t1_i,j+torid_i-0.5*dealt_time_i,j

否则：

t_start_i,j＝torid_i

次声信号检索结束时间为：

t_end_i,j＝torid_i+t3_i,j*dealt_time_i,j

其中torid_i为地震事件发生时间；

S6、确定次声台站记录的P相和S相信号，得到其次声信号搜索起止时间：

tmin_p_i,j＝dist_i,j·pslow_min-ptime_delta+torid_i

tmax_p_i,j＝dist_i,j·pslow_max+ptime_delta+torid_i

tmin_S_i,j＝dist_i,j·Sslow_min+Stime_delta+torid_i

tmax_S_i,j＝dist_i,j·Sslow_max+Stime_delta+torid_i

其中，pslow_min为P相信号慢度最小值、pslow_max为P相信号慢度最大值、ptime_delta为P相信号走时误差限值、Sslow_min为S相信号慢度最小值、Sslow_max为S相信号慢度最大值、Stime_delta为S相信号走时误差限值；

S7、对各次声台站的方位角范围内的次声信号进行搜索，找出符合条件的次声信号作为由地震事件诱发的次声信号，并与地震事件关联，形成声震联合关联事件。