《单程和双程波动方程叠前深度偏移方法》.pdfVIP

东 北 石 油 大 学 学 报 第38卷 第 4期 2014年 8月 JOURNALOFNORTHEAST PETROlEUM UNIVERSITY Vo1．38 No．4 Aug． 2014 单程和双程波动方程叠前深度偏移方法 田东升 ，王云专 ，李义鹏 ，石 颖 ，柯 璇 ，李婷婷 ，刘淑芬 (1．东北石油大学 地球科学学院，黑龙江 大庆 163318； 2．北华航天工业学院 电子工程系 ，河北 廊坊 065000) 摘 要：叠前深度偏移是获得地下构造映像的有效手段 ，而基于波动方程的叠前深度偏移方法对速度横 向变化剧烈 的地层有更好的适应性．分析基于单程波方程的相移法、相移加插值法、频率一空间域有限差分法、傅里叶有限差分法和 基于双程波方程 的逆时偏移方法 ，借助于地堑模型与盐丘模型 ，测试 5种逆 时偏移方法成像复杂构造 的精度和适应性． 结果表明，基于波动方程的叠前深度偏移方法可实现横向变速地下构造成像 ，相 比于基于双程波方程 的逆时偏移方法 ， 单程波方程方法对垂直断层等高陡倾角构造成像有局限性；逆时偏移方法对垂直断层、盐丘下边界等复杂构造可 以清晰 成像 ，辅 以精确 的地层速度 ，逆时偏移方法在地震资料成像领域中有广阔的发展和应用前景． 关 键 词 ：叠前深度偏移 ；相移法；相移加插值 ；频率～空间域有限差分 ；傅里叶有限差分 ；逆时偏移 中图分类号 ：TE132．1 文献标识码 ：A 文章编号：2095—4107(2014)04—0039—06 0 引言 随着油气勘探 目标 日趋复杂，具体表现在断块小、倾角陡、纵横向速度变化剧烈等方面，常规偏移方法 很难达到地震资料处理要求，加之巨大的数据处理量，亟需研究高精度、高效率的偏移算法，为精细地质构 造解释及储层识别提供重要依据．叠前深度偏移是获得精确地下构造 的有效途径．在数学解法上 ，叠前深 度偏移分为两类 ：基于波动方程积分解 的射线偏移和基于波动方程微分解的波动方程叠前深度偏移．波动 方程叠前深度偏移较射线偏移在处理横向变速问题上具有更强的适应性 ，不存在射线偏移法中成像点的 多值走时问题．其中波动方程叠前深度偏移分为单程波法和双程波法． 2O世纪 7O年代 ，ClaerboutJF_l】提 出应用波动方程进行偏移，采用有限差分法求解得到单程波动方 程 15。近似公式 ，该方法在主传播方向小范围内具有较好的效果 ，对宽方位地震波传播模拟并不理想 ，尤 其对于地震数据 中包含的水平和陡倾角反射信息的偏移效果不明显．GazdagJ最初提出的频率一波数域 相移法 仅适用于对地下横 向速度不发生变化 的介质进行成像 ；进而提出的相移加插值方法E。可以适应 存在横向速度变化 的介质 ，但需要频繁计算参考波场 ，效率较低Ⅲ．随后，人们将有限差分法与傅里叶法结 合研究，提出混合偏移方法[5]，对横向速度变化强烈的介质成像效果 良好．王玉学等 推导上行波方程 的两种高阶近似表达式．冯凤萍等_8]研究加吸收层的三维45。上行波方程的隐式差分格式．然而 ，所有单程 波偏移方法存在局限性，即单程波算子在成像大角度传播的波时将发生相位改变和振幅减弱的现象，无法 对陡倾角进行成像 ；另外 ，单程波法也无法成像 回转波． 采用双程波动方程进行偏移能很好地适应剧烈的横 向速度变化 ，可以有效解决复杂地质体成像问题 ， 最典型的方法就是逆时偏移方法．该方法最早 由WhitmoreDN等 。提出，最初应用于处理叠后资料．逆 时偏移方法简单 、易于实现，不对波动方程做任何近似，从而不存在倾角的限制 ，可以对透射波 、多次波 、绕 射波等进行成像．近些年 ，随着计算机硬件技术的迅猛发展及勘探要求 的 日益提高 ，逆时偏移方法的研究 也从叠后走向叠前口 “]，从二维走 向三维．GPU加速计算技术的引入 和噪音压制策略的研究E 推动 逆时偏移技术的发展．同时，如全波形反演 的精准速度建模方法 的研究 ¨也加速波动方程叠前深度偏移 的研究进程 ． 收稿 日期 ：2013—02—12；编辑 ：任志平 基金项 目：黑龙江省教育厅科学技术研究项 目 作者简介 ：田东升(1989一)，男 ，硕士研究生 ，主要从事地震资料处理方面的研究 通讯作者：石 颖，E—mail：shiyin