三维粘弹介质地震波场有限差分并行模拟
|
摘要
有限差分法在三维粘弹性复杂介质正演模拟地震波的传播中对计算机内存和计算速度要求比较高,单个PC机或工作站只能计算较少网格内短时间的波场。本文介绍一种基于MPI的并行有限差分法,可在PCCluster上模拟较大规模三维粘弹性复杂介质中地震波传播时的波场;可预测地震波在此类条件下传播时的运动学和动力学性质。对于更好地理解波动传播现象,解释实际地震资料及反问题的解决等均具有重要的理论与实际意义。
When finite difference(FD) method is used in modeling the propagation of seismic wave in 3-D viscoelastic complex media,it consumes vast quantities of computational resources.So on single PC or workstation,3-D calculations are still limited to small grid sizes and short seismic wave traveltimes.In this paper,the parallel FD algorithm which based on a message passing interface(MPI) is introduced to solve above problem properly.Using PCCluster we can calculate the wavefield of the large 3-D viscoelastic complex models,furthermore predict and understand the kinematic and dynamic properties of seismic waves propagating through the models of the crust.It helps us in every stage of a seismic investigation.
When finite difference(FD) method is used in modeling the propagation of seismic wave in 3-D viscoelastic complex media,it consumes vast quantities of computational resources.So on single PC or workstation,3-D calculations are still limited to small grid sizes and short seismic wave traveltimes.In this paper,the parallel FD algorithm which based on a message passing interface(MPI) is introduced to solve above problem properly.Using PCCluster we can calculate the wavefield of the large 3-D viscoelastic complex models,furthermore predict and understand the kinematic and dynamic properties of seismic waves propagating through the models of the crust.It helps us in every stage of a seismic investigation.
引文
[1]Day S M,J B Minster.Numerical simulation of attenuatedwavefield using a pade approcimant method[J].Geophys.J.R.astr.Soc.,1984,78:105-118.
[2]Alterman Z S.Finite difference solutions to geophysical prob-lems[J].Journal of Physics of the Earth,1968,16:113-128.
[3]Johan O A Robertsson,Joakim O Blanch.Viscoelastic finite-difference modeling[J].Geophysics,1994,59(9):1444-1456.
[4]Liu H P,Anderson D L,Kanmori H.Velocity dispersion dueto anelasticity implication for seismology and mantle composi-tion[J].Geophys.J.R.astr.Soc..1976,47:41-58.
[5]陈国良.并行计算——结构.算法.编程[M].北京:高等教育出版社,1999.
[6]陈国良.并行算法实践[M].北京:高等教育出版社,2004.
[7]J O Blanch,J O A Robertsson,William W Symes.Modelingof a constant:methodology and algorithm for an efficient andoptimally inexpensive viscoelastic technique[J].Geophysics,1995,60(1):176-184
[8]Carcione J M,Kosloff D,Koslff R.Wave propagation simula-tion in a viscoelastic medium[J].Geophysical Jourmal of theRoyal Astronomical Society,1988,93:597-611.
[9]A R Levander.Fourth-order finite-difference P-SV seismo-grams[J].Geophysics,1988,53:1425-1436.
[10]Thomas Bohlen.Parallel 3-D viscoelastic finite differenceseismic modeling[J].Computers&Geosciences.2002,28:887-899.
[11]J O Blanch,A C Cheng,G L Varsamis.Attainable computa-tional speed for large-scale seismic modeling on PC-basedcluster[A]∥SEG Annual Meeting[C].2005.
[12]Roth,Korn Michael.Seismic waves in random media[J].Journal of Applied Geophysics.1993,29:247-269.
[2]Alterman Z S.Finite difference solutions to geophysical prob-lems[J].Journal of Physics of the Earth,1968,16:113-128.
[3]Johan O A Robertsson,Joakim O Blanch.Viscoelastic finite-difference modeling[J].Geophysics,1994,59(9):1444-1456.
[4]Liu H P,Anderson D L,Kanmori H.Velocity dispersion dueto anelasticity implication for seismology and mantle composi-tion[J].Geophys.J.R.astr.Soc..1976,47:41-58.
[5]陈国良.并行计算——结构.算法.编程[M].北京:高等教育出版社,1999.
[6]陈国良.并行算法实践[M].北京:高等教育出版社,2004.
[7]J O Blanch,J O A Robertsson,William W Symes.Modelingof a constant:methodology and algorithm for an efficient andoptimally inexpensive viscoelastic technique[J].Geophysics,1995,60(1):176-184
[8]Carcione J M,Kosloff D,Koslff R.Wave propagation simula-tion in a viscoelastic medium[J].Geophysical Jourmal of theRoyal Astronomical Society,1988,93:597-611.
[9]A R Levander.Fourth-order finite-difference P-SV seismo-grams[J].Geophysics,1988,53:1425-1436.
[10]Thomas Bohlen.Parallel 3-D viscoelastic finite differenceseismic modeling[J].Computers&Geosciences.2002,28:887-899.
[11]J O Blanch,A C Cheng,G L Varsamis.Attainable computa-tional speed for large-scale seismic modeling on PC-basedcluster[A]∥SEG Annual Meeting[C].2005.
[12]Roth,Korn Michael.Seismic waves in random media[J].Journal of Applied Geophysics.1993,29:247-269.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心