井中微地震P/S波事件识别方法和装置与流程

技术特征：

1.一种井中微地震P/S波事件识别方法，该方法包括：

根据射孔数据P波事件水平分量偏振分析，获得检波器方位角，利用该检波器方位角对压裂数据水平分量进行方位一致性校正；

基于校正结果，根据给定的极性判断准则对压裂数据Z分量进行分离，获得分别含有P波、S波微地震信息的两个分量Zp、Zs；

利用射孔P/S波事件相对旅行时对上述两个分量Zp、Zs进行叠加处理，获得两个叠加道Zp_stk、Zs_stk；以及

利用改进能量比法，分别计算与两个叠加道相应的事件识别曲线，并基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置。

2.根据权利要求1所述的方法，其中，利用射孔P/S波事件相对旅行时对上述两个分量进行叠加处理包括：拾取射孔数据P波、S波事件同相轴相对第一个检波器旅行时时差tp、ts，并且利用tp、ts，对含有P波、S波微地震信息的两个分量Zp、Zs进行叠加处理，获得两个叠加单道Zp_stk、Zs_stk。

3.根据权利要求1所述的方法，其中，利用射孔P/S波事件相对旅行时对上述两个分量进行叠加处理包括：对两个分量Zp、Zs的每个样点振幅值Zp_i、Zs_i进行如下叠加处理，得到两个叠加道Zp_stk、Zs_stk：

Zp_stk(i)＝∑[Zp_i+(Zp_i+Δtp₂)+(Zp_i+Δtp₃)+......+(Zp_i+Δtp_N)]/N

Zs_stk(i)＝∑[Zs_i+(Zs_i+Δts₂)+(Zs_i+Δts₃)+......+(Zs_i+Δts_N)]/N、

其中，N为检波器个数，Zp_stk(i)、Zs_stk(i)为两个叠加道第i个样点的振幅值(0,Δtp₂,Δtp₃,......,Δtp_N)、(0,Δts₂,Δts₃,......,Δts_N)分别为拾取的P波、S波相对旅行时，i为样点序号。

4.根据权利要求1所述的方法，其中，基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置包括：索搜出事件识别曲线上大于门槛值且为区域极大值的值，将该值对应的时间作为井中微地震P/S波事件在第一个检波器出现的时间。

5.根据权利要求1所述的方法，其中，基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置包括：

1)针对P波和S波，分别在给定时间范围内在事件识别曲线上搜索第一个大于或等于门槛值的样点位置T_P、T_S：

2)针对P波和S波，在事件识别曲线上分别搜索给定时窗W_P、W_S范围内的区域极大值F_P1、F_S1，极大值F_P1、F_S1对应的时间样点Time(P1)、Time(S1)即为第一个可能的P/S波事件时间位置；

3)针对P波和S波，分别对事件识别曲线在上述时间样点Time(P1)、Time(S1)附近的振幅值进行归零处理；

4)以T_p以后的时间范围作为给定时间范围，重复执行1)-3)，直至识别出所有大于门槛值的P波事件时间以及S波事件时间。

6.一种井中微地震P/S波事件识别装置，该装置包括：

用于根据射孔数据P波事件水平分量偏振分析，获得检波器方位角，利用该检波器方位角对压裂数据水平分量进行方位一致性校正的部件；

用于基于校正结果，根据给定的极性判断准则对压裂数据Z分量进行分离，获得分别含有P波、S波微地震信息的两个分量Zp、Zs的部件；

用于利用射孔P/S波事件相对旅行时对上述两个分量Zp、Zs进行叠加处理，获得两个叠加道Zp_stk、Zs_stk的部件；以及

用于利用改进能量比法，分别计算与两个叠加道相应的事件识别曲线，并基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置的部 件。

7.根据权利要求6所述的装置，其中，利用射孔P/S波事件相对旅行时对上述两个分量进行叠加处理包括：拾取射孔数据P波、S波事件同相轴相对第一个检波器旅行时时差tp、ts，并且利用tp、ts，对含有P波、S波微地震信息的两个分量Zp、Zs进行叠加处理，获得两个叠加单道Zp_stk、Zs_stk。

8.根据权利要求6所述的装置，其中，利用射孔P/S波事件相对旅行时对上述两个分量进行叠加处理包括：对两个分量Zp、Zs的每个样点振幅值Zp_i、Zs_i进行如下叠加处理，得到两个叠加道Zp_stk、Zs_stk：

Zp_stk(i)＝∑[Zp_i+(Zp_i+Δtp₂)+(Zp_i+Δtp₃)+......+(Zp_i+Δtp_N)]/N

Zs_stk(i)＝∑[Zs_i+(Zs_i+Δts₂)+(Zs_i+Δts₃)+......+(Zs_i+Δts_N)]/N、

其中，N为检波器个数，Zp_stk(i)、Zs_stk(i)为两个叠加道第i个样点的振幅值(0,Δtp₂,Δtp₃,......,Δtp_N)、(0,Δts₂,Δts₃,......,Δts_N)分别为拾取的P波、S波相对旅行时，i为样点序号。

9.根据权利要求6所述的装置，其中，基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置包括：索搜出事件识别曲线上大于门槛值且为区域极大值的值，将该值对应的时间作为井中微地震P/S波事件在第一个检波器出现的时间。

10.根据权利要求6所述的装置，其中，基于给定门槛值在事件识别曲线上搜索出所有可能的P/S波事件时间位置包括：

1)针对P波和S波，分别在给定时间范围内在事件识别曲线上搜索第一个大于或等于门槛值的样点位置T_P、T_S：

2)针对P波和S波，在事件识别曲线上分别搜索给定时窗W_P、W_S范围内 的区域极大值F_P1、F_S1，极大值F_P1、F_S1对应的时间样点Time(P1)、Time(S1)即为第一个可能的P/S波事件时间位置；

3)针对P波和S波，分别对事件识别曲线在上述时间样点Time(P1)、Time(S1)附近的振幅值进行归零处理；

4)以T_p以后的时间范围作为给定时间范围，重复执行1)-3)，直至识别出所有大于门槛值的P波事件时间以及S波事件时间。