地震反演原理及其应用ppt课件.ppt

《地震反演原理及其应用ppt课件.ppt》由会员分享，可在线阅读，更多相关《地震反演原理及其应用ppt课件.ppt（90页珍藏版）》请在三一办公上搜索。

1、地震反演原理及其应用,引言：,地震波阻抗反演是地震数字处理的终结结果之一，是油气藏勘探、开发中储层横向预测的重要手段。而且随着勘探、开发程度的提高，对反演资料的精度要求也越来越高。但是，由于现有的地震资料受地震波激发和采集系统本身条件的限制，就决定了所记录的地震波仅有相似的振幅信息，分辨率有限且包含各种干扰。因此，波阻抗反演无疑是一个多解的非线性问题。,油气田开发的工作多是针对储层进行的。而地震勘探长期以来只是利用岩层的声学特征确定岩性的分界面。这就使地震资料与油田地质的结合发生困难。为了使地震资料能与钻井资料直接连接对比，就要把界面型的反射剖面转换成岩层型的测井剖面，把地震资料变成可与钻井直

2、接对比的形式。实现这种转换的处理过程就是地震反演技术。,引言,地震反演原理及其在油田开发中的应用,主要内容地震反演概述与地震反演技术相关的几个概念地震反演方法,所谓反演就是正演模型处理的逆过程。反演处理的数学算法实际上只能和正演模型相同。正演模型是比较简单的。当前关于正演模型的处理技术没有争议，对于任一给定的地质模型只能有一个地震响应与之对应。反演过程要复杂得多：某些正演模型的过程可能没有逆。对于一个给定的地震响应也可能会有多个地质模型与之对应。,地质模型,地震响应,模型计算,反演计算,地质模型,地震响应,输入,处理,输出,正演模型（模拟）,反演模型（反演）,地震反演原理及其在油田开发中的应用

3、,地震反演概述：,所谓正演过程就是人工模拟地震记录的形成。即 而反演过程则是估算一个子波的逆反子波，用反子波和地震道进行褶积运算，得到反射系数 r(t),然后由上式导出的递推公式逐层递推计算出每一层的波阻抗，即 用Gardner公式从波阻抗中分离出速度和密度。,地震反演概述,地震记录的形成和反演原理,正演合成记录,反演合成声波测井,地震反演方法,叠前反演,旅行时间反演（层析X射线成象）,振幅反演(AVO),叠后反演,振幅反演,递归反演,稀疏脉冲反演,基于模型的反演,波动方程反演,地震反演概述,目前在生产中用的较多的是叠后振幅反演。较为流行的应用软件有STRATA、Jason、ISIS、VanG

4、urad 等。其以一维褶积模型为理论基础，仅使用指定的地质界面的振幅信息（而不是它的反射时间）。其基本反演方法包括有递归反演、稀疏脉冲反演和基于模型的反演等。 另外，近年来由于地震属性分析技术、神经网络识别技术及弹性物理模量参数的引入，使传统的反演方法在原有的基础上有了长足的发展。例如：以多元地震属性分析为基础的储层测井特性的模拟技术（非阻抗反演方法）；弹性反演及以神经网络识别为特征的数值逼近技术等。使储层识别及其物性预测的精度和可靠性都有较大地提高。,地震反演概述,波动方程反演技术是近几年发展起来的一种非常前沿的非线性叠后反演技术。通过地震波的数值模拟来研究地下地层分布及储层物性、含油气性。

5、 其基本思想认为：含油、气储层是非均匀各向异性的具有固体状态与流体状态的三维三相介质。建立在均匀各向同性纯固体基础上的岩石弹性理论与波动传播理论，难以逼真描绘地震波在含油、气地层中传播的复杂物理现象。为此，从实际应用出发，研究横向各向同性的非均匀的粘弹介质与双相弹性介质中地震波传播的数值模拟方法，以研究地震波在地层中传播的物理机制及其与含油、气地层的定量关系。 目前，以物探局范祯祥先生开发的非线性波动方程地震波反演技术软件ANGEL2001为代表的地震反演软件包已投入生产应用。本技术的基本技术思路是：采用有限元波动方程对地震波剖面数值模拟，借助于最优选择理论对地震波反映的物性参数与几何形态进行

6、逐次逼近，以实现物性参数的反演。在此基础上借助于神经网络分析对所反演的物性参数进行非线性标定。小波技术，模糊识别技术，分形技术穿插应用其中。,地震反演概述,AVO反演技术是一种使用叠前动校后的角度道集资料来识别储层岩性及其含油气性的反演方法，尤其以识别气层而见长。 此方法的基本思想是：地震资料采集时，检波器是按一定的炮检距布置在远离炮点的位置上，地震射线是以非零入射角入射到岩层界面上。而并非在一维褶积模型中关于入射波垂直入射的假设。反射系数不仅与纵波速度和密度有关，还是横波速度的函数。反射系数随炮检距的变化而变化。所谓 AVO (Amplitude Versus Offset)一词就是“振幅随

7、炮检距的变化”。由于振幅随炮检距的变化是纵、横波速度和密度的函数；因此就可以对常规的叠加前共反射点道集上的振幅随炮检距的变化进行反演，提取纵、横波速度和密度。 在实际应用中，由于砂岩储层含气后会造成砂岩速度的降低。利用低速度标志可以圈定气藏的边界。但由于流体受力不产生剪切形变，故横波传播速度不受砂岩孔隙中所含流体成分的影响，只反映砂岩骨架的传播速度。因此，联合应用反演得到的纵、横波剖面就可以达到检测气层的目的。,所有的反演方法都不可避免地带有“非唯一性”的问题。这意味着对于一特定的地震响应，可能有多个地质模型与之相对应。在各种可能的地质模型之间来决定那一个与实际地质模型相符合的唯一途径是使用其

8、他信息（且是不包含在地震数据中的）。这种附加信息通常是由下面两种途径来提供的：初始猜测模型最终结果偏离初始猜测程度的约束条件。这意味最终反演结果除地震资料之外还要依赖于这种“附加信息”。,地震数据,最佳剖面,约束条件,地质模型,合并及转换,最终的反演,地震反演概述,与地震反演技术相关的几个概念,褶积模型反射系数地震子波分辨率与地层反褶积的应用剩余相位校正,地震反演原理及其在油田开发中的应用,地震道可由下式表示: 地震记录 =子波 反射系数 + 噪声即 （1）这里 意思是褶积或 （2）,与地震反演技术相关的几个概念,褶积模型：,关于褶积模型的假设条件:叠后数据 地震道是零炮间距的；没有多次波；没

9、有 AVO 效应；噪声是随机的, 即是白噪的, 与地震不相关； 没有相干噪声。子波是恒定的 不随时间变化。地震数据已经是偏移过的 每一道地震记录仅仅取决于地震道位置垂直向下的反射系数序列。,w(t),R(t),S(t),在时间域，地震道可由下列方程表示: 地震记录 =子波 反射系数 + 噪声,与地震反演技术相关的几个概念,在频率域，褶积就是反射系数的频谱与子波的频谱的乘积。 其结果意味我们所能得到的地震记录的频谱仅能保留部分与子波频谱相对应的成分，这样就丢失了对应于反射系数谱中高频和低频的内容。 反演就是要试图恢复这些丢失的频谱成分即根据反演算法来充填丢失成分的方法。,与地震反演技术相关的几个

10、概念,反演可认为是，给定地震道 s(t) ，确定其反射系数的处理过程。在方程(1)中，反射系数与地层阻抗有关，其关系可用下式表示 （3）这里反演的目的是为了从地震道本身估算地层速度。显然，首先需要从褶积模型中提取反射系数的一个估算值，因此，我们可进行反褶积的相关处理。在讨论反褶积和反演之前，让我们更详细地讨论一下褶积模型的两个主要组成部分： 反射系数和地震子波。,与地震反演技术相关的几个概念,反射系数,无论是炸药震源还是一个与地面接触的平板的突然撞击（可控震源），当地震震源的能量被释放时，能量被作为弹性波通过地球传播。我们知道的最简单类型的波是压缩波（即纵波或P波）。一种岩石允许声波通过的能力

11、由它的声阻抗决定，它是压缩波速度和岩石密度的乘积。正如我们期望的，岩石越致密，它的声阻抗就越高；如压实砂岩的声阻抗通常比泥岩的声阻抗高；另外，孔隙度对速度也有影响，孔隙度越高，速度、密度越低。层状介质中，每当声阻抗有变化时，就会在具有不同阻抗的两层之间的界面处产生一个地震反射。假如声波传播到两个具有不同声阻抗层之间界面，反射系数可被写为方程（3）。方程（3）表明，反射系数可以为正或负，取决于I(j-1)或I(j)哪一个较大，且其绝对值不能大于1。前面做的分析是针对一个下行波而言的，它从上方传播到界面。上行波的反射系数将被简单地作为下行波反射系数的负数。,与地震反演技术相关的几个概念,显而易见，

12、不是所有的入射能量都能被反射（只有海上的记录在空气与水的界面处几乎全被反射，它的反射系数接近于1）。透射能量的值为入射振幅与反射能量的差，或是 （4） 注意：方程（3）和方程（4）只适用于P波垂直入射一个界面的情况。对于非垂直入射，将发生类型转换，产生反射和透射横波。在进行地震岩性精确测量时，了解这些影响是很重要的。对于叠后情况，我们将假设方程（3）能准确地预测反射系数。,与地震反演技术相关的几个概念,在多数反演算法中对地层反射系数做如下假设：反射系数是平稳随机的白噪序列有了这个统计学假设，许多反褶积方法才得以实现。但事实上，地层反射系数不可能是白噪的。地层沉积总是具有旋回性，一个岩层的顶界反

13、射系数的后面通常跟着这个岩层底面的反射系数。为了使反射系数尽量满足白噪假设，在资料处理中通常采用一些方法对地震道自相关函数进行修正，以使地震道自相关函数中的反射系数部分逼近白噪。在反褶积、反演计算过程中通常是在零时间自相关值上加一个常数，抬高零时间上的自相关值达到提高反射系数谱的白噪化程度白噪系数。,与地震反演技术相关的几个概念,地震子波是接受到的来自一个单个界面的反射波（不是一个单位脉冲），一个在时间上具有一定延续长度的波动单元。地震子波的波形决定于震源激发波形、表层面波干涉、近地表条件变化、大地滤波特性、检波器和记录仪器滤波特性、资料处理对波形的滤波作用等，因此，要确切地知道子波是很难的。

14、在频率域，子波的确定由两部分组成：确定振幅谱确定相位谱（相位谱的确定迄今仍然很难，且是反演的主要误差源）相位谱类型有： 零相位、常相位、最小相位或非最小相位等。,与地震反演技术相关的几个概念,地震子波,与地震反演技术相关的几个概念,最小相位子波,零相位子波,与地震反演技术相关的几个概念,常相位子波,非最小相位子波,对零相位和常相位子波而言，可将其看作是一系列不同振幅和频率正弦波的集合，所有的正弦波都是零相位或常相位(如900)。 最小相位不是一个简单而易理解的概念，可以有许多不同的方法描述它。最简单直观的方法是，最小相位在零时间之前不存在(数学上，我们称这种子波是物理可实现的)，并且能量主要分

15、集中在前部。事实上，最小相位的术语是指在所有的有相同振幅谱和不同相位谱的因果子波中，相位谱最接近于零的那一个。,与地震反演技术相关的几个概念,子波的确定方法大体上可分为类：直接确定性的完全统计性的使用测井的直接确定法指直接用地面接收器和其他方法测量子波。 完全统计法只用地震数据确定子波。这种方法很难准确 地确定子波的相位谱。 测井曲线法指除使用地震数据外，还使用测井曲线信息。 子波是空变且时变的，对地震剖面应该确定大量的子波。但实际反演的应用上，试图确定可变的子波可能引起更大的不确定性。因此，实际使用的是对整个剖面提取一个“平均”子波。,与地震反演技术相关的几个概念,统计方法提取子波 使用自相

16、关统计原理，从地震数据中提取地震道的振幅谱来作为子波的振幅谱。对子波的相位谱则做最小相位或零相位的假设。这样就可以唯一地确定子波。在此方法中，子波长度是一个关键的参数，该参数决定了地震道振幅谱的平滑程度，当子波长度增加时，子波谱就接近于分析时窗内地震数据的谱。测井曲线法提取子波 指除使用地震数据外，还使用测井曲线信息。从原理上讲，这样可以在井点位置提取准确的相位信息。问题是，这种方法关键依赖于测井曲线与地震数据之间有一个好的连接关系标定准确，因为将深度采样的测井曲线转换为双程旅行时采样的测井曲线时，可能会出现降低结果质量的闭合差。,与地震反演技术相关的几个概念,影响提取子波的主要因素：,使用测

17、井曲线提取子波的方法极其依赖井和地震道之间的匹配关系，由于深度采样的测井曲线在转换成双程旅行时间的测井曲线时，会出现误差，从而降低了子波的质量。地震记录本身是非零相位的，而用于标定的合成记录确是零相位的，在地震道和合成记录之间必然存在一个相位差。因此，要想获得一个好的子波，除认真分析测井与地震道的对应关系做好标定外，还应根据抽取子波的相位谱对地震道进行相位校正处理。,与地震反演技术相关的几个概念,分辨率与地层反褶积的应用,分辨率是地震反演的一个很关键的问题。 在油气田勘探中，地震反演多是应用于 大套的储层识别。但在开发中，不仅要求能识别储层及描述储层分布（横向分辨率），而且还要求能划分出小层（

18、垂直分辨率）。因此，如何能提高地震分辨率，成为开发地震技术中重要的一环。,与地震反演技术相关的几个概念,常规地震记录由于是一个有限带宽的信 号，故其分辨率是有限的。既没有10Hz 以下的低频信息，又缺60-70Hz以上的 高频成分。频谱多表现为一个近似钟形 的振幅曲线。 地层反射系数是由测井声阻抗计算的， 频带很宽。频谱上的大峰值反映了具有 周期性的优势地层分量。尤其是在陆相 地层中薄互层占有相当的比重，致使高 频成分很丰富。,a. 地层反射系数频谱,b. 地震反射频谱,如果能使有限的地震带宽展成类似于地层反射系数谱的特征，就能使地震分辨能力的得到提高。,与地震反演技术相关的几个概念,不同频带

19、对薄层的分辨能力,使地震频带变窄的原因：,震源是带限的；大地是个低通滤波器，越深的反射频率越低；检波器是高通的，不能记录10Hz以下的低频；震源组合和检波器组合是低通的；记录仪器也是带限的；资料处理中的所有叠加和偏移等都是低通的。 尤其是对于高频端，绝大多数原因都限制了对高频的接收。因此提高地震记录的分辨率是非常困难的。对以上任何方面的改进，都将有助于改善地震分辨率。,与地震反演技术相关的几个概念,提高地震分辨率的手段：,高分辨率采集技术包括使用小药量多井组合激发；使用高频检波器接收；小道距、小组合基距接收，并提高采样率。在处理中采用高分辨分频处理流程。反褶积处理 反褶积处理技术是建立在褶积模

20、型的基础上的。反褶积就是计算一个反子波，与地震道褶积（相当于对地震道进行滤波），消除子波的褶积效应，从而获得一个近似的地层反射系数系列，以便用于反演高分辨率的岩层型剖面。其方法根据对子波的处理方式分为两类：统计反褶积子波反褶积（地层反褶积）,与地震反演技术相关的几个概念,常规剖面与高分辨率剖面对比图,与地震反演技术相关的几个概念,与地震反演技术相关的几个概念,统计性反褶积： 对褶积模型的基本假设为：地震道的自相关等于子波的自相关。由于在得到的自相关函数中没有子波的相位信息，故必须对子波做一个最小相位或零相位的假设。只要知道振幅谱，就可以唯一地确定子波。这样就可以计算出反褶积滤波器。如脉冲反褶积

21、。最小相位反褶积零相位反褶积,零相位统计反褶积 + 带通滤波,原始地震记录,与地震反演技术相关的几个概念,原始地震记录,最小相位统计反褶积 + 带通滤波,与地震反演技术相关的几个概念,与地震反演技术相关的几个概念,统计反褶积的使用多是为了回避确定子波的困难。实际上，上述假设不是完全正确的。而且只有在输入满足子波相位假设的条件下，其结果才是最佳的。因此基于此褶积模型和子波假设的反褶积，效果不会太好。,子波反褶积 在已知子波 的前提下，按反褶积滤波器计算出反子波 f（t）：即 (5) 用反子波 f(t) 对地震道进行确定性子波处理，从而得到反射系数。如在Strata反演软件中就用了一种叫做双边脉冲

22、反褶积的方法做确定性子波反褶积。另外，为了控制噪声水平和分辨率，还可以加入白噪因子和算子长度等参数来优化反褶积。,与地震反演技术相关的几个概念,原始地震记录,反褶积后的地震记录,与地震反演技术相关的几个概念,原始地震记录,子波反褶积+带通滤波,与地震反演技术相关的几个概念,与地震反演技术相关的几个概念,合成地震记录,反射系数,无噪声记录子波处理,含10%地震噪声记录子波处理,子波反褶积也有一个较大的弱点抗噪能力差。如果要获得一个理想的结果，原始输入的噪声水平必须有所控制。,与地震反演技术相关的几个概念,反褶积的应用在一定程度上消除了子波的褶积效应，进而获得一近似的反射系数剖面。事实上，要想完全

23、消除子波的褶积效应是不可能的。因此只能是在一定程度上展宽了反射波的频谱，改善反射剖面的分辨率。并最终使反演结果的分辨率有所提高。地震反射剖面中仍残留有剩余子波。而且在多数情况下，反褶积之后剖面的信噪比有所下降，且在有效频宽之外常有假频出现。一般而言，信噪比的降低并不影响反演的最终结果，因为在模型反演中是没有噪声引入的。但信噪比也不能降的太低，以至于无法解释、追踪有效波的同相轴，使初始模型更加偏离实际的地层模型。至于假频可用带通滤波加以切除。,剩余相位校正,在地震反演中，要把地震反射剖面反演成岩层型剖面，从理论上讲地震剖面必须是零相位的。一般而言，应用零相位化子波反褶积原则上能把反射剖面处理成零

24、相位剖面。但是，若子波估算有误差，也就难做到真正的零相位化。相位校正通常是通过纯相位整形滤波来实现的。而相位校正量则是通过使用测井资料计算子波的相位谱而得到的。因而根据使用测井资料提取子波的方法不同，其相位校正的结果亦不相同。,与地震反演技术相关的几个概念,常相位在提取子波的的包络内，最大波峰若偏离零时刻的时间为Dt，提取子波即存在一初始相位角 f0 ，有f0 相当于地震道与合成记录之间相位差的平均值，是一个与频率无关的常数相位角。常相位的相位校正，就是把这个常相位f0 校正为零值，使地震道零相位化。,f,f,f,t,f,f0,f,f,与地震反演技术相关的几个概念,非线形相位严格说来，实际子波

25、的波形是变频的。子波相位谱也不是一个常数，多表现为一个有斜率的曲线。如果去掉斜率，就是一个沿截距f0 摆动的曲线。用这个曲线做相位校正，就可以同时把常相位和非线形相位都校正掉了。,f0,t,f,f,f,f,f,f,与地震反演技术相关的几个概念,一般而言，常相位校正作为一个常规处理手段，已能使反演的对井效果得到基本保证。若要做得更好一些，就应在正确标定的基础上，用实际提取子波的相位谱做非线形相位校正。 子波提取与相位校正这两项处理工作应当结合起来进行，但前提必须是建立在正确、合理的井震标定的基础之上。相位校正不但解决了子波估算中相位差的问题，同时还解决了地震剖面的极性问题。不论输入的 反射剖面是

26、什么极性，经相位校正后，剖面的极性都将趋于和合成记录一致即，接近于零相位剖面，至少在井旁道附近是这样。,与地震反演技术相关的几个概念,地震反演方法：,叠后地震反演也是一种处理过程，通过地震反演来分析地震道，并且试图重建地层的速度和声阻抗。反演基于的基础模型是一维褶积模型。反演可以看成是给定地震道s(i)时，确定反射系数r(i)的过程。 依据对子波处理及处理与反演有关的非唯一解问题的不同，一般地可分为以下三种类型：递归反演基于模型的反演稀疏脉冲反演,地震反演原理及其在油田开发中的应用,递归反演,递归反演（也称带限反演）是反演方法中最简单的一种。其基本思想是：根据反射系数公式 (3) 第 j 层的

27、阻抗 I(i) 可以由第 j-1 层的阻抗 I(i-1) 求得，即 (6) 若取 I(0) 作为波阻抗的初始值，则上式(6)变为 (7) 根据此式就可以连续地求出每一层的阻抗值。,地震反演方法,但在实际应用中，是将上式取对数后做泰勒级数展开，整理得 (8) 上式右边是反射系数 r(j) 的变限求和。在实际资料处理中，把反褶积后的地震道 s(t) 当作反射系数 r(j) 使用。因此，对地震道s(t) 做变限求和，而作为直流分量的 I(0) 已被滤掉，得到的是有限带宽的对数相对波阻抗，公式(8)将改写为 (9)这就是应用于递归反演处理的最主要的数学过程。,地震反演方法,地震反演方法,这是在理想条件

28、下递归反演的应用。,一般地，递归反演（带限反演）包含三个步骤：使用声波测井、RMS速度或两者结合建立一个低频速度模型。使用（9）式的递归(Recursive)反演流程反演地震道，给出声阻抗或速度的中间频带（1060Hz）部分。把低频和中频信息结合起来，计算反演的输出。,地震反演方法,地震数据,测井数据,井震标定,反褶积,提取子波,地震解释,低频模型,变限求和,反演结果,递归反演处理流程,递归反演方法存在的问题：,用于反演的地震道的样点是实际的反射系数序列，就是说在地震记录中没有子波的褶积效用。被视为反射系数的地震道的样点做了适当的归一化处理，其数据范围介于-1，1之间。地震数据的样点可以是任何

29、一个振幅值。且递归反演公式假定振幅的绝对值是正确的。由于递归反演公式是应用以递推的方式从道的顶端推向底部，其可能产生的误差的影响是累积的。这就意味着在底部的误差可能比顶部要大。误差累计的最大影响还在于最终反演结果的低频分量或低频趋势。通常这种低频趋势很难定义，而且非常贫乏。所以习惯上常常从解中去掉这种趋势，而代之以初始猜测模型的低频分量趋势（如STRATA、JASON）。递归反演之所以也称之为有限带宽反演是因为最终反演阻抗的频带与地震数据的频带相同。,地震反演方法,影响因素：,地震反演方法,递归反演中子波的影响： 三点子波被递归反演看作是三个分离的反射系数。,递归反演中限定带宽的影响： 反演结

30、果显示递归反演存在两个问题：低频成分丢失，并且在反演结果中没有恢复出来。其结果只能是一个相对阻抗。子波被忽视，从而导致反演结果中产生错误的旁瓣。,地震反演方法,地震反演方法,递归反演中噪声的影响： 地震道中的任何噪声都被递归反演算法解释为一个反射系数，并能产生一个假的同相轴。,地震反演方法,带限反演原始地震数据,地震反演方法,带限反演合成声阻抗,小结：,递归反演是一种早期开发的叠后反演过程。在原理上，这种技术是很简单的。如果我们假设地震道代表地层反射系数的一种近似，那么这种反射系数可以被反演来给出声阻抗。然而地震道是带限的，反射系数的低频端（010Hz）和高频端（60Hz以上）均丢失。另外，地

31、震道的相位误差和噪声污染也是使这一问题复杂化的因素。 因此，其结果的分辨率是很有限的，对于识别较薄的含油气储层是不可行的。但是递归反演也有其不可替代的优点，因为在反演的整个过程中没有人为约束条件的加入，且多解性小。因此，其结果具有很高的真实性。对于识别和了解大套储层及其宏观展布特征也还是十分有用的。,地震反演方法,基于模型的反演,模型反演是从地震正演模型技术发展而来的一种反演技术。其基本原理是：首先使用测井和地震解释成果建立对油藏地质模型的初始猜测初始模型（用速度、密度表示）。通过正演算法制作合成地震剖面地震模型。然后，将地震模型与实际地震剖面比较，根据比较结果，反复修正初始模型，计算新的地震

32、模型，以最佳地符合地震资料。当地震模型与地震资料对比误差足够小时即满足期望的边界条件，此时的油藏地质模型就是反演的最终结果。 模型反演的算法:基于褶积模型的基于弹性波波动方程理论,地震反演方法,合成剖面,子波,最佳拟合？,模型修改,地震剖面,初始模型,最终模型,基于模型反演过程的示意图,目前较为流行的模型反演方法：,模型叠代反演具有代表性的软件有：ISIS、SLIM等测井约束反演具有代表性的软件有：Strata、Jason、VanGurad,地震反演方法,模型迭代反演：,模型叠代反演的基本思想是在充分肯定原始地震资料的基础上，对初始油藏地质模型只做一个很粗略、随机的猜测（只是沿模型剖面在若干控

33、制点上用一系列不同深度、厚度相对较大的速度和密度层来定义）。在剖面的每一个道位置上，计算模型的反射系数序列，与一个估算子波或给定的子波褶积，便得到合成地震剖面。将合成地震剖面与实际地震剖面的对应道逐道进行对比，计算它们最小平方误差和（或者是）相似系数。最初的初始模型与地震道的差别较大，但经过反复迭代修改模型，并使用地震道作为约束条件，最终达到合成地震剖面与实际地震剖面每一道都有最佳拟合为止，此时的初始模型就是反演的最终结果。模型叠代反演的关键是迭代算法的优化，例如，ISIS软件中就采用了全局优化和局域优化相结合的非线性模拟退火方法，可使误差函数能量达到全局最小。,地震反演方法,地震反演方法,I

34、SIS反演处理流程,这种反演方法由于是做正演模型，因而不会受地震频宽的限制，即分辨率不受限制。本方法虽然使用了测井数据，但只是用于提取子波，而并没有作为约束条件参与反演的过程。从原理上来讲，此方法更加强调地震数据，而没有引入测井数据中的高频成分并体现在初始模型中，其结果也必然更加依赖地震数据和经优化的数学模型（反演目标函数）。虽然反演结果的分辨率较原始地震数据有所提高，但距分辨小于10米以下的薄储层的目标仍有一定的距离。而且，反演结果也不是唯一的。,地震反演方法,小结：,测井约束反演:,测井约束反演的基本思想与模型叠代反演大体相同，所不同的是：初始模型的建立是由使用地震解释层位的控制把测井声阻

35、抗曲线外推内插获得，初始模型较模型迭代法的要细致得多。使用井孔测井资料作为约束条件，而模型迭代法则是用粗略的地震道做约束条件。,地震反演方法,基本原理： 首先，使用地震解释层位的控制将测井声阻抗曲线沿层位横向逐道外推。外推的同时，按层序的厚度变化对声阻抗曲线进行拉伸和压缩。多井约束时，井间则按井的距离加权内插。这样，初始模型就保留了测井曲线的高分辨率。但井的外推和线性内插仍不能解决油藏的横向非均质性问题，因此就要用地震资料来修正井外推和线性内插的声阻抗曲线，实现非线性外推和内插。具体做法就是在给定的约束条件下不断地在初始模型与地震道之间求解一个反射系数校正值 ，并用此来修正初始模型的阻抗曲线，

36、制作合成地震道，使其尽可能地与地震道吻合。最终修正的模型就是反演得到的波阻抗。在这一过程中，克服多解性的约束条件则是指独立于地震数据之外的测井信息。,地震反演方法,地震反演方法,STRATA软件中测井约束反演的处理流程,反演约束条件（以Strata软件为例）,约束反演 约束反演使用初始猜测约束作为反演的起点，并用一个最大阻抗变化参数（初始猜测平均阻抗的百分比）作为限定反演计算的阻抗偏离初始猜测的硬边界。在反演计算中，阻抗参数可以自由的改变，但不能越过固定的边界。即标准正则方程满足 其中：L0样点处的初始猜测阻抗 LAV 输入约束L0的平均阻抗 当不存在约束或约束很宽时，由目标函数的最小平方解系

37、统，可以得到与地震道最佳拟合的期望输出。且其低频趋势由初始模型来实现而不是由数据解出。反之，最大阻抗变化参数减小时，约束变紧。而当其趋于零时，则引起期望输出无限地逼近初始模型。,地震反演方法,地震反演方法,测井约束反演地震数据,地震反演方法,初始猜测模型,地震反演方法,测井约束反演（硬约束）合成声阻抗,随机反演：随机反演把地震道和初始猜测阻抗视为两种数据，把初始猜测视为独立信息与地震道加权相加。虽然不象约束反演中设定一个硬边界来约束反演阻抗值的变化，但计算出一个随着计算阻抗偏离初始猜测而增加的补偿。由其目标函数 其中，w1为区间0，1的加权因子模型约束参数，表示我们对两种数据的相对信任程度。即

38、，随机反演同时平衡两种信息地震道和初始猜测阻抗。,地震反演方法,地震反演方法,测井约束反演（随机匹配）合成声阻抗,测井约束反演与模型迭代法相同，由于不受地震频带的限制，因而其结果具有较高的分辨率，也没有噪声的引入。由于使用了测井资料作为约束条件，提高了反演结果的唯一性。至少在钻井位置及其附近是如此。但这也就决定了此方法对测井资料的品质以及井震匹配相关性的好坏有着非常高的要求，反演对此非常敏感。因此，用于反演的测井数据一定要做环境校正、剔除野值；而且操作者不但要熟悉反演方法，还要了解目标区的地质资料，并认真作好井震标定。,地震反演方法,小结：,反演处理实例,资料背景输入地震资料：L9625 AV

39、O零炮检距的正极性剖面； 采样率：1 ms； CDP范围：230940。输入测井资料：夏盐3井声波、密度测井曲线； 深度范围：21003770 m目 标 处 理 层：盆地腹部夏盐北背斜侏罗系三工河组气层， 夏盐3井气层井下深度为 26862692 m。,反演,地震反演原理及其在油田开发中的应用,测井数据（AC、Density）,地震数据,地震地质层位标定（井的拉伸/压缩）,提取子波,子波是否零相位？,反褶积+带通滤波,使用提取子波的相位谱进行相位校正。,N,Y,带限反演,储层解释及成果图件输出,地震解释,建立初始模型(测井阻抗曲线的横向内插),测井约束反演,反演处理实例,反演处理流程,经反复标

40、定和相位校正处理后最终抽取的子波，并使用这个子波对反射剖面做子波反褶积。,反演处理实例,反演处理实例,统计性反褶积与子波反褶积的比较 a 子波反褶积 b 统计性反褶积,反演处理实例,反褶积后最终抽取的用于反演计算的子波。,反演处理实例,初始模型（阻抗曲线）,反演处理实例,测井约束反演 （最大允许阻抗变化=25%）,反演处理实例,测井约束反演（放大）,反演处理实例,约束反演的孔隙度剖面（放大）,反演处理实例,约束反演的合成剖面,反演处理实例,约束反演的误差剖面,反演处理实例,随机反演 （模型约束参数=0.5）,反演处理实例,随机反演（放大）,反演处理实例,随机反演的孔隙度剖面（放大）,反演处理实例,随机反演的合成剖面,反演处理实例,随机反演的误差剖面,反演处理实例,