基于改进逆透视变换的智能车测距技术_油气地球物理技术进展（一）：偏移、反演及定量解释...

　 

欢迎广大石油界的朋友加入到石油物资大数据当中来，点击左下角蓝色字体阅读原文进行石油物资查询

油气地球物理技术进展(一)：高效地震采集技术

---

编者按：油气地球物理技术的新进展主要体现在高效地震采集、偏移反演技术、定量地震解释、物探理论研究及岩石物理实验、自动化智能化物探技术等方面。《石油物探技术信息》微信公众号将对近期部分油气地球物理技术进展进行梳理，并分四期进行连续报道。本期专题作为第二篇，将对偏移和反演技术、定量地震解释技术进展进行报道。

---

1.偏移和反演技术

地震资料处理中，更加注意流程规划、方法选择、参数选取时的地质、地球物理意义，为后续解释和重新采集处理过程提供可靠的资料，体现了采集、处理、解释一体化的趋势。地震数据处理技术研究应用的主要内容包括高密度、宽方位资料的插值与规则化、静校正、去噪、建模、偏移成像、反演、OVT域处理、针对性属性计算等，技术进展中最为重点的，集中于偏移技术和全波形反演技术。

偏移和反演是求解波动方程的不同途径。一般而言，偏移给出地下构造情况，反演能揭示地下物性分布，偏移与反演的融合是发展大方向。联合偏移反演是朝着这个方向的努力之一。

1.1 偏移技术

地震偏移成像是地震资料处理中的关键一步，偏移成像结果直接决定了构造解释、储层预测的精度。随着油气勘探开发程度的不断加深，地震勘探逐渐走向更加复杂的地下构造和复杂储层区域。复杂几何形态的盐体、不同喷发形式形成的火山岩、古潜山等构造体的存在使得地震波场十分复杂，建模和偏移难度加大，这些地质体的底边界、侧翼、下覆构造的成像更为困难。偏移技术的研究聚焦于提高解决复杂地质问题的能力(复杂地表、复杂目标体的成像能力、成像结果的精度和分辨率等)、偏移成像的计算效率。叠前深度偏移已经成为常规流程。逆时偏移、束偏移的研究不断深入，并取得很多成功应用。最小二乘偏移重新受到重视。

1.1.1 最小二乘偏移

标准的偏移方法(包括了逆时偏移)是基于正演模拟算子的共轭转置(而不是它的逆)，使得成像结果受假像、不均匀的照明的影响，特别是在复杂地区。最小二乘偏移(LSM)则是对正演模拟算子的逆的近似，因而能减轻上述问题。

LSM的优势有：提高成像分辨率；减少观测不规则、采样稀疏、孔径不足时的偏移假象，改善成像聚焦；补偿吸收衰减、几何扩散、各向异性等引起的振幅问题，提高成像保幅性。

声波、弹性波最小二乘成像都有新的成果。最小二乘全波场偏移技术(LS-FWM)，考虑一次波和高阶反射能量，显著增强了成像的照明度和分辨率，成像较少受串音干扰，并改善了波数成分，振幅更加均衡，取得较好的断层构造落实效果。使用基于双域相移的偏移和反偏移算法，将LSM和稀疏反演相结合，成像分辨率、频谱带宽都得到改善，还可以重新设计观测系统，降低每次迭代偏移和反偏移的震源数目，从而降低计算成本。最小二乘Q偏移(LSQKir)方法，结合了LSM和QPSDM的优点，能改善照明度、振幅保真度、成像分辨率、进行Q补偿，断层成像得到改进。单次迭代(非迭代)最小二乘深度偏移能改进偏移结果的保幅性，提高分辨率和照明补偿效果，也能够进行考虑衰减问题的偏移。单次迭代最小二乘深度偏移计算量小于目前常规的迭代最小二乘偏移方法(常规迭代LSM的计算量约为常规偏移的2N倍(N为迭代次数)。基于反射率的反射地震波传播方程(反射波动方程)的LSM方法能够得到地下反射层界面的物性参数，而常规LSM得到的是地下反射层的物性参数)。

LSM也存在一些不足。一是计算量太大。二是LSM偏移约束条件的选取困难，目前仍然没有统一的结论。

LSM未来发展方向是寻求效果与效率的平衡。主要有两条路径：利用高效求解Hessian矩阵逆方法代替多次迭代法；LSM偏移与高效偏移算子相结合，在实现高精度成像的同时兼顾效率实用性。而最小二乘偏移与其它偏移方法的结合，成为近年来最引人注目的进展，如最小二乘逆时偏移成像(LSRTM)。

LSRTM是近年发展比较快的一种成像方法。通常地震成像会受到不规则采样及照明不足的影响，大多都采用加权方法尽量去弥补采集脚印或成像失真。使用最小二乘可以克服采集假频、补偿照明不足。

基于线性化弹性波全波形反演与伪Hessian预处理的弹性波LSRTM方法，使用梯度最小二乘来解决优化问题，基于伪Hessian对角化的预处理使成像精度更高，人为因素更小，相比弹性波RTM，具有更好的振幅均衡性，而且能够去除P、S波阻抗扰动的影响。基于能量范数的弹性波LSRTM能够压制人为因素造成的背景散射，加速收敛，对任意非均匀介质的准确成像，并能改善模型边缘的振幅表现。

虽然LSRTM相比常规偏移具有较大的优势，但仍然面临着速度场的精度、震源子波的估计等问题。为此发展出了不依赖子波的LSRTM(SILSRTM)等算法。SILSRTM用观测波场卷积模拟波场的参考道、模拟波场卷积观测波场的参考道形成目标泛函，目标泛函的两项中同时存在观测子波和模拟子波，从而去除子波影响。

1.1.2 无需速度模型的偏移成像技术

对于速度模型的依赖，是偏移技术中一直未能摆脱的困境之一(其它如：漫长的时间代价、对计算和存储能力的无止境需求等)。偏移结果的好坏在很大程度上依赖于速度模型建立得恰当与否。而且，越是偏移成像能力强的方法，对输入速度模型的准确性要求就越高。从另外一个角度来看，地震波的传播速度是地震勘探中一个至关重要的参数，如果已知了传播速度，也就没必要再进行偏移了。速度与偏移的关系，陷入了“先有鸡还是先有蛋”式的怪圈。

解决上述问题，至少有三条途径可以同时努力：改善速度估计方法，建立尽量准确的初始模型；研究反演方法(比如全波形反演等)，提供接近实际的物性分布情况；寻求弱速度依赖性甚至无需速度模型的偏移技术。在第三条路径中，逆散射子级数(Inverse scatting subseries，ISS)法、小曲线(Curvelet)法是两项研究较多的无需速度模型的偏移技术。

ISS深度成像技术考虑了逆散射序列中的高阶项，突破了基于Born线性化近似的逆散射成像仅对小扰动速度模型有效的局限性，并且无需水平层状介质的假设，可以实现有横向速度变化情况下的叠前深度成像，已经经过实际资料的验证。ISS给无速度模型偏移成像带来了希望，此外，还在地表、层间多次波消除，非线性AVO，Q补偿，地震资料插值等领域得到应用并取得较好成果。

地震数据经Curvelet变换后，分解为一系列的系数。Curvelet系数含有方向信息，可以进行各向异性的偏移，其优点是只需要处理为数不多的小曲线系数，降低了偏移的计算量。

小曲线域成像方法具有如下特点：Curvelet分解提供了关于局部倾斜的内在信息；Curvelet系数含有方向信息，可以进行各向异性偏移；能处理多路径的情况；算法直接使用压缩后数据为数不多的Curvelet系数，降低了偏移的计算量；除了压缩过程，其它步骤的完成都非常快；成像较清晰，没有偏移“微笑”；成像过程中无需建立速度模型。除此之外，Curvelet变换在地震中的应用同样非常广泛。

1.1.3Marchenko成像

Marchenko成像是一种不含有多次反射引起的假象的成像技术，它是基于Marchenko基准面重建的(只用地面地震记录、背景速度模型，将地面地震震源和检波器延拓到地下任意深度)，重建出反射响应的基准面。重建出的资料只包含新基准面深度以下的反射响应，基准面以上的介质与反射记录无关。

墨西哥湾实际海洋资料的应用结果表明，Marchenko重建基准面、成像更适应于野外资料的有效成像，特别是它压制了层间多次波引起的成像假象。

Marchenko方法在波场分解、震源子波估计等方面也有应用。

1.2 全波形反演技术

全波形反演(FWI)利用地震资料的运动学和动力学特征，对提高成像精度和储层预测可靠性具有重要意义，获得广泛的研究，并在海上、陆上取得较好应用成果。全波形反演是一种利用迭代的线性化反演方法进行求解的强非线性反演，具有计算效率低与局部极值高风险性的特点。其成功应用主要依赖于长偏移距、低频信息丰富、采集系统规则、高信噪比的优质地震数据，精确的初始速度，准确的子波估计等关键信息。

1.2.1 反射波全波形反演

基于反射的全波形反演(RFWI)方法能用于补偿背景速度模型的长波长，并能提供使用潜波时难以得到的速度更新。RFWI用分解的全波形反演(FWI)的斜率高波数成分首先更新密度，然后使用低波数成分更新速度。RFWI能够改善微型盆地中难以确定的盐上的变化的页岩速度，改善盐内和盐下速度，减轻盐体解释的模糊性，成像道集连贯平滑，构造连续性更好。在宽方位地震数据中使用自动校正盐体几何形态的RFWI方法，不仅能够拾取盐上的浅层慢速页岩体、高速碳酸盐，并且能够改善基底盐反射的一致性和连续性。使用数据域微分相似系数最优化方法来建立运动学校正模型，特别是在不能被透射能量照明的深层区域，用可调节的目标函数代替数据域微分相似系数目标函数，能够改进RFWI流程的稳健性。使用扩展源时间域RFWI有利于避免周期跳跃问题，克服了局部最小值问题，放松了常规FWI受局部最小值限制在应用时所需的苛刻要求，扩展了解空间。RFWI还适用于存在强反射界面的地区，以及初始模型不准、原始数据缺失低频等情况。各向异性多参数全波形反演可获得信噪比高、叠加和道集品质好的成像结果。

1.2.2 参数化FWI

进行多参数FWI时，由于模型空间内参数众多，参数之间不得不进行折衷，从而影响反演结果。在正交各向异性介质中，弹性近似甚至声波近似所需要的参数更多，因此上述折衷的影响也被放大得更严重。可以将FWI问题转化为求最影响数据(即数据对其最敏感)的那些参数，这些参数的类型依模型的不同详细程度、尺度和所反映的地质特征而不同。

TariqAlkhalifah等的研究指出，对ORT介质，数据将主要对四个参数的散射特性比较敏感：x1方向的水平速度、ε(主要是x-z平面的近零偏移距处的散射)、εd(x1、x2方向的水平速度的平方的比)、δ3(水平面上的非椭圆率的参数)。采用这种参数组合，水平速度和ε的散射特性都与方位无关了，可以对这两个参数进行VTI反演求得初始值，用上述另外的两个参数来拟合资料的方位变化。

1.2.3  自动化FWI

基于机器学习引入FWI实现的自动时窗拾取FWI方法，降低了全波形反演中的循环次数。而FWI方法的性质，暗示着FWI方法的自动化是非常有希望实现的。

FWI已经在海上地震资料中实用化，陆上资料的应用主要还以浅、中层为主。FWI还需解决近地表问题，断层成像不准确等问题。需要继续研究的技术难点有：初始模型建立、子波估计、缺失低频和长偏移距数据的FWI、避免局部最小陷阱、快速全局寻优、提高计算效率等。

2.定量地震解释技术

地震解释技术由定性半定量发展为定量化解释。定量解释技术的进展主要是发展了新的反演方法，例如一维随机反演方法、储层物性反演、基于贝叶斯随机地震反演的薄储层预测新方法等；深度域反演方法从深度偏移数据体反演输出反射率和声阻抗；提高复杂环境下地震构造解释和定量地震解释的可靠性与一致性，有助于降低复杂储层环境下的不确定性，改进储层定量解释。

目前，叠前地震反演已经成为定量地震解释的常规流程。近两年在定量地震解释方面又涌现出新的方法研究和应用。(1)一维随机反演方法(ODiSI)用于估算储层物性和相，并能估算相关的不确定性。ODiSI已成功应用于多个地质条件相对简单的硅质碎屑岩地区，例如安哥拉近海，北海和尼罗河三角洲。反演得到的净毛比和孔隙度数据体已经成功的用于井位部署，同时，也可以更好地了解储层地质特征、直接输入到储层模型中。(2)随着技术的发展，叠前同步反演已扩展到同时方位反演，如：压力和裂缝方向的定量方位反演方法，基于HTI、VTI介质的Rüger方程定义了弹性参数变换，通过各向同性正、反演准确模拟地震数据的各向异性反射率。合成数据的可行性研究表明，基于HTI介质的各向同性弹性参数和Thomsen各向异性参数可以准确预测宽方位角地震反演实验结果。该技术已成功应用于多个实际数据体。(3)基于确定性模型的叠前地震反演。传统叠前地震反演的直接结果通常是纵波阻抗(PI)、横波阻抗(SI)和密度。新的基于确定性模型的叠前地震反演是一种储层物性反演方法，它基于确定性模型，将弹性参数(纵波阻抗(PI)、横波阻抗(SI)和密度)转化为储层物性(岩性、孔隙度和含水饱和度)，用于地震反演结果的定量解释。模型数据和实例研究的反演结果表明该反演方法是可靠的。(4)贝叶斯随机反演。利用了子波、调谐和频谱分析、测井和局部变化的各向异性(LVA)等信息作为约束，贝叶斯地震反演的过程用于薄储层预测，能够准确反演比传统意义上地震分辨率还薄的储层。实例表明该方法的弹性阻抗反演结果非常接近实际测井，并可给出定量的误差或不确定性估计，减小开发风险。

深度域地震处理解释技术也继续向前发展。斯伦贝谢公司推出提高复杂环境下地震构造解释和定量地震解释的可靠性和一致性的深度域反演技术，能够给出反射率图像和照明补偿后的声波阻抗，从而由深度偏移数据体得到更一致的、可靠的成像结果和岩性属性。斯伦贝谢深度域反演服务已在复杂地质环境下得到成功应用。上覆复杂地质构造，例如固结的砂岩侵入体，能够影响储层的成像质量进而引起定量解释的错误或不确定性，深度域声阻抗反演改进了北海渐新世砂岩侵入体的振幅保真度和分辨率，得到更可靠的岩石声学参数，提高了地震振幅的可靠性，进而改进储层定量解释。今后，地震解释技术逐渐向深度域地震解释和智能化解释方向发展。

据《石油物探科技动态》2019年第1期 《油气地球物理技术进展》

撰稿/陈伟 李霞 李红敬 杨瑞娟

黎娜 周小慧 杨江峰 汪忠德

免责声明：海贝能源公众号旨在分享和传播优秀文化思想，无任何商业目的，版权归原作者所有，如有侵权请及时告之删除！欢迎广大朋友投稿。

点击左下角，关注石油物资大数据↓↓↓