弹性波叠前逆时深度偏移与GPU加速.pdfVIP

中国地球物理2011 ·735· (20)地震波传播与成像探查 弹性波叠前逆时深度偏移与GPU加速 Elastic reversetime and acceleration GPU pre-stack depthmigration 马庆珍’李振春黄建平 Ma Zhenchun Qingzhen’Li HuangJianping 中国石油大学(华东)地球资源与信息学院青岛266555 近年来，多波多分量地震勘探由于能够比纵波勘探得到更多的地下介质信息，有效减少地震勘探的多 解性，因而受到了广泛的关注。弹性波叠前逆时深度偏移是当前最为精确的一种多波多分量资料处理方法， 相对于其他方法，如纵、横波分离结合常规的纵波资料处理的方法，由于它考虑了的波的矢量特征，成像 噪音少且可以同时对转换波等非常规波动进行成像。但弹性波叠前逆时深度偏移庞大的计算量阻碍了其应 用推广。本文利用GPU对弹性波逆时偏移进行加速，有效地提高了运算效率，降低了偏移成本。 1．弹性波叠前逆时偏移原理 一般来说，逆时偏移主要由三部分组成：①成像条件的求取，②接收波场的逆时延拓，③成像条件的应 用。弹性矢量波叠前逆时深度偏移采用矢量波动方程对上下行波场进行延拓，即通过有限差分等数值方法对 矢量波动方程进行数值离散，得到波场延拓算子。其中下行波场(震源波场)的正向延拓即成像条件的求取 等价于一次弹性矢量波数值模拟，而上行波场(检波波场)的反向延拓即逆时延拓可以看做数值模拟的逆过 程。成像条件的应用是影响逆时深度偏移效果的核心，本文采用现在应用最为广泛的零延迟互相关成像条件。 逆时偏移低频噪音的去除是逆时偏移后期需要处理的一个主要问题，对比现有的各种滤波方法，本文 采用的拉普拉斯滤波在去除低频假象的同时能够很好地保留有效信号，且实现较为简便。另外，震源归一 化的零延迟互相关成像条件能够有效地提高有效信号的能量，尤其是深层信号能量的增强效果明显。 2．GPU高性能运算 GPu高性能运算从2003年提出到现在发展的已经比较完善。一系列GPU高性能运算设备的推出以及 加速在各个学科领域都取得了很好的效果。在地震波传播与成像方面，李振春等人利用GPU对弹性波有 限差分数值模拟进行了加速，刘红伟(2010)等人应用GPU对声波叠前逆时偏移进行了实现…。 3．程序实现及加速效果分析 通过对GPU高性能运算的分析研究，本文将弹性矢量波叠前逆时深度偏移算法进行了GPU上的编程 实现。由于矢量波波场延拓算子本身具有良好的并行性，因此可在整个二维模型区域完全并行展开进行计 算。由于逆时偏移中互相关的两个波场延拓方向不同，因此常规方法中通过现正向延拓并保存正演波场， 应用成像条件时再读取的办法来解决。为避免存取波场带来的巨量存储和读取延迟问题，本文采用 Eric(2008)中的策略四来解决逆时偏移的存储问题，即通过计算换存储【2】。 为验证程序的精确性和测试加速效果，本文设计了一系列模型进行测试，并同已有的CPU程序进行 对比分析。通过模型测试及对比分析可以发现：首先，利用GPU进行弹性矢量波叠前逆时深度偏移加速 运算对各种构造都能够取得同CPU程序同样的成像效果。其次，利用GPU对偏移算法进行加速运算，无 需复杂的算法设计和优化即可得到相当可观的加速效果，对降低弹性矢量波叠前逆时偏移成本，提高弹性 矢量波叠前逆时偏移在工业上的应用将会有很好的推动作用。 学引进人才项目(Y090104)支持资助。 参考文献 【1】刘红伟，等．地震叠前逆时偏移商阶有限差分算法及GPU实现【J】．地球物理学报，2010，53(7)：1725一i733 Dussaud，et forreverse-time Annual 【2】Eric a1．Computationalstrategies migration[C]．SEGMeeting，2008．