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

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

SSR与DSR组合的波动方程速度建模方法及应用

 波动方程叠前深度偏移是目前解决复杂构造成像最有效的方法，而建立复杂构造的速度模型又是叠前深度偏移的关键。本文结合单平方根（SSR）和双平方根（DSR）波动方程偏移方法，提出并实现了一套行之有效的且适应纵、横向变速的波动方程偏移速度建模方法，将并行批量处理与交互式可视化技术相配合，开发了交互式波动方程偏移速度建模软件。该方法的特点是用单平方根（SSR）相移法偏移成像，用双平方根（DSR）方法做横向变速处理，所建模型既可是时间坐标的层速度，又可是深度坐标的层速度。经对理论模型及实际资料处理，证明了方法的正确性和有效性。     

倾斜叠加三维叠前波动方程偏移

本文提出的倾斜叠加三维叠前波动方程偏移的理论和方法，是对不同的共中心点三维地震数据进行τ－ｐ变换，形成共中心点的三维τ－ｐ数据体；从不同的共中心点三维τ－ｐ数据体中抽选出共Ｐ值的三维数据体，并对各个共Ｐ值三维数据在Ｆ－Ｋ域中进行偏移延拓；然后，将偏移后不同Ｐ值的三维数据进行叠加，即得倾斜叠前三维偏移的三维数据体。通过点脉冲响应试验和点绕射偏移试验，说明了本理论方法的正确性和有效性。     

波动方程共炮检距道集叠前深度偏移

目前所用的几种波场外推算法（分步傅里叶、傅里叶有限差分、广义屏和频率--空间域有限差分）可以在共炮道集上获得较好的成像效果。在上述算法中,共炮道集为了兼顾各种波传播角,每炮偏移都要考虑相当多的边道,并且炮点和检波点要分别向下外推,因而计算量很大。针对上述问题,本文提出了基于双平方根算子的波动方程共炮检距道集的深度偏移方法,它适用于二维和三维叠前深度成像,其步骤为：①在高频假设条件下,把共炮检距道集波场延拓公式中的积分运算进行稳相近似,得到波场延拓的相移公式;②把速度场分裂为层内常速背景和变速扰动,可求得整个均匀层波场深度延拓的偏移时移量、各层的偏移时移量及振幅校正系数,进而可得到最终波场延拓值。脉冲响应测试和理论模型试算表明,该方法具有较高的计算效率和成像精度,可适用于复杂地质体成像。     

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

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

中国石化石油勘探开发研究院南京石油物探研究所 波动方程叠前偏移成像技术研究达到国际先进水平

2004年9月29日，由中国石化股份有限公司科技开发部组织的验收委员会和鉴定委员会，对南京石油物探研究所承担的中国石化科技攻关项目“复杂条件下深度成像方法技术研究”项目进行验收和技术成果鉴定。     

基于双平方根方程的双域保幅地震偏移方法

从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到保幅的单程波动方程;然后按照叠前深度偏移中的波场双重向下延拓理论实现共炮域公式到CMP域公式的转换,推导出由双平方根方程定义的保幅单程波动方程;基于反问题求解中常用的摄动理论,把速度场分裂为层内常速背景和变速扰动,分别在频率波数域和频率空间域求得波场深度延拓的偏移时移量及振幅校正系数,得到最终的双域保幅波场延拓算子。理论模型试算和实际资料处理表明,该方法不但具有较高成像精度,而且可以为后续的地震属性分析提供更真实的振幅信息。     

山前带地震数据的波动方程叠前深度偏移方法

作为复杂介质中最精确的地震波成像技术，波动方程叠前深度偏移在山前带复杂构造成像中大有可为，但该方法对速度误差非常敏感，抗噪能力较差，其成败取决于叠前去噪、表层速度结构反演和偏移速度分析的有效性，要使山前带地震勘探取得更大突破，就需改善单程波动方程深度偏移对高陡、倒转界面的成像精度，研究出可提高初至层析分辨率的新方法以及自地表开始的偏移速度建模方法。基于广角隐式有限差分单程波传播算子与波场“逐步累加”技术，研制了直接自起伏地表进行波场延拓与成像的叠前深度偏移算法。SEG山前逆掩推覆构造模型偏移试验表明，该算法具有很高的精度；川西龙门山实际宽线地震资料成像的试处理结果，也展示了该方法相对于叠后深度偏移的优势。研究还表明，面向地质目标优化地震波照明与接收效果也是需要重视的问题。     

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

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

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

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

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

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

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

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

基于GPU计算平台的三维波动方程叠前深度偏移

GPU/CPU协同并行计算是一项全新的、具有革命性的技术,将对地震资料处理行业产生深远的影响。本文在充分发挥GPU存储带宽宽、多寄存器和多处理器结构特点的基础上,通过改进地震成像并行计算模块的核心函数,形成一套波动方程地震成像GPU高效并行算法,并取得理想的加速比和可扩展性。所形成的一套GPU/CPU平台下波动方程叠前深度偏移软件,已投入地震资料处理实际应用。理论模型试算和实际资料的处理验证了本方法的有效性和实用性。     

三维转换波叠前方位各向异性校正技术

转换波成像与常规纵波成像不同的是转换波受方位各向异性的影响更为严重,相同反射点的快、慢波时差可达到40ms,严重影响转换波成像质量和分辨率。本文以横波分裂理论为基础,提出了一种叠前方位各向异性校正技术。该技术是基于共检波点叠加的层剥离转换波方位各向异性分析方法,并在叠前共检波点道集上进行方位各向异性校正。通过该技术在喇嘛甸工区的应用,得出T2层方位各向异性的主方位在67°～84°之间,快慢波时延在10～40ms之间。经过方位各向异性校正后,径向分量能量增强,切向分量能量减弱,提高了转换波成像质量和垂向分辨率。     

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

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

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

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

DMO和叠前时间偏移的共同起点

DMO和叠前时间偏移对于改善速度分析及叠加效果,最终提高剖面成像质量都是很有用的工具,一般在共偏移距道集中进行。其实现方法有多种,如等价偏移距方法（equivalent offset migration,Bancroft et al,1994)、DMO PSI(prestack imaging)方法（Gardner,1986）和非成像的炮集偏移（non-imaging shot-record migration,Berryhill,1996）等。这些方法的一个共同特点是用一个新定义的偏移距替代原来的偏移距。在每一种方法中,成像过程都归结为以新偏移距定义的共反射点双曲线型时－距关系上的动校正。本文中,我们证明了不同方法的DMO和叠前时间偏移有其共同起点－－双平方根算子,说明它们在原理上是等价的;最后我们提出了基于波动方程的DMO和叠前时间偏移方法。     

转换波叠前偏移技术新进展

转换波（C波）资料的叠前偏移是C波资料处理的重要环节，它能够解决纵波（P波）与C波最终剖面的波组对比、C波波场正确归位及提高C波成像质量等。介绍了C波叠前偏移的新技术。内容涉及到横波（S波）偏移速度分析、各向异性偏移、剩余偏移以及矢量（张量）偏移等领域，展示了C波叠前偏移技术目前发展的趋势。