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

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

高精度屏算子地震偏移成像方法研究

从波场延拓的非稳态相移公式出发,基于反问题求解中常用的摄动理论,利用单平方根算子的渐进展开,推导出了波动方程广义屏叠前深度偏移算子方程的高阶形式;针对散射波场计算项对于横向变速介质的不稳定性问题,通过数学近似提出了有效提高稳定性的策略并应用到波场递归外推过程中,从而得到一种稳定的高精度屏算子地震偏移成像方法。数值试验和实际资料处理表明,该方法具有更高的精度,对宽角度信息成像更好,是对传统广义屏算子的有效补充和发展。     

波方程傅氏有限差分法叠前深度偏移

通过对几种叠前深度偏移算法进行对比，分析了各自的优缺点。着重阐述了近年来发展的波动方程傅氏有限差分算法的原理及实现过程，该方法将波场延拓算子分解成频率波数域和频率空间域的三个算子进行运算。结合了相移法和有限差分算法的优点，克服了两种算法的不足。与目前广泛应用的波动方程Kirchhoff积分法叠前深度偏移相比，具有成像精度高，保持地震波动力学特征等优点。应用FFD偏移成像原理，研制了波动方程傅氏有限差分法叠前深度偏移软件，在Marmousi模型上成功地进行了FFD叠前深度偏移处理，取得了理想的成像效果。     

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

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

带误差补偿的有限差分叠前深度偏移方法

基于波动方程有限差分法深度偏移对介质速度纵横向变化的强适应性,本文介绍了一种优化系数的单程波方程叠前深度偏移算子。它在限差分波场延拓计算的基础上,增加了针对算子误差的补偿校正,从而提高了算子的成像精度;同时也可以基本消除差分频散,提高成像剖面信噪比。该算子具有方程阶数低且能对陡倾角成像的特点,能适应速度场的任意变化。文中一切计算均在频率域进行,与时间域有限差分算法相比,具有计算效率高、成像方便的优点。脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏称进强横向变速情况下“三高”（高分辨率、高信噪比与高保真）数据处理的有效手段。     

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

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

基于波动方程的聚焦点控制照明叠前深度偏移

 基于波动方程的聚焦点控制照明叠前深度偏移技术借助于差分计算，把速度、密度等介质参数的影响体现在差分计算的矩阵方程中，能够自动适应速度场的任意变化，快速傅里叶变换的使用也加速了波场延拓的计算速度。因此此法兼具有限差分偏移方法和傅里叶偏移方法的优点，既可适应速度场的剧烈变化，又可保证对陡倾地层的成像效果，是目前针对复杂构造最有效的成像方法之一。对于单个聚焦点及其周围的成像步骤为：①采用矩形网格情况下绕射走时的有限差分计算方法生成合成算子；②应用合成算子来合成面炮震源和面波记录；③对合成的面炮震源和面波记录做傅里叶有限差分法波动方程叠前深度偏移，得到该聚焦点及其附近区域的成像结果。按照上述成像步骤，将震源波场和炮集记录依据相应的外推公式进行延拓，最终应用成像条件求取成像值。在地质目标处选取多个聚焦点，可以得到面向目标的控制照明偏移成像，在多个层位上选取多个聚焦点进行控制照明叠前深度偏移，可以得到整块区域的成像。通过对Marmousi模型的试算，取得了较好的效果。     

共偏移距道集波动方程叠前深度偏移的Green函数法

在基于波动方程的偏移速度分析中，共偏移距道集数据的波动方程叠前深度偏移是十分重要的，因为它能提供地质上可解释的偏移图像（如不同偏移距的共成像点道集）。本文根据共偏多距道集数据的具体物理特性，采用Green函数法实现共偏移距道集数据的波动方程叠前深度偏移。在Green函数算子的具体构造中使用了稳定的Rytov近似计算慢度横向扰动引起的散射波场。用Marmousi模型数据进行了试验，结果表明共偏移距道集波动方程叠前深度偏移Green函数法的偏移结果不仅与常规的利用双平方根算子的共偏移距道集波动方程叠前深度偏移方法的结果相当，而且还能为偏移速度分析提供了不同偏移距的共成像点道集，对野外各种观测系统的适应性也很强。但本文的Green函数法的计算量较常规共偏移距波动方程叠前深度偏移法有明显增加。     

分步傅立叶波动方程三维叠前深度偏移算法及应用

把分步傅立叶波场外推算子推广到三维叠前深度偏移成象,得出频率波数域与空间频率域混合域中运算的三维波场延拓算子,提出一种基于波动方程波场延拓的炮域三维叠前深度偏移算法.用SEG/EAGA盐丘模型标准数据试算,获得了良好的成象效果.对胜利探区埕北地区的海上实际资料进行了试处理,取得了良好的效果.     

խ��λ����˫ƽ��������ƫ��Ӧ���о�

双平方根（DSR）方程为波动方程偏移提供了一种新的理论框架。基于“沉降观测”概念的DSR方程叠前深度偏移已经成为了一类重要的地震波成像方法。与单程波方程偏移方法相比，其计算效率要高很多，便于输出角度域的共成像道集，有利于开展波动方程偏移速度分析。DSR方程全偏移方法对计算机条件要求太高，且对地震勘探中大量存在的窄方位三维地震数据不太适用。文章结合DSR方程叠前深度偏移的波场传播算子与成像条件，讨论了它在窄方位三维地震数据成像中的应用问题，对比分析了几种稳相近似DSR方程偏移方法的特点与适用条件。数值试例表明，基于DSR方程的窄（共）方位角叠前深度偏移方法对盐丘（或礁体）、古潜山等复杂构造具有很强的成像能力，值得在石油天然气勘探中推广应用。     

平面波偏移、分角度成像与AVA道集生成

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的可靠方法，但存在着计算量大、对观测系统适应性差等缺点。平面波偏移是利用波动方程实现精确叠前成像的另一类方法，其基本原理是：通过地表延迟放炮的方式生成平面波震源，利用下行波方程进行波场正向延拓得到下行波场；对地表采集的炮集记录，以组合延迟放炮的方式叠加，得到地表平面波记录，利用上行波方程进行波场反向延拓得到上行波场；二者互相关求和，实现平面波地下波场成像。分析表明：平面波成像技术的精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当，但计算速度要快得多，且易于并行计算。二维Marmousi模型数值计算表明，射线参数的范围和间隔是影响平面波成像质量的主要因素；不同角度入射的平面波对最终成像结果的贡献是不同的，据此可以有针对性地选择射线参数进行平面波成像。     

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

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

波动方程叠前深度偏移技术及其应用

简要介绍了波动方程叠前深度偏移技术的发展趋势,讨论了P道集波动方程速度分析技术,推导出了共方位角波动方程叠前深度偏移方法的计算公式。利用上述技术和方法,借助国产神威I型机,完成了埕北地区130km2三维地震资料的叠前深度偏移,实现了共方位角波动方程偏移算子的并行运算,并行效率达到了线性加速。偏移结果表明,此方法能有效改善地震数据的成象精度,适用于构造复杂地区。     

起伏地表波动方程叠前深度偏移技术——以川东复杂地区应用为例

在我国西部或南方地区的山地、山前带地震勘探中,起伏地表是影响复杂构造成像的重要因素。为此,开展了直接从起伏地表进行波场延拓与成像的波动方程叠前深度偏移方法研究。包括：基于Reshef提出的“逐步-累加延拓”思想,采用带误差补偿的频率-空间域有限差分单程波延拓算子实现单炮记录的波场延拓;基于时间一致性成像原理提取地下反射界面的成像值;利用层析反演技术和剥层技术建立起伏地表深度域层速度场。理论模型和实际资料的偏移结果表明,该方法能够对复杂近地表和复杂地下地质构造很好地成像。     

共方位角波动方程偏移速度分析方法及其应用

叠前深度偏移是复杂地质构造成象的有力工具,其成象的质量主要取决于所用速度—深度模型的精度。本文根据波场外推理论,给出了层剥离波动方程偏移速度分析方法。该方法利用沿层剩余速度分析直接修改层速度,适用于3D偏移速度分析,已在复杂潜山构造成象中取得较好的效果。     

复杂地质构造叠前地震成像方法在东部M区块的应用

复杂地质构造成像面临的主要难题是地表及地下地质条件复杂,速度纵横向变化剧烈。叠前深度偏移能够使复杂构造准确成像,通过对比Kirchhoff积分法和波动方程法叠前深度偏移优缺点,给出叠前深度偏移实现过程中的数据准备、时间域模型建立、速度模型建立及叠前深度偏移方法选择。利用Kirchoff积分法能够对目标线进行多次迭代优化速度的优势,获得准确的速度—深度模型,用2种偏移方法对研究区的不同资料进行偏移。应用实例表明,2种方法对复杂构造都能够准确成像,但对信噪比的依赖程度不同,在信噪比相对较高的区域波动方程具有明显的优势,而在信噪比极低的区域,Kirchoff积分法更具有优势。经分析得出不同偏移方法对采用速度模型的精度要求不同,在具有不同信噪比的地区得到的成像结果表现为不同的效果。     

双平方根叠前深度偏移的广义高阶屏方法

叠前深度偏移技术既可在炮点-接收点域实现,也可在共中心点一炮检距域实现。前人将裂步延拓算子推广到共中心点一炮检距域相移法双平方根叠前偏移中。在波场向下延拓的每一步长内,仅通过一次时移量来校正常速相移延拓产生的误差,得出了简单、高效的裂步双平方根叠前深度偏移方法,但精度较低。本文基于波场延拓的非稳态相移公式,通过引入参考速度,并对双平方根项中的两个平方根项作泰勒级数展开,经过适当的数学推导,得出了共中心点-炮检距域波场延拓的双平方根非稳态相移新的高阶屏近似公式。该公式可直接在共中心点-炮检距域高效地延拓叠前波场。通过增加屏的阶数,提高了剧烈横向变速条件下叠前深度偏移的精度。它是裂步双平方根叠前深度偏移方法的推广。理论模型和实际资料的试算结果表明,本文方法是有效而实用的。