[發明專利]一種以振幅為特征函數的自動拾取P波震相方法

有效

• 申請號: 201811415093.9
• 申請日: 2018-11-26
• 公開（公告）號: CN109283576B
• 公開（公告）日: 2020-06-09
• 發明（設計）人: 李啟成;郭雷;何書耕;閔也
• 申請（專利權）人: 遼寧工程技術大學
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 沈陽東大知識產權代理有限公司 21109
• 代理人: 劉曉嵐

1.一種以振幅為特征函數的自動拾取P波震相方法，其特征在于，具體流程如下：

步驟1：輸入地震豎直方向記錄x_i(t)，i＝1...N，其中，x_i(t)為第i個豎直方向加速度記錄，其單位是cm/s²，N為地震數據的記錄長度；

步驟2：用公式(1)計算特征函數P_i；

P_i＝|x_i(t)ⁿ| (1)

其中，P_i為第i時刻x_i(t)對應的特征函數，n是大于或等于4的整數；

步驟3：用公式(2)計算長短時均值比序列{P_i}；

其中，P_i和P_j都是與地震記錄相關的特征函數，{P_i}是經過長短時均值比變換得到的序列，i是地震記錄點序號，N₄為長短時均值比方法的短尺度，N₃為長短時均值比方法的長尺度；j是計算過程中疊加運算需要的地震記錄序號；

步驟4：在長短時均值比序列{P_i}中，判斷長短時均值比計算結果，即判斷序列{P_i}中第i時刻長短時均值比是否大于或等于閾值，如果大于或等于閾值，繼續步驟5，如果小于閾值，則i＝i+1，判斷下一時刻長短時均值比結果是否大于或等于閾值，所述閾值為閾值下限的2倍，閾值上下限定義為：定義用長短時均值比拾取有效波即P波的第一個峰值平均值為拾取到時閾值的上限，拾取干擾波最大峰值平均值為閾值下限；

步驟5：在長短時均值比序列{P_i}中，從長短時均值比結果大于或等于閾值的對應時刻開始，按時間順序，尋找長短時均值比結果的極大值，該極大值對應的時刻就是P波到時；具體方法是：在長短時均值比結果達到閾值后，標記該時刻，把該時刻的長短時均值比結果與下一時刻的長短時均值比結果做比較，如果下一時刻的長短時均值比結果小于或等于該時刻的長短時均值比結果，則該時刻的長短時均值比結果對應的時刻就是P波到時；如果下一時刻長短時均值比結果大于該時刻的長短時均值比結果，則i＝i+1，重復步驟5，直到找到極大的長短時均值比結果，該極大的長短時均值比結果對應的時刻就是P波震相到時。

2.根據權利要求1所述一種以振幅為特征函數的自動拾取P波震相方法，其特征在于，所述特征函數僅用豎直方向地震記錄參與計算，東西、南北向地震記錄不參與計算。