三维交错网格有限差分地震波模拟的GPU集群实现三维交错网格有限差分地震波模拟的GPU集群实现.pdf

第26卷第6期 地球物理学进展 V01．26，No．6 2011年12月(页码：1938--1949)PR(ⅪRESSINGE()PHYSICS Dec．，2011 2903．2011．06．007． doi：10．3969／j．issn．1004 G Y X 3D Finite—difference LongH，ZhaoB，I。i F，eta1．Accelerating SeismicWave onGPUcluster． Staggeredgrid Modeling ProgressinGeophys．(in 三维交错网格有限差分地震波模拟的GPU集群实现 龙桂华1’2…， 赵宇波1， 李小凡2， 高 琴1， 王 周1 (1．中国科学院深圳先进技术研究院，深圳518055；2．中国科学院地质与地球物理研究所，北京100029 3．浙江大学地球科学系，杭州310027) 摘要有限差分实现简单、速度快，作为地震波场模拟一种有效数值方法，被广泛用于正演计算密集的波形反演 和逆时偏移中．三维地震波正演模拟计算量大，一直以来制约着三维叠前逆时偏移和反演的工业化应用，GPU通 用计算技术的产生及其内在的数据并行性有望改变这一现状．本文通过分析三维交错网格有限差分方法在GPU 上的实施，利用片内共享存储器实现了三维地震波数值模拟的高效算法，取得了较单核CPU快79x～108x的加速 比；通过区域分解技术将单GPU上不能计算的地质体模型沿z轴方向进行粗粒度分解，采用消息传递接口交换边 界数据，运用MPI+CUDA的方式实现了大尺度三维地震波场模拟，并着重分析了影响GPU并行计算效率的一些 关键因素．大尺度三维地震波场模拟的加速实现，为促进叠前逆时偏移和波形反演技术的工业化转化提供了可能， 因此具有重要的研究意义． 关键词GPU，交错网格，有限差分，图形处理器，CUDA 2903．2011．06．007 A doi：10．3969／j．issn．1004 中图分类号P315 文献标识码 3D Finite。_differenceSeismic AcceleratingStaggered_。grid Wave onGPUcluster Modeling LONGGui—hual2…，ZHA()Yu GAO WANGZhoul B01，LIXiao—fan2， Qinl， (1．ShenzhenInstitutes 518055，China； ofAdvancedTechnology，ChineseAcademyofSciences，Shenzhen 2．Institute and ofGeologyGeophysics，ChineseAcademyofSciences，Beijing100029，China 310027，China) 3．DepartmentofEarthSciences，ZhejiangUniversity，Hangzhou AbstractAsanefficientnumerical di