[发明专利]利用纵波、转换波数据联合反演储层纵横波速度比的方法

权利要求书

1.一种利用纵波、转换波数据联合反演储层纵横波速度比的方法，其特征在于，所述方法包括下述步骤：

输入地震数据，从而基于相同成像点位置的纵波成像道集和转换波成像道集，根据目的储层所在时间范围来选取纵波反演时窗[T_a,T_b]和转换波反演时窗[T_c,T_d]；

根据选取的时窗的每一个纵波成像道集和转换波成像道集，建立纵波反射系数表达式和转换波反射系数表达式：

Rpp(θ)=12(ΔVpVp+Δρρ)+(12ΔVpVp-4γ2ΔVsVs-2γ2Δρρ)sin2θ]]>

Rps(θ)=(-2γΔVsVs-2+γ2γΔρρ)sinθ]]>

其中，R_pp为纵波的反射系数，R_ps为转换波的反射系数，θ为纵波入射角，V_p为纵波速度，V_s为横波速度，ρ为介质的密度，γ为纵横波速度比；

根据建立的纵波反射系数表达式和转换波反射系数表达式，得到与成像道集地震反射信息相关的γ变化率表达式：

Δγγ|a=2γan+1-γanγan+1-γan=(2+γ2)A-γ2B+(1+γ2)C]]>

其中，为与成像道集地震反射信息相关的γ变化率，

A为纵波成像道集在选定时窗内给定样点的截距，其计算式为：

A=12(ΔVpVp+Δρρ)]]>

B为纵波成像道集在选定时窗内给定样点的梯度，其计算式为：

B=(12ΔVpVp-4γ2ΔVsVs-2γ2Δρρ)]]>

C为转换波成像道集在选定时窗内给定样点的梯度，其计算式为；

C=(-2γΔVsVs-2+γ2γΔρρ);]]>

在选定时窗内纵波成像道集和转换波成像道集的给定样点处，根据旅行时信息得到纵横波速度比：

γt=2ΔTpsΔTpp-1]]>

其中，γ_t为纵横波速度比，ΔT_pp和ΔT_ps分别为给定样点的上下层旅行时差；

根据得到的纵横波速度比，建立与旅行时信息相关的γ变化率表达式：

Δγγ|t=2γtn+1-γtnγtn+1-γtn]]>

其中，为与旅行时信息相关的γ变化率；

对于与成像道集地震反射信息相关的γ变化率和与旅行时信息相关的γ变化率，在选定时窗内所有样点建立目标函数：

E=Σi=1i=N[Δγγ|ai-Δγγ|ti]2]]>

其中，i＝1,2,…,N；

采用基于阻尼最小二乘解的广义线性反演方法，求解如下方程组：

∂E∂(Δγ/γ)1=0,∂E∂(Δγ/γ)2=0,......,∂E∂(Δγ/γ)N=0]]>

从而反演得到所述成像点位置每个样点的纵横波速度比γ_i。