GPU加速的二维地震波场模拟研究.PDFVIP

第 37卷 第 9期 计 算 机 工 程 2011年 5月 V0l1I37 No．9 ComputerEngineering M av 2011 · 博士论文 · 文章编号：10Oo_．3428(20n)09_一l009-__03 文献标识码：A 中图分类号：TP391 GPU加速的二维地震波场模拟研 究 王占刚，苑春方 (中国矿业大学(北京)煤炭资源与安全开采国家重点实验室，北京 100083) 摘 要：采用交错网格有限差分方法模拟二维地震弹性／粘弹性波场要花费大量的计算时问，为此利用GPU并行处理特点和绘制管道 ，将 计算区域划分为内部区域和PML边界处理区域，整个计算过程 由顶点编程和片段编程处理，采用FBO技术实现差分迭代结果的纹理转换。 实验结果表明，与CPU实现相比，GPU方法提高了模拟效率 ，并且随着网格规模的增加，其效率不断提升，可以实现大规模的高效模拟。 关健诃：有限差分；图形处理器 ；粘弹介质；波场传播；二维 Research0nGPUAccelerated2D SeismicW aveM odelling WANG Zhan·gang，YUANChun-fang (StateKeyLaboratoryofCoalResourcesandMineSa~ty，ChinaUniversityofMiningandTechnology，Beijing100083，China) [Abstract]ThestaggeredgridfinitedifferencemethodisinefficientforseismicwavemodellinginTwo—Dimension(2D)elastic／viscoelasticmedia onasinglePC．Toimprovecomputationspeedup，aGraphicsProcessingUnit(GPU)acceleratedmethodproposedformodem GPU offersan excellentcost—to—performance—ratioparallelism ．ThegeophysicalmodelisdecomposedintosubdomainsforPMLabsorbingconditions．Thevegex andfragmentprocessingareusedtosolveFD schemesandtheupdatedframesareswappedinFramebufferObject(FBO)attachmentsforthenext simulationstep．Simulationprogram runningonmodernPC providessignificantspeedupoveraCPU implementation，whichmakesitpossibleto simulaterealtimecomplexseismicpropagationinhighresolutiongridsizesonlow。COStPC． IKeywords]finitedifference；GraphicsProcessingUnit(GPU)；viscoelasticmedia；wavepropagation；Two—Dimension(2D) DoI：10．3969／j．issn．1000—3428．2011．09．004 1 概述 交错网格有限差分方法是 目前研究二维弹性 ／粘弹性波 孥：妻c++ 场模拟最主流方法之一。但是当模型精度高、网格规模很大 时，有限差分方法要花费大量的计算时问，因此，研究复杂 dt = 石等cdx+dy一2罢誓dy+ 地质模型高精度波场模拟往往需要机群实现” 。 这些方法采用显示差分格式，通过合理地计算区域划分 dt =刀 等cdx+dy一z号t，