起伏地表叠前逆时深度偏移与并行实现 - 石油地球物理勘探.pdf

年 月 第 卷 第 期 ２０１２ ４ ４７ ２ 　 　 　 ·处理技术 · 文章编号： （ ） １０００７２１０２０１２０２０２６６０６ 起伏地表叠前逆时深度偏移与并行实现  王延光 匡 斌 　 　 （胜利油田分公司物探研究院，山东东营 ２５７０２２） 王延光，匡斌 起伏地表叠前逆时深度偏移与并行实现 石油地球物理勘探， ， （）： ． ． ２０１２４７２ ２６６ ２７１ ～ 摘要 本文系统讨论了逆时偏移中常用的波动方程求解方法，并简述了其各 自的特点和优势。文中提出用逆 　 时偏移实现基于起伏地表的波动方程叠前深度偏移，理论模型数据的试算结果与单程波方程和 Ｋｉｒｃｈｈｏｆｆ积分 方法结果进行对比，展示了逆时偏移的优越性。自行研发的基于 ／ 平台的并行软件，推动了逆时偏移 ＣＰＵ ＧＰＵ 算法实用化的运行，在应用中取得了良好效果。 关键词 逆时偏移 起伏地表 图形处理器（ ） ＧＰＵ 　 　 　 中图分类号： 文献标识码： Ｐ６３１ Ａ 　　 解的单程波方程是一个拟微分方程，对它的求解需 １ 引言 要进行有理化近似，如运用波动方程的旁轴近似， 　 从一个深度到另一个深度，在频率域外推波场。这 现有的偏移方法可以划分为两个主要的类型： 些方法虽然能对接近波场外推方向（通常为垂直方 基于射线的Ｋｉｒｃｈｈｏｆｆ积分方法和基于波动方程的 向）的波很好地成像，但对于大角度波场传播，如接 方法。 世纪 年代末期， 积分方法在 近和超过 倾角传播的地震波不能有效地成像， ２０ ９０ Ｋｉｒｃｈｈｏｆｆ ９０° 地震成像处理中占据主导地位。该方法适用于陡倾 从而使单程波方程偏移算法对陡倾角的构造成像 角和复杂地质构造的成像，并且运算速度快、灵活性 受到限制。 强，可以很容易地实现基于地表的叠前深度偏移。 全程波偏移算法通过直接求解原始的全声波方 然而，由于其旅行时计算时采用了对波动方程的高 程来实现，能模拟所有方向传播的地震波，包括反射 频近似，使得 Ｋｉｒｃｈｈｏｆｆ方法不能有效地对复杂地 波和透射波。也就是说，全程波方程的逆时偏移算 质构造成像。传统的积分偏移方法仅考虑单一的射 法不受倾角限制，甚至能对复杂构造和柱面波成像。