基于图形处理器的伪谱和高阶有限差分混合方法地震波数值模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
伪谱和高阶有限差分混合方法,在垂直方向采用交错网格有限差分算子,利用其并行程度高的特点,在水平方向采用伪谱算子,保留其高精度的优势,是计算地震波场的有效方法.图形处理器(graphic processing unit,简写为GPU)由于其高度并行性,在计算此类问题中有显著的优势.由英伟达(NVIDIA)公司推出的统一计算设备架构(compute unified device architecture,简写为CUDA)平台极大地简化了GPU编程的难度.为提高计算效率,本文实现了基于CUDA平台的混合方法二维地震波场模拟.然后基于二维均匀介质模型将CPU与GPU版本的运行时间进行对比.实际测试结果表明,基于CUDA的并行模拟方法在保证计算精度的同时显著地提高了计算速度,为开展大规模非均匀地球介质地震波传播数值模拟提供了一种可选的方法.
Hybrid PSM/FDM method,approximating the horizontal spatial derivative with pseudospectral method(PSM)while using finite difference method(FDM)to approximate the vertical spatial derivative,is an efficient technique in seismic wave simulation because it makes usage of both the high level of parallelism in FDM and the great accuracy in PSM.With natural parallelism,graphic processing unit(GPU)is very suitable for calculating such kind of problem,and compute unified device architecture(CUDA)drastically simplifies GPU programming.In order to improve the computing efficiency,in this paper,we implemented the 2D hybrid method for seismic wave simulation based on CUDA.The GPU-based hybrid method brings about significant speedup according to numerical tests performed in this paper and proves to be a reliable way for simulating large scale seismic wave propagation.
引文
龙桂华,李小凡,江东辉.2010.基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案[J].地球物理学报,53(12):2964--2971.Long G H,Li X F,Jiang D H.2010.Accelerating seismic modeling with staggered-grid Fourier pseudo-spectral differentiation matrix operator method on graphics processing unit[J].Chinese Journal of Geophysics,53(12):2964--2971(in Chinese).
    秦艳芳,王彦宾.2012.地震波传播的三维伪谱和高阶有限差分混合方法并行模拟[J].地震学报,34(2):147--156.Qin Y F,Wang Y B.2012.Three-dimensional parallel hybrid PSM/FDM simulation of seismic wave propagation[J].Acta Seismologica Sinica,34(2):147--156(in Chinese).
    单蕊,朱伟.2013.基于GPU并行的谱元法地震波数值模拟[J].中国煤炭地质,25(3):52-54,62.Shan R,Zhu W.2013.Parallel SEM seismic wave simulation based on GPU[J].Coal Geology of China,25(3):52--54,62(in Chinese).
    王德利,雍运动,韩立国,廉玉广.2007.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报,29(1):30--34.Wang D L,Yong Y D,Han L G,Lian Y G.2007.Parallel simulation of finite difference for seismic wavefield modelingin 3-D viscoelastic media[J].Northwestern Seismological Journal,29(1):30--34(in Chinese).
    王月英,孙成禹.2006.弹性波动方程数值解的有限元并行算法[J].中国石油大学学报:自然科学版,30(5):27--30.Wang Y Y,Sun C Y.2006.Finite element parallel algorithm for numerical solution of elastic wave equation[J].Journalof China University of Petroleum,30(5):27--30(in Chinese).
    魏星,王彦宾,陈晓非.2010.模拟地震波场的伪谱和高阶有限差分混合方法[J].地震学报,32(4):392--400.Wei X,Wang Y B,Chen X F.2010.Hybrid PSM/FDM method for seismic wavefield simulation[J].Acta SeismologicaSinica,32(4):392--400(in Chinese).
    严九鹏,王彦宾.2008.重叠区域伪谱法计算非均匀地球介质地震波传播[J].地震学报,30(1):47--58.Yan J P,Wang Y B.2008.Modeling seismic wave propagation in heterogeneous medium using overlap domainpseudospectral method[J].Acta Seismologica Sinica,30(1):47--58(in Chinese).
    严武建,王彦宾,石玉成.2013.基于伪谱和有限差分混合方法的兰州盆地强地面运动二维数值模拟[J].地震工程学报,35(2):302--310.Yan W J,Wang Y B,Shi Y C.2013.2Dsimulation of the strong ground motion in Lanzhou basin with hybrid PSM/FDM method[J].China Earthquake Engineering Journal,35(2):302--310(in Chinese).
    张明财,熊章强,张大洲.2013.基于MPI的三维瑞雷面波有限差分并行模拟[J].石油物探,52(4):354--362.Zhang M C,Xiong Z Q,Zhang D Z.2013.3Dfinite difference parallel simulation of Rayleigh wave based on messagepassing interface[J].Geophysical Prospecting for Petroleum,52(4):354--362(in Chinese).
    郑亮,张怀,尹凤玲,石耀霖.2013.基于GPU的大规模三维地震波传播有限元程序实现[C]∥中国地球物理2013:第七分会场论文集.昆明,2013年10月13日.Zheng L,Zhang H,Yin F L,Shi Y L.2013.Finite element program implementation of large scale seismic wave propagation modeling based on GPU[C]∥Chinese Geophysics 2013:Proceedings of the 7th Branch.Kunming,October 13,2013.
    Aoi S,Nishizawa N,Aoki T.2012.Large scale simulation of seismic wave propagation using GPGPU[C]∥15th WorldConference on Earthquake Engineering.Lisbon,Portugal,September 24to 28,2012.
    Cai C,Chen H Q,Deng Z,Chen D,Khan S U,Zeng K,Wu M X.2012.GPGPU-aided 3Dstaggered-grid finite-difference seismic wave modeling[C/OL]∥12th International Conference on Scalable Computing and Communications(ScalCom).[2014-05-12].http:∥newtheme.hgpu.org/p=8610.
    Cerjan C,Kosloff D,Kosloff R,Reshef M.1985.A nonreflecting boundary condition for discrete acoustic and elasticwave equations[J].Geophysics,50(4):705--708.
    Daudt C R,Braille L W,Nowack R L,Chiang C S.1989.A comparison of finite-difference and Fourier method calculations of synthetic seismograms[J].Bull Seismol Soc Am,79(4):1210--1230.
    Furumura T,Kennett B L N,Takenaka H.1998.Parallel 3-D pseudospectral simulation of seismic wave propagation[J].Geophysics,63(1):279--288.
    Kosloff D,Baysal E.1982.Forward modeling by a Fourier method[J].Geophysics,47(10):1402--1412.
    NVIDIA.2014.CUFFTLibraryUser’s Guide[M/OL].[2014-04-26].http:∥docs.nvidia.com/cuda/cufft/.
    Okamoto T,Takenaka H,Nakamura T,Aoki T.2010.Accelerating large-scale simulation of seismic wave propagationby multi-GPUs and three-dimensional domain decomposition[J].Earth Planets Space,62(12):939--942.
    Okamoto T,Takenaka H,Nakamura T,Aoki T.2012.Large-scale simulation of seismic wave propagation of the 2011Tohoku-Oki M9earthquake[C]∥Proceedings of the International Symposium on Engineering Lessons Learnedfrom the 2011 Great East Japan Earthquake.Tokyo:Japan Association for Earthquake Engineering:349--360.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心