波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口(MPI)实现了弹性波方程的并行求解,极大地提高了数值求解弹性波方程的计算效率.
Wave equation forward modelling plays a significant role in seismic data acquisition,processing,interpretation and inversion.At present,however,forward modelling algorithms based on the seismic wave-equation are restrained in large-scale industrial application because of their massive computational cost.This article studies the basic methods for solving the 2D elastic wave-equation with finite difference and parallel algorithms,and gives the computation domain division method and communication scheme that are suitable for the parallel algorithm.Also,the article analyzes the difference of the computation time cost,the speedup ratio and efficiency of the parallel programs under different parameters.Numerical experiments realized the parallel solution of elastic wave-equation with Message Passing Interface(MPI),and show that the introduction of MPI greatly improves the computation efficiency.