本发明专利技术提供了一种保幅平面波叠前深度偏移方法。该方法包括：对原始炮集数据进行预处理，以得到统一基准面的炮集数据；针对预定的不同入射角度平面波，基于保幅偏移算子，获得对应于各个平面波入射角度的保幅偏移结果；以及对所获得的对应于各个平面波入射角度的保幅偏移结果进行叠加，获得最终的成像结果。根据本发明专利技术的方法，在获得高精度平面波偏移成像结果的同时，还可以具有相对振幅保持性，实现了平面波保幅叠前深度偏移成像，为后续的振幅随偏移距和振幅随角度变化分析提供更精确的道集。
【技术实现步骤摘要】
本专利技术涉及石油勘探中地震资料处理领域的叠前深度偏移技术，具体地涉及一种。
技术介绍
地震偏移理论与技术的发展是地震勘探形势、应用地球物理学以及计算机技术等相关领域飞速发展推动的结果。地震偏移成像提供着精度越来越高的成像结果，更好的服务于油气藏勘探。基于波动方程的偏移方法被公认为是具有较好偏移效果的方法，其中包括了基于射线、单程波方程的和全声波方程的偏移方法。基于射线理论的KirchhofT积分法具有较高的计算效率，但对于复杂介质中的焦散、多重路径和干涉等现象，该方法显得无能为力。 基于全声波方程的叠前逆时偏移方法具有较高的计算精度，可以适应速度场的强横向变换，但是计算量大，计算效率比较低。相比较而言，基于单程波方程的偏移方法在可以适应速度场强横向变化的同时具有较高的计算效率。在单程波叠前深度偏移方面，从早期开始就出现了相移加内插、分步傅里叶(文献)、傅里叶有限差分法(文献)以及相屏类的双域偏移方法。数学和物理领域新兴的理论和工具为波动方程偏移成像技术提供新的思路和方法，出现了相空间小波分析地震波场偏移成像的算法和Hamilton体系下地震波场延拓的辛群算法和李群算法。文献将小波变换方法引入到单程波偏移方法中，提出了不同基函数下的小波束叠前深度偏移方法。 上述的单程波叠前深度偏移方法只能够保证走时的正确性，保证构造成像的准确性。随着勘探的深入，工业界希望在得到构造成像的同时，能够得到地下反射系数信息，也就是保幅偏移技术。基于文献提出的单程波分裂方法，文献将常规偏移方法改造为保幅的偏移算法并且证明其在高频渐近的意义下结果等价于KirchhofT反演的结果，在波场延拓过程中通过有限差分算法求的振幅补偿项补偿振幅损失。文献、和采用了有限差分法保幅叠前深度偏移算子，但是，由于单程波保幅偏移中引入了振幅补偿项使得计算效率也随之降低，尤其是在处理大面积区块资料的时候更加制约了其应用。引用文献 Stoffa P. L . Split-step Fourier migration . Geophysics, 1990，55(4):410-421 Ristow D. and Riihl T. Fourier f inite-difference migration . Geophy sics,1994，59(12):1882-1893ffu, R. S. , Y. Wang and M. Q. Luo. Beamlet migration using local cosine basis · Geophysics, 2008，73(5):S207-S217张关泉.波动方程的上行波方程和下行波的耦合方程组.应用数学学报，1993，16(2) :251 263Zhang, Y. , G. Q. Zhang, and N. Bleistein. True amplitude wave equationmigrationarising from true amplitude one-way wave equations:1nverse Problems, 2003b, 19，1113-1138张宇.振幅保真的单程波方程偏移理论.地球物理学报，2006，49(5): 1410-1430刘定进，印兴耀.傅立叶有限差分法保幅叠前深度偏移方法.地球物理学 报，2007，42(1) 11-16崔兴福，张关泉，吴雅丽.三位非均匀介质中真振幅地震偏移算子研究.地 球物理学报，2004，47 (3) =509-51
技术实现思路
本专利技术的目的在于提高单程波保幅叠前深度偏移的计算精度和效率，更好的服务 于岩性油气藏的勘探。本专利技术的一个方面提供了一种，该方法包括预处理步骤对原始炮集数据进行预处理，以得到统一基准面的炮集数据；偏移步骤针对预定平面波入射角度范围中的各个平面波入射角度，基于保幅偏 移算子，获得对应于各个平面波入射角度的保幅偏移结果；以及叠加步骤对所获得的对应于各个平面波入射角度的保幅偏移结果进行叠加，以 获得最终的偏移结果。根据上述，可以提高偏移成像保幅程度，同时具有 较高的计算效率。在速度场准确的情况下，该方法在得到高精度构造成像的同时，还可具有相对振 幅保持成像的能力，提供不同入射角平面波的成像结果，为后续的振幅随偏移距(AVO)和振 幅随角度变化(AVA)等提供更精确的道集和结果。附图说明所包括的附图用来提供对本专利技术的进一步的理解，其构成了说明书的一部分，例 示了本专利技术的优选实施方式，并与文字说明一起用来解释本专利技术的原理，其中图1例示了根据本专利技术一个实施方式的的流程图2例示了根据本专利技术一个实施方式的的偏移步 骤的具体处理；图3是根据一个示例的平层模型偏移试算效果图，其中(a)为平层模型速度场， (b)为单炮记录；图4是基于传统分步傅里叶算子和保幅算子不同入射角度平面波偏移结果，其中 Ca)基于传统分步傅里叶算子，-30度入射平面波偏移；(b)基于传统分步傅里叶算子，O度 入射平面波偏移；(C)基于传统分步傅里叶算子，30度入射平面波偏移；(d)基于保幅算 子，-30度入射平面波偏移；(e)基于保幅算子，O度入射平面波偏移；(f)基于保幅算子，30 度入射平面波偏移；图5是从基于分步傅里叶偏移算子和保幅算子的平面波偏移结果上拾取的层位 成像值对比，其中(a)对应于-30度入射角，(b)对应于O度入射角，(C)对应于30入射角；图6是Marmousi模型速度场以及记录,其中(a)为速度模型，(b)为单炮记录,(C) 为垂直入射平面波合成记录；以及图7是基于分步傅里叶偏移算子和保幅算子的对Marmousi模型的平面波偏移结果，选取了 21个射线参数，从-20度到20度，角度间隔是2度，其中(a)基于分步傅里叶偏移算子，(b)基于保幅算子。具体实施方式下面参照附图对本专利技术的实施方式进行例示性描述。 参照图1，本专利技术的一个实施方式提供了一种，该方法包括预处理步骤(SllO):对原始炮集数据进行预处理，以得到统一基准面的炮集数据；偏移步骤(S120):针对预定平面波入射角度范围中的各个平面波入射角度，基于保幅偏移算子，获得对应于各个平面波入射角度的保幅偏移结果；以及叠加步骤(S130):对所获得的对应于各个平面波入射角度的保幅偏移结果进行叠加，以获得最终的偏移结果。在预处理步骤SllO中，优选地，还可以对所得到的炮集数据进行规则化，使得每炮具有相同道数。下面，参照图2，对偏移步骤S120进行例示性描述。在偏移步骤S120中，针对所述预定平面波入射角度范围(例如，-90度到90度)中的各个平面波入射 角度，执行以下处理合成处理(S210):在地表构建各平面波入射角度的平面波合成算子，应用该平面波合成算子将多个单炮震源和单炮记录合成为平面波震源和平面波记录；相位校正处理(S220):对合成的平面震源执行相位校正处理，使其能够满足保幅叠前深度偏移对炮集数据输入的要求；波场延拓处理(S230):从基准面开始，基于保幅波场延拓算子，沿深度域方向对合成的所述平面波震源和所述平面波记录进行波场延拓；以及成像处理(S240):在各个延拓深度层上应用反褶积成像条件成像，直至最大深度层。合成处理S210可以通过以下方式来执行通过将地表合成算子作用于本文档来自技高网...
【技术保护点】
一种保幅平面波叠前深度偏移方法，该方法包括：预处理步骤：对原始炮集数据进行预处理，以得到统一基准面的炮集数据；偏移步骤：针对预定平面波入射角度范围中的各个平面波入射角度，基于保幅偏移算子，获得对应于各个平面波入射角度的保幅偏移结果；以及叠加步骤：对所获得的对应于各个平面波入射角度的保幅偏移结果进行叠加，以获得最终的偏移结果。
【技术特征摘要】
【专利技术属性】
技术研发人员：叶月明，庄锡进，胡冰，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。