基于GPU三维精细积分法正演模拟.pdfVIP

·484· 中国地球物理2011 (15)信息技术与地球物理 基于GPU的三维精细积分法正演模拟 ofthe3D forseismic method Applicationpreciseintegration basedon modeling GPU 段玉婷‘李靖宇 胡天跃 Duan Li Hu YutingJingyuTianyue 北京大学地球与空间科学学地球与空间科学学院地球物理系 北京100871 1．引言 地震波正演模拟技术是模拟复杂近地表低速带和地质构造深部精细结构的有效方法。随着油气资源的 勘探程度不断提高，尤其是陆上油气勘探到达的深度和勘探工作要求的精细度也逐步加大，三维勘探技术 也趋向成熟。三维勘探带来地球物理数据规模的急剧增加，庞大的计算规模会导致巨大的计算量，因而计 算的效率是必须要考虑的因素之一。本文针对最新发展的图形处理器(GPU)并行计算技术，对精细积分 方法的实现过程进行优化，并利用该方法实现三维地震正演模拟。 2．三维精细积分法原理 精细积分法的基本思想是对波动方程的空间域用任意差分方式离散，形成一组关于时间的差分方程， 然后再对时间域用子域精细积分的方法解出该差分方程组。对于三维声波方程可以表示为 V。HV+f 位移矩阵。采用有限差分网格离散后，进行时间域积分，得到 V”“一TV“+(T—I)H-1f“ 其中，T(At)=e胁。在解出V后，取其前一半，可以得到波场U的显示差分格式。 在模拟过程中，采用伪6函数法，将爆炸点源模拟成一个以点源为中心的扩展区域。并利用最佳匹配 吸收层(PML)的方法来消除边界反射波。 理论上，精细积分法在在时间域得到的解是精确的，解决了有限差分法频散严重的问题，从而提高了 整个计算格式的精度。本文采用改进型精细积分法，对四阶声波方程进行求解(Tang，2007)。 3．GPU高性能并行计算 计算过程中采用CPU+GPU异构模式，由CPU负责执行复杂的逻辑运算和事务管理等不适合数据并 行的计算，由GPU上的每个计算单元来负责循环中每个网格点上的迭代计算同时进行。这种利用GPU强 大处理能力和高带宽来弥补CPU性能不足的计算方式在发掘计算机潜在性能方面，以及成本控制与性价 比方面具有显著优势，大大的提高了计算效率。 4．三位正演模拟结果 本文采用有斜坡的四层地层模型进行三维正演模拟分析，得到不同平台上正演模拟得到的叠加剖面和 时间偏移剖面，通过与模型剖面的对比验证了计算的正确性。精细积分程序的GPU并行版本和CPU串行 版本对比，加速比达到了40倍。 参考文献 ofa methodinforwardseismic SEG 77th Annual 【1】Tang,G，Hu,T．，Yang,J．Applicationspreciseintegration modeling[C]．The International 130--2134． Abstract,2007，2 Meeting：Expanded seismicwave within 【21 complex Yang,J．，Liu,T．，Tang,Ct，Hu,T．Modelingpropagation structure