反演技术原理

1、反演技术刖言反演的概念、目的反演的发展历史及趋势反演的基本方法四 地震反演难题的解决方案反演的实质反演的基本流程AVO反演处理简介地震、测井、钻井是石油工作者认识地下地质构造、地层、岩性、物 性、含油气性的最重要的信息来源。虽然测井、钻井仅能提供井孔附 近的有关信息，尤其是有关岩性、物性、含油气性的信息，但是这些 信息往往具有很高的分辨率，可信度、准确性，能确切地指出含油气 层的位置，定量化分析与储层、油藏有关的参数。然而一个油气田勘 探、开发方案的设计、实施、调整仅靠测井、钻井资料是远远不够的， 必须与地震资料相结合进行综合分析才能取得良好效果。地震资料的分辨率虽然远远不及测井、钻井，但是随

2、着地震勘 探技术的发展，从光电记录、模拟记录到数字记录，从二维到三维， 地震资料的信噪比、分辨率、成像的准确性都获得了极大的提高，由 于地震资料包含大量地下地质信息，覆盖面积广，具有三维特性，所 以这项技术的使用越来越受到石油工作者的重视，如何利用地震资料 研究地下地质构造、地层？如何进行储层预测、油藏描述？如何进行 油藏、含油气层的预测？这些问题促使地球物理学家、地质学家开发应用了一系列地震 资料特殊处理技术，如地震资料反演技术、地震属性分析技术、AVO 分析技术，这些技术充分利用测井、钻井、地震的长处，使人们对地 下储层、油藏的研究从点到面、从二维到三维、从三维可视化研究到 油藏动态监测、

3、从定性研究到定量化研究，大大提高了钻探成功率， 有效地指导了油田开发，为提高油田最终采收率起到了积极的作用， 因此地震技术被列为二十一世纪石油工业发展的首要技术，相信地震 资料特殊处理技术（地震资料反演技术、地震属性分析技术、AVO分 析技术）也必将在我国油田勘探、开发中起到越来越重要的作用。一.反演的概念、目的地震资料反演技术就是充分利用测井、钻井、地质资料提供的丰 富的构造、层位、岩性等信息，从常规的地震剖面推导出地下地层的 波阻抗、密度、速度、孔隙度、渗透率、沙泥岩百分比、压力等信息。 那么如何理解这个概念？还是让我们看看什么是正演吧！正演的概念如果我们已知地下的地质模型，它的地震响应如

4、何？通过模拟野 外地震采集，得到单炮记录，再通过速度分析、动校正、叠加、偏移 得到合成剖面这一过程就是正演。0J5600.005400.005200.005000.004800.004600.004400.004200.004000.00-5(-10地质模型mhpD-侣-20-25-3000 111111111111111110 5001000150C200C250O300O3500W0050C500C550O600O65007000750OB00OB50C000054005000460042003800340030002600220018001400射线追踪Shot Gather单炮记录CM

5、P道集记录速度分析Gornrnon Offset叠加剖面反演的概念了解了正演的概念，反演就好理解了。地震剖面代表了什么？ 地震剖面的同相轴实质上代表的是反射系数，同相轴追踪着反射系数 而不是砂岩地层，只有转换成波阻抗，才能真实地反应砂层的变化。地震资料反演可分为两部分：（1）.通过有井（绝对）、无井（相对）波阻抗反演得到波阻抗、速度 数据体。（2）利用测井、测试资料结合波阻抗、速度数据进行岩性反演，得到孔隙度、渗透率、砂泥百分比、压力等数据。反演就是由地震数据得到地质模型，进行储层、油藏研究。常规地震剖面构造、岩性解释特殊岩性体地震剖面上图是包括火山锥、火山溢流相的特殊岩性体地震剖面，从常规

6、剖面上地质工作者很难进行岩性体的解释和追踪。下图是常规剖面和 反演剖面的对比，从图中可以看出反演剖面更有利于地质体的解释、 识别、和横向追踪。常规剖面与反演剖面对比3反演的目的反演提供各种岩性剖面，目的就是将已知井点信息与地震资料 相结合，为油田工作者提供更多的地下地质信息，建立储层、油藏的 概念模型、静态模型、预测模型，提高油田采收率。下面几张图反映 了井点与反演结果的空间结构关系。体现了只有通过地震资料反演， 才使储层、油藏的空间定量化研究成为可能。钻井+反演储层、油藏研究利用反演结果进行油藏可视化研究二.反演发展历史及趋势反演的发展历史地震勘探技术经历了模拟技术、数字技术、多次覆盖技术、

7、三维 地震技术一次又一次地飞跃，在油气勘探中的作用越来越大。目前正 在向开发地震技术发展，在油气田开发中也将起到重要作用。特别是 “三高一准”（高信噪比、高保真度、高分辨率和准确成像）地震技术 的发展使它的精度大幅度提高，已经逐步从做简单的构造图，扩展到 进行地层、沉积、构造、储层物性、生油层评价和超压预测等研究， 进而到油气藏流体研究和油藏动态监测。如何利用地震资料进行油藏定量化研究？地震资料所包含的信 息（振幅、频率、相位、速度、吸收、反射结构等）与岩性、物性、 含油气性有何关系？这些问题一直是物探工作者努力探索、研究的焦点，主要可以分 为以下几个阶段：六十年代后期：用地震速度计算砂泥岩含

8、量。七十年代早期：亮点技术的应用。七十年代中期：声阻抗技术应用初期。主要方法 是反射系数求 和法、反射函数积分法和波阻抗递推法等。七十年代后期：纵横波勘探兴起，亮点、声阻抗技术广泛应用。八十年代早期：AVO分析技术开始应用，声阻抗技术中的稀疏脉冲 法、以模型为基础的方法也已出现并得到应用和迅速推广。九十年代：多波勘探、层析成像、井间地震、高分辨率地震采集、 处理技术等新技术不断涌现，地震反演技术也得到迅速发展，声阻抗、 岩性反演与地质统计学、模型正反演相互融合使反演多解性问题得到 较为妥善解决。以沉积模式控制建立三维地质模型多井约束构建参数模型020406080100sand none反演资料

9、对比分析测井、反演的发展趋势地震资料反演技术目前正由叠后到叠前、叠前、叠后相结合，由 单一的波阻抗反演到利用地质统计学、分形分维、神经网络等技术与 测井、测试、钻井、地质综合研究相结合，由单一的资料反演到正、反演相结合，储层建模、约束反演、油藏数值模拟相互验证，其目的 就是要通过多约束条件解决反演多解性，提供准确结果为油田的勘 探、开发服务。三.反演的基本方法当今的地震反演方法，一是建立在波动理论基础之上的。由其固 有的不适定性，反演的效果尚不如以褶积模型为基础的方法，加之这 种方法尚不够完善而未能得到普遍推广。二是建立在褶积模型基础上 的，即叠后资料反演,Brain.h.Russell将其分

10、为三类，即：带限反演，或称常规递推法。地震资料实际上是地震子波与反射系数的褶积。S（t）=R（t）*W（t）S（t）地震纪录R（t）反射系数W（t）子波反射系数：当波阻抗反差不大时，则有:p V -p V TOC o 1-5 h z R = 221-p V +p V2 21 1p V -p V = ApV2 21 1ApVR qpV因此反射系数取积分便近似有：所以反“出牝2 誉=2ln py射系数的积分正比于波阻抗的自然对数，这是一种最简单的波阻抗概念，我们通常称之为道积分（相对 波阻抗方法）。美国某家石油公司在墨西哥湾海上寻找非背斜油藏， 一直取得成功。其成功诀窍有三：亮点、积分地震道和AV

11、O技术。积分地震道的优点是：（1）递推时累计误差小；（2）计算简单，不需要反射系数的标定；（3）没有井的控制也能作。缺点是它不知道波阻抗的绝对值。常规递推法的绝对波阻抗反演通常分两步:（1）令时刻tn时波阻抗为Zn，时刻tn+St时的波阻抗为Zn+1递推公式为：显然，如果利用测井资料已知波阻抗初始值Z0,则可利用上式依 次推算出任意时刻的波阻抗。大家知道，在野外记录时，由于受接收 仪器的频率特性影响，在记录中基本不含低于6-8HZ的频率分量。此 时的波阻抗为相对波阻抗。（2）利用地震速度场或测井资料建模速度场通过滤波仅保留其低频 成分，把所得到的低频成分加到相对波阻抗曲线上，就可得到可用于 储

12、层定量计算和油藏描述的绝对波阻抗。稀疏脉冲法，包括最大似然反褶积、L1范数反褶积、最小熵反褶 积、最大熵反褶积、同态反褶积等，稀疏脉冲反演是基于脉冲反褶积 基础上的递推反演方法，其基本假设是地层的强反射系数是稀疏分布 的。从地震道中根据稀疏的原则提取反射系数,与子波褶积后生成合 成地震记录;利用合成地震记录与原始地震道残差的大小修改参与褶 积的反射系数个数,再作合成地震记录;如此迭代，最终得到一个能最 佳逼近原始地震道的反射系数序列。该方法适用于井数较少的地区， 其主要优点是能够获得宽频带的反射系数，较好地解决地震反演的多 解性问题，从而使反演结果更趋于真实。约束稀疏脉冲反演采用一个快速的趋势

13、约束脉冲反演算法，用解释层 位和井约束控制波阻抗的趋势和幅值范围，脉冲算法产生了宽带结果, 恢复了缺失的低频和高频成分；同时,再加入根据井的波阻抗的趋势 约束。约束稀疏脉冲反演最小误差函数是J=S (ri)p+入 q (di-si)q+a 2S (ti-Zi)2(1)式中：ri为样点的反射系数；zi为样点的波阻抗；di是原始地震道;si 是合成地震道；Zi介于井约束的最大和最小波阻抗之间；ti是用户提 供的波阻抗趋势；a为趋势最小匹配加权因子；p,q为L模因子；i是地 震道样点序号;入为数据不匹配加权因子。如果从最大似然反褶积中求反射系数r(t),则在上述过程中为了 得到可靠的反射系数估计值，

14、可以单独输入波阻抗信息作为约束条件, 从而求得最合理的波阻抗模型Z(t)=Z(t-1)(1+r(t)/(1-r(t)(2)稀疏脉冲法假设反射系数是稀疏的、离散的，利用测井资料可以 得到井旁道的准确反射系数，通过上述反褶积方法，在测井资料、地 质模型的约束下，逐道递推子波、反射系数，从而反演出波阻抗、速 度等数据。常规递推法与稀疏脉冲反演法主要是利用反褶积方法来恢复反 射系数序列，由经过标定的反射系数序列递推出相对波阻抗，然后加 上从声波测井和地质模型中得到的低频分量，最终得到反演波阻抗。 这两类方法的主要缺陷是选择可靠低频信息较为困难，由反射系数递 推波阻抗过程中误差积累快，当反射系数存在较大

15、误差时，递推出来 的波阻抗剖面会面貌全非。此外，经过反褶积处理的结果，并不代表真正的反射系数序列， 稀疏脉冲法在地质结构复杂的条件下使用效果很差。其精度也难以满 足储层预测、油藏描述的需要。2.以模型为基础的反演方法。如地震岩性模拟、广义线性反演、宽 带约束反演等。基于模型反演的基本原理是:假设一个N层地层模型,各层厚度、速度、密度参数分别为l(i)、v(i)、p (i)(i=1,2,3，,N),地震 波在各层垂直传播时间为t(i)=2d(i)/v(i),则第i层底部的反射时间 为t(i) =St(j) i=1,2,N (3)其地震褶积模型为M(i) = E r(j)wi-T (j)+1 i=

16、1,2，Ns(4)定义目标函数:o（m）= D- F +WIIM -Mpri+W|V M xMprj第一项：记录残差，要使反演结果的模型响应（F）尽可能逼近实际记录（D） 第二项：先验约束，即反演结果不能偏离先验值太远第三项：是要保证反演结果具有一定的横向连续性，使解更合理D:实际地震记录F:合成地震记录MI :波阻抗模型参数M pri :波阻抗模型参数的先验值V * :横向梯度WJ :波阻抗模型先验值的约束权系数WC :波阻抗横向连续性的约束权系数IIp : L这种反演方法采用最优化算法，迭代速度与稳定性都很好，克服 了波阻抗的相对和绝对标度在递推反演中的缺陷，改善了波阻抗界面 的分辨率，消

17、除了子波的剩余效应所造成的畸变，受地震资料带限性 质的影响小，提高了反演结果的可信度。以模型为基础的反演方法以测井资料为约束条件，采用正、反演结 合进行迭代，求取地下波阻抗将反演方法推向非线性问题。这种新方 法利用了测井资料的高频和低频信息，大幅度拓宽了地震信号的频 带，可以更好地获得薄层、薄互层的波阻抗信息。因而表现出强劲的 发展势头，是目前国内外各软件公司重点发展的技术，也是油田开发 阶段进行储层预测、油藏描述的主要应用技术。以模型为基础的反演主要分三个阶段进行：（1）首先应综合地震、测井和地质等资料得到的波阻抗曲线、层位 解释结果和岩性信息，确定一个初始波阻抗模型。这个初始模型把应 用地

18、质知识解释的层位、断层和岩性信息反馈到反演中去。（2）把地震道的估计结果与实际地震道相比，得到剩余误差值。利 用这个误差，通过随机算法（或模拟退火、神经网络、遗传算法等非 线性全局最优化方法），在噪声和模型协方差估计值得约束下，迭代 修改模型，直到获得一个可以接受的剩余误差为止。最终控制反演过 程的稳定性与分辨率，进而处理出高质量的波阻抗剖面。（3）利用测井岩性、物性、波阻抗反演结果采用拟和、地质统计学（克里金、协克里金）方法求取相关关系，进行岩性反演。目前国内外反演软件众多，多以国外软件为主。在我国较有影响 的软件有以下几家：（1）jason反演软件，该软件从软件设计上充分考虑了勘探、开发不

19、 同阶段的资料情况，有较完善的建模软件（SeiRIS、EarthModel、 VelMod），子波分析软件（Wavelets），反演方法主要有Invertrace（基于道反演技术，包括常规递推反演、约束稀疏脉冲反演、道合并 和孔隙度估算），InverMod（基于模型参数反演技术可以反演在时间 和深度域里与测井分辨率相同的模型），StatMod（地质统计学和随机模拟、随机反演技术可 以精细描述油藏非均质性特征）。（2）EPS的优势：全三维建模和反演，确保反演结果空间闭合。融合地质、测井、地震等多元地学信息到同一模型。采用遗传算法（GA）和模拟退火（SA）全局寻优的非线性反演方 法比较成功的解决了

20、地震反演所面临的数据和模型之间高度非线 性及目标函数具有多个极小值问题。提高反演的信息使用量、信息匹配精度和反演结果的置信度。（3）Strata反演软件，该软件使用简便但在模型控制方面不如jason 软件，反演方法主要有常规递推反演、宽带约束反演、稀疏脉冲反演。（4）RM反演软件，是Schlumberger （斯伦贝谢）GeoQuest软件的储 层预测及油藏描述软件，该软件主要以常规递推反演为主，今年可能 推出以模型为基础的反演。（5）ISIS反演软件，属于模型技术的反演方法，针对基于模型技术 的反演方法所存在的多解性问题，它利用波阻抗垂向变化、水平连续 性、信噪比、反演结果和低频模型的偏差，

21、用一个反演的概率空间来 限制多解性。该方法是正、反演相结合的方法，类似于西方地球物理 公司的SLIM （地震岩性反演），缺点是分辨率难以满足开发的需要，但在勘探阶段还是较好的方法。InterWell综合测井与地质约束三维多道地震地层岩性反演，该 软件在我国还未多见，方法上类似于jason的InverMod。西方地球物理公司的PAIT、PIVT、SLIM、LithoSeis等反演软件， 该软件这些年基本没多大发展,基本还处于带限反演、稀疏脉冲反演 阶段。其他的反演软件还有很多，如很多地震资料处理软件都含有反演 方面的程序，但相比前面所说的几家软件从使用、效果、方法以及商 业化程度上都无明显优势，

22、因此在生产中还应用较少。反演的方法、软件众多，正如前文所述最终结果必须能解决储层、 油藏研究的问题，从下面几张图可以看出，测井、地质、物探、钻井、 油藏工程人员必须紧密结合，合理、准确、适当的应用反演资料。A=cat.，_l R2E1:S-利用反演结果预测含油气性图JASON与STRATA反演软件对比表表5-1JASONSTRATA优点1、具有合理的单井地震地质模型建立方法和多种质量控制手段2、具有强大的初始模型建摸功能（较精确），对断层处理效果好3、在横向上具有较好的可靠性1、在纵向上具有较高的分辨率和信噪比2、和井的吻合程度较高3、识别薄砂层的能力较强4、反演处理速度较快缺点1、在井数较少

23、时纵向上的分辨率较低2、反演处理速度较慢1、标定过程中可随意修改原始测井曲线2、横向上分辨率的可靠程度不一致3、对断层处理效果差四、地震反演难题的解决方案计算量和数据量非常庞大“基于延迟脉冲模型的快速算法。数据和模型之间高度非线性“模拟退火非线性反演方法。目标函数具有多个极小值模拟退火全局最优的反演方法。具有多解性复杂地质建模技术，实现多元地学信息在模型空间的有机融 合，提高反演的信息使用量。五、反演的实质从地震道中消去子波的影响而达到反射系数序列，势必存在多解性问题，目 前的方法除了递推反演之外在某种程度上都是基于模型的约束反演，区别只在于 如何平衡非地震信息的约束作用和如何拓宽频率。任何反

24、演都存在多解性的困扰，对地球物理反演来说，增加合理的地质约束是多解性中 的最优解，是地质含义上的惟一解。六、反演的基本流程目前的常见地震资料反演处理流程主要可以分为两部分：1.前处理用于反演的地震资料要求是高保真的地震资料，通常需要在提供 的保幅叠偏地震数据基础上进行必要的提高分辨率、信噪比处理，流 程如下：叠后地震资料处理流程2.反演处理通过前处理，认真分析反演工区的地质特点、地震资料品质、 进行测井曲线的归一化处理之后，反演处理流程如下：原始地震成果数据叠后提高信噪比和分辨率保幅处理解释成果输入1合成记录标定1子波分析及提低频模型构造主要过井剖面预反演是全区波阻抗反演岩性模型岩性反演成果输

25、出常规反演流程五、反演的质量控制因素地震资料反演结果的精度从反演的方法上来看不难想象，不论 是带限反演、稀疏脉冲反演、基于模型的反演，还是目前还未广泛推 广的波动方程反演大致主要受以下几个方面的限制：（1）地震资料本身的限制，用于反演的地震资料必须是高保真 度资料，如果资料本身不是保持振幅处理的，那么在进行反演时，所得到的波阻抗就反映不了地下地层的信息，也就根本就不能用于储层 预测和油藏描述。此外，地震资料还应尽量高信噪比、高分辨率、准 确成像，这些都将影响反演的分辨率，一句话，地震资料的品质将最 终影响反演的分辨率。Viei.-j 1 Spectr-um Trace range: 170 to 280 Time range： 1200 to 2700 ms255075Frequency (Hz)100资料品质分析（2）用于反演的测井资料品质，在进行反演前必须对工区内参加反 演的测井曲线进行环境校正，进行全区曲线的归一化处理。如果测井 曲线没有进行环境校正，归一化处理，在反演时同样也反映不了地下 地层的地质信息。（3）子波的提取，以褶积模型为基础的反演方法需要子波来进行反 射系数的提取，如何估算子波？怎样提取的子波才更适合于反演？才 能得到准确的