炮域波动方程叠前深度偏移的偏移距域和角度域共成像点道集计算

炮域波动方程叠前深度偏移是一种实用化的成像算法,特别在稀疏炮集宽方位角观测系统采集的数据成像中非常有效。本文将偏移距定义为叠前偏移中下行震源波场和上行接收点波场之间的位移,在偏移过程中提取偏移距域共成像点道集,利用最早应用于炮—检偏移的径向记录道映射方法,将偏移距域道集转换到角度域。计算过程包括多偏移距成像和炮—检点坐标系向中心点坐标系的转换。该方法能够在炮集偏移过程中生成随角度变化的反射波成像结果,但在稀疏炮集观测系统下,共成像道集的质量仍会受到炮点假频的影响。     

波动方程叠前深度偏移角度道集求取方法

针对常规直接转换法和直射线叠前时间偏移(PSTM)法提取的角度道集无法满足地质构造复杂地区的AVA分析与反演精度的问题,研究基于共炮检距道集波动方程保幅叠前深度偏移方法,提取角度域共成像点道集。首先,对地震数据进行叠前深度偏移得到炮检距域共成像点道集;其次,在频率-波数域对炮检距域共成像点道集利用快速插值映射法抽取角度域共成像点道集。通过与直接转换法和直射线PSTM法提取的角度道集比较,分析不同方法求得的角度道集信息的准确性和AVA特征。数值试算和实际资料处理结果表明:共炮检距道集波动方程保幅叠前深度偏移方法求取的角度道集的振幅相对保真,角度范围更广,更适用于AVA分析。     

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

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

弹性波Kirchhoff叠前深度偏移速度分析

本文提出的弹性波偏移速度分析方法是基于Kirchhoff叠前深度偏移形成纵波和转换波炮检距域共成像点道集拉平准则,分别对纵波和横波偏移速度进行更新。当两分量地震数据成像深度不-致时,通过调整横波偏移速度进行深度匹配,完成高精度的弹性波场偏移速度分析。文中分别给出速度更新及深度匹配方法。模型数据和实际资料试算结果表明了该方法的可行性和有效性。     

波动方程偏移角度域共成像道集

对于复杂介质，炮域共成像道集和偏移距域共成像道集遇到了麻烦：对于相同到达旅行时和相同的水平慢度，反射同相轴可能来自地下界面不同的点，即地下界面的点存在多解性。对于地下界面的每一点，角度域共成像道集却惟一地定义了射线对，因此，数据体中的每一个同相轴只对应地下界面一个点。利用波动方程偏移能够生成可以用于偏移速度分析和振幅随入射角变化分析的角度域共成像道集。     

角度域共成像点道集的提取与叠加成像

基于共成像点道集的速度分析是速度分析方法中常用的一种方法。共成像点道集有偏移距域共成像点道集、炮域共成像点道集和角度域共成像点道集等。在复杂地区，偏移距域共成像点道集和炮域共成像点道集存在运动学和动力学假象，而角度域共成像点道集则有效地克服了这个问题。首先分析了角度域共成像点道集的特性，讨论了如何较准确地提取角度域共成像点道集；然后分别利用单平方根算子和双平方根算子波动方程偏移方法提取了偏移距域共成像点道集；接下来应用插值法将偏移距域共成像点道集投影到角度域，生成角度域共成像点道集；最后进行了角度域叠加成像。模型试算结果表明：与常规波动方程偏移方法相比，角度域叠加成像能够改善局部成像质量。     

基于起伏地表的合成平面波叠前深度偏移

实践表明,波场延拓叠前深度偏移比Kirchhoff积分法成像精度高,但计算量大;而平面波叠前深度偏移比炮域叠前深度偏移运算效率高,且成像精度相当,但只适应水平地表.为此,介绍了一种适应起伏地表的合成平面波叠前深度偏移方法,其基本思路是:以地震排列的最高点所在平面为波场延拓起始面,将起伏地表的地震排列观测数据(检波点或炮点)向下延拓到地表最低点所在的水平基准面,实现波动方程基准面校正;在此平面上应用p变换将全部炮点合成为平面波震源,从而使全部炮记录分解成平面波记录;运用下行波方程、上行波方程分别将平面波震源波场、平面波记录波场沿深度方向外推,在每个深度进行波场相关并累加,获得该深度的成像波场值,得到共分角度的平面波偏移道集;将所有不同共分角度的平面波偏移道集按坐标叠加,得到基于起伏地表的合成平面波叠前深度偏移成像结果.四川龙驹坝地质模型的理论试算及四川实际山地资料HNT12线的处理结果表明:基于起伏地表的合成平面波叠前深度偏移成像结果的质量与传统炮域叠前深度偏移的结果相当,但运算效率显著提高,且适应起伏地表.     

双平方根算子波场延拓道集速度分析

在中心点一半偏移距域，利用波动方程的叠前双平方根算子．实现炮点及检波点波场向下延拓。通过速度扫描，在波场延拓后的深度域共中心点道集和经分选后的共成像点道集上，分析速度变化对成像点道集的影响。利用分偏移距成像，修正速度场，以达到提高成像精度的目的。     

三维VSP数据的波动方程偏移成像

常规的三维VSP地震数据偏移成像大都使用基于射线理论的Kirchhoff积分方法在共炮点道集中进行，因此对地下速度横向变化大的复杂构造区适应性差、效率很低。基于三维VSP地面激发井中接收观测系统的特殊性，利用地震波场的互易性原理，提出了一种在共接收点道集中进行三维VSP地震数据波动方程偏移成像的方法，不仅极大地提高了三维VSP偏移成像的计算效率，而且还适用于速度横向变化大的复杂构造区。为了利用时间域偏移成像方法对速度模型的相对不敏感特性，在三维VSP偏移成像中，借助于波动方程叠前深度偏移成像算法的概念，还提出了一种可以满足时间域速度模型横向变化要求的波动方程叠前时间偏移成像方法，拓展了波动方程叠前时间偏移成像的应用范围。     

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

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

应用散射时距关系的三维三分量VSP成像

由于观测系统的特殊性,无法采用地面地震技术处理三维三分量（3D3C） VSP数据。现有的VSP成像技术中,可借鉴的高精度成像方法较少,缺乏相匹配的处理软件,成像速度分析多需要迭代,成像难度大。为此,提出一种适用于3D3C-VSP数据的无需迭代的速度建模与基于散射时距关系的成像方法。首先,建立了适用于3D3C-VSP观测系统的散射时距曲线方程,将地震道上各时刻的振幅分别搬至对应的反射面上,通过干涉叠加,抽取共成像点道集;然后,针对上、下行波射线路径不对称以及偏移速度不一致的情形,把共成像点道集中消除正常时差的上行波波场校正至具有双曲时差特性的虚拟炮检距道集中;最后,对虚拟炮检距道集进行精细速度分析、动校正和叠加,获得高精度的PP波与PS波成像剖面。该方法用于四川盆地A井3D3C-VSP地震数据的偏移,获得了优于常规VSP-CDP叠加的成像效果。     

基于地质层位约束的2D波动方程偏移速度分析方法

叠前深度偏移和偏移速度分析是复杂构造成像处理中两个相互依赖且不可分离的过程。一方面,叠前深度偏移是使地下复杂构造精确成像的惟一途径,而叠前深度偏移需要一精确的速度模型;另一方面,因为叠前深度偏移对速度非常敏感,偏移后资料的剩余曲率提供了指示偏移速度正确与否的信息,可以作为偏移速度分析的工具。利用2D波动方程叠前深度偏移抽取共成像点的偏移距道集,通过分析共成像点道集的剩余时差与偏移速度误差之间的关系,给出了基于地质层位约束的迭代剖面偏移剩余速度分析方法,并用一个合成数据实例验证了该方法的有效性。     

高精度深度域层析速度反演方法

深度域层速度场是叠前深度偏移最重要的参数,直接决定了偏移成像的精度。传统层析速度建模方法由三部分组成,首先利用Dix公式或基于层位约束的Dix公式生成初始层速度模型,然后利用炮检距域道集对初始速度模型进行层析迭代,最后利用网格层析方法局部修饰速度模型。本文改进了传统层析速度建模流程：首先利用炮域波动方程相干反演方法建立初始速度模型,再利用角度域共成像点道集进行层析迭代,最后利用基于层位约束的网格层析方法局部修饰速度模型。渤海M区块三维实际数据试算证明本文对层析速度建模流程的改进是有效的,改进后的新流程可以大幅提高层速度场的精度。     

三维共散射点道集叠前时间偏移速度分析

基于共散射点(CSP)道集的叠前时间偏移速度分析方法的理论基础是散射理论和Kirchhoff积分偏移方法。等价炮检距为研究叠前地震偏移和速度估计方法提供了新的思路。本文首先阐述了二维情况下单平方根旅行时方程的推导,以及共散射点道集的提取与构建,总结了共散射点道集的优势;在此基础上研究了三维情况下共散射点道集映射原理、构建方法及速度分析方法。实际三维资料的试算表明,共散射点道集对应的速度谱和叠前时间偏移剖面的精度高于常规CMP道集对应的速度谱和叠前时间偏移剖面。     

一种MZO的新算法及其计算机上实现

文章提出了一种实现零偏移距化 (MZO)的新算法。其方法原理简明 ,类同于倾斜时差校正(DMO)。新算法巧妙的应用波场外推的技术及波动方程偏移的有限差分方法。新算法是在共偏移道集上实现的。对于三维地震资料需应用同方位角共偏移距道集。由于此方法采用了快速算法 ,比同类的DMO方法快一倍 ,其处理效果相当     

三维Born近似波动方程炮域叠前深度偏移

基于射线理论的Kirchoff三维偏移方法，在解决复杂构造的成像时遇到了困难。基于波动方程的三维叠前深度偏移对横向剧烈变速具有较强的适应性，且成像精度较高。Born近似傅里叶偏移方法，是基于小扰动假设，引入参考慢度，可以适应横向变速。利用频率波数域与空间频率域混合域偏移算子，构造出炮域波动方程三维叠前深度偏移算法，形成相应的处理流程。对SEG/EAGE盐丘模型数据集C的宽方位角三维数据进行了测试，实现真正全三维波动方程叠前深度偏移处理；对胜利油田CB地区部分海上实际资料进行了试处理。