消息传递接口在声波方程正演中的应用
|
摘要
声波方程正演在地震资料采集、处理、解释与反演中均发挥重要作用,但现有的基于求解地震波动方程的正演算法由于受庞大计算量的制约而难于大规模应用于工业生产。从声波方程出发,研究了利用有限差分法并行求解该方程的基本思路与方法,给出了适于并行求解的计算空间划分方法,分析了不同参数条件下并行程序的运行时间、加速比与效率,引入消息传递接口(MPI)实现了声波方程的并行求解,极大地提高了数值求解声波方程的计算效率。
Forward modeling of seismic wave equations play a significant role in seismic data acquisition,processing,interpretation and inversion.At present however,forward modeling algorithms based on seismic wave-equation are restrained from large-scale industrial application because of their massive computational cost.This paper studies the basic methods for solving acoustic equation with finite difference and parallel algorithm,and proposes a computational domain division method and proper communication scheme for parallel algorithm.We also analyzed the difference of computational time cost,the speedup and efficiency of parallel programs under different parameters.Numerical experiments of parallel solution of acoustic equation with message passing interface (MPI)showed that the introduction of MPI greatly improves the computation efficiency.
Forward modeling of seismic wave equations play a significant role in seismic data acquisition,processing,interpretation and inversion.At present however,forward modeling algorithms based on seismic wave-equation are restrained from large-scale industrial application because of their massive computational cost.This paper studies the basic methods for solving acoustic equation with finite difference and parallel algorithm,and proposes a computational domain division method and proper communication scheme for parallel algorithm.We also analyzed the difference of computational time cost,the speedup and efficiency of parallel programs under different parameters.Numerical experiments of parallel solution of acoustic equation with message passing interface (MPI)showed that the introduction of MPI greatly improves the computation efficiency.
引文
1孙卫涛,杨慧珠.各向异性介质弹性波传播的三维不规则网格有限差分方法[J].地球物理学报,2004,47(2):332~337
2张会星,何兵寿,张晶,等.复杂各向异性介质中的地震波场有限差分模拟[J].煤炭学报,2008,33(11):1 257~1 262
3董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411~419
4牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005.215
5陈伟.起伏地表条件下二维地震波场的数值模拟[J].勘探地球物理进展,2005,28(1):25~31
6佘德平.波场数值模拟技术[J].勘探地球物理进展,2004,27(1):16~21
7 Berfenger J P.A perfectly matched layer for the ab-sorption of electromagnetics waves[J].Journal Compu-tation Physics,1994,114(2):185~200
8都志辉.高性能计算并行编程技术——MPI并行程序设计[M].北京:清华大学出版社,2001.352
2张会星,何兵寿,张晶,等.复杂各向异性介质中的地震波场有限差分模拟[J].煤炭学报,2008,33(11):1 257~1 262
3董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411~419
4牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005.215
5陈伟.起伏地表条件下二维地震波场的数值模拟[J].勘探地球物理进展,2005,28(1):25~31
6佘德平.波场数值模拟技术[J].勘探地球物理进展,2004,27(1):16~21
7 Berfenger J P.A perfectly matched layer for the ab-sorption of electromagnetics waves[J].Journal Compu-tation Physics,1994,114(2):185~200
8都志辉.高性能计算并行编程技术——MPI并行程序设计[M].北京:清华大学出版社,2001.352
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心