基于GPU%2fCPU叠前逆时偏移研究与适应性的分析.pdfVIP

130 SPG/SEG 深圳 2011 国际地球物理会议论文集 基于 GPU/CPU 叠前逆时偏移研究 及适应性分析 ① ①② ① 王宇超 刘文卿 王西文 （①中国石油勘探开发研究院西北分院；②成都理工大学“油气藏地质及开发工程”国家重点实验室） 摘要 基于双程波的逆时偏移是当前精度最高的深度域成像方法,可以对全波场逆时外推，使高角度反 射界面的反射波，甚至超过 90o反射界面的反射波进行精确成像。通过研究，基于有限差分方法提出了 精确且高效的差分格式逆时波场外推算法。针对多次反射波严重和构造十分复杂时，由互相关成像条 件引入的低频噪音通过优化Laplace算子进行压制。以CPU/GPU作为计算核心，将计算量最大的波场逆 时外推通过GPU实现。并利用随机速度边界的思想提高波场外推算法的并行性，解决了大规模存储的I/O 问题。通过NWGI模型试算验证算法及方法的正确性，分析了几种成像方法存在的问题及其影响因素， 对比研究了当前流行的几种叠前成像方法及其适用性，最后指出了逆时偏移方法未来的发展趋势和研 究重点。 关键词 逆时偏移 波动方程 成像条件 GPU/CPU 协同计算 随机边界 1 引言 常用的叠前深度偏移方法主要有单程波波 方程逆时延拓算法的研究主要包括波场延拓算 动方程偏移和克希霍夫积分偏移。单程波波动 子构造、数值频散压制和边界条件构建等内容。 方程偏移是基于双程波方程的单向波分解。只 有在速度变化不大的情况下成立。单程波方程 波动方程逆时外推算法主要是基于有限差分方 在描述大角度传播时，相位和振幅都存在问题， 法求解波动方程的算法。逆时偏移的瓶颈是计 导致成像误差较大，不能对陡倾角反射精确成 算量大，在逆时偏移实现策略方面，国内外学 像。此外，单程波方程不能对回转波等特殊波 [5] 者做了大量研究。2007 年，Symes 提出在 场进行成像。克希霍夫积分偏移能适应陡倾角 成像。但克希霍夫积分偏移不能处理单成像点 Checking-Points上记录两个时间层的波场。节 多值走时的问题。另外，克希霍夫积分法只能 省存储外推的震源波场的硬盘空间。2009 年同 描述波在光滑介质中的传播。因此，单程波波 济大学王华忠老师、中科院李博博士基于GPU 动方程和克希霍夫积分偏移都无法对复杂地下 [6~7] 开展了逆时偏移研究 。应用随机边界思想， 构造特别是盐下构造精确成像，地球物理勘探 家将目光重新转向逆时偏移。 使波场外推过程不用考虑吸收边界条件，在 1978 年Hereonl 阐述了有限差分法求解 GPU/CPU协同平台上有很高的加速比。本文由 波动方程实现逆时偏移的思想 [1] ，1982 年 [2] 声波方程的高阶有限差分近似入手，在前人研 Whitemore 真正提出逆时偏移实现方法。此 后，Baysal、McMechan等[3，4]将其应用于叠后 究基础上基于GPU/CPU异构平台实现了逆时 偏移，因存在艰难的数值计算效率和存储等问 偏移，分析了高阶有限差分格式的稳定性与频