探地雷达衰减补偿逆时偏移成像方法

探地雷达信号在地下介质中传播时易受到电导率所产生的衰减影响,从而使得传统偏移成像结果在高衰减区域变得模糊.本文提出了衰减补偿的逆时偏移方法来消除电导率的影响.该方法基于麦克斯韦方程组实现电磁波的正演模拟和逆时传播.通过改变衰减项的正负号,保证了逆时传播的时间对称性,从而能够重构出原始波场,实现衰减补偿.数值实验比较了传统逆时偏移方法和衰减补偿逆时偏移方法在存在高导异常区域情况下的成像效果,结果证明了衰减补偿逆时偏移方法能够很好地恢复由电导率造成的信号衰减,从而提高探地雷达剖面的分辨率.     

基于干涉成像条件的探地雷达衰减补偿逆时偏移成像

探地雷达(GPR)衰减补偿逆时偏移方法采用互相关成像条件对实际地下结构成像时,由介质小尺度不均匀所产生的散射波和噪声能量在电磁波衰减补偿逆时外推中会急剧增大,严重影响成像分辨率.为此,本文提出一种基于干涉成像条件的GPR衰减补偿逆时偏移方法.该方法通过在逆时外推过程中,人为改变电磁波动方程中衰减项前的正负号,以保持逆时外推的时间对称性,补偿衰减的电磁波能量;干涉成像条件利用Wigner分布函数(WDF)对所有时刻的逆时外推电磁波场进行滤波,以有效压制逆时偏移剖面中的低频噪声及成像伪影.数值试验结果表明：相比于互相关成像条件,干涉成像条件可有效压制衰减补偿逆时偏移过程中由于背景小尺度不均匀性所产生的散射波和低频噪声,且几乎不影响目标体的成像能量;在计算效率相当的同时,成像分辨率更高.     

一种基于逆时偏移的高效照明补偿成像方法

对地震波偏移成像而言,由于观测采集系统不规则、不完备采样以及地下复杂构造对波场传播的影响,在成像点处通常存在照明不均匀的现象,对偏移成像结果进行照明补偿是提高地质体成像精度的有效方法.本文给出了一种基于逆时偏移实现的高效照明补偿成像方法.在高频近似下,本文首先定义了衡量地下照明均匀程度的雅克比矩阵,并以此为基础给出照明...     

基于电磁波衰减补偿的探地雷达三维逆时偏移成像

常规探地雷达（GPR）逆时偏移大都未考虑高频电磁波在地下高电导率介质中传播的强衰减特性,致使高衰减区域的成像质量低;实际地下结构呈三维空间分布,二维GPR逆时偏移难以将反射波准确归位和绕射波完全收敛于真实位置.为此,本文构建了一种基于电磁波衰减补偿的三维GPR逆时偏移算法.该算法采用三维时域有限差分法计算正传和反传电磁波场,并通过改变反传电磁波方程中包含电导率的衰减项的正负号,人为保持反传时的时间对称性和反转不变性,以补偿正传时衰减的能量;零时刻成像条件用于获取地下三维结构的成像结果.数值试验结果表明：相比于常规GPR逆时偏移算法,基于电磁波衰减补偿的三维GPR逆时偏移可精确补偿电磁波在地下高电导率介质中正传衰减的能量,高电导率区域的成像精度更高,分辨率更好,抗干扰能力更强,其结果更有利于指导后续雷达剖面的解译.

     

一种基于正则化策略的稳定衰减补偿逆时偏移

地下介质通常具有黏滞性,地震波在地下介质的传播过程中将不可避免地伴随着振幅衰减和速度频散,进而影响地震成像的准确性和分辨率.基于黏滞介质的衰减补偿逆时偏移能沿地震波的传播路径恢复其所经历的振幅衰减和相位畸变,有效提升成像效果.然而,由于地层对地震波的吸收衰减效应呈指数变化,衰减补偿过程中高低频分量的非同步增长易导致补偿算法数值不稳定.为此,本文提出了一种基于正则化策略的稳定衰减补偿逆时偏移方法.该方法基于解耦的常Q分数阶拉普拉斯算子黏滞声波方程描述地震波在地下介质中的传播效应,将振幅正则化因子引入该方程的时间-波数域的解析解中,以确保在补偿过程中地震波场能稳定延拓.二维、三维合成数据以及实际资料的偏移算例均证实了该方法的可行性和有效性,所提出的方法能有效地处理衰减补偿中的不稳定问题,明显提升地震资料的偏移成像质量.

     

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

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

添加正则化项的黏声逆时偏移成像方法研究

在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。     

基于时空域自适应高阶有限差分的声波叠前逆时偏移

叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率.     

衰减补偿的黏弹性最小二乘逆时偏移

地震勘探是发现和认识油气藏的重要方法,随着国内油气勘探领域、勘探对象的日趋复杂化,油气勘探迫切需求地震成像的精准度、分辨率和保真度的大幅度提升.相比于常规逆时偏移方法,最小二乘逆时偏移成像有较少的伪像,成像振幅均衡,成像分辨率高.目前常用弹性最小二乘逆时偏移方法来对地下介质进行成像.但由于地下介质通常表现出黏弹性性质,...     

基于起伏地形平化策略的弹性波逆时偏移成像方法

随着能源和资源勘查开采工作的深入,地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一,能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓,由于传统的地震波场延拓技术往往基于水平地表条件,相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法：基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法,采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造,获得了相应起伏地形条件下的复杂几何模型,开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性,可适应不同类型起伏地表采集的地震资料,显示出该方法在地震勘探领域的良好应用前景.     

Strategies for stable attenuation compensation in reverse‐time migration

Attenuation in seismic wave propagation is a common cause for poor illumination of subsurface structures. Attempts to compensate for amplitude loss in seismic images by amplifying the wavefield may boost high‐frequency components, such as noise, and create undesirable imaging artefacts. In this paper, rather than amplifying the wavefield directly, we develop a stable compensation operator using stable division. The operator relies on a constant‐Q wave equation with decoupled fractional Laplacians and compensates for the full attenuation phenomena by performing wave extrapolation twice. This leads to two new imaging conditions to compensate for attenuation in reverse‐time migration. A time‐dependent imaging condition is derived by applying Q‐compensation in the frequency domain, whereas a time‐independent imaging condition is formed in the image space by calculating image normalisation weights. We demonstrate the feasibility and robustness of the proposed methods using three synthetic examples. We found that the proposed methods are capable of properly compensating for attenuation without amplifying high‐frequency noise in the data.     

黏声介质最小平方逆时偏移

介质的黏滞性是普遍存在的.黏滞性介质中的真振幅成像需要校正由介质的黏滞性引起的振幅衰减与速度频散,然而常规的反Q偏移方法存在不稳定问题.本文在反演的框架下求解黏声介质成像问题,在有效避开不稳定的同时实现真振幅成像.首先将波动方程线性化,并依此建立黏声介质最小平方逆时偏移（LSRTM）的目标函数;然后推导波动方程伴随算子,并在此基础上借助伴随状态法推导迭代求解的具体算法;最后通过引入动态相位编码技术将计算量降至与常规逆时偏移相同的数量级.该方法在真振幅成像过程中考虑了介质黏滞性的影响,更接近实际情况,具有更好的振幅保持性.相对于常规逆时偏移,该方法能够自动压制成像噪声,具有更高的成像分辨率及精度.通过模型试算验证了方法的正确性.     

Implementation aspects of attenuation compensation in reverse‐time migration

Attenuation compensation in reverse‐time migration has been shown to improve the resolution of the seismic image. In this paper, three essential aspects of implementing attenuation compensation in reverse‐time migration are studied: the physical justification of attenuation compensation, the choice of imaging condition, and the choice of a low‐pass filter. The physical illustration of attenuation compensation supports the mathematical implementation by reversing the sign of the absorption operator and leaving the sign of the dispersion operator unchanged in the decoupled viscoacoustic wave equation. Further theoretical analysis shows that attenuation compensation in reverse‐time migration using the two imaging conditions (cross‐correlation and source‐normalized cross‐correlation) is able to effectively mitigate attenuation effects. In numerical experiments using a simple‐layered model, the source‐normalized cross‐correlation imaging condition may be preferable based on the criteria of amplitude corrections. The amplitude and phase recovery to some degree depend on the choice of a low‐pass filter. In an application to a realistic Marmousi model with added Q, high‐resolution seismic images with correct amplitude and kinematic phase are obtained by compensating for both absorption and dispersion effects. Compensating for absorption only can amplify the image amplitude but with a shifted phase.     

井中地震三维高效黏声逆时偏移成像方法及应用

井中地震是井中激发地面接收的一种新型地震采集方法, 与地面地震相比, 具有采集成本低, 作业效率高的优点.此外, 井中地震资料信噪比高、频带宽、波场丰富, 可实现井周小尺度构造精细成像.在油藏开发中, 井中地震可实现油藏描述与监测.鉴于上述原因, 近年来, 井中地震受到业内越来越多的关注.然而, 井中地震震源能量弱, 地层吸收衰减效应的影响强于地面地震.因此, 需要发展针对性的井中地震衰减补偿偏移成像方法.本文基于Kelvin-Voigt模型推导新的振幅衰减和相位频散解耦的黏声波动方程, 通过改进激发振幅成像条件, 实现高效的井中三维地震黏声逆时偏移成像方法.此外, 本文采用时变低通滤波器压制高波数噪声, 提升黏声逆时偏移方法的稳定性.数值算例及实际资料应用表明, 本文提出的井中三维黏声逆时偏移成像方法计算效率高、稳定性好, 实用化潜力大.

     

带正则化校正的TTI介质qP波方程及其逆时偏移方法和应用

各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性.     

TTI介质有限差分逆时偏移的稳定性探讨

在沉积学中,可假设在相同时期的沉积层具有相近的物理性质和演化过程.因此,沿层传播的地震波和垂直于地层传播的地震波具有各向异性的特点.在纵波资料的处理中,考虑各向异性对逆时偏移的影响,通常假设介质的横波速度为零,这样可以得到纵波在TTI介质中的传播方程,但是该方程在实际计算中仍存在数值稳定性问题.本文加入横波分量可有效解决数值稳定性问题,并选取适当的横波速度减小对纵波成像的影响,实现地震波在TTI介质中的逆时偏移.实际测算表明,P-SV波的方程中包含横波分量,若假设SV的速度为零,则会导致方程的差分格式不稳定;若加入SV波,选择合适的SV波速度可以使SV波的全区各向异性和反射系数达到极小,并可有效的抑制SV波对纵波勘探的影响.本文的方法是一种稳定的TTI介质中的逆时偏移方法.