[發明專利]一種基于Khinchin-Einstein定理的微震事件檢測方法和系統

在審

• 申請號: 201910920368.2
• 申請日: 2019-09-26
• 公開（公告）號: CN110703319A
• 公開（公告）日: 2020-01-17
• 發明（設計）人: 翟明岳
• 申請（專利權）人: 廣東石油化工學院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/24
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于Khinchin-Einstein定理的微震事件檢測方法，其特征在于，包括：

步驟1，輸入實測的微震信號序列S；

步驟2，生成微震信號差分原始序列，具體為：

ΔS^O＝[ΔS₁,ΔS₂,…,ΔS_N-1,ΔS_N]＝[S₂-S₁,S₃-S₂,…,S_N-S_N-1,S_N-S₁]；

其中:

ΔS^O：所述信號差分原始序列，長度為N；

ΔS_n：所述信號差分原始序列ΔS^O中第n個元素[n＝1,2,…,N]；

S_n：所述微震信號序列S中的第n個元素；

S＝[S₁，S₂，…，S_N]:所述微震信號序列，長度為N；

步驟3，生成微震信號差分次生序列，具體為：

ΔS^K＝[ΔS₁,ΔS₂,…,ΔS_K]；

其中:

ΔS^K：K點處的所述微震信號差分次生序列[K＝1,2,…,N]；

步驟4，計算頻域窗口函數H(f)，具體為：

其中：

Δf＝f_MAX-f_MIN：頻帶寬度；

f_MAX：所述信號序列S的最高頻率；

f_MIN：所述信號序列S的最低頻率；

步驟5，求取K點處的對數化窗口能量E_K，具體為：

其中：

P_K(f)：序列ΔS^O+ΔS^K的功率譜密度；

步驟6，根據Khinchin-Einstein定理檢測微震事件；具體為：如果K點處的對數化窗口能量E_K≤2.595log[1+f_MIN]，則在所述信號序列S的第K點處檢測到微震事件；否則，未檢測到微震事件。