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

摘要

本发明提供一种基于正则化的非线性地震叠前弹性参数反演方法，包括以下步骤：根据已知的测井资料与射线弹性阻抗得出多角度的射线弹性阻抗井曲线REIi；引入已知的层位、断层资料以及上述多角度的射线弹性阻抗井曲线REIi，以建立整个工区不同角度的射线弹性阻抗模型；使用叠前偏移距道集，计算出上述不同入射角的部分叠加数据体；通过稀疏脉冲反演算法进行反演，得出上述入射角对应的射线弹性阻抗数据体；将上述入射角对应的射线弹性阻抗数据体，导入非线性Landweber-Fridman迭代算法，以得出整个工区的纵、横波速度和岩石密度，进而计算出其它岩性参数。本发明利用多个角度的REIi数据，并通过非线性Landweber-Fridman迭代法进行反演计算，以获得准确的岩性参数。

主权项

一种基于正则化的非线性地震叠前弹性参数反演方法，包括以下步骤：根据已知的测井资料与射线弹性阻抗REI利用下式得出多角度的射线弹性阻抗井曲线REIi，其中，所述测井资料中包括纵波速度、横波速度、密度与入射角， REI i = REI ( α , β , ρ , θ i ) = ρα cos θ i ( 1 - β 2 α 2 sin 2 θ i ) 4 , - - - ( 23 ) 在上式中，α为纵波速度，β为横波速度，ρ为岩石密度，θ为入射角；引入已知的层位、断层资料以及上述多角度的射线弹性阻抗井曲线REIi，以建立整个工区不同角度的射线弹性阻抗模型；使用叠前偏移距道集，计算出上述不同入射角的部分叠加数据体；将多角度的部分叠加数据体、连同多角度的射线弹性阻抗模型，通过稀疏脉冲反演算法进行反演，得出上述入射角对应的射线弹性阻抗数据体；将上述入射角对应的射线弹性阻抗数据体，导入非线性Landweber‑Fridman迭代算法，以得出整个工区的纵、横波速度和岩石密度；根据纵波速度α、横波速度β和岩石密度ρ，利用下式得出整个工区的岩性物理参数， σ = 0.5 * ( α / β ) 2 - 1 ( α / β ) 2 - 1 , - - - ( 24 ) λρ = I P 2 - 2 * I s 2 , - - - ( 25 ) λ/μ＝(α/β)2‑2，                                （26）在上式中，Ip为纵波阻抗ρα，Is为横波阻抗ρβ，σ为泊松比，λ为拉梅常数，μ为岩石的剪切模量。