基于gpu的逆时偏移方法在地震数据处理中的应用.pdfVIP

基于 GPU 的逆时偏移方法在地震数据处理中的应用 邢怡* 天津地质调查中心，天津 300170 摘要：采用双程波求解微分波动方程的逆时偏移算法，从理论上弥补了常规地震偏移的成像缺陷，但是存在运算时间长、储存 量大和低频假象等缺点。本文利用图形处理器（GPU）加速技术、随机边界算法等方法对该算法进行改进，有效克服了上述缺点， 使该方法能快速有效地应用于叠前深度偏移处理。用基于 GPU的逆时偏移算法对实际地震数据进行叠前深度偏移处理，并与 Kirchhoff算法进行比较，结果表明：此算法用时短、噪声小、轴向清晰，具有较好的成像效果。 关键字：逆时偏移；叠前深度偏移；GPU加速；随机边界算法 1、引言 [1] 逆 时 偏 移 的 思 想 最 早 由 Hemon 于 1978 年 提 出 的 ， 直 到 1983 年 通 过 Baysal、Whitmore、Loewenthal、McMechan等几位科学家分别将之用于纵波资料的叠后偏移处理中，才真正将逆 时偏移引入地震勘探这一领域。与此同时，Gazdag等还对纵波资料的叠前逆时偏移方法进行了研究，并处理 实际资料。 逆时偏移方法的基本核心思想是通过双程波波动方程在时间域上求解波动方程（对人工给予的震源子波 正向传播和接收到的地震资料进行反向传播），然后利用成像条件实现偏移。在此处理过程中使用了双程波波 动方程，避免了上下行波的分离，是一种以矢量波理论为基础的深度域偏移方法，它对地表记录到的波场在 时间方向上进行逆时延拓，然后利用成像条件对地下各点进行成像。这种方法允许波向各个方向传播，没有倾 角限制，精确度较高，并能使回转波和多次波较好的成像，能够解决复杂的地质成像问题，所以对该方法的 研究和改进具有很重要的实际应用价值。但是这样需要先利用人工震源进行正向传播，然后保存其波场，再利 用地震资料进行反向传播，反向传播的同时读出正向传播的波场，并利用成像条件成像，这就导致了处理时 需要巨大的存储量。此外，地震波场模拟还会产生巨大的计算量。这些导致逆时偏移处理存在着最致命的缺点 就是计算效率低。为使其能更好地应用于实际生产中，就要寻找新的技术将逆时偏移方法进行改进。 本文从地震波场分析入手，为提高逆时偏移方法效率、 小储存量和克服低频假象缺陷，提出 相应的解 决对策。 2、算法原理介绍 [2] 对于地震波场模拟所需要的巨大计算量，本文利用最新的GPU加速技术 予以解决。对于存储的问题，利 [3] 用Robert G.Clapp 提出的随机边界方法解决。 ⑴、应用GPU技术对逆时偏移算法的实现 在GPU上实现逆时偏移方法的成像步骤是： 首先，要完成一次波动方程正演计算，并保存每一时间步的震源波场信息。然后，逆时外推记录波场，保 存每一时间步的记录波场信息。最后，分别读取保存的同一时刻的震源波场和记录波场，做互相关成像。具体 算法如下： ①计算人工震源的正向传播（从0 到T），保存波场。 1 2 stt − ∆s = δ (x − a)f (t) (1) v s a f (t) 其中 是待求解的位移方程，是空间和时间的函数； 表示震源的位置； 为合适的子波。 ②计算逆时波场成像(从T 到0)。  1 r − ∆r = ∑δ (x − r )g (t)  2 tt i i v i 