CAJ | 학술논문

波动方程正演在地震资料采集、处理、解释与反演中均能发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产.本文从二维弹性波动方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法与通信方案,分析了不同参数条件下并行程序的运行时间、加速比与效率.引入消息传递接口(MPI)实现了弹性波方程的并行求解,极大地提高了数值求解弹性波方程的计算效率.
파동방정정연재지진자료채집、처리、해석여반연중균능발휘중요작용,단현유적기우구해지진파동방정적정연산법유우수방대계산량적제약이난우대규모응용우공업생산.본문종이유탄성파동방정출발,연구료이용유한차분법병행구해해방정적기본사로여방법,급출료괄우병행구해적계산공간화분방법여통신방안,분석료불동삼수조건하병행정서적운행시간、가속비여효솔.인입소식전체접구(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.