[发明专利]一种特定地区地震事件的报警方法

权利要求书

1.一种特定地区地震事件的报警方法，其特征在于：包括以下步骤：

S1、在给定时间段范围内和特定地区范围内选取合适的报警台阵；

S2、读取报警台阵各子台的实时数据，进行实时F-K计算，并得到数据中能量最大信号的方位角和慢度；具体包括以下步骤：

S2.1、对报警台阵的实时数据进行F-K计算，得到在不同慢度条件下的信号谱P：

式中S_n和S_e是慢度在南北向和东西向上的分量，F_i(f)是报警台阵第i个子台的实时数据与自身作卷积后的序列，f₁和f₂分别是滤波频带的起止频率，J为报警台阵的子台个数，dnorth_i和deast_i是报警台阵第i个子台相对于中心子台在南北向和东西向的偏移量；

S2.2、记录信号谱P在最大值时对应的慢度向量S_nm、S_em；

S2.3、计算慢度向量S_nm、S_em对应的方位角数值Azi：

S2.4、计算慢度向量S_nm、S_em对应的慢度数值Slo：

S3、根据S2获得的方位角和慢度对报警台阵数据进行聚束得到聚束通道结果，然后根据聚束通道结果计算信号检测结果，即短时绝对平均值和长时绝对平均值的比值；

S4、将计算得到的方位角、慢度和信号检测结果与预先设置的阈值进行比较，对符合条件的信号进行事件报警，具体如下：

S4.1、由选取的报警台阵的坐标和特定地区的范围，得到该地区事件到报警台阵的方位角范围[Azi_min Azi_max]；

S4.2、由选取的报警台阵的坐标和特定地区的范围，得到该地区事件到报警台阵的初至波震相类型，给出该震相的慢度范围[Slo_min Slo_max]；

S4.3、根据计算得到方位角、慢度和信号检测结果的最大值SNR_max同时满足下列条件时，给出报警：

条件1：Azi_min≤Azi≤Azi_max

条件2：Slo_min≤Slo≤Slo_max

条件3：SNR_max≥Th，其中Th为预设的信号检测阈值。

2.根据权利要求1所述的一种特定地区地震事件的报警方法，其特征在于：S1具体包括以下步骤：

S1.1、计算第i个台阵到特定地区范围内的最大距离Di：

D_i＝max(distance(sta_i,range))

式中sta_i为第i个台阵中心子台的经纬度坐标，range为特定地区经纬度范围，distance为台阵和特定地区之间距离计算函数，max为取最大值函数；

S1.2、选取距离在给定距离D内的台阵：

D_i＜D

S1.3、统计给定时间段范围内特定地区范围内发生的事件总数为N，S1.2中选取的台阵中每个台阵检测到的事件数为K_i，其中1<K_i≤N，计算每个台阵的监测灵敏度Ei：

S1.4、选取灵敏度大于给定阈值E的一个或多个台阵作为报警台阵：

E_i>E。

3.根据权利要求2所述的一种特定地区地震事件的报警方法，其特征在于：S3具体包括以下步骤：

S3.1、根据步骤S2.2中计算得到的慢度向量S_nm、S_em对台阵数据进行聚束，得到聚束通道时间序列数据BeamData(t)：

S3.2、根据获取的台阵聚束通道时间序列数据BeamData(t)，计算信号检测时刻短时绝对平均值和长时绝对平均值：

式中A_i表示是采样点序号i的数值，j表示信号检测时刻采样点序号，STA_j表示短时绝对平均值，LTA_j表示长时绝对平均值，SN表示短时绝对平均值时间窗口数据点数，LN表示长时绝对平均值时间窗口数据点数；

S3.3、计算STA_j/LTA_j最大值SNR_max。