逆时偏移方法技术进展综述_张春燕探究.pdfVIP

文章编号:1671-8585(2010)05-0309-09 收稿日期:2010-06-13;改回日期:2010-06-28。 第一作者简介:张春燕(1986—), 女 ,中国石油大学(华东)在读硕 士 ,主要从事逆时偏移成像方面的研究工作。 　 逆时偏移方法技术进展综述 张春燕 ,李振春 ,孙小东 (中国石油大学(华东)地球资源与信息学院 ,山东青岛 266555) 摘要:逆时偏移(RTM)利用双程波动方程 , 基于时间逆时外推对波场进行重构 , 理论上无倾角限制 , 可以实现对 回转波 、棱镜波以及多次波的正确成像 , 获得精确的动力学信息 ,具有良好的保幅性。 然而 , 庞大的存储空间需 求和计算量以及固有的低频成像噪声 ,使得 RTM无法得到广泛的应用。 针对上述问题 , 分析 、对比了各种逆时 偏移方法技术及其应用效果 ,测试了偏移去噪的几种方法 , 最后对逆时偏移的发展方向作了展望。 关键词:叠前逆时偏移;计算策略;低频成像噪声;Laplacian滤波;成像条件 中图分类号:P631.4 文献标识码:A 　　随着地震勘探对象日趋复杂 ,适应起伏地表 、 高陡构造 、复杂速度分布和复杂储层的地震成像技 术 ———逆时偏移(RTM)得到了不断发展 。 RTM基 于双程波动方程进行波场延拓 ,避免了传统单程波 偏移中的倾角限制 ,理论上可以实现对回转波 、棱 镜波以及多次反射波的成像 ,获得更为精确的振幅 等动力学信息 ,实现保幅成像[ 1, 2] 。实际应用中 , RTM可以适应介质速度的横向剧烈变化 ,能够对 常规地震成像方法所无法进行准确成像的复杂区 域(如盐丘及强断裂带)进行较准确的成像 。 由于 RTM能够对回转波以及棱镜波进行成像 , 因此可以用其对复杂速度场进行更细化 、更精确的 估计 [ 3 ～ 5] ,提高复杂背景中的速度建模精度[ 6 ～ 9] 。 RTM也可以应用于探地雷达信号处理或隧道超前 预报等 ,为高精度工程地质勘探提供有效途径[ 10] 。 然而 ,庞大的内存需求和计算量以及成像过程 中固有的低频成像噪声制约了 RTM的实用性。很 多学者已就这一难题进行了研究 ,本文在同行研究 的基础上 ,重点分析了提高 RTM计算效率的方法 技术研究进展 ,并就计算硬件的选择作了相应分 析;对低频成像噪声衰减的各种方法 ,例如波场外 推过程中去噪的方法 、应用新的成像条件去噪的方 法以及对叠加成像结果应用滤波去噪的方法等 ,作 了详细的叙述和对比分析 ,展示了相关方法的偏移 去噪效果;对 RTM在实际研究工作中的各种应用 作了介绍 ,并就 RTM的发展方向作了展望 。 1　方法原理 RTM的基本原理如下:从记录波场的最大时 刻起 ,将每个采样时间的炮记录作为边界条件 ,在 时间上反向外推波场至零时刻 ,然后在每个网格点 上提取满足成像条件的成像值。为了消除边界反 射 ,在计算过程中必须引入边界条件。 三维逆时深度偏移的定解问题可描述为: 2u x2 + 2u y2 + 2u z2 = 1 v2 2u t2 u(x, y, z, t)tTmax =0 u(x, y, z, t)z=0 =ψ(x, y, t) (1) 式中:u(x, y, z, t是声波波场 , v(x, y, z为介质速 度 , ψ(x, y, t)为地表接收的三维地震记录。利用 高阶有限差分求解方程(1),即可解决三维逆时深 度偏移成像问题 。 叠前 RTM是对单炮记录数据进行 RTM,在此 过程中首先要正向外推震源波场到最大时刻并将 中间结果保存下来;然后将记录波场数据沿时间反 方向外推至零时刻 ,记录中间结果;最后在每个网 格点应用合适的成像条件 ,将所有炮的偏移结果叠 加得到最终偏移结果 。 2　RTM的计算策略 在单炮记录数据的叠前 RTM计算过程中 ,为 了得到每个成像点的成像值(以互相关成像条件 为例),需要得到该点正向外推波场和逆向外推波 场每个时刻的波场值 ,这样就必须先正向外推震源 波场到最大时刻 , 然后逆向外推记录波场至零时 刻 ,并且每次都要将中间的波场快照保存下来 。因 此 ,整个外推过程需要庞大的计算量和存储量 ,使 得 RTM方法的实用性遇到了很大挑战。 第 33卷第 5期 2010年 10月 勘探地球物理进展 ProgressinExplorationGeophysics Vol.33, No.5 Oct., 2010 2.1　计算方法的选择 Dussaud等提出在标量场(压力场)中基于跳 步时间步进的 RTM,它有几种不同的计算策略 ,即 在空间和时间上采用不同的插值算法来求解每个 时间步长的震源波场和记录波场 ,比如在空间上使 用双线性内插 ,在时间上使用三次内插 ,通过不同 的算法组合