单程弹性波逆时偏移和相移偏移方法

本文从单分量位移弹性方程出发,统一论述了单程弹性波相移偏移和单程弹性波逆时偏移方法。并从传播特性和计算效率的角度,对单、双程弹性波逆时偏移方法作了比较,结果表明,单程方法虽可保证波场单向传播并具有较好的成像效果,但计算效率相对较低。     

单步波场延拓法单炮时间偏移

偏移算法可以分为波场和射线两种类型,其中叠前时间偏移主要使用射线法,而射线法难以达到最佳振幅响应特性,因此,波场法一直都是研究的目标。现有的波场法叠前时间偏移主要是共炮检距和共平面波两种F-K域实现方式。一般而言,单炮总是最适合的选择,它具有方程简单、操作方便和容易实现的特点。根据这个认识,本文通过相移加内插的方法实现单炮时间偏移,采用均方根速度进行单步波场延拓,从而取得良好的成像质量和振幅特性。     

快速相移插值偏移及在复杂波场成像中的应用

在海拉尔地震工区,由于构造复杂、断裂带多,常规的偏移方法很难完全实现复杂波场的偏移成像,现有的相移插值偏移方法也因计算效率较低、受速度因素的影响大而未得到广泛应用,因此对相移偏移方法进行了改进。经改进后的相移插值偏移方法提高了计算效率,能适应大倾角的偏移归位,该方法在处理高陡构造波场成像方面比45°有限差分偏移方法具有明显的优势。     

一种解耦的起伏地表弹性波高斯束偏移方法

随着勘探开发逐步深入, 成像已逐步从单纯利用纵波数据发展为弹性波数据。针对双重复杂构造探区, 在弹性波高斯束偏移的基本原理中, 引入高程及地表倾角信息, 推导了由高斯束位移矢量表示的解耦弹性波场反向延拓公式, 利用互相关成像条件且通过转换波极性校正得到起伏地表条件下的P-P波和P-S波高斯束偏移算子, 发展了一种解耦的起伏地表弹性波高斯束偏移方法, 此方法无须对地震数据进行纵横波分离。在实现方法的基础上, 通过典型起伏地表模型的弹性波现象分析、成像结果分析及成像角道集提取, 证实了解耦的起伏地表弹性波高斯束偏移成像方法的正确性、适用性及良好的成像效果。     

波动方程法井间地震成像

井间地震波场丰富，含有反射波、直达波和管波等多种波型。为实现井间地震成像，首先对井间地震资料进行波场分离，提取反射波信息；其次通过直达波层析成像估算偏移速度场；最后利用Ｇａｚｄａｇ相移方法对实际资料进行偏移处理，获取井间高分辨率成像剖面。     

三维VSP资料波动方程叠前深度偏移研究

已有的一些三维VSP资料处理方法，存在保幅性较差、对复杂构造不能正确归位、成像分辨率低、噪声强和计算效率低等问题，为此，提出了三维VSP资料波动方程叠前深度偏移成像方法，其实质是把地表数据的三维单程波动方程叠前深度偏移成像方法推广应用到三维VSP数据。对方法进行了阐述，其基本原理是将三维VSP数据按井中检波器的位置抽成共检波点道集，对共检波点道集进行三维叠前深度偏移。实现过程为将VSP波场分解为背景波场和扰动波场，背景波场采用相移法偏移，扰动波场采用频率空间域有限差分法偏移，然后对偏移后的上下波场进行相关处理，获得成像结果。数值模拟证明方法正确可行，实际应用表明方法实用。     

叠前多分量偏移

用位移势能外推法进行叠前弹性偏移是一种混合的、系统的、功能块的矢量波场偏移算法。根据位移→势能→势能外推→位移这一顺序,本文提出一种新的不均匀介质波场外推法,它相对比较精确,且不费机时。利用一种改进的适用于复杂构造的褶积声波正演模拟程序计算二种直接下行P波和S波的旅行时,这种P波和S波在弹性偏移中是必需的。本文提出一种基于时间一致原理的新的成像条件 把外推和分解后的PP波和SS波位移振幅与成像条件部面相乘就可得到二种图像(PP波和SS波)。成像条件剖面中各种波之间的串扰最小。与以前的方法相比,新提出的成像条件公式更加优越,这是因为它能用于纵横波联合震源产生的多分量地震数据的偏移。数值试验的结果是非常鼓舞人心的,它清楚地表明这种技术是可靠的。为了解决成像条件中射线的角度和极性问题,还需开展进一步的工作     

叠前多分量偏移

用位移势能外推法进行叠前弹性偏移量种混合的、系统的、功能块的矢量波场偏移算法。根据位移→势能→势能外推→位移定顺序，本文提出一种新的不均匀介质波场外推法，它相对比较精确，且不费机。利用一种改进的适用于复杂构造的褶积声波正演模拟程序计算二种直接下行Ｐ波和Ｓ波的旅行时，这种Ｐ波和Ｓ波在弹性偏移中是必需的。本文提出一种基于时间一致原理的新的成像条件。反外推和分解后的ＰＰ波和ＳＳＩ皮位移振幅与成像条件剖面     

地震偏移中的时移成像条件

通过对成像点处的源点波场与接收波场匹配分析可实现基于单散射近似的地震成像。互相关是一种波场匹配分析的常用方法。经Kirchhoff积分法偏移或波动方程偏移后,地表地震数据成为关于空间和时间变量的深度域波场。简单的成像条件即利用零时间延迟的源点波场与接收波场互相关提取成像值。互相关成像条件可以在时间域和空间域实现应用。基于延迟波场的互相关成像可用于成像准确性分析以及实现角度域成像。介绍了一种时移(时延)波场的互相关成像条件。由该成像条件得到的成像结果是关于源点波场与接收波场之间时移量的函数,可用于积分法偏移、波动方程偏移或逆时偏移的时移道集和反射角度道集成像。利用模型数据数值试验展示该方法的主要特点。     

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

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

各向异性介质弹性波高斯束偏移方法

弹性波高斯束偏移是一种兼顾成像精度与计算效率的多分量地震成像方法,但当前的研究多集中在各向同性介质声波、弹性波成像和各向异性介质声波成像,针对各向异性介质弹性波偏移的文献相对较少。为此,在前人的研究基础上,发展了各向异性介质弹性波射线追踪方程,以二维各向异性弹性波Kirchhoff-Helmholtz积分为基础,利用弹性动力学表征的格林张量推导了各向异性介质中弹性波波场正、反向延拓公式,提出了一种针对各向异性介质多分量地震数据的成像方法。通过在成像公式中引入权函数,很好地压制各向异性介质成像中不同波型引起的串扰。通过引入符号函数,解决了转换波成像过程中的极性反转问题。TTI介质洼陷模型和TTI介质多层构造模型试算结果表明：与各向同性介质弹性波高斯束偏移方法相比,各向异性介质弹性波高斯束偏移成像结果的信噪比较高、同相轴连续性较好、构造位置更准确;使用多分量地震记录,利用纵、横波波场信息成像,有效提高了成像分辨率;与传统的基于标量波场的成像方法相比,所提方法的适应性更好。     

基于波场分离的逆时偏移成像条件研究及在准噶尔盆地的应用

与传统深度偏移成像方法相比,基于双程波的逆时偏移方法针对地下陡倾角地区和复杂地下地质条件地区的成像有着无可比拟的优势。逆时偏移方法的实现通常有3个步骤:①在时间域实现震源波场的正向外推;②在时间域实现检波点波场的反向外推;③应用合理的成像条件构建成像结果。在实际应用过程中,通常是采用互相关成像条件来构建成像,然而,互相关成像条件构建的逆时偏移成像结果会产生强振幅的低频噪音,这些噪音会极大地降低资料的信噪比,容易导致浅层构造被噪音掩盖而无法识别,仅采用单纯的带通滤波并不能很好地滤除这种噪音。提供了一种新的基于波场分离理论的逆时偏移成像条件,此方法能有效地消除这些特定的噪音,其核心思路就是将震源及检波点波场分离成它们的单程波传播分量,然后仅在反射点处利用互相关的成像条件构建成像结果,从而实现逆时偏移成像。Marm-ousi模型数据验证了该方法的准确性;同时,采用文中提到的成像条件在准噶尔盆地百口泉地区进行了逆时偏移测试,测试结果充分证明新成像条件对野外资料的适应性。     

VTI介质中波场分离算子特征研究

垂直轴对称横向各向同性(VTI)介质弹性波逆时偏移成像在每一步波场外推过程中都需要对纵、横波进行分离,以获得各自独立的qP波、qSV波成像剖面,因此纵、横波的分离效果决定了最终成像剖面的质量。本文针对纵横波波场分离过程中最为关键的分离算子进行了研究,分析了VTI介质中基于交错网格有限差分法波场外推过程中qP波、qSV波分离算法及波场分离算子的形态特征和影响因素。结果表明,VTI介质波场分离算子的中间部分取决于算子的阶数,两侧部分则体现了VTI介质的各向异性特征。     

各向异性弹性波高阶有限差分法叠前逆时深度偏移

采用高阶交错网格有限差分法算子构建了各向异性介质多分量弹性波动方程的高精度正演模拟和叠前逆时深度偏移的数值离散方程,在正演过程中采用最大绝对振幅能量法实现各向异性介质多分量初至旅行时的计算,并以此作为多分量双程弹性波叠前逆时成像条件,同时给出计算所需的稳定性条件、多分量叠前逆时成像的基本原理、吸收边界条件等,采用凹陷模型合成了共炮点道集,并进行多分量双程各向异性弹性波叠前逆时成像研究,并与对正演炮集进行多分量各向同性叠前逆时成像结果进行对比。计算结果表明,考虑了各向异性弹性波场的矢量特性,能够更为准确地实现叠前多分量弹性波场的成像问题,使地质层位中的断层、断点等复杂目标成像更加清晰准确,并且偏移成像精度较高,因此开展各向异性叠前逆时成像可为当前高精度岩性地震勘探提供方法指导。     

地震记录的Q值正演模拟与反滤波

地震记录的正演模拟与反演方法是地震勘探中重要的处理和解释手段之一。目前普遍采用的正演方法都是已知地层速度、密度或波阻抗正演一定观测系统条件下的合成记录，并没有考虑地下介质非弹性吸收的影响。正演获得的合成记录具有较多的假设条件，与实际地震记录差距较大。本文提出的模拟方法是利用波场成像原理，在递推上行波时应用Ｆｕｔｔｅｒｍａｎ模型同时考虑介质的吸收作用，推导出频波域类似于相移正演算法的Ｑ值正演模拟合成     

各向异性介质中弹性波逆时偏移技术的研究现状与展望

基于矢量波动理论的多波地震技术在解决各向异性地层的精确勘探方面具有理论基础优势,纵、横波深度域成像是多波地震资料处理的关键之一,多分量联合逆时偏移是实现纵、横波深度域成像的重要方法之一。总结了横向各向同性(TI)介质中弹性波逆时偏移技术的研究现状,重点分析了弹性波延拓和纵、横波成像中各个环节的实现思路与存在的主要问题,在此基础上探讨了该领域未来的研究方向。TI介质中弹性波方程逆时偏移领域重点技术主要包括:三分量地震数据之间的频谱一致性处理技术、各向异性随机边界构建技术、各向异性逆时偏移的噪声压制技术、数据驱动的TI介质中的纵、横波解耦技术、更准确的纵、横波传播方向求取技术和横波三叉区的处理与成像技术等。     

单程波外推高陡倾角反射成像技术的应用研究

基于单程波外推方法的二次反射波成像技术基本原理是：在将地震波场进行炮点下行波和检波点上行波常规波场延拓的基础上,进一步对检波点波场进行下行波延拓,对延拓的3个波场通过上、下行波两两交叉相关,进行上行和下行波场成像,该技术具有类似于逆时偏移的功能。以一个垂直构造模型为例,简要分析了高陡构造的波场特征;基于一个二维和一个三维高陡构造模型数据,进行了偏移成像处理试验,结果表明,可以对倾角大于90°的高陡构造精确成像。将该技术应用于某山前带实际地震资料的偏移处理,获得了较好的效果,常规偏移剖面中成像模糊的层位边界在该偏移剖面上成像清晰。     

弹性波高斯束叠前深度偏移

本文提出一种新的多分量地震成像方法--弹性波高斯束偏移,同以往的弹性波偏移（Kirchhoff偏移,逆时偏移）方法相比,兼具较高的成像精度与计算效率。该方法通过局部倾斜叠加将多分量地震记录分解为不同波型、不同初始方向的局部平面波,并利用弹性动力学高斯束的走时、振幅以及极性信息进行矢量波场的延拓成像。本文以二维弹性波Kirchhoff-Helmholtz积分为基础,推导了基于高斯束表示的弹性波波场反向延拓公式,然后利用互相关成像条件求取反射波以及转换波的成像值。针对转换波成像剖面中的极性反转问题,本文根据反射界面处转换波的传播以及偏振特点,给出了一种校正方法。最后,通过数值模型试算,验证了该弹性波成像方法的正确性和有效性。     

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

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

波动方程叠前深度偏移成像软件系统的研制及应用

基于波动方程的叠前深度偏移成像技术可以将叠前深度偏移技术的应用领域从复杂构造成像扩大到复杂地质条件下的岩性地层成像。研制开发了波动方程叠前深度偏移成像软件系统，该系统包含工区管理、数据管理、偏移速度分析、构造建模、2D／3D地震速度建模、地震偏移成像、三维可视化、辅助计算等一系列功能模块；具有独特的速度建模、叠前深度偏移成像、高效并行计算、三维可视化和性能优化等特色技术。对该软件系统进行了SEG／EAGE3D盐丘模型数据测试，在所获得的成像剖面上盐丘边界和断层清晰。将该软件应用于实际地震资料的处理，河南油田泌阳凹陷的高陡构造、胜利油田的古潜山内幕都得到了很好的成像。