基于MPI的OBC三维多波多分量地震观测正演模拟并行算法实现
|
摘要
针对海上OBC三维地震观测正演数值模拟的海量计算需求,研究并实现了基于信息传递接口(Message Passing Interface,MPI)并行平台的交错网格有限差分数值模拟并行算法。该算法能够以三维任意偶阶空间精度从X,Y和Z3个方向同时实现并行运算,可以任意选择空间并行计算方向以及并行计算的进程个数。在算法实现过程中,引入了进程坐标系的概念,使得每一个参与计算的进程对应一个进程坐标,通过进程坐标来判断边界条件的使用;同时引入虚拟进程以方便模型边界处的数据交换。通过层状模型采用不同计算进程个数的正演计算验证了该并行算法的可行性和运行效率;最后对海上某靶区三维地质模型进行了Patch观测系统下的多波多分量地震观测正演并行模拟。
For massive computation on forward numerical modeling of 3Dseismic survey in OBC observation,a parallel staggered-grid finite-difference modeling algorithm which based on message passing interface(MPI)is studied and implemented.This algorithm can achieve parallel computing from the X,Yand Zthree directions simultaneously with any order spatial accuracy and choose any direction of space and the number of process in parallel computing.During the realization of this parallel algorithm,we introduce the concept of the process coordinate system.Therefore each process has a process coordinate which can judge the use of boundary conditions as well.Virtual process is introduced to facilitate the exchange of data at the boundary of the model.The different numbers of the total processes are tested in forward modeling of a simple model to verify the feasibility and efficiency of this algorithm.Then this algorithm is used for the parallel modeling of 3D multi-wave and multi-component seismic data acquisition in a marine target region with OBC Patch observation.
For massive computation on forward numerical modeling of 3Dseismic survey in OBC observation,a parallel staggered-grid finite-difference modeling algorithm which based on message passing interface(MPI)is studied and implemented.This algorithm can achieve parallel computing from the X,Yand Zthree directions simultaneously with any order spatial accuracy and choose any direction of space and the number of process in parallel computing.During the realization of this parallel algorithm,we introduce the concept of the process coordinate system.Therefore each process has a process coordinate which can judge the use of boundary conditions as well.Virtual process is introduced to facilitate the exchange of data at the boundary of the model.The different numbers of the total processes are tested in forward modeling of a simple model to verify the feasibility and efficiency of this algorithm.Then this algorithm is used for the parallel modeling of 3D multi-wave and multi-component seismic data acquisition in a marine target region with OBC Patch observation.
引文
[1]Caldwell J.Marine multicomponent seismology[J].The Leading Edge,1999,18(11):1274-1282
[2]张树林,夏斌,何家雄.海上多波多分量地震采集技术的应用[J].天然气地球科学,2005,16(1):103-107Zhang S L,Xia B,He J X.The offshore multi-waves and multi-components seismic acquisition technique:a case study of Yinggehai Basin[J].Natural Gas Geoscience,2005,16(1):103-107
[3]张树林,李绪宣,姜立红.海上多波多分量地震技术新进展与发展方向[J].物探化探计算技术,2000,22(2):97-107Zhang S L,Li X X,Jiang L H.Improvement and development of China offshore multiwave and multicomponent seismic technique[J].Computing Techniques for Geophysical and Geochemical Exploration,2000,22(2):97-107
[4]刘海波,全海燕,陈浩林,等.海上多波多分量地震采集综述[J].中国石油勘探,2007,12(3):52-57Liu H B,Quan H Y,Chen H L,et al.Overview of marine multi-wave and multi-component seismic acquisition[J].China Petroleum Exploration,2007,12(3):52-57
[5]He H Y,Zhou H Z,Fu D D.Multicomponent seismic-exploration in Yinggehai Basin[J].The Leading Edge,2002,21(9):914-920
[6]Virieux J.SH wave propagation in heterogenous media:velocity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1957
[7]Virieux J.P-SV wave propagation in heterogenous media:velocity-stress finite-difference method[J].Geophysics,1986,51(4):889-901
[8]Levander A R.Fourth-order finite-difference P-SV seismograms[J].Geophysics,1988,53(11):1425-1436
[9]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419Dong L G,Ma Z T,Cao J Z,et al.A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(3):411-419
[10]裴正林.三维各向同性介质弹性波方程交错网格高阶有限差分模拟[J].石油物探,2005,44(4):308-315Pei Z L.Numerical simulation of elastic wave equation in 3-D isotropic media with staggered-grid highorder difference method[J].Geophysical Prospecting for Petroleum,2005,44(4):308-315
[11]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5):23-29Pei Z L.Three-dimensional numerical simulation of elastic wave propagation in 3-D anisotropic media with staggered-grid high-order difference method[J].Journal of the University of Petroleum,China,2004,28(5):23-29
[12]张军华,臧胜涛,单联瑜,等.高性能计算的发展现状及趋势[J].石油地球物理勘探,2010,45(6):918-925Zhang J H,Zang S T,Shan L Y,et al.Development status and trends for high performance computing[J].Oil Geophysical Prospecting for Petroleum,2010,45(6):918-925
[13]都志辉.高性能计算并行编程技术—MPI并行程序设计[M].北京:清华大学出版社,2001:13-15Yu Z H.Parallel program technique of high performance computation—MPI parallel program design[M].Beijing:Tsinghua University Press,2001:13-15
[14]Furumura T,Koketsu K,Wen K L.Parallel PSM/FDM hybrid simulation of ground motions from the1999Chi-chi,Taiwan,earthquake[J].Pure and Applied Geophysics,2002,159(9):2133-2146
[15]Bohlen T.Parallel 3-D viscoelastic finite difference seismic modeling[J].Computers&Geosciences,2002,28(8):887-899
[16]Micikevicius P.3Dfinite difference computation on GPUs using CUDA[J].Expanded Abstracts of 2nd Workshop on General Purpose Processing on Graphics Processing Units,2009,79-84
[17]王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1):30-34Wang D L,Yong Y D,Han L G,et al.Parallel simulation of finite difference for seismic wavefiled modeling in 3-D viscoelastic media[J].Northwestern Seismological Journal,2007,29(1):30-34
[18]王月英.基于MPI的三维波动方程有限元法并行正演模拟[J].石油物探,2009,48(3):221-243Wang Y Y.3Dwave-equation finite-element forward modeling based on message passing interface parallel computing[J].Geophysical Prospecting for Petroleum,2009,48(3):221-243
[19]何兵寿,张会星,韩令贺.弹性波方程正演的粗粒度并行算法[J].地球物理学进展,2010,25(2):650-656He B S,Zhang H X,Han L H.Forward modeling of elastic wave equation with coarse-grained algorithem[J].Progress in Geophysic,2010,25(2):650-656
[20]魏星,王彦宾,陈晓非.模拟地震波场的伪谱和高阶有限差分混合方法[J].地震学报,2010,32(4):392-400Wei X,Wang Y B,Chen X F.Hybrid PSM/FDM method for seismic wavefield simulation[J].Acta Seismologica Sinica,2010,32(4):392-400
[21]秦艳芳,王彦宾.地震波传播的三维伪谱和高阶有限差分混合方法并行模拟[J].地震学报,2012,34(2):147-156Qin Y F,Wang Y B.Three-dimensional parallel hybrid PSM/FDM simulation of seismic wave propagation[J].Acta Seismologica Sinica,2012,34(2):147-156
[22]公续飞,杜启振.三维弹性波有限差分数值模拟及其GPU并行实现[C]∥中国地球物理学会第二十七届年会论文集.长沙:中国地球物理学会,2011:485Gong X F,Du Q Z.Finite-difference simulation of3Delastic wavefield and its implementation on GPU[C]∥Chinese Geophysical Society 27th Annual Meeting Proceeding.Changsha:Chinese Geophysical Society,2011:485
[23]段玉婷,李靖宇,胡天悦.基于GPU的三维精细积分法正演模拟[C]∥中国地球物理学会第二十七届年会论文集.长沙:中国地球物理学会,2011:484Duan Y T,Li J Y,Hu T Y.Application of the 3D precise integration method for seismic modeling based on GPU[C]∥Chinese Geophysical Society27th Annual Meeting Proceeding.Changsha:Chinese Geophysical Society,2011:484
[24]张明财,熊章强,张大洲.基于MPI的三维瑞雷面波有限差分并行模拟[J].石油物探,2013,52(4):354-362Zhang M C,Xiong Z Q,Zhang D Z.3Dfinite difference parallel simulation of Rayleigh wave based on message passing interface[J].Geophysical Prospecting for Petroleum,2005,44(4):308-315
[25]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005:16-18Mou Y G,Pei Z L.Seismic numerical modeling for 3-D complex media[M].Beijing:Petroleum Industry Press,2005:16-18
[2]张树林,夏斌,何家雄.海上多波多分量地震采集技术的应用[J].天然气地球科学,2005,16(1):103-107Zhang S L,Xia B,He J X.The offshore multi-waves and multi-components seismic acquisition technique:a case study of Yinggehai Basin[J].Natural Gas Geoscience,2005,16(1):103-107
[3]张树林,李绪宣,姜立红.海上多波多分量地震技术新进展与发展方向[J].物探化探计算技术,2000,22(2):97-107Zhang S L,Li X X,Jiang L H.Improvement and development of China offshore multiwave and multicomponent seismic technique[J].Computing Techniques for Geophysical and Geochemical Exploration,2000,22(2):97-107
[4]刘海波,全海燕,陈浩林,等.海上多波多分量地震采集综述[J].中国石油勘探,2007,12(3):52-57Liu H B,Quan H Y,Chen H L,et al.Overview of marine multi-wave and multi-component seismic acquisition[J].China Petroleum Exploration,2007,12(3):52-57
[5]He H Y,Zhou H Z,Fu D D.Multicomponent seismic-exploration in Yinggehai Basin[J].The Leading Edge,2002,21(9):914-920
[6]Virieux J.SH wave propagation in heterogenous media:velocity-stress finite-difference method[J].Geophysics,1984,49(11):1933-1957
[7]Virieux J.P-SV wave propagation in heterogenous media:velocity-stress finite-difference method[J].Geophysics,1986,51(4):889-901
[8]Levander A R.Fourth-order finite-difference P-SV seismograms[J].Geophysics,1988,53(11):1425-1436
[9]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419Dong L G,Ma Z T,Cao J Z,et al.A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(3):411-419
[10]裴正林.三维各向同性介质弹性波方程交错网格高阶有限差分模拟[J].石油物探,2005,44(4):308-315Pei Z L.Numerical simulation of elastic wave equation in 3-D isotropic media with staggered-grid highorder difference method[J].Geophysical Prospecting for Petroleum,2005,44(4):308-315
[11]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5):23-29Pei Z L.Three-dimensional numerical simulation of elastic wave propagation in 3-D anisotropic media with staggered-grid high-order difference method[J].Journal of the University of Petroleum,China,2004,28(5):23-29
[12]张军华,臧胜涛,单联瑜,等.高性能计算的发展现状及趋势[J].石油地球物理勘探,2010,45(6):918-925Zhang J H,Zang S T,Shan L Y,et al.Development status and trends for high performance computing[J].Oil Geophysical Prospecting for Petroleum,2010,45(6):918-925
[13]都志辉.高性能计算并行编程技术—MPI并行程序设计[M].北京:清华大学出版社,2001:13-15Yu Z H.Parallel program technique of high performance computation—MPI parallel program design[M].Beijing:Tsinghua University Press,2001:13-15
[14]Furumura T,Koketsu K,Wen K L.Parallel PSM/FDM hybrid simulation of ground motions from the1999Chi-chi,Taiwan,earthquake[J].Pure and Applied Geophysics,2002,159(9):2133-2146
[15]Bohlen T.Parallel 3-D viscoelastic finite difference seismic modeling[J].Computers&Geosciences,2002,28(8):887-899
[16]Micikevicius P.3Dfinite difference computation on GPUs using CUDA[J].Expanded Abstracts of 2nd Workshop on General Purpose Processing on Graphics Processing Units,2009,79-84
[17]王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1):30-34Wang D L,Yong Y D,Han L G,et al.Parallel simulation of finite difference for seismic wavefiled modeling in 3-D viscoelastic media[J].Northwestern Seismological Journal,2007,29(1):30-34
[18]王月英.基于MPI的三维波动方程有限元法并行正演模拟[J].石油物探,2009,48(3):221-243Wang Y Y.3Dwave-equation finite-element forward modeling based on message passing interface parallel computing[J].Geophysical Prospecting for Petroleum,2009,48(3):221-243
[19]何兵寿,张会星,韩令贺.弹性波方程正演的粗粒度并行算法[J].地球物理学进展,2010,25(2):650-656He B S,Zhang H X,Han L H.Forward modeling of elastic wave equation with coarse-grained algorithem[J].Progress in Geophysic,2010,25(2):650-656
[20]魏星,王彦宾,陈晓非.模拟地震波场的伪谱和高阶有限差分混合方法[J].地震学报,2010,32(4):392-400Wei X,Wang Y B,Chen X F.Hybrid PSM/FDM method for seismic wavefield simulation[J].Acta Seismologica Sinica,2010,32(4):392-400
[21]秦艳芳,王彦宾.地震波传播的三维伪谱和高阶有限差分混合方法并行模拟[J].地震学报,2012,34(2):147-156Qin Y F,Wang Y B.Three-dimensional parallel hybrid PSM/FDM simulation of seismic wave propagation[J].Acta Seismologica Sinica,2012,34(2):147-156
[22]公续飞,杜启振.三维弹性波有限差分数值模拟及其GPU并行实现[C]∥中国地球物理学会第二十七届年会论文集.长沙:中国地球物理学会,2011:485Gong X F,Du Q Z.Finite-difference simulation of3Delastic wavefield and its implementation on GPU[C]∥Chinese Geophysical Society 27th Annual Meeting Proceeding.Changsha:Chinese Geophysical Society,2011:485
[23]段玉婷,李靖宇,胡天悦.基于GPU的三维精细积分法正演模拟[C]∥中国地球物理学会第二十七届年会论文集.长沙:中国地球物理学会,2011:484Duan Y T,Li J Y,Hu T Y.Application of the 3D precise integration method for seismic modeling based on GPU[C]∥Chinese Geophysical Society27th Annual Meeting Proceeding.Changsha:Chinese Geophysical Society,2011:484
[24]张明财,熊章强,张大洲.基于MPI的三维瑞雷面波有限差分并行模拟[J].石油物探,2013,52(4):354-362Zhang M C,Xiong Z Q,Zhang D Z.3Dfinite difference parallel simulation of Rayleigh wave based on message passing interface[J].Geophysical Prospecting for Petroleum,2005,44(4):308-315
[25]牟永光,裴正林.三维复杂介质地震数值模拟[M].北京:石油工业出版社,2005:16-18Mou Y G,Pei Z L.Seismic numerical modeling for 3-D complex media[M].Beijing:Petroleum Industry Press,2005:16-18
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心