[發明專利]一種微地震剪切源弱事件S波初至自動拾取方法

有效

• 申請號: 201910497390.0
• 申請日: 2019-06-10
• 公開（公告）號: CN112068193B
• 公開（公告）日: 2022-05-03
• 發明（設計）人: 余波
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種微地震剪切源弱事件S波初至自動拾取方法，其特征在于，包括：

S100，針對微地震P波信號，從原始微地震三分量矢量波場中分解出沿波傳播方向P波分量；

S200，針對微地震S波信號，從原始微地震三分量矢量波場中分解出垂直傳播方向S波分量和切向S波剩余分量T；

S300，根據P波分量計算P波長短時窗能量比，并基于P波長短時窗能量比構建P波特征曲線；

S400，在P波特征曲線上搜索出P波長短時窗能量比大于等于給定的P波初至識別門檻值的峰值，拾取所述峰值對應的時間樣點作為P波初至時間t_P；

S500，根據S波分量和切向S波剩余分量T，構建S波瞬時協方差矩陣，計算其最大特征值，通過計算最大特征值前后時窗能量比，構建S波特征值瞬時能量比曲線；

S600，確定與P波初至時間t_P相匹配的S波初至時間t_S的可能的時間范圍，在此時間范圍內，從S波特征值瞬時能量比曲線上搜索出與P波初至時間t_P相匹配的S波初至時間t_S；

S700，在P波初至時間t_P和S波初至時間t_S位置附近，對S波特征值瞬時能量比曲線的數據進行置零處理，以去除S波特征值瞬時能量比曲線上已經識別出的P波初至和S波初至；

S800，在經過置零處理的S波特征值瞬時能量比曲線上，搜索出S波特征值瞬時能量比大于等于預設的弱事件S波初至識別門檻值的極值，拾取所述極值對應的時間樣點作為所求的微地震剪切源弱事件S波初至時間t^*_S。

2.如權利要求1所述的微地震剪切源弱事件S波初至自動拾取方法，其特征在于，所述步驟S100包括以下步驟：

S110，針對微地震P波信號，通過矢端曲線分析，對原始微地震三分量中的X分量和Y分量進行偏振分析與旋轉處理，獲得波傳播徑向分量R和水平切向分量；

S120，再通過矢端曲線分析，對波傳播徑向分量R和原始微地震三分量中的Z分量進行偏振分析與旋轉處理，獲得沿波傳播方向P波分量和垂直傳播方向N分量。

3.如權利要求1所述的微地震剪切源弱事件S波初至自動拾取方法，其特征在于，所述步驟S200包括以下步驟：

S210，針對微地震S波信號，通過矢端曲線分析，對原始微地震三分量中的X分量和Y分量進行偏振分析與旋轉處理，獲得波傳播徑向分量R和切向S波剩余分量T；

S220，再通過矢端曲線分析，對波傳播徑向分量R和原始微地震三分量中的Z分量進行偏振分析與旋轉處理，獲得垂直傳播方向S波分量和切向S波剩余分量T。

4.如權利要求1所述的微地震剪切源弱事件S波初至自動拾取方法，其特征在于，在所述步驟S300中，按照下式計算P波長短時窗能量比，以構建P波特征曲線：

其中，ER_Pi為時間樣點i對應的P波長短時窗能量比，P_i為時間樣點i對應的P波分量，L₁和L₂分別為給定的長時窗和短時窗。