[发明专利]基于正则化的非线性地震叠前弹性参数反演方法

权利要求书

1.一种基于正则化的非线性地震叠前弹性参数反演方法，包括以下步骤：

根据已知的测井资料与射线弹性阻抗REI利用下式得出多角度的射线弹性阻抗井曲线REI_i，其中，所述测井资料中包括纵波速度、横波速度、密度与入射角，

REIi=REI(α,β,ρ,θi)=ραcosθi(1-β2α2sin2θi)4,---(23)]]>

在上式中，α为纵波速度，β为横波速度，ρ为岩石密度，θ为入射角；

引入已知的层位、断层资料以及上述多角度的射线弹性阻抗井曲线REI_i，以建立整个工区不同角度的射线弹性阻抗模型；

使用叠前偏移距道集，计算出上述不同入射角的部分叠加数据体；

将多角度的部分叠加数据体、连同多角度的射线弹性阻抗模型，通过稀疏脉冲反演算法进行反演，得出上述入射角对应的射线弹性阻抗数据体；

将上述入射角对应的射线弹性阻抗数据体，导入非线性Landweber-Fridman迭代算法，以得出整个工区的纵、横波速度和岩石密度；

根据纵波速度α、横波速度β和岩石密度ρ，利用下式得出整个工区的岩性物理参数，

σ=0.5*(α/β)2-1(α/β)2-1,---(24)]]>

λρ=IP2-2*Is2,---(25)]]>

λ/μ＝(α/β)²-2，                                （26）

在上式中，I_p为纵波阻抗ρα，I_s为横波阻抗ρβ，σ为泊松比，λ为拉梅常数，μ为岩石的剪切模量。