逆时偏移角道集构建及速度分析方法

合集下载

逆时偏移中用Poynting矢量高效地提取角道集

逆时偏移中用Poynting矢量高效地提取角道集王保利;高静怀;陈文超;张唤兰【摘要】Reverse time migration has shown its superior performance for imaging the complex subsurface structures, but it needs a more accurate velocity model. Common-image gathers are an important output of prestack depth migration. They not only can provide velocity model for depth migration, but also can provide amplitude and phase informationfor subsequent subsurface attribute interpretation. In this paper, we have developed a method for extracting angle domain common-image gathers during reverse time migration. In our method, vector wave equation is used to calculate the wavefield, energy flux density vector (Poynting vector) is used to calculate the reflection angle, and finally cross-correlation imaging condition is used to obtain the angle domain common-image gathers. Finally, a model and a field marine data are done to show the effectiveness and advantages of this new method.%逆时偏移在提高复杂介质的成像质量方面表现出了优越的性能,但逆时偏移对速度精度的要求比较高.共成像点道集是一种非常重要的叠前深度偏移输出,它除了能为深度偏移处理提供速度信息外,还能够提供振幅和相位等信息,为后续的属性解释提供依据.本文提出一种在逆时偏移成像过程中提取角度域共成像点道集的方法,该方法采用矢量波动方程进行波场传播,并用能流密度矢量(Poynting vector)计算反射角,最后应用互相关成像条件输出角度域共成像点道集.该方法简单易于实现,且几乎不需要额外的计算量和存储量,非常适合于进行逆时偏移速度分析,同时提出的角道集也能用于进行AVA等分析.最后通过模型算例和实际数据检验了方法的有效性和优越性.【期刊名称】《地球物理学报》【年(卷),期】2013(056)001【总页数】7页(P262-268)【关键词】逆时偏移;角道集;速度分析;Poynting矢量;AVA【作者】王保利;高静怀;陈文超;张唤兰【作者单位】西安交通大学电子与信息工程学院,西安710049;中煤科工集团西安研究院,西安 710077;西安交通大学电子与信息工程学院,西安710049;西安交通大学电子与信息工程学院,西安710049;西安科技大学地质与环境学院,西安 710054【正文语种】中文【中图分类】P6311 引言20世纪80年代就提出的逆时偏移［1－2］在对复杂地下结构进行成像时体现出了较强的优势，因此越来越受到人们的重视［3－10］.多路径、多次波、大倾角、回转波、棱镜波以及保幅等关键问题上，逆时偏移几乎都可以较好地解决.但准确的逆时偏移成像比其它成像方法（如单程波偏移、Kirchhoff偏移和beam偏移等）对速度模型的准确度更加敏感.为了给逆时偏移建立一个较准确的速度模型，需要一种迭代的方法来不断地更新速度.而共成像点道集（CIGs）是叠前深度偏移中的很重要的输出道集，它可以提供速度建模所需要的信息，以及反映地下岩性的振幅和相位等信息.正确的速度模型得到的共成像点道集中的同相轴是平的，而不正确的速度得到的是弯曲的，且弯曲的程度与速度校正量是有关系的，通过这种关系就可以不断地更新速度模型［11］.通常，有两种共成像点道集，一种是偏移距域的（common－offset CIGs，简称为COCIGs），一种是角度域的（angle－domain CIGs，简称为 ADCIGs）.在地下反射体是水平层状的情况，这两种道集是可互换的.Kirchhoff偏移和beam偏移方法可以很方便快速地输出COCIGs［12］.但多路径问题会影响到COCIGs的质量，进而会影响速度分析和AVO分析的精度.而单程波偏移和逆时偏移中的COCIGs是很难提取的.Save和 Formel（2003）［13］利用单程波偏移在二维情况下首先提取了地下局部偏移距域的CIGs，然后再转换为ADCIGs.Fomel（2004）［14］将这种方法推广到三维情况，并指出三维情况下该方法非常复杂.Save和Fomel（2006）［15］又提出一种新的方法，该方法需要先生成延时CIGs，然后再转换为 ADCIGs.Xu等（2010，2011）［16－17］提出一种直接从三维逆时偏移中提取ADCIGs的方法，该方法在频率波数域求取传播方向，进而计算出反射角和方位角信息，最后使用给定的成像条件直接输出ADCIGs.这种方法可以得到质量较好的ADCIGs，但需要保存所有时刻的炮点波场值和接收点波场值，并进行四维的傅里叶变换，计算量和存储量都非常得大，不适合进行逆时偏移速度分析等处理.Yoon和 Marfurt（2003）［5］用坡印廷矢量来确定传播方向并用此修改成像条件来去除逆时偏移中的低频噪声，取得了较好的效果.本文借助于这种思想，用一阶应力速度方程（与常规标量应力波动方程相比，增加了额外的存储量，但减少了计算坡印廷矢量的耗时）进行波场延拓，并用坡印廷矢量求取成像点处的传播方向，最后使用给定的成像条件输出逆时偏移中的ADCIGs.2 方法原理2.1 波场传播逆时偏移是通过波动方程在时间上对采集到的地震数据进行反向传播来实现的.其偏移过程使用的是双程波波动方程，因此不受倾角限制，且能正确地对回转波、多次波等进行成像.逆时偏移的原理本文不再赘述.通常，逆时偏移采用的声波波动方程为标量应力波动方程［18－22，10］，其形式如下：其中，υ为介质纵波速度，P为质点应力波场.该波动方程不利于使用坡印廷矢量（原因见下节）.因此本文采用一阶应力－速度声波波动方程：式中：ρ为介质密度，υ为介质纵波速度，v和w 分别为x和z轴上的质点振动速度.本文对公式（2）采用交错网格技术［18］来构建波场传播算子.考虑到逆时偏移的计算量和为避免数值频散问题对逆时偏移结果的影响，在时间方向采用二阶而在空间方向采用高阶差分格式（≥8阶）.2.2 Poynting矢量求取反射角坡印廷矢量也称为能流密度，1884年由坡印廷提出用来测量单位时间内穿过垂直于此矢量方向的单位表面的电磁场能量.Poynting矢量值等于电场强度与磁场强度的叉积，矢量方向代表电磁波的传播方向.Yoon和 Marfurt（2003）［5］给出了地震波场中的坡印廷矢量计算式：图1 坡印廷矢量快照图（a）水平分量；（b）垂直分量.Fig.1 The snapshots of Poynting vector（a）Horizontal component；（b）Vertical component.其中，v表示位移速度矢量，v＝（v，w），P为应力波场.通过（3）式可知，如果使用（1）式的标量声波波动方程，需计算应力波场值的时间和空间导数，计算比较复杂.而如果采用（2）式的一阶应力速度方程，则应力P和位移速度矢量v都是逆时偏移过程中已有的，不需任何额外计算.图1是均匀介质中的某时刻计算的坡印廷矢量快照图，可以很明显地判断出波的大致传播方向：在水平分量图（a）中蓝色和红色分别代表向右传和向左传，而在垂直分量图（b）中蓝色和红色则分别代表向下传和向上传.用SS和SR分别表示炮点波场的坡印廷矢量和接收点波场的坡印廷矢量，则根据余弦定理，可以求出两矢量的夹角θ，由此可以得到波场的反射角为β＝.为了增加抗噪性能，夹角θ可在一小时窗内进行计算，如2.3 提取ADCIGs获得反射角β后，便可以进行ADCIGs的提取.常规的逆时偏移成像条件是对所有角度道集的叠加，因此不能直接用来提取ADCIGs，需要进行修改.考虑到计算得到的反射角通常与我们给定的采样角度（一般是等间隔的）不一致，我们用高斯采样函数将常规的成像条件其中S和R分别表示炮点波场和接收点波场（本文中为应力波场值），而β表示当前时刻该成像点处计算得到的反射角，βk表示ADCIGs中的离散角，σ表示高斯函数的方差.根据（7）便能容易地提取出ADCGIs.需特别说明的是，（5）式中，分母在很小甚至为零情况下，计算得到的角度β误差较大，但此时，其能量贡献也非常小，因此不会影响最终角道集的质量.具体的计算流程如图2.3 模型试验3.1 Marmousi模型的ADCIGs用本文提出的方法提取出的Marmousi模型的ADCIGs如图3所示.角度从0°～90°，角度间隔为2°，每100个CMP显示一个ADCIG，一共23个ADCIGs.从图上可以看出每个ADCIG上的同相轴都是拉平的，且低频噪声主要分布在大角度的道上.图4为第500个CMP位置处的ADCIG.3.2 低频噪声逆时偏移的低频噪声主要是由透射波成像造成的，而入射波和透射波的夹角在模型比较光滑情况下接近于180°（传播路径一微小段范围内），因此可通过在角道集上对大角度的道进行衰减或直接置零来对逆时偏移成像剖面进行去噪.理论上这种去噪方法等同于拉普拉斯滤波［7，11］.图5、图6和图7分别为用0°～40°、40°～70°和70°～90°等角道集叠加后的结果，可明显看出在0°～40°内几乎没有任何低频噪声，而在70°～90°的偏移剖面几乎全部是低频噪声.图2 Poynting矢量提取角道集的流程图Fig.2 The flow chart of extracting ADCIGs using Poynting vector图4 Marmousi模型第500个CMP位置处的ADCIGFig.4 A ADCIG located at 500th CMP of Marmousi model3.3 AVA角道集的另一个用途就是用于AVA 分析［23－25］，本文使用单层模型分析了角道集上振幅的特征.图8为用振幅随角度递增的模型参数（上下层速度分别为1500m／s和1800m／s，密度均为2.1g／cm3）得到的角道集AVA曲线与理论AVA曲线的对比.图9为用振幅随角度递减的模型参数（上下层速度分别为1500m／s和1300m／s，密度均为2.1g／cm3）得到的角道集AVA曲线与理论AVA曲线的对比.对比可以看出，两种模型得到的角道集的AVA曲线和理论AVA曲线基本一致.4 实际数据试验图10为从某地区的一条测线中提取出的部分ADCIGs.可以看出1.3～1.8km段中同相轴基本是平的，而在1.8～2.3km断层同相轴是向下弯曲的，说明偏移过程中所用的速度偏高.通过不断迭代地进行速度调整即可实现逆时偏移速度分析.图10 实际数据中提取的部分角道集Fig.10 Part of the ADCIGs from field data 5 结论本文在逆时偏移过程中采用一阶应力速度波动方程进行波场传播，并用Poynting 矢量来确定传播方向并计算反射角，进而通过给定的成像条件进行逆时偏移角道集的提取.整个计算过程中几乎不额外增加计算量和存储量，提取的共角度道集可用于进行速度分析和AVA分析等.该方法可用于三维逆时偏移中角道集等的提取（包括方位角道集），实现简单，易于进行并行加速（如GPU等）.后续工作将在此工作基础上，使用提取的角道集进行逆时偏移速度分析，以得到较好的偏移速度剖面.参考文献（References）［1］Whitmore N D.Iterative depth migration by backward time propagation.53rd Annual International Meeting，SEG Expanded Abstracts，1983，2：827－830.［2］Baysal E，Kosloff D D，Sherwood J W C.Reverse timemigration.Geophysics，1983，48（11）：1514－1524.［3］张美根，王妙月.各向异性弹性波有限元叠前逆时偏移.地球物理学报，2001，44（5）：711－719.Zhang M G Wang M Y.Prestack finite element reverse－time migration for anisotropic elastic waves.Chinese J.Geophys.（in Chinese），2001，44（5）：711－719.［4］Yoon K，Shin C，Hong S.3Dreverse－time migration using acoustic wave equation.An experience with the SEG／EAGE data set.The Leading Edge，2003，22（1）：38－41.［5］Yoon K，Marfurt Kurt J，Starr W.Challenges in reversetime migration.74th Annual International Meeting，SEG Expanded Abstracts，2003：1057－1060.［6］Bednar J B，Bednar C J，Shin C.Two－way versus one－way：A case study style comparison.76th Annual International Meeting，SEG Expanded Abstracts，2006，25：2343－2347.［7］Zhang Yu，Xu S，Bleistein N，et al.True－amplitude，angledomain，common－image gathers from one－way wave－equationmigrations.Geophysics，2007，72（1）：S49－S58.［8］Zhang Y，James S. Practical issues of reverse time migration：True amplitude gathers，noise removal and harmonic－source encoding.First Break，2009，27（1）：53－60.［9］Liu F Q，Zhang G Q，Morton S A，et al.Reverse time migration using one way wavefield imaging condition.77th Annual International Meeting，SEG Expanded Abstracts，2007，26（1）：2170－2174［10］Liu F Q，Zhang G Q，Morton S A，et al.An effective imaging condition for reverse－time migration using wavefielddecomposition.Geophysics，2011，76（1）：S29－S39.［11］Zhou H B，Gray S H，Young J，et al.Tomographic residual curvature analysis：The process and its components.73rd Annual International Meeting，SEG Expanded Abstracts，2003，22：666－669. ［12］Xu S，Chauris H，Lambaré G，et mon－angle migration：A strategy for imaging complex media.Geophysics，2001，66（6）：1877－1894.［13］Sava P C，Formel S.Angle－domain common－image gathers by wavefield continuation methods.Geophysics，2003，68（3）：1065－1074. ［14］Fomel S.Theory of 3－D angle gathers in wave－equation imaging.74th Annual International Meeting，SEG Expanded Abstracts，2004，1（1）：1053－1056.［15］Sava P，Formel S.Time－shift imaging condition in seismic migration.Geophysics，2006，71（6）：S209－S217.［16］Xu S，Zhang Y，Tang B.3Dcommon image gathers from Reverse time migration.80th Annual International Meeting，SEG Expanded Abstracts，2010：3257－3260.［17］Xu S，Zhang Y，Tang B.3Dangle gathers from Reverse time migration.Geophysics，2011，76（2）：S77－S92.［18］李博，刘国峰，刘洪.地震叠前时间偏移的一种图形处理器提速实现方法.地球物理学报，2009，52（1）：245－252.Li B，Liu G F，Liu H.A method of using GPU to accelerate seismic pre－stack time migration.ChineseJ.Geophys.（in Chinese），2009，52（1）：245－252.［19］李博，刘红伟，刘国峰等.地震叠前逆时偏移算法的CPU／GPU实施对策.地球物理学报，2010，53（12）：2938－2943.Li B，Liu H W，Liu G F，et putational strategy of seismic pre－stack reverse time migration on CPU／GPU.Chinese J.Geophys.（in Chinese），2010，53（12）：2938－2943.［20］刘红伟，李博，刘洪等.地震叠前逆时偏移高阶有限差分算法及 GPU实现.地球物理学报，2010，53（7）：1725－1733.Liu H W，Li B，Liu H，etal.The algorithm of high order finite difference pre－stack reverse time migration and GPU implementation.Chinese J.Geophys.（in Chinese），2010，53（7）：1725－1733.［21］刘红伟，刘洪，邹振等.地震叠前逆时偏移中的去噪与存储.地球物理学报，2010，53（9）：2171－2180.Liu H W，Liu H，Zou Z，et al.The problems of denoise and storage in seismic reverse time migration.ChineseJ.Geophys.（in Chinese），2010，53（7）：2171－2180.［22］刘红伟，刘洪，李博等.起伏地表叠前逆时偏移理论及GPU加速技术.地球物理学报，2011，54（7）：1883－1892.Liu H W，Liu H，Li B，et al.Pre－stack reverse time migration for rugged topography and GPU acceleration technology.Chinese J.Geophys.（in Chinese），2011，54（7）：1883－1892.［23］杨午阳.模拟退火叠前AVA同步反演方法.地球物理学进展，2010，25（1）：219－224.Yang W Y.Pre－stack simultaneous AVA inversion algorithm with simulated annealing.Progress in Geophys.（in Chinese），2010，25（1）：219－224.［24］王保利.弹性阻抗反演及应用研究.地球物理学进展，2005，20（1）：89－92.Wang B L.Elastic impedance inversion and its application.Progress inGeophys.（in Chinese），2005，20（1）：89－92.［25］喻岳钰，杨长春，王彦飞等.叠前弹性阻抗反演及其在含气储层预测中的应用.地球物理学进展，2009，24（2）：574－580.Yu Y Y，Yang C C，Wang Y F，et al.Application of prestack seismic elastic impedance inversion to gas reservoir.Progress in Geophys.（in Chinese），2009，24（2）：574－580.。

偏微分方程中的反问题和逆时偏移算法

偏微分方程中的反问题和逆时偏移算法在偏微分方程领域中，反问题是一个具有挑战性的研究领域，它涉及从一些已知的观测数据中恢复出未知的模型参数或边界条件。

而逆时偏移算法则是一种常用于地震勘探领域的反问题求解方法。

本文将介绍偏微分方程中的反问题以及逆时偏移算法的基本原理和应用。

一、偏微分方程中的反问题1. 反问题的定义及挑战性在偏微分方程的求解过程中，通常需要确定一些未知量，如模型参数或边界条件。

反问题即是对于已知的观测数据，求解这些未知量的问题。

而由于观测数据的不完备性或噪声的存在，反问题的求解往往具有困难和挑战性。

2. 常见的反问题求解方法在求解反问题过程中，常见的方法包括正则化方法、Bayes方法、最小二乘法以及变分方法等。

利用这些方法可以对反问题进行数值求解，从而恢复出未知的模型参数或边界条件。

二、逆时偏移算法1. 逆时偏移算法的基本原理逆时偏移算法是一种常用于地震勘探中的反问题求解方法。

其基本思想是利用地震波场数据和Born近似模型来恢复地下介质信息。

逆时偏移算法可以分为两个步骤：正演模拟和逆时偏移。

2. 逆时偏移算法的步骤逆时偏移算法的第一步是正演模拟，即通过假设一组地下模型参数和边界条件，计算出地震波场数据。

这一步骤可以使用传统的有限差分或有限元方法进行计算。

逆时偏移算法的第二步是逆时偏移，即通过将正演得到的波场数据与观测到的数据进行匹配，估计出地下模型的参数和边界条件。

这一步骤可以使用最小二乘法或变分方法等进行求解。

3. 逆时偏移算法的应用逆时偏移算法广泛应用于地震勘探领域，可以用于识别地下构造和油气储层等。

此外，逆时偏移算法还可以用于医学成像、非损检测、声学波传播等领域。

总结本文介绍了偏微分方程中的反问题和逆时偏移算法。

反问题是一个具有挑战性的研究领域，需要从已知的观测数据中恢复出未知的模型参数或边界条件。

逆时偏移算法是一种常用的反问题求解方法，通过正演模拟和逆时偏移两个步骤可以恢复地下介质的信息。

基于三角网格的有限差分法叠后逆时偏移

基于三角网格的有限差分法叠后逆时偏移
三角网格有限差分法叠后逆时偏移（Triangular-mesh Finite Difference Method Applied after Inversion of Tme Migration）是一种基于有限差分法的走时偏移修正技术，由于该技术具有计算效率高，抗粒子带宽增大的优势，近年来得到了广泛的应用。

三角网格有限差分法叠后逆时偏移是由三角网格有限差分（TFD）方法结合叠后合成算法（FIP）和振幅位移补偿（AGC）技术等而构建的，它的基本思想是将在水平面上有效模拟的有限差分正演波累积，通过一维正演和逆时走时，有效地降低走时偏移带来的影响。

技术优势除了高效率，FIP和AGC技术使得三角网格有限差分法叠后逆时偏移具有长程和宽粒子带的优势，比如FIP技术可以有效地克服孔隙效应，AGC技术可以提高振幅精度，而基于三角网格有限差分的正演模拟计算效率也能达到较高水平。

应用优势三角网格有限差分法叠后逆时偏移技术在测井、地震勘探、地球物理以及石油和天然气资源勘查等应用中得到了广泛应用，特别是在大地深处地震测量和目标地层分析方面，对比传统走时偏移技术而言，整个方案能够更好地提高处理的精度。

总的来说，三角网格有限差分法叠后逆时偏移技术由高效率、抗噪声强度和长程宽粒子带等优势，近年来被石油勘探界广泛应用。

它不仅是实现精确、长距离测量的有力方法，还为勘探新技术奠定了坚实的基础，对推动勘探技术的发展起到了重要作用。

地震波逆时偏移方法研究

地震波逆时偏移方法研究1. 引言1.1 背景介绍地震是地球内部产生的突发性能量释放现象，往往会引发破坏性的地质灾害，对人类生命财产安全造成威胁。

为了更加准确地研究地震发生的机制和预测地震的可能性，地震勘探技术逐渐成为地震学领域中的重要工具。

地震波逆时偏移方法是一种重要的地震勘探技术，通过逆时偏移处理地震波数据，可以有效地提高地震成像的分辨率和准确性。

地震波逆时偏移方法基于地震波在介质中的传播原理，利用地震波数据中的反射波和折射波信息，通过逆时偏移算法将地下物质的分布情况反演出来。

这一方法在地震勘探领域有着广泛的应用，可以用于石油勘探、地质灾害监测等方面。

由于地震波逆时偏移方法具有较高的分辨率和准确性，因此在地震学研究和勘探实践中得到了广泛的应用。

1.2 研究意义地震波逆时偏移方法是地震勘探领域中一种重要的地震成像技术，其具有极为重要的研究意义。

地震波逆时偏移方法可以有效地提高地震成像的分辨率和准确性，帮助地质学家更准确地解释地下结构，为油气勘探和地质灾害预测提供重要的技术支持。

地震波逆时偏移方法还可以在地震勘探过程中减少噪音干扰，提高数据处理效率，从而降低勘探成本。

随着地震波逆时偏移方法的不断发展和完善，其在地震监测、地质构造研究、矿产勘探等领域的应用也将得到进一步拓展，为地质科研和资源勘探带来更多的可能性。

深入研究地震波逆时偏移方法的原理和应用，对于推动地震勘探技术的进步，促进地球科学研究的发展具有重要的意义。

1.3 研究目的地震波逆时偏移方法是一种重要的地震成像技术，通过对地震波传播路径进行逆向推断，可以获取地下结构的高分辨率成像。

本文旨在探讨地震波逆时偏移方法在地震勘探领域的应用及其优势，并对其存在的问题进行深入分析。

具体目的包括以下几点：1.深入理解地震波逆时偏移方法的原理及步骤，掌握其核心技术，为进一步研究和实践提供坚实的基础。

2.探讨地震波逆时偏移方法在地震勘探中的优势，比较其与传统方法的差异，为选择适合的成像技术提供参考。

波场分解算法与逆时偏移角道集

第57卷第4期 2018年7月石油物探GEOPHYSICAL PROSPECTING FOR PETROLEUMVol .57,No .4Jul .，5018王美霞，张晓慧，张钋，等.波场分解算法与逆时偏移角道集*+石油物探，2018,07 (4): 050-069WANG Meixia#ZHANG Xiaohui # ZHANG Po,et al.'Wavefield decomposition and RTM angle gathers [J+Cjeophysical Prospecting for Petroleum , 2018,07(4) : 000-069(l.S ta to il G u lf Services ,H ou sto n 77〇42;2.前斯塔托尔北京技术服务有限公司，北京100022)摘要:逆时偏移生成的角度域共成像点道集(ADCIGs )有利于提高成像质量，可用于偏移速度分析和A V A /A V O 研究等。

使用Poyntmg 矢量估计波传播方向可以快速生成角道集，但当模型和波场比较复杂时，不同方向传播的波场分量相互重叠，无法准确计算波的传播方向。

波场分解可以提高Poyntmg 矢量的准确性。

为此，研究了波场分解算法，在计算Poyntmg 矢量之前先将波场分解成不同方向传播的分量，用于生成角道集。

为解决传统的傅里叶变换在频率-波数域进行波场分解时计算量大的问题，对波场分解算法进行了优化，并提出了 SSE (单指令多数据流扩展）向量化和多线程并行计算实现算法，从而提高了计算效率。

理论和数值计算结果表明:使用波场分解可以显著提高逆时偏移角道集的质量，此优化算法能够得到与传统方法一致的结果并且提高了计算效率。

二维S gsbedA 模型和三维SEAM T T I 模型的数值计算也证明了该算法用于生成高质量角道集和逆时偏移图像的有效性和稳定性。

优化后的波场分解算法并不局限于上、下行波分解，可容易地拓展到其它方向的波场分解，也可应用于最小二乘偏移和全波形反演，以消除成像噪声，提高成像质量。

逆时偏移

Hemon, C., 1978 Whitmore，1983；McMechan， 1983；Baysal等，1983；Loewenthal等，1983
逆时偏移技术原理、算法及应用
• 国际上，借助CPU实现的三维逆时偏移需要数千个CPU核，使得我国大部分勘探单位望尘莫及。不仅如此，由于规模化实验困难，国内研发三维叠前逆时偏移、前全波形反演、三维波动方程消除多次波曾经长期不能工业化，使得国际大型地球物理服务公司得以对我国技术封锁，成为提高地震数据定量描述能力的主要技术瓶颈。
速度模型
逆时偏移技术原理、算法及应用
“海量数据高端成像技术”
• 塔里木盆地在8000米深，仍存在液态烃；今年五月，华北油田在霸县凹陷雾迷山组6067米碳酸盐岩潜山，打出了油气当量日产1300吨的探井； • 勘探深度不断加深，造成了勘探面积不断增加，地震数据量呈几何级数增加,中石油中亚滨里海盐下勘探区数据达到1.4T，川东北一个勘探区块达到7个TB。 • 开发区高密度地震同样面临海量数据量问题。大庆油田、胜利油田和准噶尔油田分别采集了15TB、 15TB 和5TB的高密度资料。 • 成像处理方法购买、运行、维护、配套成本大大增加，研发和推广十分困难， “海量地震数据成像处理技术”长期受制于国外
逆时偏移技术原理、算法及应用
图形处理器（GPU）超算为高精度地震成像处理技术带来新的机会（加速5-10倍）优势：多核小粒度大规模并行计算，速度快，显存与GPU带宽大限制：缓存小、内存小、PCIE(显卡与CPU）数据交换带宽有限，编程难度较大
•GPU的全称是图形处理器(Graphic Processing Unit)，是NVIDIA公司推出的一款用于计算机显示的设备，俗称“显卡”。 2007年

地震波逆时偏移方法研究

地震波逆时偏移方法研究地震波逆时偏移方法（Reverse Time Migration，RTM）是一种新型的地震成像方法，具有较高的精度和分辨率，广泛应用于油气勘探、地震地质研究等领域。

本文介绍了地震波逆时偏移方法的基础原理、算法和应用研究现状。

地震波逆时偏移方法是利用地震波在地下传播与反射的特性实现对地下结构的成像。

其基本原理是以地震波源为中心，将地震记录数据在时间轴上倒序反演到地震波源处，然后进行反射成像。

具体来说，地震波逆时偏移方法主要包括以下步骤：1、前向传播：在地震波源处施加指定波形的地震震源，将地震波信号传播到每个模型单元。

2、反演求解：根据反演方程，利用上一时刻网格单元中的压力场信息和速度模型，计算当前时刻的速度场和压力场。

同时，计算观测数据的残差，通过残差的逆时中心分散源分布对速度模型进行校正。

3、反向传播：反推每个时刻的波场信息，得到在地震波源处反射回来的应力波形，从而实现成像。

在地震波逆时偏移方法的实现中，需要采用适当的算法来计算速度模型和波场信息。

下面分别介绍常用的有限元方法、有限差分方法和偏移算法。

1、有限元方法有限元方法是一种数值方法，通过将地下结构离散化成有限个结构单元，采用形函数法和单元刚度矩阵计算波场信息。

有限元方法的优点是可以很好地处理波传播和反射现象，但计算量较大，需要较高的计算效率和处理力。

有限差分方法是一种数值离散方法，采用差分算子计算相邻单元间的差分，采用传播规则更新波场信息。

有限差分方法计算速度模型较为简单，但需要大量的内存和计算资源。

3、偏移算法偏移算法是一种基于波动方程的成像算法，具有较高的成像精度和分辨率。

偏移算法主要由反演、卷积、积分三个部分组成。

通过从地震数据中提取反射信息，根据波动方程求解反传波场信息，再与传播波场信息卷积运算得到成像结果。

地震波逆时偏移方法已经成为研究地下结构、油气田勘探等领域的重要工具。

目前，该方法在地震资料处理、反演成像、油气勘探等方面得到广泛应用。

一种提取逆时偏移角道集的方法[发明专利]

专利名称：一种提取逆时偏移角道集的方法
专利类型：发明专利
发明人：田坤,于海铖,王常波,张学涛,步长城,李国磊,王彦军,张敬仁,揭景荣
申请号：CN201710448348.0
申请日：20170919
公开号：CN107479091A
公开日：
20171215
专利内容由知识产权出版社提供
摘要：本发明涉及油气勘探地震资料处理领域，特别是涉及一种提取逆时偏移角道集的方法，该提取方法提高了复杂区域角道集的提取精度，并且很好地解决了同相轴交叉的问题，对复杂区域的成像效果好，有利于进行偏移速度分析。

该提取方法包括有对预设声波模型进行有限次差分正演计算得到输入炮记录；自震源点正向外推得到炮点波场；自检波点反向外推得到检波点波场；傅里叶变换得到炮点局部慢度域、检波点局部慢度域，分解得到炮点局部角度分量、检波点局部角度分量；确定角度域成像条件提取得到逆时偏移角道集等步骤。

申请人：中国石油化工股份有限公司,中国石油化工股份有限公司胜利油田分公司物探研究院地址：100728 北京市朝阳区朝阳门北大街22号
国籍：CN
代理机构：东营双桥专利代理有限责任公司
代理人：方圆
更多信息请下载全文后查看。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

解决多路径问题，减少偏移假象，提高道集成像质量；其次，推导了新的纵波角度域速度更新公式，并且实现了逆时
偏移角度域速度分析方法，模型数据的数值实验结果表明了该方法的正确性和有效性。
关键词：逆时偏移；波前矢量；角道集，速度分析
中图分类号：Ｐ６３１。４
Ａｂｓｔｒａｃｔ：Ａｔｐｒｅｓｅｎｔｔｈｅｐｒｅｃｉｓｉｏｎｏｆｍｉｇｒａｔｉｏｎｖｅｌｏｃｉｔｙｍｏｄｅｌｉｓｈｅｔｋｅｙｐｒｏｂｌｅｍｆｏｒｔｈｅａｐｐｌｉｃａｔｉｏｎｅｆｅｃｔｏｆｒｅｖｅｒｓｅｔｉｍｅｍｉｇｒａｔｉｏｎ（ＲｒＭ）．
中国石油勘探
第１９卷第３期
ＣＨＩＮＡＰＥＴＲＯＬＥＵＭＥＸＰＬＯＲＡＴＩＯＮ２０１４年５月
ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７２－７７０３．２０１４．０３．００８
ｄｏｍａｉｎ．Ｔｈｅｎｕｍｅｉｃｒａｌｔｅｓｔｒｅｓｕｌｔｓｏｆｍｏｄｅｌｄａｔａｓｈｏｗｔｈｅｃｏｒｒｅｃｔｎｅｓｓｏｆｔｈｅｍｅｈｏｔｄｉｎｔｈｉｓｐａｐｅｒ．
文献标识码：Ａ
ＭｅｔｈｏｄｆｏｒＲＴＭＡｎｇｌｅＤｏｍａｉｎａｎｄＶｅｌｏｃｉｔｙＡｎａｌｙｓｉｓ
ＬＢｉｎ
（ＮｏｒｔｈｗｅｓｔＢｒａｎｃｈｏｆＰｅｔｒｏＣｈｉｎａＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＰｅｔｍ￣ｕｍＥｘｐｌｏｒａｔｉｏｎ＆Ｄｅｖｅｌｏｐｍｅｎｔ）
逆时偏移角道集构建及速度分析方法
吕彬
（中国石油勘探开发研究院西北分院）
摘
要：速度模型精度问题已经成为制约逆时偏移实际应用效果的瓶颈，尤其是在复杂构造区，受制于目前积分
法偏移速度分析固有的多路径问题缺陷，得到的成像道集质量和偏移速度模型精度往往难以满足逆时偏移的成像要求。基于此，开展逆时偏移角道集构建及速度分析方法的研究。首先，开展了波前矢量逆时偏移角度域共成像点道集构建方法的研究，该方法物理意义明确，计算效率较高。与传统偏移距域及炮域成像道集相比，逆时偏移角道集能够有效
Ｂａｓｅｄｏｎｔｈｉｓ，ｗｅｐ￣ｓｅｎｔＲＴＭ — ｂａｓｅｄｍｉｒａｇｔｉｏｎｖｅｌｏｃｉｔｙａｎａｌｙｓｉｓｍｅｈｏｔｄｉｎｈｅｔｎｇａｌｅ－ｄｏｍａｉｎｉｎｈｉｔｓｐａｐｅｒ．Ｆｉｒｓｔ，ｗｅｄｅｖｅｌｏｐｈｅｔｒｅｓｅａｒｃｈｏｆｗａｖｅｆｒｏｎｔｐｏｌａｒｉｚａｔｉｏｎｍｅｈｏｔｄｔｏｃａｌｃｕ１ａｔｅＲＴＭ－ｂａｓｅｄｎｇａｌｅｄｏｍａｎｉｃｏｍｍｏｎ — ｉｍａｇｅｇａｔｈｅｒｓ．Ｔｈｉｓｍｅｔｈｏｄｈａｓｄｅｉｎｆｉｔｅｐｈｙｓｉｃａｌｍｅａｎｉｎｇａｎｄｈｉｇｈｃｏｍｐｕｔａｔｉｏｎａｌｅｉｃｆｉｅｎｃｙ．Ｃｏｍｐｒｅａｄｗｉｔｈｃｏｎｖｅｎｉｏｔｎａｌｏｆｓｅｔｄｏｍａｉｎｎｄａｓｈｏｔｄｏｍａｉｎｃｏｍｍｏｎ－ｉｍａｇｅｇａｔｈｅｒｓ，ｎｇａｌｅｄｏｍａｉｎｃｏｍｍｏｎ・ｉｍａｇｅｇａｔｈｅｒｓｃｏｕｌｄｅｆｅｃｔｉｖｅｌｙｓｏｌｖｅｍｕｌｔｉ — ｐａｔｈｐｒｏｂｌｅｍ，ｒｅｄｕｃｅｍｉｇｒａｔｉｏｎａｒｔｉｆａｃｔｓｎｄａｉｍｐｒｏｖｅｉｍａｇｉｎｇｑｕａｌｉｔｙｏｆｇａｔｈｅｒｓ．Ｓｅｃｏｎｄ，ｗｅｄｅｉｖｒｅｎｅｗａｎｇｌｅ－ｄｏｍａｉｎｍｉｒａｇｔｉｏｎｖｅｌｏｃｉｔｙｕｐｄａｔｅｅｑｕａｔｉｏｎａｎｄｈａｖｅｆｏｒｍｅｄＲＴＭ－ｂａｓｅｄｍｉｇｒａｔｉｏｎｖｅｌｏｃｉｔｙａｎａｌｙｓｉｓｍｅｈｏｔｄｉｎｈｅｔａｎｇｌｅ —