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

摘要

本发明提供一种利用纵波、转换波数据联合反演储层纵横波速度比的方法，所述方法采用道集地震反射信息获取的γ变化率与旅行时信息获取的γ变化率无限逼近这一思想，构建反演目标函数，目标函数的关键参数包括纵波截距、纵波梯度、转换波梯度及旅行时信息计算得到的纵横波速度比，这些参数均可由实际输入的纵波、转换波道集的AVO信息和旅行时信息获取，因此目标函数的构建更加合理，最终反演的纵横波速度比弹性参数也更加精确。

主权项

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。