GPU加速的二维地震波场模拟研究
|
摘要
采用交错网格有限差分方法模拟二维地震弹性/粘弹性波场要花费大量的计算时间,为此利用GPU并行处理特点和绘制管道,将计算区域划分为内部区域和PML边界处理区域,整个计算过程由顶点编程和片段编程处理,采用FBO技术实现差分迭代结果的纹理转换。实验结果表明,与CPU实现相比,GPU方法提高了模拟效率,并且随着网格规模的增加,其效率不断提升,可以实现大规模的高效模拟。
The staggered-grid finite difference method is inefficient for seismic wave modelling in Two-Dimension(2D) elastic/viscoelastic media on a single PC.To improve computation speedup,a Graphics Processing Unit(GPU) accelerated method proposed for modern GPU offers an excellent cost-to-performance-ratio parallelism.The geophysical model is decomposed into subdomains for PML absorbing conditions.The vertex and fragment processing are used to solve FD schemes and the updated frames are swapped in Framebuffer Object(FBO) attachments for the next simulation step.Simulation program running on modern PC provides significant speedup over a CPU implementation,which makes it possible to simulate realtime complex seismic propagation in high resolution grid sizes on low-cost PC.
The staggered-grid finite difference method is inefficient for seismic wave modelling in Two-Dimension(2D) elastic/viscoelastic media on a single PC.To improve computation speedup,a Graphics Processing Unit(GPU) accelerated method proposed for modern GPU offers an excellent cost-to-performance-ratio parallelism.The geophysical model is decomposed into subdomains for PML absorbing conditions.The vertex and fragment processing are used to solve FD schemes and the updated frames are swapped in Framebuffer Object(FBO) attachments for the next simulation step.Simulation program running on modern PC provides significant speedup over a CPU implementation,which makes it possible to simulate realtime complex seismic propagation in high resolution grid sizes on low-cost PC.
引文
[1]Bohlen T.Parallel 3-D Visco-elastic Finite Difference Seismic Modeling[J].Computers&Geosciences,2002,28(8):887-899.
[2]Sun Weitao,Shu Jiwu,Zheng Weimin.Parallel Seismic Propagation Simulation in Anisotropic Media by Irregular Grids Finite Difference Method on PC Cluster[C]//Proc.of International Conference for Computational Science and Its Applications.Singapore:[s.n.],2005:762-771.
[3]Kelly K R.Numerical Study of Love Wave Propagation[J].Geophysics,1983,48(7):833-853.
[4]Virieux J P.P-SV Wave Propagation in Heterogeneous Media:Velocity-stress Finite-difference Method[J].Geophysics,1986,51(4):889-901.
[5]Graves R W.Simulating Seismic Wave Propagation in 3D Elastic Media Using Staggered Grid Finite Differences[J].Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
[6]Carcione J M,Kosloff D,Kosloff R.Wave Propagation Simulation in a Viscoelastic Medium[J].Geophysical Journal of the Royal Astronomical Society,1988,93(2):597-611.
[7]Robertsson J O,Blanch J O,Symes W W.Viscoelastic Finite-difference Modeling[J].Geophysics,1994,59(9):1444-1456.
[8]Robertsson J A.A Numerical Free-surface Condition for Elastic/Viscoelastic Finite-difference Modeling in the Presence of Topography[J].Geophysics,1996,61(6):1444-1456.
[9]Luebke O,Harris G,Lefohn K A,et al.A Survey of General-purpose Computation on Graphics Hardware[C]//Proc.of Annual Conference of the European Association for Computer Graphics.Dublin,UK:[s.n.],2005:21-51.
[10]卢俊,张保明,黄薇,等.基于GPU的遥感影像数据融合IHS变换算法[J].计算机工程,2009,35(7):261-263.
[11]Mark H.Mapping Computational Concepts to CPUs[M].[S.l.]:Addison-Wesley,2005.
[12]Berenger J P.A Perfectly Matched Layer for The Absorption of Electromagnetic Waves[J].Journal of Computational Physics,1994,114(2):185-200.
[2]Sun Weitao,Shu Jiwu,Zheng Weimin.Parallel Seismic Propagation Simulation in Anisotropic Media by Irregular Grids Finite Difference Method on PC Cluster[C]//Proc.of International Conference for Computational Science and Its Applications.Singapore:[s.n.],2005:762-771.
[3]Kelly K R.Numerical Study of Love Wave Propagation[J].Geophysics,1983,48(7):833-853.
[4]Virieux J P.P-SV Wave Propagation in Heterogeneous Media:Velocity-stress Finite-difference Method[J].Geophysics,1986,51(4):889-901.
[5]Graves R W.Simulating Seismic Wave Propagation in 3D Elastic Media Using Staggered Grid Finite Differences[J].Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
[6]Carcione J M,Kosloff D,Kosloff R.Wave Propagation Simulation in a Viscoelastic Medium[J].Geophysical Journal of the Royal Astronomical Society,1988,93(2):597-611.
[7]Robertsson J O,Blanch J O,Symes W W.Viscoelastic Finite-difference Modeling[J].Geophysics,1994,59(9):1444-1456.
[8]Robertsson J A.A Numerical Free-surface Condition for Elastic/Viscoelastic Finite-difference Modeling in the Presence of Topography[J].Geophysics,1996,61(6):1444-1456.
[9]Luebke O,Harris G,Lefohn K A,et al.A Survey of General-purpose Computation on Graphics Hardware[C]//Proc.of Annual Conference of the European Association for Computer Graphics.Dublin,UK:[s.n.],2005:21-51.
[10]卢俊,张保明,黄薇,等.基于GPU的遥感影像数据融合IHS变换算法[J].计算机工程,2009,35(7):261-263.
[11]Mark H.Mapping Computational Concepts to CPUs[M].[S.l.]:Addison-Wesley,2005.
[12]Berenger J P.A Perfectly Matched Layer for The Absorption of Electromagnetic Waves[J].Journal of Computational Physics,1994,114(2):185-200.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心