基于GPU并行的谱元法地震波数值模拟.pdfVIP

第 25卷 3期 中 国 煤 炭 地 质 Vo1．25No．O3 20l3年 3月 C0AL GEOLOGY OFCHINA MaT． 2O13 doi：10．3969~．issn．1674—1803．2013．03．13 文章编号：1674—1803(2013)3-0052—03 基于GPU并行的谱元法地震波数值模拟 单 蕊 。朱 伟 (1．中国煤炭科工集 团西安研究院，陕西 西安 710077； 2．中国地质大学地球物理与空间信息学院，湖北 武汉 430074) 摘 要：谱元法是有限元法 的一个重要分支，具有对模型适应能力强、计算精度较高和易于实现等优点。有别于有限 元 的近似计算 ，谱元法的质量矩阵采用数值积分得到的对角阵，既简化了计算 ，又可避免计算精度的下降。以用于通 用计算领域 的GPU为例 ，介绍了在 CUDA编程平台下实现谱元法地震波数值模拟并行化的方法。在实例应用中为 避免浮点数据的原子操作 ，提高计算效率 ，模型被采样为串行执行的四个并行单元集合体。 关键词 ：谱元法 ，GPU，地震波数值模拟 中国分类号 ： P631．4 文献识别码 ：A ParallelSEM SeismicW aveSimulatiOilBasedonGPU ShanRui，ZhuWei (1．XianResearchInstituteofChinaCoalTechnology EngineeringGroupCorp，Xian，Shaanxi710077； 2．InstituteofGeophysicsandGeomatics，CUG，Wuhan，Hubei43o074) Abstract：Spectralelementmethod (SEM)isoneofimportantbranchesoffiniteelementmethod．Itcanappliedtocomplexmodels， withhighprecision andeasytorealize．Withtherecentdeveloped generalcomputingdeviceGPU and programmingmodelCUDA， SEM seismic wave simulation isparallelimplemented．In orderto avoid floatatomic operation and to improve efficiency，the geometricalmodelisresampledtofourserialcomputedpartsinwhichelementsareexecutedparalle1． Keywords：spectralelementmethod；GPU；seismicwavesimulation 0 前言 时间轴方向可采用显式差分递推格式[61。 庞大的计算和存储需求限制了有限元法的应 地震波数值模拟是勘探地震学的重要研究 内容 用。有限元程序一般需要大型并行机做为硬件基础。 之一。波场数值模拟的实现方法主要是有限差分法 随着 图形处理器 (GPU)在通用计算领域 的实现 ，计 和有限元法 有限差分方法将波动方程直接离散 ，易 算机的科学计算能力成倍增长。微机工作站配置数 于实现 ，计算和储存量都相对较小。因此 ，有限差分 块 GPU计算卡即能实现大规模并行计算 。本文利用 法在勘探地震领