真振幅全倾角单程波方程偏移方法

常规单程波方程偏移方法不能保持声波波场的传播振幅,大倾角的波场振幅被削弱,同时,波场的传播受到90°倾角的限制。针对这些问题.提出了真振幅全倾角单程波方程偏移方法。该方法首先向下延拓激发波场和接收波场,然后再将回转波场向上传播,从而将波场的传播角度范围拓宽到0°～180°,可以对任意倾角的构造进行成像,并通过引入真振幅校正项加强了大倾角结构的成像。将常规单程波方程偏移方法和真振幅全倾角单程波方程偏移方法应用于墨西哥湾地区复杂盐丘模型和实际资料,结果显示,盐丘的陡倾角侧翼边界在常规单程波偏移中不能得到很好成像,经过全倾角单程波方程偏移处理后被清晰地显示出来。真振幅全倾角单程波方程偏移方法的计算量大约是常规方法的两倍。     

真振幅偏移方法综述

在传统偏移方法基础上发展起来的真振幅偏移在偏移后的ＡＶＡ／ＡＶＯ／ＡＶＰ分析及属性预测中起到了重 要作用。简要回顾了真振幅偏移的发展历程，主要介绍了在Ｋｉｒｃｈｈｏｆｆ型真振幅偏移中起关键作用的积分权函 数的几种不同形式，以及基于声波波动方程分裂的单程波真振幅偏移和粘弹性波动方程真振幅偏移的理论与方 法，并分析了它们实际应用中存在的一些问题及可能的解决途径。     

裂步傅里叶真振幅叠前深度偏移

 由于传统裂步傅里叶叠前深度偏移方法是一种基于单程波波动方程的偏移方法，只能反映波传播的运动学特征，并不能反映波传播的动力学特征，因此不能保持地震数据的真振幅特性。为此本文在传统裂步法的基础上引入了裂步傅里叶真振幅叠前深度偏移算法及其实现流程，即在实现传统裂步法偏移的两个步骤中分别引入不同的真振幅校正算子项对传统的偏移算子进行校正，实现真振幅成像。文中给出此法对Marmousi模型进行试算的结果及实际地震数据的偏移结果，表明该算法具有较好的振幅保持特性。     

SEG第七十四届年会论文题录（三）

九、偏移 1线束和子波 （1）基于照度分析的波动方程（2）波动方程偏移中角域和真振幅成像的采集孔径校正（3）可控孔径波径-波动方程偏移（4）应用具有偏移窗口的局部余弦基进行波束偏移（5）真振幅波束偏移     

考虑薄层效应的真振幅偏移

基于波动方程的深度偏移中主要的处理是波场反向外推。常用的波场反向外推匹配滤波器法忽略了与薄层有关的角相关振幅和相位畸变。这将产生弥散的图像和错误的AVA(振幅随入射角变化)效应。根据“能量互易性”建立发射波场和反射波场之     

Kirchhoff型真振幅偏移与反偏移

对Kirchhoff型真振幅偏移和反偏移理论中的基本概念、基本假设和基本公式进行了简要的回顾和综合评述，并探讨了在今后发展中所面临的问题。从研究问题的出发点来看，已发表的真振幅成像理论与波动方程没有直接的关系。在数学上，加权绕射叠加和等时线叠加的表达方式具有多样性，真振幅加权函数具有多解性,这为真振幅偏移和反偏移的实现提供了多个途径。如果速度模型过于复杂，已发表的真振幅加权函数将在波场的焦散区失效。因此，在今后的研究中要探讨焦散区的加权函数计算问题。此外，要找寻最佳偏移孔径的快速估算方法，建立相应的预处理体系，并将现有的标量理论推广到矢量波场和各向异性介质。     

波动方程叠前深度偏移方法综述

:对波动方程叠前深度偏移的基本概念、实现方法进行了分类和阐述。波动方程叠前深度偏移方法主要分为2类：一类是单平方根方程偏移,在偏移过程中,上、下行波分别向下延拓,并通过互相关成像条件来提取成像值；另一类是基于“沉降观测”概念的双平方根偏移,在偏移过程中,炮点和检波点同时向下外推,当两者重合时（零偏移距）,零时间的波场值就作为该空间点的成像值。对共炮真振幅偏移进行了阐述,并指出也可以在角度域道集实现该算法。理论模型的处理效果证明，波动方程叠前深度偏移成像技术是解决强横向变速情况下复杂构造成像的一种有效手段。     

基于单程波方程的合成平面波保幅偏移

基于单程波方程的保幅偏移方法不仅能够提供波场的运动学信息,还能够提供动力学信息,实现-定程度上的真振幅偏移。相对于传统偏移算子的单程波偏移方法,保幅偏移加入了振幅补偿项,降低了计算效率。为了提高单程波方程保幅偏移的计算效率,本文提出了基于单程波方程的合成平面波保幅偏移方法,该方法在成倍提高计算效率的同时,可以提供地下目标区域的方向性照明。应用平层模型,对比分析了传统偏移算子与保幅平面波偏移算子的成像效果,应用复杂的Marmousi模型验证了该方法对复杂构造的适应性和有效性。     

相对振幅保持的转换波叠前时间偏移方法研究

 振幅保持的转换波叠前偏移不但能使构造准确成像，而且还可为AVO分析提供更多可靠信息。本文以纵波的单程波方程为基础，给出了转换波的单程波方程。根据反褶积成像条件，利用Kirchhoff Helmholtz积分推导出转换波的共炮保幅叠前时间偏移公式。由于最后反演积分的结果是体现反射系数的信息，所以又导出了转换波乘积型成像条件的共角度道集的偏移公式，这样使偏移后道集的振幅峰值成为入射角的函数，同时也大大提高了计算的稳定性。理论模型试算证明，无论是水平层还是倾斜层，应用振幅保持的转换波叠前时间偏移都能使偏移后的振幅峰值比较真实地反映原始反射系数的大小。     

波动方程最优分裂步相移(OSP)偏移方法研究及应用

波动方程叠前深度偏移是复杂构造准确成像的最有效手段,它用可描述波在复杂介质中传播过程的算子来实现波场外推,更适用于复杂介质中波的成像,此外,它还具有保振幅性的优点。波动方程叠前深度偏移的核心问题就是波场外推算子,一般要求算子能适应陡倾角反射的成像和较剧烈的横向速度变化,同时还要有较高的计算效率。介绍了OSP算子的基本原理及特点,并与其他几种偏移方法进行了比较;同时用OSP方法研制的软件在大庆油田进行了实际应用,取得较好的效果,表明该方法具有良好的应用前景。     

基于优化系数的混合域Fourier有限差分叠前深度偏移

为了提高复杂高陡构造区成像精度,提出了基于系数优化的混合域Fourier有限差分叠前深度偏移方法,该方法利用padé近似的有理函数对波场外推算子进行展开,然后利用切比雪夫多项式优化展开式系数,推导得到新的波场外推算子,降低了与波动方程精确波场外推算子的相对误差,提高了对波场外推算子的逼近程度,且在保证计算效率的同时提高了高陡构造区地震偏移成像的精度。对比改进后的混合域Fourier有限差分偏移方法与常规的傅里叶有限差分偏移方法(FFD)对Marmousi模型偏移剖面的成像效果,系数优化后的混合域Fourier有限差分偏移方法具有更高的成像精度。     

基于改进的时移成像条件的保幅叠前深度偏移研究

为了实现保幅叠前深度偏移,在波动方程保幅偏移方程的基础上详细推导了傅里叶有限差分波动方程延拓算子。由于波动方程叠前深度偏移技术的应用受计算效率的制约,因而,提出了利用大延拓步长进行波场延拓,在延拓层位,利用相关成像条件进行成像。对延拓步长层间的层位,结合含有绕射聚焦项的时移映射函数,对延拓层位上未基于频率叠加的上、下波场的互相关值进行傅里叶逆变换,并结合零时刻成像原理求出延拓层间的成像值。脉冲响应测试和Marmousi模型试算表明该方法可以实现保幅叠前深度偏移,并可以大幅度提高偏移成像效率。     

几种地震波叠后深度偏移方法的比较

地震波深度偏移方法是解决复杂地质构造成像的有力手段。为了研究复杂构造和速度分布条件下地震资料的叠后精确成像方法,分别采用相移加内插法、分裂步傅里叶法、傅里叶有限差分法和逆时偏移法对复杂的Marmousi模型进行叠后深度域偏移试验,计算结果表明,逆时偏移法具有最佳的成像效果,使得复杂构造成像清晰,绕射波场等完全收敛,计算精度和计算效率较高,是一种高精度的叠后深度偏移方法。     

转换波叠前深度偏移

叠前深度偏移是转换波最佳成像方法。在研究Ｐ波叠前深度偏移的基础上，进一步研究了转换波叠前深度偏移方法，该方法涉及多波速度建模、震源函数求取以及延拓成像，其中偏移成像采用有较高精度的傅里叶有限差分（ＦＦＤ）方法实现。经理论模型及实际资料处理证明了方法及软件的正确性。     

混合法VSP共炮记录叠前深度偏移

VSP处理方法反演地下结构直接、精确，是寻找复杂构造的、隐蔽的油气藏的最重要的方法系列之一。文章提出的混合波场延拓法VSP共炮记录单程波动方程叠前深度偏移，就是一种VSP对复杂构造精确成像的方法。首先把波动方程分解为上行和下行单程波，把波动方程解析解法(傅氏变换法)和数值解法(有限差分法)相结合，提出VSP偏移的混合波场延拓法，以使波场延拓适应速度的纵横向变化；再结合VSP的特点，在波场外推过程中，当在某深度处继续下延波场时，把从上部延拓至此的波场加上此处埋置的检波器的记录结果作为此处的波场；借鉴地面地震叠前深度偏移原理，提出了可适应速度纵横向变化的、可对复杂构造精确成像的单程波动方程混合波场延拓VSP共炮记录叠前深度偏移方法。纵、横向变速介质的数值模拟结果显示，该方法精确、有效。     

二维Offset平面波有限差分法叠前时间偏移

在岩性油气藏勘探中,保振幅或保波形成像是进行储层特征分析和参数估计的基础.由于Kirchhoff积分叠前时间偏移、F-K域波动方程叠前时间偏移不适应弱横向变速介质情形,因此,基于偏移距(Offset)域平面波偏移思想,提出了Offset平面波有限差分法叠前时间偏移方法.讨论了CMP道集平面波分解原理,推导了Offset平面波方程及其有限差分解法,给出了角度道集生成方法.利用大庆油田的一个横向缓变的速度模型,对Offset平面波有限差分法叠前时间偏移方法的有效性进行了验证.结果表明,该方法适应缓横向变速介质情形,波场外推算子保幅,计算效率高,可以为AVP/AVA分析提供较好的基础数据.     

平面波射线追踪与面炮波场时间滞后深度成像

面炮(Arealshot)偏移技术是对实际炮集数据按照不同方向传播的平面波进行时间延迟分解并合成组合炮记录,分别进行组合震源下行波的向下波场外推和组合记录上行波的向下波场延拓,在相遇深度得到该方向平面波的反射波成像。本文基于面炮深度偏移技术,充分利用射线追踪计算量小和波动方程波场延拓的良好振幅特性,提出了用平面波射线追踪法计算平面波组合震源下行波前走时和面炮波场滞后时间成像的深度偏移方法,并着重讨论了程函方程有限差分法计算平面波传播时间和组合记录上行波场时间滞后的深度成像方法。通过对Marmousi模型数据的成像检验,表明该法要比积分法偏移效果优越。     

应用GPU的傅里叶有限差分逆时偏移

逆时偏移是针对复杂构造的主流叠前深度偏移成像技术,常规逆时偏移技术大多采用高阶有限差分和成像后滤波算法。由于实际偏移的计算网格较大造成空间频散现象和成像后的去低频噪声滤波会损害有效低频信息,因此成像频宽受到一定影响。为此,提出基于GPU的傅里叶有限差分逆时偏移方法。即采用傅里叶有限差分算法进行波场传播计算,波场外推中直接在时间—波数域做显式波场分解和互相关成像;同时采用基于GPU的算法,显著提高了傅里叶有限差分逆时偏移的计算速度。实际资料应用结果表明,该技术能实现对高、低频有效信息的充分保护和利用,从而提高了地震成像分辨率。     

相移插值法偏移中一种新的变延拓步长方法

在叙述相移插值法波动方程偏移原理、对该偏移过程中的波场外推公式及最终成像方程进行适当的改进变形后,重新标定了向下外推算子,使得成像过程可借助于 FFT 来完成,并给出了一种不同于时间域内插的新的变延拓步长的方法。本方法与传统的相移插值法偏移中时间域内插值法变延拓步长相比,向下外推步长可远远大干地震波视波长的四分之一,从而满意地解决了相位插值偏移算法效率低的问题。文中给出的理论模型计算与实际煤田及陡倾角地震数据偏移处理的结果表明:本方法是一种快速经济的,适合于陡倾角、复杂地质构造归位的 f-k 域偏移成像技术。