[发明专利]一种微震信号P波初至点综合拾取方法

权利要求书

1.一种微震信号P波初至点综合拾取方法，其特征在于，包括以下步骤：

步骤1：导入微震时间序列x(n)

导入微震时间序列x(n)(n＝1,2,…,N)，其中N为微震信号的总采样点数，取 N＝4000～7000，微震信号采样频率f＝4000～7000Hz；

步骤2：拾取STA/LTA法首个阀值触发点k₁

由公式计算x(n)的长短时窗均值比值STA(k)/LTA(k)，k为第 k个采样点，k＝W_LTA,W_LTA+1,…,N，W_STA、W_LTA分别为短、长时窗采样点个数，W_STA＝50～70， W_LTA＝500～700，取STA/LTA首个大于阀值的触发点作为第一P波初至点k₁；

其中，STA/LTA阀值为λ取2～3；

步骤3：拾取PAI-K法峰度最大值点k₂

由公式计算x(n)的峰度值K(k)，M为滑动时窗 的长度，M＝200～300，k＝M,M+1,…,N，取峰度最大值点作为第二P 波初至点k₂；

步骤4：拾取AIC法自相关最小值点k₃

由公式AIC(k)＝k·log{var(x[1,k])}+(N-k-1)·log{var(x[k+1,N])}计算x(n)的自相 关值AIC(k)，k＝1,2,…,N，var(x[1,k])指第1点至第k点的微震时间序列方差，var(x[k+1,N]) 指第k+1点至第N点的微震时间序列方差，取自相关最小值点作为第三P波初至点k₃。

步骤5：综合法确定P波初至点k

按照以下准则确定P波初至点k：

(1)k₁无拾取时，取k＝k₃；

(2)k₁有拾取时：①|k₁-k₃|≤a，取k＝k₃；②|k₁-k₂|≤a，且|k₁-k₃|>a，取k＝k₂；③|k₁-k₂|>a， 且|k₁-k₃|>a，取k＝k₁。

2.根据权利要求1所述的微震信号P波初至点综合拾取方法，其特征在于，取 N＝4000～7000，f＝4000～7000，W_STA＝50～70，W_LTA＝500～700，λ＝2～3，M＝200～300， a＝180～300。