[发明专利]一种基于时间慢度域的偶极反射横波三维成像方法

权利要求书

1.一种基于时间慢度域的偶极反射横波三维成像方法，其特征在于，包括：

步骤(1)：对偶极阵列声波进行低通滤波，所述低通滤波的截止频率为4.7kHz；

步骤(2)：判断地层与井轴是否平行；

若地层与井轴平行，对偶极阵列声波进行组合，得到共发射器组合以及共接收器组合；

使用频率波数域滤波方法将共发射器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成下行反射横波；

使用频率波数域滤波方法将共接收器组合中不同慢度波场进行开窗滤波，保留慢度范围为的波场，形成上行反射横波；

若地层与井轴不平行，则对偶极阵列声波进行组合，得到共源距组合；

使用频率波数域滤波方法将共源距组合中不同慢度波场进行开窗滤波，保留慢度范围为-2s～0的波场，形成上行反射横波；保留慢度范围为0～2s的波场，形成下行反射横波；

步骤(3)：利用波形相干叠加法或慢度时间相关法对上行反射横波以及下行反射横波进行近似度分析，得到慢度-时间近似度关系图；

步骤(4)：将慢度-时间近似度关系图转换为深度-到井轴距离近似度关系图；慢度、时间、深度、到井轴距离之间满足：以及

其中，s表示慢度，Z表示发射器与界面之间的距离，v表示地层中传播的横波速度，T表示时间，z表示发射器与第一个接收器之间的距离，α表示界面与井轴夹角；

步骤(5)：对相同深度的数据进行叠加，形成深度与到井轴距离上的二维矩阵。

2.根据权利要求1所述的一种基于时间慢度域的偶极反射横波三维成像方法，其特征在于，利用波形相干叠加法对上行反射横波以及下行反射横波进行近似度分析，所述近似度满足：

其中：f_k为上行反射横波以及下行反射横波的波形数据，IW为窗口长度。

3.根据权利要求1所述的一种基于时间慢度域的偶极反射横波三维成像方法，其特征在于，利用慢度时间相关法对上行反射横波以及下行反射横波进行近似度分析，所述近似度满足：

其中：s_i(t)是第i个接收器的波形，i＝1,...,m；m是接收器个数；δ是接收间距；T_W是窗口长度。