[發明專利]一種實現地震預警的方法

有效

• 申請號: 201710842194.3
• 申請日: 2017-09-18
• 公開（公告）號: CN107479094B
• 公開（公告）日: 2018-11-30
• 發明（設計）人: 李啟成
• 申請（專利權）人: 遼寧工程技術大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 沈陽東大知識產權代理有限公司 21109
• 代理人: 梁焱

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。