Abstract: 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.

Key words: finite difference, Graphics Processing Unit(GPU), viscoelastic media, wave propagation, Two-Dimension(2D)

摘要： 采用交错网格有限差分方法模拟二维地震弹性/粘弹性波场要花费大量的计算时间，为此利用GPU并行处理特点和绘制管道，将计算区域划分为内部区域和PML边界处理区域，整个计算过程由顶点编程和片段编程处理，采用FBO技术实现差分迭代结果的纹理转换。实验结果表明，与CPU实现相比，GPU方法提高了模拟效率，并且随着网格规模的增加，其效率不断提升，可以实现大规模的高效模拟。

关键词: 有限差分, 图形处理器, 粘弹介质, 波场传播, 二维

CLC Number: 

• TP391