基于MPI的三维波动方程有限元法并行正演模拟
|
摘要
在三维空间进行地震波动方程有限元正演模拟时,采用基于消息传递界面(MPI)的并行算法可以克服基于单PC机串行算法对数据容量和计算速度的局限。其基本原理是:将模拟区域剖分成多个小区域,每一个进程处理其中一个小区域;在运算过程中,各个进程之间互通相邻区域节点的物性参数和前一时刻的位移值,以备计算当前时刻区域内的位移值,共同实现整个模拟空间的正演模拟。在正演模拟时,采用质量矩阵近似方法来提高程序的并行性,压缩数据量和运算量,最佳并行效率的进程个数则依据并行系统能够满足程序内数据容量来确定。通过3层水平层状介质模型的数值模拟,对方法的可行性和有效性进行了验证。
In 3D spatial seismic wave-equation finite-element forward modeling,the parallel algorithm based on message passing interface (MPI) can overcome the limitation of data capacity and computation velocity in single PC with serial algorithm.The basic principles are as follows:the simulated area is divided into several small subregions, and each processor computes one small region;in computation, the physical property parameters and the displacement at the former moment of the nodes from neighboring regions are communicated in processors each other,to prepare for calculating the displacement of the region at current moment and realize forward modeling of the whole simulated space together.In forward modeling, mass matrix approximation method was adopted to improve the parallelism of the program,compress the data volume and calculated amount.Processor number for optimal parallel efficiency was set in terms of parallel system satisfying data capacity in program. Through the numerical simulation of three-layer horizontal layered medium model,the method was proved feasible and effective.
In 3D spatial seismic wave-equation finite-element forward modeling,the parallel algorithm based on message passing interface (MPI) can overcome the limitation of data capacity and computation velocity in single PC with serial algorithm.The basic principles are as follows:the simulated area is divided into several small subregions, and each processor computes one small region;in computation, the physical property parameters and the displacement at the former moment of the nodes from neighboring regions are communicated in processors each other,to prepare for calculating the displacement of the region at current moment and realize forward modeling of the whole simulated space together.In forward modeling, mass matrix approximation method was adopted to improve the parallelism of the program,compress the data volume and calculated amount.Processor number for optimal parallel efficiency was set in terms of parallel system satisfying data capacity in program. Through the numerical simulation of three-layer horizontal layered medium model,the method was proved feasible and effective.
引文
1孙成禹,万学娟,倪长宽.弱连接界面地层地震反射波场的正演方法[J].石油物探,2007,46(5):425-432
2郭智奇,刘财,杨宝俊,等.粘弹各向异性介质中地震波场模拟与特征[J].地球物理学进展,2007,22(3): 804-810
3李胜军,孙成禹,高建虎,等.地震波数值模拟中频散压制方法分析[J].石油物探,2008,47(5):444-448
4杜世通.变速不均运介质中波动方程的有限元法数值解[J].华东石油学院学报,1982,6(2):1-20
5杜世通.地震波动力学理论与方法[M].山东:中国石油大学出版社,2008.189-251
6 Padovani E,Priolo E,Seriani G.Low-and high-order finite element method:experience in seismic modeling ??[J].Journal of Computational Acoustics,1994(2): 371-422
7周辉,刘斌,徐世浙.二维各向异性介质中三分量波动方程有限元法模拟[J].青岛海洋大学学报,2000, 30(4):668-672
8杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报,2002,45(4):575-583
9孙卫涛,杨慧珠.各向异性介质弹性波传播的三维不规则网格有限差分方法[J].地球物理学报,2004, 47(2):332-337
10王德利,雍运动,韩立国,等.三维粘弹性介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1): 30-34
11陶果,何峰江,王兵,等.声反射成像测井在地层中的三??维波场模拟[J].中国科学D辑:地球科学,2008,38(增刊Ⅰ):166-173
12殷文.基于频率域高阶有限差分的正演模拟及并行算法[J].吉林大学学报(地球科学版),2008,38(1): 144-151
13 Sarma G S,Mallick K,Gadhinglajkar V R.Nonreflecting boundary condition in finite-element formulation for an elastic wave equation[J].Geophysics,1998, 63(3):1006-1016
14王月英,宋建国.波场正演模拟中Sarma边界条件的改进[J].石油物探,2007,46(4):353-362
15薛东川,王尚旭.波动方程有限元叠前逆时偏移[J].石油地球物理勘探,2008,43(1):17-21
2郭智奇,刘财,杨宝俊,等.粘弹各向异性介质中地震波场模拟与特征[J].地球物理学进展,2007,22(3): 804-810
3李胜军,孙成禹,高建虎,等.地震波数值模拟中频散压制方法分析[J].石油物探,2008,47(5):444-448
4杜世通.变速不均运介质中波动方程的有限元法数值解[J].华东石油学院学报,1982,6(2):1-20
5杜世通.地震波动力学理论与方法[M].山东:中国石油大学出版社,2008.189-251
6 Padovani E,Priolo E,Seriani G.Low-and high-order finite element method:experience in seismic modeling ??[J].Journal of Computational Acoustics,1994(2): 371-422
7周辉,刘斌,徐世浙.二维各向异性介质中三分量波动方程有限元法模拟[J].青岛海洋大学学报,2000, 30(4):668-672
8杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报,2002,45(4):575-583
9孙卫涛,杨慧珠.各向异性介质弹性波传播的三维不规则网格有限差分方法[J].地球物理学报,2004, 47(2):332-337
10王德利,雍运动,韩立国,等.三维粘弹性介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1): 30-34
11陶果,何峰江,王兵,等.声反射成像测井在地层中的三??维波场模拟[J].中国科学D辑:地球科学,2008,38(增刊Ⅰ):166-173
12殷文.基于频率域高阶有限差分的正演模拟及并行算法[J].吉林大学学报(地球科学版),2008,38(1): 144-151
13 Sarma G S,Mallick K,Gadhinglajkar V R.Nonreflecting boundary condition in finite-element formulation for an elastic wave equation[J].Geophysics,1998, 63(3):1006-1016
14王月英,宋建国.波场正演模拟中Sarma边界条件的改进[J].石油物探,2007,46(4):353-362
15薛东川,王尚旭.波动方程有限元叠前逆时偏移[J].石油地球物理勘探,2008,43(1):17-21
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心