一种微震信号P波初至点综合拾取方法

摘要

本发明公开了一种微震信号P波初至点综合拾取方法，包括如下步骤：导入微震时间序列数据；拾取STA/LTA法首个阀值触发点k₁；拾取PAI-K法峰度最大值点k₂；拾取AIC法自相关最小值点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₁。通过综合拾取法有效解决了STA/LTA法无拾取事件较多、拾取精度较低，PAI-K法和AIC法拾取稳定性较差的技术问题。此方法具有处理简便、适用性强、准确性高等特点。

主权项

一种微震信号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₁。