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

摘要

本发明涉及一种次声信号与地震事件的自动关联方法，首先确定待关联次声台站清单；对于每一个待关联次声台站，分别计算其与待关联地震事件间的震中距，台站至地震事件的方位角，根据地震事件发生时间确定次声信号的检索起止时间，确定次声台站记录的P相和S相信号，得到其次声信号搜索起止时间；对各次声台站的方位角范围内的次声信号进行搜索，找出符合条件的次声信号作为由地震事件诱发的次声信号，并与地震事件关联，形成声震联合关联事件。本发明综合多种信号关联约束条件和事件所在区域历史信息，可以识别出由地震事件诱发的次声信号，提高事件分析准确性；还可将次声监测数据序列多源信号中确定部分排除，降低信号误判风险。

主权项

1.一种次声信号与地震事件的自动关联方法，其特征在于：包括以下步骤：S1、确定待关联地震事件对应的次声信号可监测震中距mag_th：首先提取待关联地震事件的体波震级mb、近震震级ML，选取其中的最大值m作为台站筛选震级，确定待关联地震事件对应的次声信号可监测震中距mag_th与m的关系如下：根据上述关系，确定次声监测台网中符合震级‑震中距的具有监测数据的可关联次声台站列表1；S2、从历史次声事件数据库中搜索与待关联地震事件距离小于特定震中距阈值的参考事件，确定该地区历史次声事件曾关联的次声台站，形成可关联次声台站列表2，综合可关联次声台站列表1和可关联次声台站列表2，形成待关联次声台站清单；S3、对于每一个待关联次声台站，分别计算其与待关联地震事件间的震中距disti,j，台站至地震事件的方位角azi,j；其中i为待关联地震事件序号，j为待关联次声台站序号；利用震中距disti,j计算信号频率上限f maxi,j、不同等效声速条件下的信号理论走时，选取走时误差delt_timei,j；fmaxi,j＝4‑0.055*disti,jt1i,j＝111.1949*disti,j/vht2i,j＝111.1949*disti,j/vmt3i,j＝111.1949*disti,j/vl其中vh、vm、vl分别表示等效声速的上限、平均值及下限；dealt_time1、dealt_time2分别为次声信号走时误差识别低误差系数及走时误差识别高误差系数；dealt_time1、dealt_time2根据需求预先设定；S4、关联次声台站信号方位角误差；S5、根据地震事件发生时间确定次声信号的检索起止时间t_starti,j和t_endi,j；对于次声信号检索起始时间，若t1i,j‑0.5*dealt_timei,j>0：t_starti,j＝t1i,j+toridi‑0.5*dealt_timei,j否则：t_starti,j＝toridi次声信号检索结束时间为：t_endi,j＝toridi+t3i,j*dealt_timei,j其中toridi为地震事件发生时间；S6、确定次声台站记录的P相和S相信号，得到其次声信号搜索起止时间：tmin_pi,j＝disti,j·pslow_min‑ptime_delta+toriditmax_pi,j＝disti,j·pslow_max+ptime_delta+toriditmin_Si,j＝disti,j·Sslow_min+Stime_delta+toriditmax_Si,j＝disti,j·Sslow_max+Stime_delta+toridi其中，pslow_min为P相信号慢度最小值、pslow_max为P相信号慢度最大值、ptime_delta为P相信号走时误差限值、Sslow_min为S相信号慢度最小值、Sslow_max为S相信号慢度最大值、Stime_delta为S相信号走时误差限值；S7、对各次声台站的方位角范围内的次声信号进行搜索，找出符合条件的次声信号作为由地震事件诱发的次声信号，并与地震事件关联，形成声震联合关联事件。