矢量黏弹性衰减补偿高斯束偏移

基于弹性波动理论的多波多分量高斯束偏移具有计算效率高和成像准确等优点.但是目前此方法没有考虑实际地下介质的黏弹性对地震波传播的影响,从而无法补偿能量衰减和校正相位畸变,这使得该方法对一些含高黏弹性地层的成像效果不佳.针对衰减区域的成像问题,本文提出一种黏弹性衰减补偿高斯束偏移方法,该方法以多波多分量矢量波场弹性高斯束偏移方法为基础,在偏移过程中沿射线路径通过引入品质因子Q来考虑黏弹性影响并进行衰减补偿.该方法能够在偏移过程中实现PP波和PS波的自动分离及分别成像.同时,本文给出了在矢量波场偏移过程中提取角度域共成像点道集的方法,以便用于成像质量控制,并为后续速度和黏弹性参数反演提供所需的数据.本文利用2D层状模型和洼陷模型进行了方法测试,其成像结果验证了本文所提出的黏弹性衰减补偿高斯束偏移方法的可行性和有效性.

     

高斯束方法及其应用

高斯束方法是在动力学射线追踪和旁轴近似方法基础上发展而来的.高斯束是波动方程在特定射线附近的高频渐近解.本文介绍了高斯束方法及其与射线方法相比的优越性,以及高斯束在地震波场正演模拟以及偏移方面的应用,并提出了实现高斯束偏移的一种新方法.     

多波多分量高斯束叠前深度偏移

本文对基于弹性波动理论的多波多分量高斯束偏移进行了完整且详细的分析和公式推导,实现了3D空间多分量(矢量)波场的直接成像.由于当前多数基于弹性波动方程的偏移方法只是假设应力边界条件为自由地表边界条件,这种假设不符合垂直地震剖面(VSP)和海底电缆(OBC)等地震数据.为此本文详细分析了实际应用中常见的三种弹性各向同性介质模型的应力边界条件:自由空间、海底和自由地表模型.在上行传播假设情况下,获得了应力边界条件与位移边界条件的关系式.在此基础上,准确推导了3D多波多分量高斯束波场延拓和偏移成像公式,并在偏移过程中实现了完整的多波型自动分离.由于常规的互相关成像条件不适用于矢量波场成像,本文引用了散度/旋度互相关成像条件.通过约定PS转换波的正向旋转方向解决了3D空间PS成像极性翻转问题.利用2D和3D模型数据偏移成像验证了我们所提出的多波多分量高斯束偏移方法的可行性.     

三维最小二乘弹性高斯束叠前深度偏移

与声波高斯束成像相比,弹性波高斯束偏移更适用于复杂油气藏多波多分量地震勘探.但是由于观测系统的局限性和深部构造的复杂性,该方法同样存在成像分辨率低、照明不均衡等问题.本文结合最小二乘偏移和弹性高斯束偏移的优势,提出了一种通过弹性高斯束叠加构建Born正演（反偏移）算子和偏移算子的三维最小二乘叠前深度偏移方法.依据最小二乘反演理论,建立基于反偏移数据与实际观测数据残差的目标函数,采用共轭梯度算法迭代更新来建立地下真实的反射率.与传统弹性高斯束偏移方法相比,该方法不仅提高了成像分辨率,而且使复杂构造特别是陡倾角地层的成像照明也得到了补偿.理论模型测试结果证明了本文方法的可行性和有效性.

     

基于吸收衰减补偿的多分量高斯束逆时偏移

高斯束逆时偏移结合了射线类偏移的高计算效率和波动方程逆时偏移的高精度,能很好地处理焦散点、大倾角成像问题,并且具有面向目标成像的能力.多分量地震资料的偏移技术可以对地下复杂构造进行更准确的成像,由于实际地下介质具有黏滞性,研究黏弹性叠前逆时偏移具有一定的现实意义.本文采用高斯束逆时偏移方法对多分量地震数据进行吸收衰减补偿,首先分别给出纵波和转换波共炮域高斯束叠前逆时偏移方法原理,在此基础上推导补偿吸收衰减的表达式,校正Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的多分量高斯束叠前逆时偏移.数值模型的测试结果显示,在考虑地下介质的黏滞性时,本文方法具有更高的成像分辨率.     

基于黏声介质的角度域高斯束逆时偏移方法研究

高斯束逆时偏移是一种兼具计算效率和成像精度的深度域成像方法,能够面向目标成像.地下介质中黏滞性普遍存在,利用传统各向同性或完全弹性的成像方法处理黏滞性探区的数据会降低分辨率,并导致成像位置不准确和振幅欠估计等问题.本文在高斯束逆时偏移的基础上,通过对震源点和检波点处的波场进行衰减补偿,并结合高斯束求解时的角度信息,实现了黏声介质角度域高斯束逆时偏移方法.最后通过模型和实际资料试算对本文方法的正确性和适用性进行了验证.试算结果表明:相比于声波高斯束逆时偏移,本文方法能够对黏滞性引起的吸收衰减进行有效补偿,同时提取的角度域共成像点道集(ADCIGs)不仅可以用于分角度叠加成像压制成像噪声,而且能够为后续的偏移速度分析提供支撑.     

TI介质局部角度域高斯束叠前深度偏移成像

各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

基于吸收衰减补偿的多分量高斯束叠前深度偏移（英文）

高斯束叠前深度偏移是对地下介质进行精确成像的方法之一,多分量地震资料的叠前深度偏移技术可以对复杂构造进行更准确的成像。由于实际地下介质具有粘滞性,研究粘弹性叠前深度偏移具有一定的现实意义。本文采用高斯束偏移方法对多分量地震数据进行吸收衰减补偿,首先利用高斯束模拟粘弹性介质中传播的地震波,并在模拟中引入与品质因子Q有关的复速度和精确的粘弹性Zoeppritz方程,合成地震记录;然后给出纵波和转换波共炮域高斯束叠前深度偏移原理,在此基础上推导补偿吸收衰减的表达式,校正Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的多分量高斯束叠前深度偏移,并分析Q误差对偏移结果的影响。数值实验表明,在考虑地下介质的粘滞性时,本文方法具有更高的成像分辨率。     

多波高斯束叠前深度偏移速度分析(英文)

多分量地震资料叠前深度偏移技术可以对地下复杂地质构造进行更准确的成像,精确成像的前提是获取准确的纵横波偏移速度。本文采用高斯束偏移方法对多波地震数据进行偏移速度分析,首先分别给出纵波和转换波共偏移距域高斯束叠前深度偏移方法原理,在此基础上抽取纵波和转换波偏移距域共成像点道集;然后根据共成像点道集拉平准则,分别对纵波和横波速度进行更新;当两种波成像深度不一致时,对纵波和转换波成像剖面进行深度匹配,完成高精度的纵横波偏移速度分析。模型数据和实际资料试算表明,该方法是一种有效的多波偏移速度分析方法。     

面向盐下成像的基于波前构建法的高斯波束偏移

高斯波束偏移是目前工业上应用最广泛的偏移成像技术之一.但是当处理盐下成像问题时,由于盐体上界面起伏剧烈且盐体与围岩速度差异较大,只有小角度入射的射线才能穿过盐体,从而导致了盐下区域照明较差.通常,采用加密初始射线的方式来缓解这个问题,但这会导致计算时间的剧烈增加.由此,本文提出了基于波前构建法的高斯波束偏移,以保持波前传播过程中射线密度近似恒定.相较于加密初始射线的方式,本文提出的方法是自适应的,对于提高射线的覆盖范围更加有效.此外,本文采用高精度的波前初始化射线追踪技术作为波前构建的方法.Sigsbee2A模型下的应用实例表明,相较于常规高斯波束偏移加密初始射线的方式,本文提出的方法不仅能够获得更好的盐下照明效果,且对计算效率的影响也更小.     

非零井源距VSP逆高斯束叠加成像

VSP-CDP叠加成像是非零井源距VSP数据成像的重要技术手段,但是当地下构造比较复杂时,基于水平层状介质假设的叠加成像方法误差较大.本文使用基于高斯射线束的叠加成像方法对复杂构造VSP数据进行成像,根据高斯束有效邻域波场近似理论提出逆高斯叠加权函数,并使用该权函数对VSP数据进行VSP-CDP转换.同常规的VSP-CDP叠加成像方法相比,本文方法不仅解决了复杂构造非零井源距VSP数据成像问题,而且可有效加密反射点覆盖次数从而提高成像的横向分辨率.模型算例及实际资料验证了本文方法的有效性和稳定性.     

复杂地表条件下高斯波束叠前深度偏移(英文)

在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。     

True amplitude elastic Gaussian beam imaging of multicomponent walkaway vertical seismic profiling data

We develop the true‐amplitude prestack migration of multicomponent data based on the use of elastic Gaussian beams for walkaway vertical seismic profile (VSP) acquisition systems. It consists in a weighted summation of multishot data with specific weights, computed by tracing elastic Gaussian beams from each imaging point of the target area towards the sources and receivers. Each pair of beams may be connected with either a pair of P‐rays (PP‐image) or the P‐ray towards sources and the S‐ray to receivers (PS‐image) and is uniquely determined by dip (the angle of the bisector between the rays and the vertical direction) and opening (the angle between the rays) angles. Shooting from the bottom towards the acquisition system helps to avoid well‐known troubles, in particular multipathing for the imaging conditions in complex velocity models. The ability to fix the dip angle and implement summation over opening angles leads to the so‐called selective images that contain mostly interfaces with desired slopes. On the other hand, a set of images computed for a range of opening angles by summation over all available dip angles is used as input of an AVO‐like inversion procedure for the recovery of elastic parameters. The feasibility of this imaging procedure is verified by synthetic data for 2D realistic elastic models.     

弹性各向同性介质一次散射波场高斯束Born正演

本文提出了一种弹性波一次散射波场的正演方法——弹性波高斯束Born正演.该方法以线性散射理论为基础,通过Born近似建立起地下散射点处不同波型的反射率同弹性波主分量波场之间的数据映射关系,利用高斯束所包含的走时、振幅和极性信息进行不同波型的局部平面波的合成,进而通过逆倾斜叠加将所合成的局部平面波转化为时空域的多分量地震记录.该方法不但保持了射线类方法高效的优点,还具备了处理多次走时波场的能力,从而保证了复杂构造的波场模拟的精度.文中两个数值模型的应用效果表明,本文所提出的弹性波高斯束Born正演算法具有近似于波动方程有限差分法的波场模拟精度以及高得多的计算效率.

     

基于高斯束与高斯波包的Gabor框架散射波模拟方法

在给出真实模型和相应光滑背景模型的情况下,如何计算扰动模型(散射体)产生的散射波场是一个有实际意义的正演问题.在Gabor变换域描述散射体,且入射波场为短时宽带信号时,散射波场可以在频率域用高斯束或时间域用高斯波包描述.相对于波动方程方法,高斯束和高斯波包的计算效率更高;背景模型光滑时,高斯束和高斯波包方法的精度也接近波动方程方法.文中导出了声波假设下应用高斯束和高斯波包计算散射波的方法.测试分析了高斯波包的计算精度.给出了一般散射体的散射波模拟策略.同时针对一个理论模型完成了本文方法计算散射波的实验,实验结果表明高斯波包散射波计算方法是有效可行的.     

共炮检距高斯波束偏移的一种快速实现算法

与共炮高斯波束偏移相比,共炮检距高斯波束偏移具有直接抽取炮检距域共成像点道集的优势.过去,共炮检距高斯波束偏移以损失成像精度的代价采用最速下降法来降低积分的维数,从而提高计算效率.但经过最速下降近似简化的偏移公式仍是频率域的,需要在每个频点进行计算.为此,本文提出一种快速实现算法来避免采用最速下降法.本文通过分析一个水平层状速度模型的偏移过程和Marmousi速度模型的成像结果来检验不同插值方法对快速实现算法的成像精度和计算效率的影响,并建议采用二维三次卷积插值方法.同时本文在Marmousi速度模型下验证了快速实现算法相对于最速下降法在成像精度和计算效率上的优势.此外,本文将采用二维三次卷积插值的快速实现算法应用于Sigsbee2A模型并获得了清晰的盐下图像.