[發明專利]一種基于臺站對EDT面概率分布函數的微震定位方法

在審

• 申請號: 202111579966.1
• 申請日: 2021-12-22
• 公開（公告）號: CN114415231A
• 公開（公告）日: 2022-04-29
• 發明（設計）人: 曾志毅;韓鵬;楊曉輝;王蕤;繆淼
• 申請（專利權）人: 南方科技大學
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京匯信合知識產權代理有限公司 11335
• 代理人: 孫民興

1.一種基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，包括：

對所有臺站的微震信號進行濾波預處理；

對濾波預處理后的所述微震信號進行拾取，得到每個所述臺站的微震信號初至數據；

根據所述初至數據，計算獲得所有所述臺站的微震信號信噪比，并對所述微震信號信噪比進行歸一化計算；

根據所述初至數據，得到所有臺站中任意兩個不同所述臺站的觀測到時差數據，并保留兩個不同所述臺站總所述微震信號信噪比小的信噪比數據；

對微震監測區域進行三維網格劃分，并計算每個網格點到每個所述臺站的旅行走時，得到計算走時表數據；

將任意兩個不同的所述計算走時表數據相減，得到計算走時差數據；

將所述觀測到時差數據、所述微震信號信噪比小的信噪比數據和所述計算走時差數據代入到EDT面概率分布函數中，得到每個網格點的概率值；

最大概率值的網格點位置為微震事件定位結果。

2.如權利要求1所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，根據所述初至數據，計算獲得所有所述臺站的微震信號信噪比包括：

信噪比公式為：

式中：為微震事件i在臺站k的信噪比；W₁為從時刻往前取長度為W的時窗內微震信號振幅值；W₂為從時刻往后取長度為W的時窗內背景噪音振幅值；和分別表示為W₁和W₂窗口內波形振幅的均方根值。

3.如權利要求2所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，對所述微震信號信噪比進行歸一化計算包括：

歸一化計算公式為：

式中：為微震事件i在臺站k的歸一化信噪比；maxⁱ(SNR)為微震事件i的信噪比的最大值。

4.如權利要求3所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，根據所述初至數據，得到所有臺站中任意兩個不同所述臺站的觀測到時差數據包括：

式中：和分別為微震事件i在臺站k和臺站l的觀測初至；為微震事件i在臺站k和臺站l的觀測到時差，即為微震事件的臺站對到時差。

5.如權利要求4所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，保留兩個不同所述臺站總所述微震信號信噪比小的信噪比數據包括：

式中：為保留微震事件i在臺站對k，l中微震信號信噪比較小的信噪比數據；為取較小值運算。

6.如權利要求5所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，將所述觀測到時差數據、所述微震信號信噪比小的信噪比數據和所述計算走時差數據代入到EDT面概率分布函數中，得到每個網格點的概率值包括：

計算公式為：

式中：L^g為網格點g的概率值；TERR為容忍度參數；為臺站對k，l的計算走時差與觀測到時差的差；為EDT面概率分布函數。

7.如權利要求6所述的基于臺站對EDT面概率分布函數的微震定位方法，其特征在于，最大概率值的網格點位置為微震事件定位結果包括：

S(G)＝max(L^g)；

式中：S(G)為概率最大值，G為微震事件定位結果。