逆时偏移技术原理及发展应用展望课件

逆时偏移技术原理及发展应用展望逆时偏移技术原理及发展应用展望五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移技术算法描述二、逆时偏移技术基本原理一、逆时偏移技术的优势主要内容：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移逆时偏移技术的优势

相对其他方法而言，逆时偏移（RTM）用全程波方程对波场延拓，避免对波动方程的近似，因此没有倾角限制，原理上可以利用转换波、棱镜波或多次反射波成像，并获得更精确的振幅等动力学信息，实现保幅成像，还可以更好地对复杂速度场进行更细化更精确的估计。成像方法不受介质速度变化的影响，能够对复杂区域进行较准确的成像。另外，RTM也可以应用于工程地质勘探，为精确工程地质勘探做出贡献。还可以为陡倾角反射层提供稳健的角度域共成像道集。逆时偏移技术的优势相对其他方法而言，逆时偏移（RTM五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移技术算法描述二、逆时偏移技术基本原理一、逆时偏移技术的优势主要内容：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移逆时偏移技术基本原理叠后逆时偏移：

叠后逆时偏移使用的是爆炸反射面成像原理，处理的是水平叠加剖面。叠后逆时计算是从时间剖面的最后一个时间采样点起，逆时外推直到零时间，此时空间所有的振幅值就组成了最终的偏移剖面。逆时偏移技术基本原理叠后逆时偏移：逆时偏移技术基本原理

叠前逆时偏移：叠前逆时偏移是对单炮记录数据进行逆时偏移，然后将各炮成像结果叠加，得到最终的成像剖面。

对单炮记录，它将炮记录的最后一个采样时刻的波场(x,z,T)作为起始平面，按时间反推，并以地震剖面资料u(x,z=0,t)作为每一步进时间的边界条件，得出时间t=0的(x,z)，应用成像条件得到最终偏移结果

u(x,z,t=0)。

逆时偏移技术基本原理叠前逆时偏移：叠前逆时偏移的成像条件：

目前普遍使用的是激发时间成像条件。激发时间成像条件

是震源到每个成像网格点的单程旅行时。

可以把地震记录看成地下不同点在不同时刻以不同强度激发的地震波在地表的叠加结果，这就是激发时间成像条件的基本原理。激发时间成像条件的求取是弹性波逆时偏移的难点之一，我们采用求解程函方程的方法求取地下各点的成像条件。逆时偏移技术基本原理叠前逆时偏移的成像条件：逆时偏移技术基本原理五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移技术算法描述二、逆时偏移技术基本原理一、逆时偏移技术的优势主要内容：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移逆时偏移技术算法描述（一）三维逆时深度偏移的定解问题可描述为：为纵横向可变的介质速度为地表接收的三维地震纪录逆时偏移技术算法描述（一）三维逆时深度偏移的定解问题可描述为逆时偏移技术算法描述上面为3个方向两端的边界条件，其中逆时偏移技术算法描述上面为3个方向两端的边界条件，其中

为导出该方程的离散差分格式，需把观测对应的地下介质分布区域或要对其进行地震波模拟的模型区域离散化，即把它们剖分成一个个的小方块。

令,为导出高阶差分方程，需把波场进行Taylor展开。

逆时偏移技术算法描述

为导出该方程的离散差分格式，需把观测对应的地下介质逆时偏移技术算法描述

然后利用微分和差分关系，得到用于正演模拟和逆时深度偏移的高阶差分方程的起始方程。另外，可以根据需要来组合不同阶次的差分格式。

利用差分法求解上面定义的定解问题可以很好地解决三维逆时深度偏移问题。逆时偏移技术算法描述然后利用微分和差分关系，得到用于正演逆时偏移技术算法描述（二）边界条件做正演和逆时偏移，很重要的一个问题就是解决人工边界的反射。若不能完全消除边界的反射，其剩余能量就会返回到计算区域，经过多次数值积分，剩余能量不断积累增大，将对记录造成一定影响。

有时使用最优系数吸收边界条件仍不能很好的解决人工边界问题。采用多次透射边界和衰减边界相结合的办法，可以较好地解决了边界反射问题。逆时偏移技术算法描述（二）边界条件五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移技术算法描述二、逆时偏移技术基本原理一、逆时偏移技术的优势主要内容：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移

基于常规的有限差分叠前逆时偏移算法，使用求解程函方程得到的激发时间成像条件，我们对凹陷模型和Marmousi模型做了相应的试算。凹陷模型60炮偏移叠加成像结果逆时偏移技术应用效果及发展前景基于常规的有限差分叠前逆时偏移算法，使用求解Marmousi模型单炮叠前偏移结果120炮偏移叠加成像剖面逆时偏移技术应用效果及发展前景Marmousi模型单炮叠前偏移结果120炮偏移叠加成像剖面

目前计算机技术已经发展到可以做3D叠前RTM。Lesage等使用混合有限差分-混谱法算法的3D逆时偏移来解双程声波方程。这一算法使得逆时偏移就可以处理3D各向同性介质、3DVTI介质和3DTTI介质。

RobertSoubaras等还提出新的求解双程波动方程的算法：两步法显式推进算法。BenD.Wards等提出了基于傅里叶变换的高保真时间推进方程，方法结合局部傅里叶变换（Gabor变换）也适用于变速介质。Faqi

Liu等在原标准波动方程的基础上引入一个频散衰减项，以此来达到抗频散的效果。逆时偏移技术应用效果及发展前景目前计算机技术已经发展到可以做3D叠前RTM。Le图1：速度模型（左）和两步法显式推进逆时偏移成像结果（右）逆时偏移技术应用效果及发展前景图1：速度模型（左）和两步法显式推进逆时偏移成像结果（右）逆

目前还发展了多步法逆时偏移，其思想是，根据速度将地下构造沿着深度大体分成2-3块，然后从上到下对各小块区域顺次应用RTM进行成像，这样就使得RTM经济有效，增加了其在实际生产中的实用性。

而且，多步法逆时偏移过程中还可以将不同偏移方法结合应用来进一步提高成像质量和成像效率。根据每一分块的速度场的复杂情况，选择合适的偏移方法，只在必要时利用RTM算法。这样，不但可以保持RTM的优势，而且充分利用其它各种成像方法的优势，降低计算时间，得到高质量的成像。逆时偏移技术应用效果及发展前景目前还发展了多步法逆时偏移，其思想是，根据速

偏移所用的速度模型。此速度模型分为三块，其中黄线定义了顶部区域的底，绿线定义了中间区域的底。逆时偏移技术应用效果及发展前景图2：速度模型分区示意图偏移所用的速度模型。此速度模型分为三块，其中黄线定义图3a：

用RTM成像结果

图3b：

混合成像结果

图3a所示为BP模型左边部分的RTM成像结果，图3b为混合成像结果，顶部和底部区域用单程波动方程偏移，中间区域用RTM进行成像。两种偏移方法所用步长都是12.5m。可以看出两种成像之间没有大的差异。两种成像结果中，盐丘体的陡倾侧翼界面都得到了很好的成像，而且复杂丘体下面的反射界面也都得到了成像。图3a：用RTM成像结果

图3c，d所示为对深源丘体区域所做的类似对比。两者大部分区域成像效果一致，但图3d所示混合成像结果对陡倾同相轴有更好的成像效果，如黄色箭头所示。和图3a所示RTM成像结果相比，图3d的陡倾同相轴更连续，振幅更均衡。图3c：

用RTM成像结果图3d：

混合成像结果图3c，d所示为对深源丘体区域所做的类似对叠前逆时偏移弹性波逆时偏移角度域弹性波逆时偏移垂直地震剖面数据的炮检距逆时偏移逆时偏移的预处理抗频散波动方程逆时偏移相移时间推进法逆时偏移利用逆时偏移为速度分析和建模打开新的领域逆时偏移技术应用效果及发展前景叠前逆弹性波逆时偏移角度域弹性波逆时偏移垂直地震剖面数据的炮RTM的研究方向：1）复杂条件下3D速度分块、相应速度分析与建模和RTM网格计算方法研究。分析整个研究区域的速度分布特点，实现研究区域沿深度方向的有效分块，基于每一分块的速度特点，有针对性地研究相应的速度分析与建模方法，建立细化精确的速度场。研究RTM网格的计算方法，推导相应的计算公式，获取相应区域的网格尺度，为复杂条件下3D多步法分区RTM成像做准备；逆时偏移技术应用效果及发展前景RTM的研究方向：逆时偏移技术应用效果及发展前景RTM的研究方向：2）复杂地质条件下高精度3D有限差分抗频散波动方程RTM算法研究。基于每一分块的速度特点和对成像精度的要求，结合现有方法对起伏地表、高陡构造和复杂储层区等复杂地质条件的适用性和方法自身存在的局限性，研究相应的高精度3D有限差分抗频散波动方程RTM算法，有效减小有限差分RTM过程中的频散；逆时偏移技术应用效果及发展前景RTM的研究方向：逆时偏移技术应用效果及发展前景RTM的研究方向：3）研究新的3D保真成像条件。分析现有成像条件的优缺点和复杂地质体地震保真偏移对成像条件的要求，研究并应用3D倾角保真校正和能量照明补偿，来消除逆散射效应及反射倾角对偏移振幅等动力学信息的影响；逆时偏移技术应用效果及发展前景RTM的研究方向：逆时偏移技术应用效果及发展前景逆时偏移技术应用效果及发展前景RTM的研究方向：4）弹性逆时深度偏移。到目前为止，人们在弹性波叠前逆时深度偏移技术研究方面取得了一定的进展，但总体来说，国内、外在这一领域还处在起步阶段,还需要投入更多的精力来研究它。

目前HuseyinDenli等提出的一种弹性成像的可选步骤是在波场重构和成像时使用全矢量场。波场重构使用弹性波方程的数值解，并且用多分量数据作为边界条件。ERTM是复杂区域真振幅偏移最有远景的工具。逆时偏移技术应用效果及发展前景RTM的研究方向：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移技术算法描述二、逆时偏移技术基本原理一、逆时偏移技术的优势主要内容：五、结论与展望四、逆时偏移技术应用效果及发展前景三、逆时偏移

原则上逆时偏移是实现全波场偏移成像的最佳方法，可以发展为真振幅逆时偏移。随着向量并行计算机及高性能集群