一种微地震剪切源弱事件S波初至自动拾取方法技术

本发明专利技术提供了一种微地震剪切源弱事件S波初至自动拾取方法及计算机存储介质。该方法通过两两分量偏振分析与旋转处理，从原始微地震三分量矢量波场中分离出P波数据和S波数据以及S波剩余分量T；然后，一方面利用长短时窗能量比法，搜索出P波数据初至时间，另一方面利用S波数据与S波剩余分量T，构建S波协方差矩阵特征值能量比曲线，并自动搜索出与P波相匹配的S波初至时间；最后，通过“零值”处理，去除S波特征值能量比曲线上已经识别出的P波和相匹配的S波数值，再通过给定较小门槛值，实现剩余的微地震弱事件S波初至时间自动拾取。本发明专利技术能够实现对微地震事件快速、准确地定位，整个定位过程简单便捷，具有很好的技术应用前景。

【技术实现步骤摘要】
一种微地震剪切源弱事件S波初至自动拾取方法
本专利技术涉及微地震监测数据处理
，尤其涉及一种微地震剪切源弱事件S波初至自动拾取方法。
技术介绍
微地震监测数据处理重点在于事件的准确定位，影响因素主要包括速度模型建立、反演算法的适用性、正演算法的精度、初至拾取等方面。其中，纵、横波初至准确拾取是震源精确定位首要条件之一。井中微地震初至拾取方法，一般采用常规地震直达波初至拾取方法，其原理主要是根据有效波与噪声在能量、偏振特性以及其他一些统计特性上存在区别，获得稳定可靠的初至时间，如能量分析法、自回归AR模型法、偏振分析法等等。每种类型方法都有各自算法特点，能量分析法是基于长短时窗能量比，当信号到达时,能量比变化快,相应的值会有一个明显的突跳，将该点时间定义为有效事件初至时间，但是信噪比较低情况容易出现误拾、漏拾；自回归AR模型法是基于信号与背景噪音属于不同AR模型，信号到达时AIC值极小值，但是不能直接判断是否是有效信号；偏振分析法是基于有效信号偏振度高、随机信号偏振度低，对应不同线性偏振系数曲线，但是该方法无法单独进行有效信号的检测本文档来自技高网...

【技术保护点】
1.一种微地震剪切源弱事件S波初至自动拾取方法，其特征在于，包括：/nS100，针对微地震P波信号，从原始微地震三分量矢量波场中分解出沿波传播方向P波分量；/nS200，针对微地震S波信号，从原始微地震三分量矢量波场中分解出垂直传播方向S波分量和切向S波剩余分量T；/nS300，根据P波分量计算P波长短时窗能量比，并基于P波长短时窗能量比构建P波特征曲线；/nS400，在P波特征曲线上搜索出P波长短时窗能量比大于给定的P波初至识别门槛值的峰值，拾取所述峰值对应的时间样点作为P波初至时间t

【技术特征摘要】
1.一种微地震剪切源弱事件S波初至自动拾取方法，其特征在于，包括：
S100，针对微地震P波信号，从原始微地震三分量矢量波场中分解出沿波传播方向P波分量；
S200，针对微地震S波信号，从原始微地震三分量矢量波场中分解出垂直传播方向S波分量和切向S波剩余分量T；
S300，根据P波分量计算P波长短时窗能量比，并基于P波长短时窗能量比构建P波特征曲线；
S400，在P波特征曲线上搜索出P波长短时窗能量比大于给定的P波初至识别门槛值的峰值，拾取所述峰值对应的时间样点作为P波初至时间tP；
S500，根据S波分量和S波剩余分量T，构建S波瞬时协方差矩阵，计算其最大特征值，通过计算最大特征值前后时窗能量比，构建S波特征值瞬时能量比曲线；
S600，确定与P波初至时间tP相匹配的S波初至时间tS的可能的时间范围，在此时间范围内，从S波特征值瞬时能量比曲线上搜索出与P波初至时间tP相匹配的S波初至时间tS；
S700，在P波初至时间tP和S波初至时间tS位置附近，对S波特征值瞬时能量比曲线的数据进行置零处理，以去除S波特征值瞬时能量比曲线上已经识别出的P波初至和S波初至；
S800，在经过置零处理的S波特征值瞬时能量比曲线上，搜索出S波特征值瞬时能量比大于预设的弱事件S波初至识别门槛值的极值，拾取所述极值对应的时间样点作为所求的微地震剪切源弱事件S波初至时间t*S。


2.如权利要求1所述的微地震剪切源弱事件S波初至自动拾取方法，其特征在于，所述步骤S100包括以下步骤：
S110，针对微地震P波信号，通过矢端曲线分析，对原始微地震三分量中的X分量和Y分量进行偏振分析与旋转处理，获得波传播径向分量R和切向分量T；
S120，再通过矢端曲线分析，对波传播径向分量R和原始微地震三分量中的Z分量进行偏振分析与旋转处理，获得沿波传播方向P波分量和垂直传播方向N分量。


3.如权利要求1所述的微地震剪切源弱事件S波初至自动拾取方法，其特征在于，所述步骤S200包括以下步骤：
S210，针对微地震S波信号，通过矢端曲线分析，对原始微地震三分量中的X分量和Y分量进行偏振分析与旋转处理，获得波传播径向分量R和切向分量T；
S220，再通过矢端曲线分析，对波传播径向分量R和原始微地震三分量中的Z分量进行偏振分析与旋转处理，获得垂直传播方向S波分量和切向S波剩余分量T。


4.如权利要求1所述的微地震剪切源弱事件S波初至自动拾取方法，其特征在于，在所述步骤S300中，按照下式计算P波长短时窗能量比，以构建P波特征曲线：



其中，ERPi为时间样点i对应的P波长短时窗能量比，Pi为时间样点i对应的P波分量，L1和L2分别为给定的长时窗和短时窗。


5.如权利要求4所述的微地震剪切源弱事件S波初至自动拾取方法，其特征在于，所述步骤S400包括以下步骤：
S410，定义P波初至识别门槛值KP，在P波特征曲线上搜索出P波长短时窗能量比大于等于P波初至识别门槛值时对应的时间样点的位置TKP，其中：
当i＜TKP且ERPi＜KP
当i＝TKP且ERPi≥KP
S420，定义P波初至搜索时窗的大小WP，搜索在P波初至搜...

【专利技术属性】
技术研发人员：余波，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11