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

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

波动方程叠前深度偏移PC集群并行计算技术

波动方程三维叠前深度偏移研究，是近年来国内外油气科研服务公司的研究焦点，运用PC集群式并行机结构的软硬件特点，进行波动方程叠前深度偏移，成为加速其庞大计算量的有效工具。本文基于地震波场成像基本算法的内在并行性，结合PC集群并行机的软硬件特点，提出2种算法处理流程对策，获得以较低的资金投入，而成功地实现了工业化生产规模的应用效果。并给出大庆油田兴城北地区60km^2三维地震资料处理的算例。结果表明，算法流程实用性能良好，具有较高的加速比，良好的扩展性及适应性。     

三维叠前波动方程共偏移距拟屏深度偏移

在中点偏移距坐标系中，横向双平方根(DSR)方程为使用屏传播算子研究三维叠前波动方程深度偏移提供了一种方便的框架。共偏移距拟屏深度偏移是偏移共偏移距、共方位角地震资     

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

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

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

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

单程波动方程叠前深度偏移并行算法

 根据单程波动方程波场外推原理，文中总结了共炮道集数据混合域波场外推与成像算法流程，该流程主要包括地震数据预处理、波动方程速度分析和波场外推成像3部分；鉴于单程波动方程波场外推算法和数据流具有非常适合粗粒度的任务级并行计算，处理器间的通讯较少，各个计算节点上可能有多个处理任务，且各节点的处理任务不等，负载极不平衡等特点，本文采用了主从模式动态负载平衡并行算法，该算法可以充分利用多个处理节点的空闲资源，提高作业的计算效率；将混合域单程波动方程波场外推算法与主从模式动态负载平衡并行算法相结合，实现了软件级断点恢复功能。应用国产“Stseis叠前深度成像处理软件系统”对胜利探区的实际资料进行了生产性应用处理，取得了较好的叠前成像结果。     

波动方程叠前深度偏移技术在苏北C工区的应用

目前解决复杂构造与速度横向变化剧烈地区地震偏移成像的最佳方法之一是波动方程叠前深度偏移技术。针对苏北地区地下构造复杂、断层多、断块小、地下速度横向变化大的特点，提出以偏移速度分析和建立正确的深度域速度模型为核心的适合于苏北C工区的叠前深度偏移处理方法。阐述了初始速度的建立、模型修正的过程及原则，通过分析苏北C工区地震资料叠前偏移处理结果，说明了波动方程叠前偏移技术的应用效果。     

单程波动方程叠前深度偏移并行计算与应用

利用单程波动方程炮道集叠前深度偏移方法,建立了相关叠前深度偏移处理流程;在研究地震资料海量数据管理特点和波场外推算子并行算法的基础上,提出了主从模式动态负载平衡并行算法;研究了大规模地震数据并行计算的断点恢复技术,实现了软件级的断点恢复策略。应用自主研究开发的地震叠前成像软件系统对胜利油田BS6探区160km2的三维实际资料进行的生产性处理取得了较好的效果。     

实用化的二维波动方程叠前深度偏移

实用化的波动方程叠前深度偏移技术正在进一步发展完善之中。本文仅比较了其中四种具有代表性的方法，即分步傅里叶法(SSF)、傅里叶有限差分法(FFD)、广义屏法(GSP)和空间一频率域有限差分法(XWFD)。这四种方法的基本思路是：将速度场分裂为背景场和扰动场；背景场的波场延拓采用相移法实现；扰动场的偏移成像采用不同的实现方式，从而构成各种方法的不同特点。文中对四种方法相应的脉冲响应进行了测试分析，同时用Marmousi模型数据和实测地震数据作了偏移成像处理，计算结果比较理想。还将波动方程偏移结果与克希霍夫积分法偏移结果进行了对比，分析表明波动方程偏移成像结果在分辨率和对弱信号的成像能力等方面明显优于克希霍夫积分法。     

共炮域波动方程叠前深度偏移技术在委内瑞拉探区的应用

委内瑞拉探区为盐丘、刺穿、逆冲构造发育的复杂构造区，共炮域波动方程叠前深度偏移技术被认为是复杂构造区地震偏移成像的一个重要而有效的工具，该技术主要分为两步：第一步利用单程波波动方程，分别从炮点和检波点对波场进行向下延拓；第二步为成像调整，即在每一个成像点上，对延拓得到的连续的检波点和炮点波场进行相关处理。共炮域波动方程叠前深度偏移技术的不足之处是高频成分损失较多，造成剖面的相干性较强。偏移背景较差。本文分析了共炮域波动方程叠前深度偏移技术的参数试验及存在的问题，认为强相干背景、CPU计算能力、对地震资料的振幅补偿精度要求更高等因素影响该技术的推广使用。文中认为可在偏移耗时可以接受的条件下，尽可能地选择较大输入频率范围，能够改善偏移背景。实例分析表明，在构造复杂地区。叠前。深度偏移成像效果要好于叠前时间偏移，共炮域波动方程叠前深度偏移成像效果要好于Kirchhoff叠前深度偏移。     

声波方程共炮记录叠前深度偏移

共炮记录叠前偏移对输入的速度非常敏感，要求偏移中所使用的延拓公式能较为精确地适应速度的纵横向变化。因此，本文在推导共炮记录偏移公式时，特别把相位移加有限差分混合法引入，使整个延拓过程既适应速度的纵、横向变化，同时也能够得到陡倾地层的精确归位，从而较好地实现共炮记录叠前偏移。正因为共炮记录叠前偏移对偏移速度非常敏感，故本文提出的方法又可用作速度分析。任何形式的叠前偏移都需大量机时，特别是对于迭代的叠前偏移速度分析，计算量更大，为此，本文针对混合法共炮记录偏移的特点，从两个方面讨论提高计算效率的方法。首先针对在频率、波数域中运算的相移法，引入F－K滤波，在波场延拓中，仅就有关的频率、波数进行运算，大大地减少了计算次数．同时又通过滤波提高了偏移的精度；然后针对整个混合法在频率域中的运算，讨论了减少频率循环次数来提高计算效率的方法；最后，作为一个特例，把共炮记录叠前偏移稍作简化，得到精确、高效的叠后混合偏移方法。总之，本文提出的基于共炮记录的算法不仅成为较为精确的、高效实用的波动方程叠前偏移方法，而且还可作为一种强有力的手段，用于偏移速度的分析之中。数值试算结果表明，本文提出的各种方法达到了预期目的。     

波动方程叠前深度偏移并行计算及其应用效果

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

直接下延法波动方程叠前深度偏移

对于地表和地下地质条件复杂地区的地震资料偏移处理来说，由于现有的偏移方法大都假设激发点和检波点在同一个水平面上，因此给偏移结果带来了一定的影响。为此，探讨了直接从起伏地表开始的波动方程叠前深度偏移方法。简述了逐步一累加法和频率一空间域有限差分叠前深度偏移方法的基本原理，在此基础上提出了复杂地表和地下地质条件下的直接下延波动方程叠前深度偏移方法，该方法不受起伏地表条件的限制，对模拟层速度的适应性强。模型试算和实际资料试处理表明，该方法直接从起伏地表开始向下偏移，将波动方程基准面校正和叠前深度偏移有机地结合起来，既能对复杂构造精确成像，又能适应任何起伏地表条件。     

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

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

波动方程共方位角三维叠前偏移方法

本文针对双平方根(DSR)单程波动方程全三维偏移方法在地震波成像实际应用中面临的计算量大,对覆盖次数很低的窄方位三维地震资料会产生很强的相干噪声,影响成像效果等难题,结合Crossline共炮检距偏移理论,首先推导出了基于双域传播算子的共方位角叠前深度偏移公式,然后通过坐标变换提出了沿双程垂直走时方向进行递归波场延拓的共方位角叠前偏移新方法。通过对SEG/EAGE盐丘模型合成数据偏移试验,证实了共方位角叠前偏移深度域和时间域成像方法的有效性。将该法和常用的Kirchhoff叠前偏移法同时用于实际地震资料偏移处理对比,展示出该法在成像精度与分辨率方面具有明显的优势。     

聚焦点控制照明波动方程叠前深度偏移初探

波动方程法叠前深度偏移能够较为精确地对地质体成像, 然而波场外推的数据量过于巨大, 影响了计算的效率。面炮组合和控制照明结合的叠前深度偏移通过对震源波场的控制以及高效的面炮偏移能够得到高质量的成像结果。进行控制照明叠前深度偏移的关键是合成算子的选择, 笔者基于有限差分的走时计算方法合成了逆时聚焦算子, 基于该合成算子利用傅里叶有限差分法进行了聚焦点控制照明叠前深度偏移, 针对M armou si 模型的试算结果表明了该方法的有效性及其在提高计算效率方面的巨大优势。     

三维合成炮叠前深度偏移

对国际上标准的SEG/EAEG模型进行了三维合成炮叠前偏移研究。与三维单炮叠前深度偏移不同 ,三维合成炮叠前偏移有一个叠前数据的合成过程 ,由于合成后数据量大大减少 ,所以合成炮偏移的计算量要比单炮偏移小得多。阐述了三维合成炮叠前偏移的实现过程 ,给出了相应公式。对SEG/EAEG盐丘模型的一个数据集进行了MPI并行计算 ,取得了良好的成像效果。所作的较多计算表明了三维合成炮叠前偏移的正确性和有效性 ,这为野外实际数据的三维叠前深度偏移提供了另一条有效途径。