一种基于杨氏模量、泊松比的两项式地震叠前反演方法技术

本发明专利技术属于地震勘探技术领域，公开了一种基于杨氏模量、泊松比的两项式地震叠前反演方法，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：叠前地震资料预处理；构造层位解释；测井资料外推与模型构建；多角度子波提取；最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。本发明专利技术通过发展两项AVO反射系数方程，将密度项以合理的方式与纵波速度建立数学关系，消除密度项，可以在缺少大角度叠前地震资料的情况下进行叠前直接反演，得到目标弹性参数，可以减小传统方法带来的累积误差。

A two item seismic pre stack inversion method based on Young's modulus and Poisson's ratio

The invention belongs to the field of seismic exploration technology, and discloses a binomial seismic pre-stack inversion method based on Young's modulus and Poisson's ratio. The binomial seismic pre-stack inversion method based on Young's modulus and Poisson's ratio includes the following steps: pre-stack seismic data processing, structural horizon interpretation, logging data extrapolation and model. Construction; Multi-angle wavelet extraction; Finally, the simulated annealing inversion method is used to obtain Young's modulus and Poisson's ratio data volume. By developing two AVO reflection coefficient equations, the density term is mathematically related to P-wave velocity in a reasonable way, and the density term is eliminated. The pre-stack direct inversion can be carried out in the absence of large-angle pre-stack seismic data, and the elastic parameters of the target can be obtained. The cumulative error caused by the traditional method can be reduced.

【技术实现步骤摘要】
一种基于杨氏模量、泊松比的两项式地震叠前反演方法
本专利技术属于地震勘探
，尤其涉及一种基于杨氏模量、泊松比的两项式地震叠前反演方法。
技术介绍
当前全球地震勘探技术发展的几个主要趋势包括：从常规油气向常规与非常规油气并重局面发展、从简单地表和浅水区向复杂地表和深水区发展、从叠后地震反演向叠前弹性反演发展。页岩油气的开采需压裂页岩生储层，因此岩区的易压裂程度成为评价优质页岩油气储层非常重要的指标。岩石物理研究表明，杨氏模量和泊松比能够较好的表征岩体的造缝难易程度，标识页岩气的“甜点”区；由于地震测线长度与储层埋深的比值通常很小，故通常很难获得超深储层的大反射倾角叠前地震数据。叠前AVO反演可以充分利用叠前地震资料包含的地质、岩性和流体信息，获取丰富的地下介质弹性参数。AVO：(Amplitudevariationwithoffset)利用地震反射波振幅随偏移距的变化特征来探讨反射系数响应随入射角的变化,进而确定反射界面上覆、下伏介质的岩性特征及物性参数的技术。主要用于石油天然气勘探。叠前AVO反演的核心是反射系数方程。通常反射系数方程由三项构成，其中之一为密度项。密度项通常不作为直接用来预测油气储层的指标，原因是具体工区的密度变化范围不大，而各种岩性的理论密度范围交集较广，且其对油气储层特征的直接相关性差，单由密度项无法进行岩性、流体的判断。在缺乏大角度地震资料或者地震资料数据信噪比较低的情况下，大角度反演的结果品质较差，提取参数的过程中，参数矩阵的稳定性较差，造成最终提取的弹性参数误差较大。另外，由于密度项的参数权值在角度有限情况下较小，反演的难度更加大。以YPD三项反射系数方程为基础，构建的YPD叠前地震反演方法，同时反演了杨氏模量、泊松比与密度。该方法为目前较新的杨氏模量、泊松比反演方法(2014年建立)。具体做法如下：假设待反演杨氏模量、泊松比及密度反射系数服从柯西分布，假设似然函数服从高斯分布，同时在反演目标函数中加入初始模型约束，并通过初始模型建立各道去相关矩阵，消除待反演参数间的互相关性，在贝叶斯反演框架下进行反演。该方法在叠前地震数据远角道集缺失的情况下，大角度入射角的引入会使参数矩阵的稳定性降低，造成最终提取的弹性参数存在较大误差。综上所述，现有技术存在的问题是：针对深层、超深层页岩气储层的地震勘探，由于远角道集数据的缺失，使得利用AVO叠前反演较难获取地下介质准确的杨氏模量、泊松比等岩石物理弹性参数。因此，提高参数矩阵稳定性，降低计算带来的累积误差成为亟需解决的问题。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种基于杨氏模量、泊松比的两项式地震叠前反演方法。本专利技术是这样实现的，一种基于杨氏模量、泊松比的两项式地震叠前反演方法，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：步骤一，叠前地震资料预处理；步骤二，构造层位解释；步骤三，测井资料外推与模型构建；步骤四，多角度子波提取；步骤五，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。进一步，所述基于杨氏模量、泊松比的两项式地震叠前反演方法具体包括以下步骤：(1)基于Aki-Richards近似方程推导得到YP两项式反射系数近似方程，作为深层杨氏模量、泊松比叠前反演方法的理论基础；Aki-Richards近似方程表述为：其中，Vp、Vs分别为纵、横波速度，ρ为密度，θ为入射角。分别是纵、横波速度反射系数与密度反射系数。(2)采用幂指数拟合的方式，引入参数可变的关系式消除密度项，由密度与纵波速度的微商关系：dρ＝FLVpL-1dVp(2)可得纵横波模量反射系数与纵横波速度、密度反射系数的关系如下所示：将(3)带入(4)(5)得令(3)整理杨氏模量、泊松比关于纵横波速度的关系，并融合以上各式，得将(9)、(10)、(11)带入Aki-Richard方程(1)，得到新的两项式近似方程其中，k为横、纵波速度比的平方，L为幂指数。进一步，所述基于杨氏模量、泊松比的两项式地震叠前反演方法在进行杨氏模量与泊松比反演计算以前，需从先验信息中求得参数L、k；在实施操作过程中，其中k为测井资料中横、纵波速度比的平方，L则通过密度与纵波速度的幂指数拟合计算得到。待反演参数的计算可简化为求解：其中，CE和Cσ与角度有关；RE和Rσ分别是杨氏模量、泊松比的反射系数时间序列反射系数与某一角度子波卷积形成一个固定角度的地震合成记录。将式(13)等号两边同时与子波卷积可推导为式SPP(t，θi)是入射角为θi的角道集，W(t，θi)是一个子波核矩阵。卷积运算之后，可得若设D＝Gm，SPP表示包含向量D的叠前角度道集，WE和Wσ构成小波核矩阵G，RE和Rσ是包含m的反射率项。反演的目标函数为f(V)＝||S-D||→min(21)式中，V为待反演弹性参数构成的参数向量；S(t)＝W*R为初始低频模型在各角度下的合成地震记录，W为一组不同角度地震子波，R为计算得到的低频地震模型的反射系数；D为实际部分叠加道集地震数据；构造f(V)的Jacobian矩阵J，通过模拟退火算法对J进行求解，得到对V的修正量ΔV0，将V0＝(V+ΔV0)作为新的反演参数向量带入式(21)，反复迭代求解f(V)的极小值，直至求得使f(V)达到允许误差的ΔVn，此时认为Vn＝(Vn-1+ΔVn)为最优解；对m的多次迭代最终值即为反演得到的最终杨氏模量、泊松比数据体。本专利技术的优点及积极效果为：通过发展两项AVO反射系数方程，将密度项以合理的方式与纵波速度建立数学关系，消除密度项，可以在缺少大角度叠前地震资料的情况下进行叠前直接反演，得到目标弹性参数，这样可以减小传统方法带来的累积误差；实际数据来自CGGHampson-Russell公司开发软件的AVO模块的测试数据。对地震资料进行预处理后抽取角道集数据，利用HRS软件依据测井资料对角道集数据提取五组不同角度的子波，并依据所建立的不同子波核矩阵与YP近似方程计算的反射系数合成出不同角度域地震剖面，带入反演流程进行迭代计算；得到最终反演结果的杨氏模量、泊松比反演剖面，反演结果与实际情况较为吻合，显示了本专利技术较好的应用潜力。附图说明图1是本专利技术实施例提供的基于杨氏模量、泊松比的两项式地震叠前反演方法流程图。图2是本专利技术实施例提供的YP近似方程地震叠前反演流程图。图3是本专利技术实施例提供的正波阻抗界面反射系数对比示意图；图中：(a)不同方程反射系数对比；(b)不同方程反射系数差值对比。图4是本专利技术实施例提供的负波阻抗界面反射系数对比示意图；图中：(a)不同方程反射系数对比；(b)不同方程反射系数差值对比。图5是本专利技术实施例提供的层状介质模型反演结果示意图；图中：(a)角度域合成地震记录；(b)杨氏模量、泊松比反演曲线对比。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。下面结合附图对本专利技术的应用原理作详细的描述。如图1所示，本专利技术实施例提供的基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：S101：叠前地震资料预处理；S102：地质层位解释；S103：测井资料外推与模型构建；S104本文档来自技高网...

【技术保护点】
1.一种基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：步骤一，叠前地震资料预处理；步骤二，构造层位解释；步骤三，测井资料外推与模型构建；步骤四，多角度子波提取；步骤五，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。

【技术特征摘要】
1.一种基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法包括以下步骤：步骤一，叠前地震资料预处理；步骤二，构造层位解释；步骤三，测井资料外推与模型构建；步骤四，多角度子波提取；步骤五，最终利用模拟退火反演方法获得杨氏模量和泊松比数据体。2.如权利要求1所述的基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法具体包括以下步骤：(1)基于Aki-Richards近似方程推导得到YP两项式反射系数近似方程，作为深层杨氏模量、泊松比叠前反演方法的理论基础；Aki-Richards近似方程表述为：其中，Vp、Vs分别为纵、横波速度，ρ为密度，θ为入射角，分别是纵、横波速度反射系数与密度反射系数；(2)采用幂指数拟合的方式，引入参数可变的关系式消除密度项，由密度与纵波速度的微商关系：dρ＝FLVpL-1dVp(2)可得纵横波模量反射系数与纵横波速度、密度反射系数的关系如下所示：将(3)带入(4)(5)得令(3)整理杨氏模量、泊松比关于纵横波速度的关系，并融合以上各式，得将(9)、(10)、(11)带入Aki-Richard方程(1)，得到新的两项式近似方程其中，k为横、纵波速度比的平方，L为幂指数。3.如权利要求1所述的基于杨氏模量、泊松比的两项式地震叠前反演方法，其特征在于，所述基于杨氏模量、泊松比的两项式地震叠前反演方法在进行杨氏模量与泊松比反演计算...

【专利技术属性】
技术研发人员：张瑞，文晓涛，刘婷，曹其壮，杨吉鑫，李雷豪，
申请(专利权)人：成都理工大学，中国石油化工股份有限公司勘探分公司，
类型：发明
国别省市：四川,51