[发明专利]一种实现地震预警的方法

权利要求书

1.一种实现地震预警的方法，其特征在于：包括以下步骤：

步骤1：实时提取三分量地震动记录{x_i}、{y_i}、{z_i}，i为地震动记录的实时记录序号；三分量地震动记录可以为位移、速度或加速度；

步骤2：对三分量地震动记录进行实时低通滤波，低通滤波器阈值根据地震预警的实际要求情况进行选择，如果考虑绝对公共安全，即可能发生破坏性地震，则p波初始周期阈值取为1s，如果考虑地震预警准确性，即一定发生破坏性地震，则p波初始周期阈值取为2s；根据p波初始周期判断是否进行地震预警计算，若p波初始周期小于p波初始周期阈值，则不能发生破坏性地震，停止地震预警计算，返回步骤1；若p波初始周期大于等于p波初始周期阈值，则能发生破坏性地震，执行步骤3；

步骤3：对经过三维低通滤波的地震记录，用公式(1)计算三维振幅相对变化率，得到P_i序列：

其中，0≤i≤N-2，N为序列长度；

步骤4：用三维振幅相对变化率P_i作为特征函数，用公式(2)对三维振幅相对变化率P_i做长短时均值比计算，得到序列{P_i}；

其中，N₄为长短时均值比方法的短尺度，N₃为长尺度，N₄＜＜N₃＜＜N，0≤i≤N-2-N₃；

步骤5：根据环境干扰、仪器干扰和地方震特点设定地震预警触发阈值P_c，根据公式(3)判断是否触发地震预警，若公式(3)满足，则触发地震预警，若公式(3)不满足，则不能触发，返回步骤1；

max{P_i}≥P_c (3) 。

2.根据权利要求1所述的一种实现地震预警的方法，其特征在于：所述步骤3中计算特征函数时为避免出现由于振幅为零使计算溢出的情况，公式(1)中所有作为分母的振幅数值都加上10^-7数量级。

3.根据权利要求1所述的一种实现地震预警的方法，其特征在于：所述步骤4中为了使p波定时准确，所用长尺度和短尺度都较短，具体如下：若地震记录的采样间隔为0.002s，则短尺度取0.1s，长尺度取0.3s；若地震记录的采样间隔为0.005s，则短尺度取0.25s，长尺度取0.75s。