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

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

Research on GPU Accelerated 2D Seismic Wave Modelling

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.