叠前地震联合反演.ppt

叠前地震联合反演  描述岩石弹性特征的主要地震参数有岩石的弹性模量、密度 、纵波速度 和横波速度 等。弹性模量反映岩石应力—应变关系的特征，密度反映岩石的比重，速度则反映地震波在岩石中的传播的特征，它是弹性模量的函数。 声阻抗和弹性波阻抗曲线对比 子波提取 在不同角度的道集上进行合成记录标定 三、弹性阻抗反演 角道集资料 四、弹性阻抗反演 应用实例 角道集资料 弹性阻抗反演 通过对单井纵横波速度进行分析可知，气层纵波速度相对降低，横波速度变化较小，单独从纵横波速度难以对储层进行识别。 从纵横波交会图上可以看出： 右上角的红的正方形表示的为含气砂岩。纵横 波速度的交会图可以明显地将气层区分开来。 提取的纵波阻抗剖面 提取的纵波速度剖面 提取的密度剖面 提取的Lambda*Rho剖面 提取的泊松比剖面 提取的纵横波速度比剖面 提取的横波速度剖面 提取的横波阻抗度剖面 纵波阻抗时间切片 纵波速度时间切片 密度时间切片 Lambda*Rho时间切片 泊松比时间切片 纵横波速度比时间切片 横波阻抗时间切片 横波速度时间切片 谢谢大家 上图为提取的横波阻抗剖面图，因横波对气层反映不敏感，气层的横波速度与其它砂岩差别不大，用横波阻抗信息难以把气层区分开来。从上面的剖面上也很难看出目的气层的范围。所插入的测井曲线为纵波速度曲线，可看出曲线与反演出的剖面吻合不好。 根据AVO理论，零炮检距（或小炮检距）剖面可以近似为声阻抗AI（Acoustic Impedance）的函数，它与岩石的密度和纵波速度有关。为了充分利用大炮检距地震振幅信息，需要利用一种与入射角有关的弹性阻抗EI（Elastic Impedance），这种弹性阻抗形式上类似于声阻抗AI，不但与岩石的密度和纵波速度，而且与横波速度、入射角有关。 弹性波阻抗(EI) 弹性波阻抗反演 寻求一个函数，使其与声阻抗有相似的性质，以便于以任意角入射时反射系数均可以用下面的公式表示： 函数 称作EI（弹性波阻抗），当地层弹性参数变化不大时，用EI的对数值来表示反射系数是准确的 弹性波阻抗反演 弹性波阻抗反演流程 用测井曲线计算井中的EI曲线，对于同一个CMP道集来说，EI是入射角的函数 将偏移距数据转化为角道集数据， 估算每个角度数据所对应的子波 利用与迭后反演相类似的反演算法计算出各个角度的弹性波阻抗剖面 声阻抗 弹性波阻抗(8度) 弹性波阻抗(16度) 子波的振幅和频率特性随角度的变化而变化，这保证了反演处理的稳定性。各角度的子波在形状上很相似，振幅和相位也都比较一致，这有利于各角度资料反演效果的对比。 5-15度的井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线 5-15度/y551井 5-15度合成记录标定 15-25度的井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线 15-25度/y551井 15-25度合成记录标定 25-35度的井的合成记录层位标定图，从左至右分别为：角度子波、地震道、合成记录和测井曲线。右边的曲线分别为：红色的波阻抗曲线、蓝色的纵波速度曲线 25-35度/y551井 25-35度合成记录标定 5-15度角道集剖面 Ang5-15度 15-25度角道集剖面 Ang15-25度 25-35度角道集剖面 Ang25-35度 Ang5-15度 Complete_impedance_5-15度 弹性波阻抗的近似算法有多种，在此我们用BP近似来计算弹性波阻抗，根据Connolly 在1999年提出的弹性波阻抗近似公式可以得到弹性波阻抗的如下表达式： 进而由上式可以得到如下表达式： 对于同一角度来说A，B，C均为常数： 从弹性波阻抗体中提取岩性参数体 对于多个不同的角度而言有： 这是一个典型的 型方程组。写成矩阵形式： 上面表达式中：左边第一项是一个n×3的常数矩阵；第二项是所要求的纵、横波速度及密度的对数形式；右边是一个n×1的常数矩阵，它是由反演得到的各个角度的弹性波阻抗体分别取对数得到的。 从弹性波阻抗体中提取岩性参数体 通过最小平方法得到各样点的 可以进一步得到：Ip、Is、泊松比、纵波模量、横波模量及拉梅参数等。 对其分别取指数就得到各采样点的岩性参数体： 从弹性波阻抗体中提取岩性参数体 参数Lamda*Rho和Mu*Rho的获取 由于Lamda*Rho和Mu*Rho能够灵敏地反映储层属性，并且物理意义明确，在岩性和流体预测方面具有重大的意义，因此我们有必要从反演中拾取这两个属性。有两种方法可