GeoEast2D-3D反演介绍——20121205.pptVIP

* 叠前弹性参数反演 不同角度EI与AI的区别 Connoly首先提出了弹性阻抗这一概念，弹性阻抗不可测量，但可由纵、横波速度、密度、入射角度参数定量计算： 根据弹性阻抗反演结果，可以根据上式反算出纵波速度、横波速度、密度等弹性参数，如果需要计算出更反映流体特征的弹性参数如 、 等，则需要通过弹性参数之间的转换公式计算过，可能会引入误差。因此，将弹性阻抗公式改造成拉梅常数、剪切模量、密度表示的形式，可以直接反演出这些弹性参数： 实际反演中采用如下的标准化形式定义弹性阻抗，可以使得弹性阻抗的数值范围与常规阻抗的数值范围一致： 叠前弹性参数反演 为得到 、 和 需三个不同角度的弹性阻抗体。将三个角度值分别代入上式可以得到线性方程组，求解线性方程组可以得到结果。 利用三个不同角度弹性阻抗体可以提取拉梅常数和密度 直接提取弹性参数λ、μ、 对标准化后的弹性阻抗方程两边取对数得到： 叠前弹性参数反演 第一项：是约束地震数据和反演合成记录间的相近程度或残差； 第二项：是Cauchy约束，用来约束反演得到的反射系数的稀疏程度。 是反射系数系数约束项系数，即Sparse Constrain Factor，取值范围【0,80】； 第三项：是模型约束，用来控制反演结果的稳定性。即Model Constrain Factor，取值范围【0,20】。 叠前弹性参数反演 反演目标函数 叠前AVA同步反演 AVA叠前同步参数反演以角道集数据作为地震数据输入，利用反射系数精确公式，进行叠前AVA弹性参数反演，生成纵波速度、横波速度、密度，进而计算出波阻抗、体积模量、剪切模量、拉梅系数、泊松比和纵横比速度比等多种岩性和属性参数，为地震资料岩性解释、储层评价、流体识别等提供依据。 * CRP道集 速度谱点数据 均方根或层速度体 角道集 纵、横波、密度曲线 叠后标定 子波 AVA道集模拟 纵、横波速度、密度模型 AVA同步反演 Vp Vs Den 层位 时域曲线 Ip K Vp/Vs 叠前AVA同步反演 叠前AVA同步反演流程 AVA叠前同步反演要求输入角道集数据体、纵波速度、横波速度、密度模型数据体。 叠前AVA同步反演 叠前AVA同步反演 密度 纵波速度 横波速度 拉梅常数 剪切模量 体模量 纵波阻抗 纵横波速度比 泊松比 *.Rho.Ava *.Vp.Ava *.Vs.Ava *.Lambda.Ava *.Mu.Ava *.Bulk.Ava *.Ai.Ava *.VpVs.Ava *.Poisson.Ava 注意事项： 1、叠前AVA同步反演是基于角道集的反演，同步反演纵、横波速度及密度，道集质量严重影响到反演结果的好坏。因此，应选择信噪比高的道集数据体。 2、如果只能获得偏移距道集，则可以使用系统提供的角道集转换模块，将偏移距道集转成角道集，这时需要提供均方根速度谱点数据，并需要事先插成速度数据体。速度数据体的质量对角道集的质量至关重要。 3、反演前，应当进行道集正演（或道集标定），主要目的是检查目的层段的叠前同相轴与测井曲线正演的叠前同相轴是否相符。如果正演道集的振幅随偏移距（或角度）变化规律与实际道集不符，甚至相反，则反演结果肯定变差。 叠前AVA同步反演 * 叠前AVA同步反演 AVA叠前同步参数反演的核心是地震波振幅正演方程和反演方法。正演方程为垂直各向异性介质（VTI）的多波振幅特征方程，该方程不仅含有三个岩性参数，而且还有纵、横波两个各向异性系数，方程直接由各向异性介质的弹性波方程得出。 叠前反演方法特点是数据量大，反演方程阶数高。反演中应用全道集数据建立反演方程，通过LSQR方法求解大型矩阵方程，实现了地震波叠前岩性参数反演。可直接反演地层纵波速度、横波速度、密度及二个各向异性系数。在此基础上还可计算速度比、泊松比、弹性模量、剪切模量、拉梅系数、声阻抗及弹性阻抗梯度等参数。 叠前AVA同步反演 对于P波入射情况，VTI介质反射、透射系数可用下列方程表示： 其中： 叠前AVA同步反演 假设模型参数向量的初始估计为M0，则地震记录Y(M,x,t)在初始模型参量M0附近用Taylor级数展开并略去高阶项得 叠前AVA同步反演 叠前反演算法流程图 输入角道集、初始模型数据 归一化地震资料和子波 计算雅可比矩阵 使用LSQR方法解方程计算模型改变量 达到最大迭代次数 是 输出反演结果 否 AVA叠前同步参数反演 * GeoEast反演子系统的限制条件： 目前所有的反演模块： 只支持时间域的数据。 只支持32位格式的数据体，因此外部加载