地震波阻抗反演和储层预测

地震波阻抗反演和储层预测第1页，课件共89页，创作于2023年2月内容波阻抗反演碳酸盐岩储层预测第2页，课件共89页，创作于2023年2月地震波阻抗反演——基本概念（1）反映岩石性质的基本参数速度、密度孔隙度、渗透率反映岩石性质间接参数声波阻抗弹性波阻抗反射系数第3页，课件共89页，创作于2023年2月声波阻抗定义：

波阻抗=速度*密度第4页，课件共89页，创作于2023年2月声波阻抗与岩性的关系声阻抗在实际的地震勘探中，有很重要的作用。如果能够从地震数据中估计出声阻抗，和与横波速度有关的参数，我们就能较方便地区分出地层的岩性以及所包含的流体成份。第5页，课件共89页，创作于2023年2月反射系数在地震水平分层的情况下，地震波的反射系数定义成：第6页，课件共89页，创作于2023年2月什么是波阻抗反演？简单地说，波阻抗反演就是把叠后地震数据集中的每一道变换为一个伪声波阻抗（Pseudo-AcousticImpedance）曲线的过程。第7页，课件共89页，创作于2023年2月地震波阻抗反演——意义地震道地震道模型地震剖面波阻抗剖面波阻抗道波阻抗剖面第8页，课件共89页，创作于2023年2月地震道模型w(t):地震子波r(t):反射系数序列s(t):地震信号地震道模型的基本假设条件——水平层状模型第9页，课件共89页，创作于2023年2月地震楔状模型低波阻抗区地震剖面高波阻抗背景波阻抗场tuningeffects第10页，课件共89页，创作于2023年2月合成记录使用的Ricker子波长度为0.24秒，采样间隔为2毫秒，主频为10Hz的Ricker地震子波

第11页，课件共89页，创作于2023年2月反演结果：波阻抗剖面彩色地震剖面地震剖面和波阻抗剖面的比较反演结果——波阻抗剖面能够正确地反映地质模型和岩性有较高的分辨率tuningeffects调谐效应消失第12页，课件共89页，创作于2023年2月实际的地震道模型(1)一般情形下反射系数序列r(t)应是宽带的，但实际的r(t)是带限的。地震道模型的基本假设条件——水平层状模型第13页，课件共89页，创作于2023年2月实际的地震道模型(2)带限的反射系数序列可以表示成<r(t)><r(t)>：平均反射系数a(t)：零相位带限平均函数（残余子波）第14页，课件共89页，创作于2023年2月根据反射系数计算波阻抗如果地震道是带限的，则由此计算的波阻抗信息也是带限的第15页，课件共89页，创作于2023年2月实际声波阻抗和地震道反演波阻抗第16页，课件共89页，创作于2023年2月地震各种信息的带宽原始反射系数序列：宽带实际地震反射数据：带限10~80Hz速度：

带限1-3Hz速度模型+<r(t)>：0~80Hz（缺失5~12Hz）第17页，课件共89页，创作于2023年2月地震反演的目的和作用地震反演的主要目的就是把地震反射数据转换成定量的岩性数据，用于储层描述(孔隙度和产层厚度)。由于反演后便于地震解释，可以补偿反演耗费的时间，提高地震解释的效率。由于测井输入的分辨率高，地震数据的井约束反演有可能得到比地震数据更高的分辨率第18页，课件共89页，创作于2023年2月反演结果——波阻抗剖面能够正确地反映地质模型和岩性子波地震道反演波阻抗道波阻抗模型第19页，课件共89页，创作于2023年2月一个反演的实际例子第20页，课件共89页，创作于2023年2月第21页，课件共89页，创作于2023年2月波阻抗反演的好处由于波阻抗反演综合了地震，测井，岩性和地质解释的多源信息，波阻抗反演结果，包含更多的信息，且不受地震子波调谐和干涉的影响声波阻抗是地层的岩性信息，而地震数据只是反映了界面的信息，是波阻抗的相对变化。因此波阻抗是联系地震数据和井数据的一个桥梁。波阻抗与岩性、孔隙度、孔隙填充物等有更直接的关系，有利于储层描述和流体分析。波阻抗数据更有利于层序地层分析，也更有利于目标解释。第22页，课件共89页，创作于2023年2月频带对反演结果的影响——带宽越宽，分辨率越高带限（Band-Limited）对反演的影响石油地震中的子波频带第23页，课件共89页，创作于2023年2月带限（Band-Limited）的含义地震频带以外的波阻抗在反演的过程中会丢失，或者原始数据资料中缺失这些频率。地震数据的带限严重影响地震资料及其波阻抗反演的分辨率第24页，课件共89页，创作于2023年2月地震波阻抗的反演方法叠后利用偏移地震资料，井约束地震反演技术。叠前CDP（CommonDepthPoint）道集上反演纵波速度、横波速度和密度的方式有好几种，但精度都不高。一种较稳定可靠的方法是，在部分叠加的基础上进行叠后反演。通过近偏移距叠加反演，可以在测井曲线标定的基础上直接计算声阻抗。第25页，课件共89页，创作于2023年2月地震波阻抗反演——算法道积分(TraceIntegration)基于地震道，地震数据是唯一输入，是一种递归算法，反演结果的带宽决定于地震带宽。层块反演(Layer-basedorblockyinversion)稀疏脉冲反演(Sparsespikeinversion)最小平方反演(Least-squaresinversion)后三种都是基于模型的宽带约束反演算法，差别在于如何处理非地震信息以补偿地震数据缺乏的高低频信息。第26页，课件共89页，创作于2023年2月道积分的原理(1)假设地震道数据经过反褶积，每一个样点都可近似看作反射系数，并记作r(t)。则根据反射系数的定义

r=(Z2-Z1)/(Z2+Z1)，进一步假设波阻抗是连续变化的，记作Z(t)，由此，在Δt时间内反射系数的变化，可以写成：

r(t)Δt

=ΔZ/(2Z+

ΔZ)假设Δt—〉0，则成立下式：

r(t)dt=dz/2Z

两边积分得，

第27页，课件共89页，创作于2023年2月道积分的原理(2)在离散情况下：第28页，课件共89页，创作于2023年2月一个气藏的道积分实例低波阻抗区圈闭的顶第29页，课件共89页，创作于2023年2月时间域递归反演基本假设从（井）低频层速度或波阻抗计算的反射系数序列+地震道的高频反射系数序列=宽带反射系数序列基本原理计算道积分，得到一个反映局部波阻抗变化的高频估计；将根据层速度或声波测井估计的低频波阻抗成分叠加到道积分结果上，得到时间域的宽带波阻抗反演结果如果把密度看作常数，则波阻抗反演结果，可以看成是虚速度曲线（Pseudo-Velocity）第30页，课件共89页，创作于2023年2月递归反演流程第31页，课件共89页，创作于2023年2月层块反演(Layer-basedorblockyinversion)这种算法把地层细分成薄层，和地震数据联系就是褶积模型；通常把地层细分和地层采样点的个数一样多，即反射系数序列密的，有时也看作白噪序列这种算法只要通过初始模型的合适选取来稳定反演过程。第32页，课件共89页，创作于2023年2月反射系数序列是密的第33页，课件共89页，创作于2023年2月合成记录使用的Ricker子波第34页，课件共89页，创作于2023年2月合成地震记录道第35页，课件共89页，创作于2023年2月稀疏脉冲反演(Sparsespikeinversion)假设地震反射只有有限的层（层比较厚），即反射系数序列以及相应的波阻抗是稀疏的；地震道数据可以模拟成少数的几个反射系数与子波的褶积，这就是稀疏脉冲的含义；为了重建地震带宽以外的频率，SSI方法使用模型来稳定和约束反演过程。第36页，课件共89页，创作于2023年2月稀疏脉冲反演流程第37页，课件共89页，创作于2023年2月最小平方反演(Least-squaresinversion)除了不考虑反射系数是稀疏的以外，其他的类似与稀疏脉冲反演最小平方反演不拓宽高频谱段，而是在初始模型的基础上，通过稳定初始模型，来恢复低频信息。第38页，课件共89页，创作于2023年2月波阻抗反演算法的评价除了道积分外，其他三种都是宽带算法；层块、稀疏脉冲和最小平方三种方法，在某种程度上都能消去子波，消除调谐效应，提高分辨率。第39页，课件共89页，创作于2023年2月波阻抗反演的目标函数R:反射系数序列D:地震道数据S:合成地震记录Z:波阻抗T:合成波阻抗测井Λλ：阻尼系数αα：阻尼系数BlockyL1SparseSpikeLpModel-BasedL2第40页，课件共89页，创作于2023年2月目标函数各项的意义反射系数约束地震道匹配地震道与合成地震记录道匹配带限低频反射系数——稀疏脉冲测井曲线匹配（先验信息，priorinformation）基于模型的波阻抗道匹配补足带限外的高低频第41页，课件共89页，创作于2023年2月第42页，课件共89页，创作于2023年2月Jason中实际使用的目标函数r:反射系数序列n=d-s:地震数据匹配Ω：过渡矩阵，开始为单位阵λ：阻尼系数第43页，课件共89页，创作于2023年2月反演过程内循环：先保持第一项不变，对地震数据匹配项作最小平方法反演，使目标函数趋于极小。当内循环反演稳定后，保持第二项不变，开始外循环反演过程。Ω的值在外循环迭代的过程中，不断修改，取前一次反射系数的结果，作为Ω的值。第44页，课件共89页，创作于2023年2月波阻抗反演的策略可以通过模型的整体控制，补偿地震信息不存在的低频成分（低频信息可来自于测井、叠前和叠后偏移速度分析的结果）通过测井曲线建立的模型，补偿地震信息欠缺的高频成分在稀疏脉冲反演初始迭代时，不考虑基于模型的目标函数项，将地震数据和模型的贡献分开。在反演的最后阶段，使用目标函数的最后一项，完成地震信息和模型之间的过渡频率反演最佳匹配。第45页，课件共89页，创作于2023年2月低频信息的反演在用测井信息获得的波阻抗测井曲线中，低频是天然存在的。通常，我们可以通过井控制的地质解释得出的波阻抗模型来获取低频信息。低频信息可以在反演的最后阶段添加进去，或者在目标函数中直接反映出来。在后一种情况下，只需要模型中非常低的那一部分频率在Jason中，上述的两种方法可以使用，在某种程度上说，这两种方法都是基于模型的。第46页，课件共89页，创作于2023年2月低频信息的来源测井曲线中的低频成分叠加速度叠前时间偏移和叠前深度偏移中使用的速度第47页，课件共89页，创作于2023年2月约束——相对反演用层位控制反演的范围(时间门，TimeGate)，中间部分插值实现。问题：上述的约束方式会不会造成反演严重依赖于地震数据的层位解释。答案：是。解决方案1：（1）在第一遍反演时放松约束，允许精度不高的层位解释；（2）在最终反演之前根据第一遍反演的结果，修正层位，给出更严格的约束。解决方案2：在第一遍反演后，修正层位时不要考虑基于模型的低频部分，同样可以达到目的。第48页，课件共89页，创作于2023年2月基于模型约束的优点可以消除多次波的影响第49页，课件共89页，创作于2023年2月子波的提取和处理在反演过程中，要确保子波的相位、频率和地震数据是匹配的。在测井曲线作过正确的时深转换，并进行过合理编辑的情况下，如果反演结果同测井曲线不匹配，则说明子波不正确。第50页，课件共89页，创作于2023年2月子波提取的原则在反演的目标区，提取子波。深层提取的子波不能用于浅层，反之亦然。具有不正确相位谱和振幅谱的子波，会导致时移，会产生多余的子波叶瓣，由此产生假的地质层位，反演结果也不可能和测井曲线匹配。第51页，课件共89页，创作于2023年2月反演结果的质量控制(QC)基于地震道的反演：应该用测井曲线（波阻抗曲线）控制反演的结果基于模型的反演：应该以合成纪录和地震记录的匹配，作为质量控制的标准第52页，课件共89页，创作于2023年2月基于地震道反演的质量控制(QC)稀疏脉冲反演在井旁得到的波阻抗曲线结果，不是算法直接计算得到的，是用户约束和目标函数设置控制的输出。很显然，这样的输出与根据测井曲线直接计算出的结果不可能一致。因此，它们的差可以作为质量控制的标准。用测井曲线计算出的波阻抗曲线，须先滤波到地震的频带，然后才能和反演出的波阻抗曲线，作为质量控制的标准。第53页，课件共89页，创作于2023年2月反演波阻抗体的质量检查（1）用没有参加反演的井，对反演出的波阻抗体进行检查（盲检：BlindTest）。同样的，和检查波阻抗曲线的方法类似，在检查比较前，须对测井曲线作合适的校正，并将测井曲线标定到地震剖面上。第54页，课件共89页，创作于2023年2月反演波阻抗体的质量检查（2）相对波阻抗检查方法对反演的波阻抗体作带限滤波，主要滤去低频成分，得到与地震数据体类似的相对波阻抗数据体；在数据体上加上层位和构造信息，相对波阻抗体应仍然能正确地反映出地质异常信息。如果全频带能反映出异常信息，而带限的相对波阻抗体不能，则说明低频信息模型有问题。第55页，课件共89页，创作于2023年2月地震波阻抗反演——数据偏移或叠后地震数据体地震解释层位测井曲线井位声波时差密度自然电位速度第56页，课件共89页，创作于2023年2月地震波阻抗反演——流程数据准备+地震数据处理+解释+检查子波提取和估计建立模型井旁反演外推得到整个剖面的反演结果第57页，课件共89页，创作于2023年2月数据准备和检查在反演前最重要的工作是检测测井曲线和由此得到的波阻抗曲线与地震剖面是否符合（匹配）测井曲线应转换到时间域，并将其滤波到地震数据的频带。不匹配的测井曲线，应作适当的编辑第58页，课件共89页，创作于2023年2月反演结果解释—波阻抗体的显示和构造解释反演波阻抗体的显示波阻抗值是正值为了在地震工作站上进行解释，可以采用以下方法：减去一个常数滤波形成相对波阻抗用常规的地震构造解释方法进行波阻抗的构造解释第59页，课件共89页，创作于2023年2月反演结果解释—属性转换波阻抗是岩性数据波阻抗也可以转换到其他属性（孔隙度，砂泥岩百分比）转换后的属性可以进行统计（体积、交会图、直方图、方差图）第60页，课件共89页，创作于2023年2月属性解释的步骤（1）利用测井数据，在给定的目标区把已知的岩性与给定频率范围内的反演波阻抗值建立起关系；第61页，课件共89页，创作于2023年2月低自然伽玛、高电阻率和低声波阻抗可以作为含碳氢化合物砂岩的标志自然伽玛和声波阻抗的交会图含碳氢化合物砂岩的岩性区第62页，课件共89页，创作于2023年2月属性解释的步骤（2）在波阻抗体上划出一个感兴趣的区域。例如一个由解释层位+时间（或深度）圈定的地层，一个感兴趣的岩性单元等。根据前一步建立的波阻抗——岩性关系，在此目标区，把不符合上述关系的数据点作透明处理，显示出的点就是感兴趣岩性的区域。第63页，课件共89页，创作于2023年2月根据前述的波阻抗和高电阻率、低自然伽玛的含碳氢化合物砂岩标志，将不符合上述关系点透明化后的显示结果第64页，课件共89页，创作于2023年2月属性解释的步骤（3）把那些有开发价值的和彼此连通的“地质体”（波阻抗数据）用不同的颜色表示出来。输出结果可包括顶底层位、层厚度、以及“地质体”内的实际属性值。第65页，课件共89页，创作于2023年2月第66页，课件共89页，创作于2023年2月属性解释的步骤（3）计算“地质体”体积、“地质体”内的孔隙度、产层厚度（netpay）和其他与声波阻抗相关的属性。因为是三维数据体，也很容易地绘出它的构造图和厚度图。第67页，课件共89页，创作于2023年2月属性解释的步骤（4）——储层模拟将单个的“地质体”转换到深度域，和由声波阻抗计算出的属性值（如孔隙度）直接输出到储层模拟的软件中。进行储层模拟和储层描述。第68页，课件共89页，创作于2023年2月第69页，课件共89页，创作于2023年2月Jason反演的数据和流程第70页，课件共89页，创作于2023年2月第71页，课件共89页，创作于2023年2月实际反演中的问题（1）测井曲线的分辨率在0.1~1m之间，频谱范围在0~125Hz之间地震信号的分辨率约10ms，频谱范围在10~80Hz之间。如果地震波的传播速度在4000m/s时，能够分辨的地层厚度约为4000m/s*10ms/2=20m第72页，课件共89页，创作于2023年2月声波阻抗的频带波阻抗的谱分布波阻抗的全谱：0~125Hz声波阻抗趋势：0~3Hz低频声波阻抗模型：0~10Hz地震+低频：0~80Hz

第73页，课件共89页，创作于2023年2月声波测井和地震频率的分布波阻抗的谱分布波阻抗的全谱：0~125Hz声波阻抗趋势：0~3Hz低频波阻抗模型：0~10Hz地震+低频：0~80Hz

第74页，课件共89页，创作于2023年2月声波时差和密度测井曲线到合成地震记录的制作过程第75页，课件共89页，创作于2023年2月Jason约束稀疏脉冲反演的频率补偿地震信息频带在15~80Hz之间，根据地震数据反演的波阻抗也在15~80Hz之间。Jason利用速度模拟(VelMod)，可以用叠加速度重建趋势模型频率(0~3Hz)。Jason�中3-15Hz的低频，在约束稀疏脉冲反演(CSSI)中通过添加趋势模型和约束来实现的。第76页，课件共89页，创作于2023年2月Jason反演低频信息补偿3~15Hz低频信息补偿的实现在CSSI反演中低频信息直接加入到目标函数中，会使反演对噪音异常敏感。因此，不在CSSI反演过程中加入，而是在反演后，用VelMod和（或）来自于EarthMod的井间插值模型联合实现第77页，课件共89页，创作于2023年2月弹性阻抗Connolly定义的