摘要
波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口(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.
出处
《地球物理学进展》
CSCD
北大核心
2010年第2期650-656,共7页
Progress in Geophysics
基金
国家自然科学基金资助项目(40804021)
国家自然科学重点基金资助项目(40839901)
国家重点基础研究发展计划(973)资助项目(2009CB219603)资助
关键词
地震波方程
正演
有限差分
并行计算
消息传递接口
seismic wave equation
forward modelling
finite difference
parallel computing
Message Passing Interface
