[发明专利]从阵列声波测井资料中提取反射波的方法

摘要

本发明公开了一种从阵列声波测井资料中提取反射波的方法，包括：根据阵列声波测井资料提供的全波列数据，提取出直达波信号和反射波信号，结合纵波速度、横波速度及斯通利波速度及仪器结构参数，对提取的直达波信号再次进行反射波提取，得到平行井轴的反射体产生的反射波，合成后得到完整的反射波。本发明充分利用了直达纵波、横波和斯通利波的速度信息且结合了仪器结构参数，对提取的直达波信号再次进行反射波提取，得到了平行井轴的反射体的反射波信号。 1

主权项

1.一种从阵列声波测井资料中提取反射波的方法，其特征在于：根据阵列声波测井资料提供的全波列数据，提取出直达波信号和反射波信号，结合纵波速度、横波速度及斯通利波速度及仪器结构参数，对提取的直达波信号再次进行反射波提取，得到平行井轴的反射体产生的反射波，合成后得到完整的反射波；

所述方法具体包括如下步骤：

a、以阵列声波原始声波波形数据作为输入，采用中值滤波进行反射波提取，令W_to_l为原始波形数据，W_zg处理得到的直达波波形数据，W_f1为处理得到的反射波波形数据，具体计算方法如下：

W_f1(i)＝W_tol(i)‑W_zg(i)

上式中，Ndep为深度行数，median为中值滤波器，L为滤波器取中值的跨度，W_zg(i)为第i个深度直达波形数据，为第i个深度下的中值滤波器中第j个波的原始波形数据，W_f1(i)为第i个深度反射波形数据；

b、对阵列声波测井数据进行处理，得到斯通利波速度V_st，横波速度V_s及纵波速度V_p；

c、利用步骤a中处理得到的直达波波形数据W_zg，结合步骤b得到的波速及仪器结构参数，分步滤除斯通利波W_st、直达横波W_s、直达纵波W_p，最终得到反映平行井轴的反射体的反射波W_f2；

d、对提取的两个反射波W_f1、W_f2进行合成，得到完整的反射波W_f：

W_f＝W_f1+W_f2。

2.根据权利要求1所述的从阵列声波测井资料中提取反射波的方法，其特征在于：所述滤除斯通利波具体为：

其中，为第i个深度t时刻的斯通利波直达波，为第i个深度下第j个接收器按△T_stij移动后的波形，△T_stij为第i个深度下第j个接收器的斯通利波时差，为第i个深度滤除斯通利波直达波后的波形，median为中值滤波器，N为接阵列声波仪器接收器个数，d为接收器距声源的距离，△T为不同接收器接收到直达波的时差。

3.根据权利要求2所述的从阵列声波测井资料中提取反射波的方法，其特征在于：所述滤除直达横波具体为：

其中，为第i个深度t时刻的横波直达波，为第i个深度下第j个接收器按△T_sij移动后的波形，△T_sij为第i个深度下第j个接收器的横波时差，为第i个深度滤横波直达波后的波形，median为中值滤波器，N为接阵列声波仪器接收器个数，d为接收器距声源的距离，△T为不同接收器接收到直达波的时差。

4.根据权利要求3所述的从阵列声波测井资料中提取反射波的方法，其特征在于：所述滤除直达纵波具体为：

其中，