[发明专利]一种从偶极横波资料中提取反射横波的方法

权利要求书

1.一种从偶极横波资料中提取反射横波的方法，其特征在于步骤如下：

第一步：对阵列声波测井采集到的原始偶极横波数据进行增益补偿和波形延时计算，得到采集到的原始横波数据Wave_dep；

第二步：根据仪器的发射频率对横波数据Wave_dep进行带通滤波，将低频噪音及高频干扰滤除，得到滤除噪音后的横波数据W_Ndep；

第三步：对横波数据W_Ndep利用时间慢度相干法进行处理，得到不同深度下慢度与时间的相关信息；将慢度转换为速度，从而得到横波速度的变化范围SV1～SV2，其中：SV1为横波最小速度，SV2为横波最大速度，单位m/s；

第四步：根据横波数据的速度变化范围SV1～SV2，建立横波数据对应的传播矩阵，并通过传播矩阵计算出相应的直达横波数据A_l(ω)exp(-iωz/v_l)；

第五步：从第二步处理得到的横波数据W_Ndep中减去第四步处理中计算出的直达横波数据A_l(ω)exp(-iωz/v_l)，得到反射横波数据。

2.根据权利要求1所述的一种从偶极横波资料中提取反射横波的方法，其特征在于：所述第一步中，具体计算方法如下：

Ori_dep＝Signal_dep/Gn_dep

Wave_dep＝Ori_dep-Tfst_dep

上式中，Gn_dep表示仪器记录到的波形增益数据；Signal_dep表示仪器通过自动增益处理过后的波形数据，Ori_dep表示通过增益恢复后的横波数据；Tfst_dep为仪器采集的波形延时数据，Wave_dep表示通过波形延时处理后得到的原始横波数据。

3.根据权利要求2所述的一种从偶极横波资料中提取反射横波的方法，其特征在于：所述第四步中，具体计算方法如下：

在横波速度变化范围(SV1～SV2)内取值，得到L个横波速度值v₁,v₂,...v_L,从而建立直达横波数据的传播模式A_l(ω)exp(-iωz/v_l),(l＝1,....,L)，其中A_l代表第l种直达横波的时间域信号，ω为角频率，A_l(ω)代表第l种直达横波的频率域信号，v_l表示第l种直达横波的传播速度，z表示深度p的位置；

当直达横波从第p个深度位置传至第q个深度位置时，可以把其传播模式表示为A_l(ω)exp(-iω(q-p)/v_l)；则在第p个深度位置可表示为下式：

ElNdep-p(ω)·Al(ω)=WNdep(ω)]]>

上式中，ω代表角频率，Ndep代表本次测井所记录的深度采样点数，表示将波形从时间域转换到频率域的转换符，W_Ndep(ω)为第二步处理后得到的横波数据的频率域形式；将上式简化为矩阵形式可得：

EgA＝W

对于上述矩阵，采用最小二乘法可求得A：

其中，:表示复共轭，T表示转置，通过上式便可求出第p个深度位置所有直达横波所对应的频率域数据A_l(ω)。