波阻抗反演和储层预测汇.pptVIP

波阻抗反演和储层预测 中国地质大学地球物理系 朱培民 内容 波阻抗反演 碳酸盐岩储层预测 地震波阻抗反演——基本概念（1） 反映岩石性质的基本参数 速度、密度 孔隙度、渗透率 反映岩石性质间接参数 声波阻抗 弹性波阻抗 反射系数 声波阻抗 定义： 波阻抗=速度*密度 声波阻抗与岩性的关系 声阻抗在实际的地震勘探中，有很重要的作用。如果能够从地震数据中估计出声阻抗，和与横波速度有关的参数，我们就能较方便地区分出地层的岩性以及所包含的流体成份。 反射系数 在地震水平分层的情况下，地震波的反射系数定义成： 什么是波阻抗反演？ 简单地说，波阻抗反演就是把叠后地震数据集中的每一道变换为一个伪声波阻抗（Pseudo-Acoustic Impedance）曲线的过程。 地震波阻抗反演——意义 地震道 地震道模型 地震剖面 波阻抗剖面 波阻抗道 波阻抗剖面 地震道模型 合成记录使用的Ricker子波 长度为0.24秒，采样间隔为2毫秒，主频为10 Hz的Ricker地震子波 地震剖面和波阻抗剖面的比较 实际的地震道模型(1) 实际的地震道模型(2) 带限的反射系数序列可以表示成 r(t) r(t)：平均反射系数 a(t)： 零相位带限平均函数（残余子波） 根据反射系数计算波阻抗 如果地震道是带限的，则由此计算的波阻抗信息也是带限的 实际声波阻抗和地震道反演波阻抗 地震各种信息的带宽 地震反演的目的和作用 地震反演的主要目的就是把地震反射数据转换成定量的岩性数据，用于储层描述(孔隙度和产层厚度)。 由于反演后便于地震解释，可以补偿反演耗费的时间，提高地震解释的效率。 由于测井输入的分辨率高，地震数据的井约束反演有可能得到比地震数据更高的分辨率 波阻抗反演的好处 由于波阻抗反演综合了地震，测井，岩性和地质解释的多源信息，波阻抗反演结果，包含更多的信息，且不受地震子波调谐和干涉的影响 声波阻抗是地层的岩性信息，而地震数据只是反映了界面的信息，是波阻抗的相对变化。因此波阻抗是联系地震数据和井数据的一个桥梁。 波阻抗与岩性、孔隙度、孔隙填充物等有更直接的关系，有利于储层描述和流体分析。 波阻抗数据更有利于层序地层分析，也更有利于目标解释。 带限（Band-Limited）的含义 地震频带以外的波阻抗在反演的过程中会丢失，或者原始数据资料中缺失这些频率。 地震数据的带限严重影响地震资料及其波阻抗反演的分辨率 地震波阻抗的反演方法 叠后 利用偏移地震资料，井约束地震反演技术。 叠前 CDP（Common Depth Point）道集上反演纵波速度、横波速度和密度的方式有好几种，但精度都不高。一种较稳定可靠的方法是，在部分叠加的基础上进行叠后反演。通过近偏移距叠加反演，可以在测井曲线标定的基础上直接计算声阻抗。 地震波阻抗反演——算法 道积分(Trace Integration) 基于地震道，地震数据是唯一输入，是一种递归算法，反演结果的带宽决定于地震带宽。 层块反演(Layer-based or blocky inversion) 稀疏脉冲反演(Sparse spike inversion) 最小平方反演(Least-squares inversion) 道积分的原理(1) 假设地震道数据经过反褶积，每一个样点都可近似看作反射系数，并记作r(t)。 则根据反射系数的定义 r = (Z2-Z1)/(Z2+Z1)，进一步假设波阻抗是连续变化的，记作Z(t)，由此，在Δt时间内反射系数的变化，可以写成： r (t) Δt = ΔZ / (2Z+ ΔZ) 假设Δt —〉0，则成立下式： r (t) dt = dz / 2Z 两边积分得， 道积分的原理(2) 在离散情况下： 一个气藏的道积分实例 时间域递归反演 基本假设 从（井）低频层速度或波阻抗计算的反射系数序列 + 地震道的高频反射系数序列=宽带反射系数序列 基本原理 计算道积分，得到一个反映局部波阻抗变化的高频估计； 将根据层速度或声波测井估计的低频波阻抗成分叠加到道积分结果上，得到时间域的宽带波阻抗反演结果 如果把密度看作常数，则波阻抗反演结果，可以看成是虚速度曲线（Pseudo-Velocity） 层块反演(Layer-based or blocky inversion) 这种算法把地层细分成薄层，和地震数据联系就是褶积模型； 通常把地层细分和地层采样点的个数一样多，即反射系数序列密的，有时也看作白噪序列 这种算法只要通过初始模型的合适选取来稳定反演过程。 反射系数序列是密的 合成记录使用的Ricker子波 合成地震记