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

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

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

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

基于束偏移的地层倾角约束网格层析建模

研究区深层构造复杂、缺少足够的地质信息及钻井资料、基于沿层层析无法建立准确速度模型。为了改善该区断块成像精度,进而落实构造圈闭,应用Tomogui叠前深度域网格层析速度建模和高斯束偏移成像。首先利用CVI(约束速度反演)方法获得平滑并符合地质特征的初始速度模型,进行高斯束体偏移;利用四维去噪技术优化处理成像道集,提高道集信噪比;在此基础上,直接进行基于道集数据驱动的网格层析反演,并利用地层倾角场约束网格层析反演过程。通过高斯束偏移与网格层析速度更新、迭代,采用深度域高斯束偏移成像明显改善了成像效果,查明了该区的断裂展布特征,断块型圈闭得到落实。     

基于三维高斯束算子解析的方位-反射角道集提取技术研究

通过三维叠前深度偏移提取地下成像点的共方位-反射角(或倾角-方位角)成像道集对于层析速度分析及振幅随角度变化分析(AVA)等至关重要。高斯束叠前深度偏移技术利用相互独立的高斯束描述波传播,并通过相邻高斯束的叠加计算波场,其实现方式解决了多路径问题,且无成像倾角限制。更重要的是,高斯束偏移提供了一种能高精度且方便高效地提取三维共方位-反射角度成像点道集的策略。通过解析求解高斯束算子,计算单个高斯束的波场传播方向矢量,进而通过不同高斯束的传播方向矢量计算得到成像点处的方位-反射角,从而实现三维高斯束叠前深度偏移方位-反射角成像道集的提取。数值计算及实际数据应用结果证明了三维高斯束叠前深度偏移共方位-反射角道集提取技术的有效性。     

弹性多波高斯束逆时偏移方法

近年来,我国部分老油区勘探开发难度日益增大,勘探目标逐渐转向裂缝油气藏和岩性油气藏等,研究弹性多波理论下的叠前深度偏移算法至关重要。根据高斯束偏移的高效性和逆时偏移的高精度,将弹性多波逆时偏移的实现思想应用于高斯束叠前深度偏移中,研究了一种弹性多波高斯束逆时偏移方法。在射线中心坐标系下,以二维弹性波KirchhoffHelmholtz积分为基础,根据多分量地震数据的矢量地震波场传播特征,利用弹性动力学高斯束表征的格林张量推导了弹性波波场正、反向延拓公式和相应的成像公式,并针对转换波的极性反转问题,提出了一种根据入射波入射角的正、负引入符号函数的校正方法。模型和实际资料试算结果验证了该方法的正确性和适用性。     

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

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

弹性波逆时偏移的角度域共成像点道集提取方法

对于多分量地震勘探,由于介质中纵、横波的耦合作用,多波角道集存在能量串扰和低波数噪声干扰等问题。针对横波Poynting矢量构建过程中存在的应力串扰问题,基于解耦延拓方程及伪横波应力构建方法,对质点振动速度和应力进行解耦,实现了无能量串扰的纵波和横波Poynting矢量构建;在基于Poynting矢量提取PP波反射角和PS波反射角的基础上,提出了基于弹性逆时偏移的PP波和PS波角道集提取方法;针对PP波角度域共成像点道集存在的低波数噪声干扰问题,提出了角度域共成像点道集的角度衰减叠加策略,提高了偏移成像精度。应用模型数据验证了方法的正确性。     

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

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

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

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

矢量Kirchhoff叠前深度偏移

由于传统的Kirchhoff积分法偏移基于标量声学波动方程,没有考虑地震波场的矢量特性,偏移成像的精度较低、保幅性较差。因此,本文从多分量地震波矢量场的特征出发,基于均匀各向同性弹性波波动方程,推导了Kirchhoff积分法矢量偏移公式;利用弹性波波动方程的矢量特性,在偏移过程中将泄漏的地震波信息还原到PP波与PS波分量中,实现了偏移过程中的波场分离。二维理论模型与实际数据的应用结果表明,基于矢量数据的深度偏移方法能够有效消除多分量地震数据的波型泄漏现象,提高偏移成像的精度。     

转换波叠前深度偏移

叠前深度偏移是转换波最佳成像方法。在研究Ｐ波叠前深度偏移的基础上，进一步研究了转换波叠前深度偏移方法，该方法涉及多波速度建模、震源函数求取以及延拓成像，其中偏移成像采用有较高精度的傅里叶有限差分（ＦＦＤ）方法实现。经理论模型及实际资料处理证明了方法及软件的正确性。     

VTI介质角度域叠前深度偏移

研究了VTI介质角度域偏移方法,以各向同性双平方根方程的角度域偏移方法为基础,从VTI介质qP波频散关系出发,推导出VTI介质角度域偏移的波场延拓算子;在频率-波数域处理以横向均匀速度传播的波场,在空间域处理具有速度扰动特征的波场以提高波场延拓精度。模型试算和实际资料处理结果表明：各向同性偏移方法由于未考虑各向异性参数的影响,绕射波不能完全收敛,波场聚焦效果差,降低了成像剖面的分辨率和信噪比,不能对地质构造精确成像;VTI介质角度域偏移可对断层、盐丘、小尺度地质体精确成像。对于角度域偏移产生的角度域共成像点道集（ADCIG）而言,各向同性偏移的ADCIG同相轴无法校平,残留断点绕射波,波场无法正确聚焦,不能正确反映局部地质特征;VTI介质偏移的ADCIG同相轴较平直,角度范围更宽,波场归位准确,精度较高。因此,利用VTI介质角度域偏移方法可对复杂构造精确成像。     

平面波射线追踪与面炮波场时间滞后深度成像

面炮(Arealshot)偏移技术是对实际炮集数据按照不同方向传播的平面波进行时间延迟分解并合成组合炮记录,分别进行组合震源下行波的向下波场外推和组合记录上行波的向下波场延拓,在相遇深度得到该方向平面波的反射波成像。本文基于面炮深度偏移技术,充分利用射线追踪计算量小和波动方程波场延拓的良好振幅特性,提出了用平面波射线追踪法计算平面波组合震源下行波前走时和面炮波场滞后时间成像的深度偏移方法,并着重讨论了程函方程有限差分法计算平面波传播时间和组合记录上行波场时间滞后的深度成像方法。通过对Marmousi模型数据的成像检验,表明该法要比积分法偏移效果优越。     

VTI介质角度叠加逆时偏移

常规各向同性逆时偏移方法在处理各向异性问题时会导致层位成像不准确、同相轴能量不聚焦等问题,并且利用拉普拉斯滤波算子等简单滤波方法不能彻底去除成像结果中的强振幅低频噪声。针对这些问题,首先应用声波近似方程模拟qP波在各向异性VTI介质中的传播,在震源激发位置设定平滑的各向同性或椭圆各向异性震源环以压制qSV波的影响;在此基础上提出了VTI介质逆时偏移的角度叠加实现策略,利用角度域共成像点道集中不同角度成像结果的叠加压制低频噪声,提高成像精度。模型试算和应用实例偏移结果表明,该方法能够对高陡复杂构造准确成像,且能量均衡、低频噪声压制效果较好。     

波动方程法井间地震成像

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

平面波最大能量叠前深度偏移

本文提出了一种新的基于最大能量震源波场的平面波叠前深度偏移方法。该方法保留了平面波偏移方法的优点。即在成像过程中,通过同时延拓频率域稀疏采样的平面波震源和最大能量旅行时,求取震源波场中的最大能量部分,从而显著地提高了平面波叠前深度偏移的效率。利用震源波场中的最大能量部分成像,还保证了较刘的成像精度。Ｍａｒｍｏｕｓｉ模型的蔗算结果表明,在复杂构造情况下,本文方法可获得同计算全波场的平面波叠前深度偏移     

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

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

叠前深度偏移速度建模方法分析

基于层位的层析反演和基于网格的层析反演是实际地震资料叠前深度偏移的两种主要速度建模方法。根据歧口凹陷典型地质模型及其正演模拟的地震波场,对两种速度建模方法进行了实验分析。实验结果表明：基于层位的层析反演方法具有较高的稳定性,但只能得到速度模型的低频分量;基于网格的层析反演方法可以获得速度场的高频分量,但受初始模型的影响较大,不容易收敛到实际速度模型。为此联合应用两种方法,即先利用基于层位的层析反演方法获得速度场的低频分量,再利用基于网格的层析反演方法获得速度场的高频分量,这样能够提高速度模型精度,改善叠前深度偏移的成像质量。     

基于角度域共成像点道集的层析速度反演方法研究

碳酸盐岩探区的复杂构造成像与速度反演是目前地震数据处理难点问题。碳酸盐岩构造复杂、目标层极深、断点发育,常规处理方法很难得到比较好的成像与速度反演结果。为此提出基于角度域共成像点道集(ADCIGs)、面向碳酸盐岩储层的层析速度分析方法,利用波动方程双平方根算子叠前深度偏移提取的角度域共成像点道集作为速度分析道集,能够更为准确地反映速度和深度的耦合关系,并且减少了假象等的干扰,使得速度分析结果精度更高。模型和实际结果表明,该方法能够较好地解决碳酸盐岩探区的速度反演问题。