间断的Galerkin方法在地震波场数值模拟中的应用概述
|
摘要
通过数值求解描述地震波传播的微分方程来模拟波的传播角度而言,在日趋繁荣的地震波场数值模拟方法的探索中,间断的Galerkin方法(Discontinuous Galerkin,DG)是继有限差分法、伪谱法、有限元法以及谱元法之后出现的另一种有效的数值模拟方法,该方法从2005年以来在地震波场模拟中得到了迅速发展,给这一学科注入了新的活力.但是从国内的发展来看,在地震波场数值模拟方面还没有得到足够的关注.尤其是现有的综述性文献都忽略了这一方法得到的有益的成果,本文将重点介绍间断的Galerkin方法的发展及其在地震波场数值模拟的最新应用研究.
In the past years,a number of different numerical methods have been developed in order to solve the seismic wave equations,such as finite difference method,pseudospectral method,finite element method and spectral element method etc.In these methods,discontinuous Galerkin(DG) method is one of the effective methods,which is of an intelligent combination of finite element and the finite volume methods,utilizing a space of basis and test functions that mimics the finite element method but satisfying the equation in as sense close to the finite volume method,appears to offer many of the desired properties as allowing for jumps of material parameters and the solutions across the element interfaces.The discontinuous Galerkin method together with a time integration approach using arbitrary high-order derivatives(ADER) has been further developed rapidly for modeling of seismic wave propagation since 2005.However,from the perspective of development in our country,it does not gained enough attention,especially in seismic modeling.In this paper,it will focus on the development of DG method and its latest application in modeling of the seismic wavefields.
In the past years,a number of different numerical methods have been developed in order to solve the seismic wave equations,such as finite difference method,pseudospectral method,finite element method and spectral element method etc.In these methods,discontinuous Galerkin(DG) method is one of the effective methods,which is of an intelligent combination of finite element and the finite volume methods,utilizing a space of basis and test functions that mimics the finite element method but satisfying the equation in as sense close to the finite volume method,appears to offer many of the desired properties as allowing for jumps of material parameters and the solutions across the element interfaces.The discontinuous Galerkin method together with a time integration approach using arbitrary high-order derivatives(ADER) has been further developed rapidly for modeling of seismic wave propagation since 2005.However,from the perspective of development in our country,it does not gained enough attention,especially in seismic modeling.In this paper,it will focus on the development of DG method and its latest application in modeling of the seismic wavefields.
引文
[1]Komatitsch D,Vilotte J P.The spectral-element method:anefficient tool to simulate the seismic response of 2Dand 3Dgeological structures[J].Bull.Seismol.Soc.Am.,1998,88:368-392.
[2]Chaljub E,Valette B.Spectral element modelling ofthreedimensional wave propagation in a self-gravitating Earthwith an arbitrarily stratified outer core[J].Geophys.J.Int.,2004,158:131-141.
[3]汪文帅,李小凡,鲁明文,张美根.基于多辛结构谱元法的保结构地震波场模拟.地球物理学报[J],2012,(10):3427-3439.WANG Wen-Shuai,LI Xiao-Fan,LU Ming-Wen,ZHANGMei-Gen.Structure-preserving modeling for seismic wavefieldsbased upon a multisymplectic spectral element method.ChineseJ.Geophys.[J](in Chinese),2012,(10):3427-3439.
[4]严红勇,刘洋.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟.地球物理学报[J],2012,55(4):1354-1365.YAN Hong-Yong,LIU Yang.Rotated staggered grid high-order finite-difference numerical modeling for wave propagationin viscoelastic TTI media.Chinese J.Geophys.[J](inChinese),2012,55(4):1354-1365.
[5]祝贺君,张伟,陈晓非.二维各向异性介质中地震波场的高阶同位网格有限差分模拟.地球物理学报[J],2009,52(6):1536-1546.ZHU He-Jun,ZHANG Wei,CHEN Xiao-Fei.Two-dimensional seismic wave simulation in anisotropic media bynon-staggered finite difference method.Chinese J.Geophys.[J](in Chinese),2009,52(6):1536-1546.
[6]程冰洁,李小凡,龙桂华.基于广义正交多项式褶积微分算子的地震波场数值模拟方法.地球物理学报[J],2008,51(2):531-537.CHENG Bing-Jie,LI Xiao-Fan,LONG Gui-Hua.Seismicwaves modeling by convolutional Forsyte polynomialdifferentiator method.Chinese J.Geophys.[J](in Chinese),2008,51(2):531-537.
[7]李一琼,李小凡,朱童.基于辛格式奇异核褶积微分算子的地震标量波场模拟.地球物理学报[J],2011,54(7):1827-1834.LI Yi-Qiong,LI Xiao-Fan,ZHU Tong.The seismic scalarwave field modeling by symplectic scheme and singular kernelconvolutional differentiator.Chinese J.Geophys.[J](inChinese),2011,54(7):1827-1834.
[8]李一琼,李小凡,张美根.基于辛格式离散奇异褶积微分算子的弹性波场模拟.地球物理学报[J],2012,(5):1725-1731.LI Yi-Qiong,LI Xiao-Fan,ZHANG Mei-Gen.The elasticwave fields modeling by symplectic discrete singularconvolution differentiator method.Chinese J.Geophys.[J](inChinese),2012,(5):1725-1731.
[9]龙桂华,李小凡,张美根.错格傅里叶伪谱微分算子在波场模拟中的应用.地球物理学报[J],2009,52(1):193-199.LONG Gui-Hua;LI Xiao-Fan;ZHANG Mei-Gen.Theapplication of staggered-grid Fourier pseudospectraldifferentiation operator in wavefield modeling.Chinese J.Geophys.[J](in Chinese),2009,52(1):193-199.
[10]张怀,周元泽,吴忠良,严珍珍,等.福州盆地强地面运动特征的有限元数值模拟.地球物理学报[J],2009,52(5):1270-1279.ZHANG Huai,ZHOU Yuan-Ze,WU Zhong-Liang,YANZhen-Zhen,et al.Finite element analysis of seismic wavepropagation characteristics in Fuzhou basin.Chinese J.Geophys.[J](in Chinese),2009,52(5):1270-1279.
[11]管西竹,符力耘,陶毅,于更新.复杂地表边界元-体积元波动方程数值模拟.地球物理学报[J],2011,54(9):2357-2367.GUAN Xi-Zhu,FU Li-Yun,TAO Yi,YU Geng-Xin.Boundary-volume integral equation numerical modeling forcomplex near surface.Chinese J.Geophys.[J].(in Chinese),2011,54(9):2357-2367.
[12]朱多林,白超英.基于波动方程理论的地震波场数值模拟方法综述[J].地球物理学进展,2011,26(5):1588-1599.Zhu D L,Bai C Y.Review on the seismic wavefield forwardmodeling[J].Progress in Geophysics(in Chinese),2011,26(5):1588-1599.
[13]李信富,李小凡,张美根.地震波数值模拟方法研究综述[J].防灾减灾工程学报,2007,27(2):241-248.Li X F,Li X F,Zhang M G.Review of seismic wavenumerical modeling methods[J].Journal of DisasterPrevention and Mitigation Engineering(in Chinese),2007,27(2):241-248.
[14]王雪秋,孙建国.地震波有限差分数值模拟框架下的起伏地表处理方法综述[J].地球物理学进展,2008,23(1):40-48.Wang X Q,Sun J G.The state-of-the-art in numericalmodeling including surface topography with finite differencemethod[J].Progress in Geophysics(in Chinese),2008,23(1):40-48.
[15]张华,李振春,韩文功.起伏地表条件下地震波数值模拟方法综述[J].勘探地球物理进展,2007,30(5):1-6.Zhang H,Li Z C,Han W G.Review of seismic wavenumerical simulation form irregular topography[J].Progressin Exploration Geophysics(in Chinese),2007,30(5):1-6.
[16]Versteeg H.,Malalasekra W.An introduction to computationalfluid dynamics:The finite volume method(2nd Edition)[M].Prentice Hall,2007.
[17]刘儒勋,舒其望.计算流体力学的若干新方法[M].科学出版社,北京,2003.Liu R X,Shu C W.Selected new numerical methods forcomputational fluid dynamics(in Chinese)[M].SciencePress,Beijing,2003.
[18]Kser M,Dumbser M.An arbitrary high order discontinuousGalerkin method for elastic waves on unstructured meshes I:The two-dimensional isotropic case with external sourceterms[J].Geophys.J.Int.,2006,166:855-877.
[19]Reed W H,Hill T R,Triangular mesh methods for theneutron transport equation[R].Los Alamos ScientificLaboratory Report,LA-UR-73-479,1973.
[20]LeSaint P,Raviart P A,On a finite element method for solvingthe Neutron Transport Equation.In mathematical aspects offinite elements in partial differential equations[C].AcademicPress,New York,1974:89-145.
[21]Johnson C,Pitkaranta J.An Analysis of the discontinuousGalerkin method for a scalar hyperbolic equation[J].Math.Comp.,1986,46:1-26.
[22]Babüska I,Zlámal M.Nonconforming elements in the finiteelement method with penalty[J].SIAM J.Num.Anal.,1973,10:863-875.
[23]Babuka I,The finite element method with penalty[J].Math.Comp.,1973,27:221-228.
[24]Baker G A,Finite element methods for elliptic equations usingnonconforming elements[J].Math.Comp.,1977,31:45-59.
[25]Wheeler M F.An elliptic collocation-finite element methodwith interior penalties[J].SIAM J.Numer.Anal.,1978,15:152-161
[26]Arnold D N.An interior penalty finite element method withdiscontinuous elements[J].SIAM J.Numer.Anal.,1982,19:742-760.
[27]Arnold D N,Berzzi F,Cockburn B,Marini L D.Unifiedanalysis of discontinuous Galerkin mehtod for ellipticproblems[J].SIAM J.Numer.Anal.,2002,39(5):1749-1779.
[28]Chavent G,Salzano G.A finite element method for the 1Dwater flooding problem with gravity[J].J.Comput.Phy.,1982,45:307-344.
[29]Chavent G,Cockburn B.The local projection p0p1-discontinuousGalerkin finite element for scalar conservation laws[J].M2AN,1989,23(4):565-592.
[30]Cockburn B,Shu C W.The Runge-Kutta local projection p1-discontinuous Galerkin method for scalar conservation laws[J].M2AN,1991,25:337-361.
[31]Shu C W,Osher S.Efficient implementation of essentiallynon-oscillatory shockcapturing schemes[J].J.Comput.Phys.,1988,77:439-471.
[32]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws II:General Framework[J].Mathematics ofComputation,1989,52(186):411-435.
[33]Cockburna B,Lina S Y,Shu C W.TVB Runge-Kutta localprojection discontinuous Galerkin finite element method forscalar conservation laws III:One-dimensional systems[J].Journal of Computational Physics,1989,84(1):90-113.
[34]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws IV:The multidimensional case[J].Mathematics of Computation,1990,54(190):545-581.
[35]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws V:Multidimensional systems[J].J.Comput.phys.,1998,141:199-224.
[36]Bassi F,Rebay S.Ahigh-or der accurate discontinuous finiteelement method for the numerical solution of the compressibleNavier-Stokes equations[J],J.Comput.Phys.,1997,131:267-279.
[37]Cockburn B,Shu C W.The local discontinuous Galerkinfinite element method for convection-diffusion systems[J].SIAM journal on numerical analysis,1998,35:2240-2463.
[38]Yan J,Shu C W.A local discontinuous Galerkin method forKDV type equations[J].SIAM J.Numer.Anal.,2002,40:769-791.
[39]于春丽.Allen-Cahn方程的局部间断Galerkin有限元方法[D],山东大学硕士学位论文,2009.Yu Chunli.Local discontinuous Galerkin finite elementmethods for Allben-Cahn equation[D].ShandongUniversity,Shandong,2009.
[40]Babuka I,Baumann C E,Oden J T.A discontinuous hpfinite element method for diffusion problems:1-D analysis[J].Comput.Math.Appl.,1999,37:103-122.
[41]Oden J T,Babuka I,Baumann C E.A discontinuous hpfinite element method for diffusion problems[J].J.Comput.Phys.,1998,146:491-519.
[42]Cockburn B,Karniadkis G E,Shu C W.The development ofDiscontinuous Galerkin Methods,in Discontinuous GalerkinMethods:Theory,Computation and Applications[M].Springer-Verlag,NY,2000.
[43]Rivière B,Wheeler M F.Discontinuous finite element methodsfor acoustic and elastic wave problems.part I:Semidiscreteerror estimates[R].TICAM report,2001,0:1-12.
[44]Qiu J,Shu C W.Hermite WENO schemes and their applicationas limiters for Runge-Kutta discontinuous Galerkin methodII:two dimensional case[J].Comput.Fluids,2005,34:642-663.
[45]Luo H,Baum J D,Lhner R.A Hermite WENO-based limiterfor discontinuous Galerkin method on unstructured grids[J].J.Comput.Phys.,2007,225:686-713.
[46]Hesthaven J S,Warburton T.High-order nodal methods onunstructured grids,I.Time-domain solution of Maxwell’sequations[J].J.Comput.Phys.,2002,181:186-221.
[47]Cockburn B,Li F,Shu C W.Local divergence-free discontinuousGalerkin methods for the Maxwell equations[J].J.Comput.Phys.,2004,194:588-610.
[48]Fezoui L,Lanteri S,Lohrengel S,Piperno S.Convergenceand stability of a discontinuous Galerkin time-domain methodfor the 3Dheterogeneous Maxwell equations on unstructuredmeshes[J].M2AN,Math.Model Numer.Anal.,2005,39:1149-1176.
[49]Grote M,Schneebeli A,Schtzau D.Discontinuous Galerkinfinite element method for the wave equation[J].SIAM J.Num.Anal.,2006,44:2408-2431.
[50]Monk P,Richter G R.A discontinuous Galerkin method forlinear symmetric hyperbolic systems in inhomogeneous media[J].J.Sci.Comput.,2005,22:443-477.
[51]Chung E T,Engquist B.Optimal discontinuous Galerkinmethods for wave propagation[J].SIAM J.Sci.Comput.2006,44:2131-2158.
[52]Etienne V,Chaljub E,Virieux J,Glinsky N.An hp-adaptivediscontinuous Galerkin finite-element method for 3-D elasticwave modelling[J].Geophys.J.Int.,2010,183:941-962
[53]Komatitsch D,Martin R.An unsplit convolutional perfectlymatched layer improved at grazing incidence for the seismicwave equation[J].Geophysics,2007,72(5):SM155-SM167.
[54]Peraire J,Persson P O.The compact discontinuous Galerkin(CDG)method for elliptic problems[J].SIAM J.Sci.Comput.,2008,30(4):1806-1824.
[55]Aagaard B T,Hall J F,Heaton T H.Characterization ofnearsource ground motion with earthquake simulations[J].Earthq.Spectra.,2001,17:177-207.
[56]Akcelik V,et al.High resolution forward and inverse earthquakemodeling on terascale computers,in SC’03:Proceedings ofthe 2003 ACM/IEEE Conference on Supercomputing[C].IEEE Computer Society,Washington,DC,USA,2003,52.
[57]Ichimura T,Hori M,Kuwamoto H.Earthquake motionsimulation with multi-scale finite element analysis on hybridgrid[J].Bull.seism.Soc.Am.,2007,97(4):1133-1143.
[58]Chin-Joe-Kong M J S,Mulder W A,Van Veldhuizen M.Higherorder triangular and tetrahedral finite elements withmass lumping for solving the wave equation[J].J.Eng.Math.,1999,35:405-426.
[59]Chaljub E,Komatitsch D,Vilotte J P,Capdeville Y,ValetteB,Festa G.Spectral-element analysis in Seismology[J].Advances in Geophysics,2007,48:365-419.
[60]Komatitsch D,Tromp J.Spectral-element simulations of globalseismic wave propagation-II.3-D models,oceans,rotation,and self-gravitation[J].Geophys.J.Int.,2002,150:303-318.
[61]Komatitsch D,Labarta J,Michéa D.A simulation of seismicwave propagation at high resolution in the inner core of theEarth on 2166 processors of MareNostrum[C],LectureNotes Comput.Sci.,2008,5336:364-377.
[62]Mercerat E D,Vilotte J P,Sanchez-Sesma F J.Triangularspectral element simulation of two-dimensional elastic wavepropagation using unstructured triangular grids[J].Geophys.J.Int.,2006,166:679-698.
[63]Pasquetti R,Rapetti F.Spectral element methods onunstructured meshes:comparisons and recent advances[J].J.Scient.Comput.,2006,27:377-387.
[64]Hughes T J R.The finite element method-linear static anddynamic finite element analysis[M].Prentice Hall,1987.
[65]Igel H.Wave propagation in three-dimensional spherical sectionsby the Chebyshev spectral method[M].Geophys.J.Int.,1996,136:559-566.
[66]Thompson L L,Pinsky P.A space-time finite element methodfor structural acoustics in infinite domains.Part 1:formulation,stability and convergence[J].Comp.Meth.Appl.Mech.Eng.,1996,132:195-227.
[67]Dumbser M,Schwartzkopff T,Munz C D.Arbitrary highorder finite volume schemes for linear wave propagation,toappear in Notes on Numerical Fluid Mechanics andMultidisciplinary Design(NNFM)[C],Springer,2005.
[68]Toro E F,Millington A C,Nejad L A.Towards very highorder Godunov schemes,in Godunov methods:Theory andapplications[M].Kluwer Academic/Plenum Publishers,Oxford,2001:907-940.
[69]Toro E F,Titarev V A.Solution of the generalized Riemannproblem for advection-reaction equations[J].Proc.Roy.Soc.London,2002:271-281.
[70]Schwartzkopff T,Munz C D,Toro E F.ADER:A high-Order approach for linear hyperbolic systems in 2D[J].J.Sci.Comput.,2002,17:231-240.
[71]Schwartzkopff T,Dumbser M,Munz C D.Fast high orderADER schemes for linear hyperbolic equations[J].J.Comput.Phys.,2004,197:532-539.
[72]Dumbser M,Munz C D.ADER discontinuous Galerkin schemesfor Aeroacoustics[J].Comptes Rendus-Mecanique,2005,333:683-687.
[73]Dumbser M,Munz C D.Arbitrary high order discontinuousGalerkin schemes,in numerical methods for hyperbolic andkinetic problems,eds.Cordier S,Goudon T,Gutnic M,Sonnendrucker E,IRMA series in mathematics andtheoretical physics[J].EMS Publishing House,2005:295-333.
[74]Butcher J C.The numerical analysis of ordinary differentialequations:Runge-Kutta and general linear methods[M].Wiley,1987.
[75]Dumbser M,Kser M.An arbitrary high order discontinuousGalerkin method for elastic waves on unstructured meshes II:The three-dimensional isotropic case[J].Geophys.J.Int.2006,167:319-336.
[76]Kser M,Dumbser M,de la Puente J,Igel H.An arbitraryhigh order discontinuous Galerkin method for elastic waves onunstructured meshes III:Viscoelastic attenuation[J].Geophys.J.Int.,2007,168:224-242.
[77]de la Puente J,Kser M,Dumbser M,Igel H.An arbitraryhigh order discontinuous Galerkin method for elastic waves onunstructured meshes IV:Anisotropy[J].Geophys.J.Int.,2007,169:1210-1228.
[78]de la Puente J,Dumbser M,Kser M,Igel H.DiscontinuousGalerkin methods for wave propagation in poroelastic media[J].Geophysics,2008,73:77-97.
[79]Dumbser M,Kser M,Toro E F.An arbitrary high orderdiscontinuous Galerkin method for elastic waves onunstructured meshes V:Local Time Stepping and p-Adaptivity[J].Geophys.J.Int.,2007,171(3):695-717.
[80]Castro C E,Kser M,Brietzke G B.Seismic waves inheterogeneous material:subcell resolution of thediscontinuous Galerkin method[J].Geophys.J.Int.2010,182:250-264.
[81]de la Puente J,Ampuero J-P,Kser M.Dynamic rupturemodeling on unstructured meshes using a discontinuousGalerkin method,J.Geophys.Res.,2009,114,B10302.
[82]Pelties C,de la Puente J,Ampuero J-P,Brietzke G B,et al.Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructuredtetrahedral meshes,J.Geophys.Res.,2012,117,B02309.
[83]Delcourte S,Fezoui L,Glinsky-Olivier N.A high-orderdiscontinuous Galerkin method for the seismic wavepropagation[J].ESAIM:Proceedings,2009,27:70-89.
[84]Piperno S,Remaki M,Fezoui L.A non-diffusive finite volumescheme for the 3DMaxwell equations on unstructured meshes[J].SIAM J.Numer.Anal.,2002,39:2089-2108.
[85]Benjemaa M,Glinsky-Olivier N,Cruz-Atienza V M,VirieuxJ,et al.Dynamic non-planar crack rupture by a finite volumemethod[J].Geophys.J.Int.,2007,171:271-285.
[86]Kser M,Dumbser M.A highly accurate discontinuousGalerkin method for complex interfaces between solids andmoving fluids[J].Geophysics,2008,73(3):23-35.
[87]Moczo P,Kristek J,Galis M,Chaljub E,Etienne V.3-Dfinite-difference,finite-element,discontinuous-Galerkin andspectral-element schemes analysed for their accuracy withrespect to P-wave to S-wave speed ratio[J].Geophys.J.Int.2011,187:1645-1667.
[88]Bernacki M,Fezoui L,Lanteri S,Piperno S.Parallel discontinuousGalerkin unstructured mesh solvers for the calculation ofthree-dimensional wave propagation problems[J].AppliedMathematical Modeling,2006,30:744-763.
[89]Kser M,Hermann V,de la Puente J.Quantitative accuracyanalysis of the discontinuous Galerkin method for seismicwave propagation[J].Geophys.J.Int.,2008,173(2):990-999.
[90]Warburton T.Rocky mountain mathematics consortium summer2008:Parallel numerical methods for partial differentialequations[OL].2008,(http://www.caam.rice.edu/timwar/RMMC/MIDG.html).
[91]Hu F,Hussaini M,Rasitarinera P.An analysis of thediscontinuous Galerkin method for wave propagationproblems[J].J.Comput.Phys.,1999,151(2):921-946.
[92]Stanescu D,Kopriva D,Hussaini M.Dispersion analysis fordiscontinuous spectral element methods[J].J.Scient.Comput.,2000,15(2):149-171.
[93]Ainsworth M.Dispersive and dissipative behaviour of high orderdiscontinuous Galerkin finite element methods[J].J.Comput.Phys.,2004,198(1):106-130.
[94]Ainsworth M,Monk P,Muniz W.Dispersive and dissipativeproperties of discontinuous Galerkin finite element methodsfor the second-order wave equation[J].J.Scient.Comput.,2006,27(1):5-40.
[95]De Basabe,J D,Sen M K.Grid dispersion and stabilitycriteria of some common finite-element methods for acousticand elastic wave equations[J].Geophysics,2007,72:T81-T95.
[96]De Basabe J D,Sen M,Wheeler M.The interior penaltydiscontinuous Galerkin method for elastic wave propagation:grid dispersion[J].Geophys.J.Int.,2008,175:83-93.
[97]龙桂华,李小凡,江东辉.基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案[J].地球物理学报,2010,53(12):2964-2971.Long G H,Li X F,Jiang D H.Accelerating seismic modelingwith staggered-grid Fourier Pseudo-spectral differentiationmatrix operator method on graphics processing unit[J].Chinese J.Geophys.(in Chinese),2010,53(12):2964-2971.
[98]刘红伟,刘洪,李博,等.起伏地表叠前逆时偏移理论及GPU加速技术[J].地球物理学报,2011,54(7):1883-1892.Liu H W,Liu H,Li B,et al.Pre-stack reverse timemigration for rugged topography and GPU accelerationtechnology[J].Chinese J.Geophys.(in Chinese),2011,54(7):1883-1892.
[2]Chaljub E,Valette B.Spectral element modelling ofthreedimensional wave propagation in a self-gravitating Earthwith an arbitrarily stratified outer core[J].Geophys.J.Int.,2004,158:131-141.
[3]汪文帅,李小凡,鲁明文,张美根.基于多辛结构谱元法的保结构地震波场模拟.地球物理学报[J],2012,(10):3427-3439.WANG Wen-Shuai,LI Xiao-Fan,LU Ming-Wen,ZHANGMei-Gen.Structure-preserving modeling for seismic wavefieldsbased upon a multisymplectic spectral element method.ChineseJ.Geophys.[J](in Chinese),2012,(10):3427-3439.
[4]严红勇,刘洋.黏弹TTI介质中旋转交错网格高阶有限差分数值模拟.地球物理学报[J],2012,55(4):1354-1365.YAN Hong-Yong,LIU Yang.Rotated staggered grid high-order finite-difference numerical modeling for wave propagationin viscoelastic TTI media.Chinese J.Geophys.[J](inChinese),2012,55(4):1354-1365.
[5]祝贺君,张伟,陈晓非.二维各向异性介质中地震波场的高阶同位网格有限差分模拟.地球物理学报[J],2009,52(6):1536-1546.ZHU He-Jun,ZHANG Wei,CHEN Xiao-Fei.Two-dimensional seismic wave simulation in anisotropic media bynon-staggered finite difference method.Chinese J.Geophys.[J](in Chinese),2009,52(6):1536-1546.
[6]程冰洁,李小凡,龙桂华.基于广义正交多项式褶积微分算子的地震波场数值模拟方法.地球物理学报[J],2008,51(2):531-537.CHENG Bing-Jie,LI Xiao-Fan,LONG Gui-Hua.Seismicwaves modeling by convolutional Forsyte polynomialdifferentiator method.Chinese J.Geophys.[J](in Chinese),2008,51(2):531-537.
[7]李一琼,李小凡,朱童.基于辛格式奇异核褶积微分算子的地震标量波场模拟.地球物理学报[J],2011,54(7):1827-1834.LI Yi-Qiong,LI Xiao-Fan,ZHU Tong.The seismic scalarwave field modeling by symplectic scheme and singular kernelconvolutional differentiator.Chinese J.Geophys.[J](inChinese),2011,54(7):1827-1834.
[8]李一琼,李小凡,张美根.基于辛格式离散奇异褶积微分算子的弹性波场模拟.地球物理学报[J],2012,(5):1725-1731.LI Yi-Qiong,LI Xiao-Fan,ZHANG Mei-Gen.The elasticwave fields modeling by symplectic discrete singularconvolution differentiator method.Chinese J.Geophys.[J](inChinese),2012,(5):1725-1731.
[9]龙桂华,李小凡,张美根.错格傅里叶伪谱微分算子在波场模拟中的应用.地球物理学报[J],2009,52(1):193-199.LONG Gui-Hua;LI Xiao-Fan;ZHANG Mei-Gen.Theapplication of staggered-grid Fourier pseudospectraldifferentiation operator in wavefield modeling.Chinese J.Geophys.[J](in Chinese),2009,52(1):193-199.
[10]张怀,周元泽,吴忠良,严珍珍,等.福州盆地强地面运动特征的有限元数值模拟.地球物理学报[J],2009,52(5):1270-1279.ZHANG Huai,ZHOU Yuan-Ze,WU Zhong-Liang,YANZhen-Zhen,et al.Finite element analysis of seismic wavepropagation characteristics in Fuzhou basin.Chinese J.Geophys.[J](in Chinese),2009,52(5):1270-1279.
[11]管西竹,符力耘,陶毅,于更新.复杂地表边界元-体积元波动方程数值模拟.地球物理学报[J],2011,54(9):2357-2367.GUAN Xi-Zhu,FU Li-Yun,TAO Yi,YU Geng-Xin.Boundary-volume integral equation numerical modeling forcomplex near surface.Chinese J.Geophys.[J].(in Chinese),2011,54(9):2357-2367.
[12]朱多林,白超英.基于波动方程理论的地震波场数值模拟方法综述[J].地球物理学进展,2011,26(5):1588-1599.Zhu D L,Bai C Y.Review on the seismic wavefield forwardmodeling[J].Progress in Geophysics(in Chinese),2011,26(5):1588-1599.
[13]李信富,李小凡,张美根.地震波数值模拟方法研究综述[J].防灾减灾工程学报,2007,27(2):241-248.Li X F,Li X F,Zhang M G.Review of seismic wavenumerical modeling methods[J].Journal of DisasterPrevention and Mitigation Engineering(in Chinese),2007,27(2):241-248.
[14]王雪秋,孙建国.地震波有限差分数值模拟框架下的起伏地表处理方法综述[J].地球物理学进展,2008,23(1):40-48.Wang X Q,Sun J G.The state-of-the-art in numericalmodeling including surface topography with finite differencemethod[J].Progress in Geophysics(in Chinese),2008,23(1):40-48.
[15]张华,李振春,韩文功.起伏地表条件下地震波数值模拟方法综述[J].勘探地球物理进展,2007,30(5):1-6.Zhang H,Li Z C,Han W G.Review of seismic wavenumerical simulation form irregular topography[J].Progressin Exploration Geophysics(in Chinese),2007,30(5):1-6.
[16]Versteeg H.,Malalasekra W.An introduction to computationalfluid dynamics:The finite volume method(2nd Edition)[M].Prentice Hall,2007.
[17]刘儒勋,舒其望.计算流体力学的若干新方法[M].科学出版社,北京,2003.Liu R X,Shu C W.Selected new numerical methods forcomputational fluid dynamics(in Chinese)[M].SciencePress,Beijing,2003.
[18]Kser M,Dumbser M.An arbitrary high order discontinuousGalerkin method for elastic waves on unstructured meshes I:The two-dimensional isotropic case with external sourceterms[J].Geophys.J.Int.,2006,166:855-877.
[19]Reed W H,Hill T R,Triangular mesh methods for theneutron transport equation[R].Los Alamos ScientificLaboratory Report,LA-UR-73-479,1973.
[20]LeSaint P,Raviart P A,On a finite element method for solvingthe Neutron Transport Equation.In mathematical aspects offinite elements in partial differential equations[C].AcademicPress,New York,1974:89-145.
[21]Johnson C,Pitkaranta J.An Analysis of the discontinuousGalerkin method for a scalar hyperbolic equation[J].Math.Comp.,1986,46:1-26.
[22]Babüska I,Zlámal M.Nonconforming elements in the finiteelement method with penalty[J].SIAM J.Num.Anal.,1973,10:863-875.
[23]Babuka I,The finite element method with penalty[J].Math.Comp.,1973,27:221-228.
[24]Baker G A,Finite element methods for elliptic equations usingnonconforming elements[J].Math.Comp.,1977,31:45-59.
[25]Wheeler M F.An elliptic collocation-finite element methodwith interior penalties[J].SIAM J.Numer.Anal.,1978,15:152-161
[26]Arnold D N.An interior penalty finite element method withdiscontinuous elements[J].SIAM J.Numer.Anal.,1982,19:742-760.
[27]Arnold D N,Berzzi F,Cockburn B,Marini L D.Unifiedanalysis of discontinuous Galerkin mehtod for ellipticproblems[J].SIAM J.Numer.Anal.,2002,39(5):1749-1779.
[28]Chavent G,Salzano G.A finite element method for the 1Dwater flooding problem with gravity[J].J.Comput.Phy.,1982,45:307-344.
[29]Chavent G,Cockburn B.The local projection p0p1-discontinuousGalerkin finite element for scalar conservation laws[J].M2AN,1989,23(4):565-592.
[30]Cockburn B,Shu C W.The Runge-Kutta local projection p1-discontinuous Galerkin method for scalar conservation laws[J].M2AN,1991,25:337-361.
[31]Shu C W,Osher S.Efficient implementation of essentiallynon-oscillatory shockcapturing schemes[J].J.Comput.Phys.,1988,77:439-471.
[32]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws II:General Framework[J].Mathematics ofComputation,1989,52(186):411-435.
[33]Cockburna B,Lina S Y,Shu C W.TVB Runge-Kutta localprojection discontinuous Galerkin finite element method forscalar conservation laws III:One-dimensional systems[J].Journal of Computational Physics,1989,84(1):90-113.
[34]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws IV:The multidimensional case[J].Mathematics of Computation,1990,54(190):545-581.
[35]Cockburn B,Shu C W.TVB Runge-Kutta local projectiondiscontinuous Galerkin finite element method for scalarconservation laws V:Multidimensional systems[J].J.Comput.phys.,1998,141:199-224.
[36]Bassi F,Rebay S.Ahigh-or der accurate discontinuous finiteelement method for the numerical solution of the compressibleNavier-Stokes equations[J],J.Comput.Phys.,1997,131:267-279.
[37]Cockburn B,Shu C W.The local discontinuous Galerkinfinite element method for convection-diffusion systems[J].SIAM journal on numerical analysis,1998,35:2240-2463.
[38]Yan J,Shu C W.A local discontinuous Galerkin method forKDV type equations[J].SIAM J.Numer.Anal.,2002,40:769-791.
[39]于春丽.Allen-Cahn方程的局部间断Galerkin有限元方法[D],山东大学硕士学位论文,2009.Yu Chunli.Local discontinuous Galerkin finite elementmethods for Allben-Cahn equation[D].ShandongUniversity,Shandong,2009.
[40]Babuka I,Baumann C E,Oden J T.A discontinuous hpfinite element method for diffusion problems:1-D analysis[J].Comput.Math.Appl.,1999,37:103-122.
[41]Oden J T,Babuka I,Baumann C E.A discontinuous hpfinite element method for diffusion problems[J].J.Comput.Phys.,1998,146:491-519.
[42]Cockburn B,Karniadkis G E,Shu C W.The development ofDiscontinuous Galerkin Methods,in Discontinuous GalerkinMethods:Theory,Computation and Applications[M].Springer-Verlag,NY,2000.
[43]Rivière B,Wheeler M F.Discontinuous finite element methodsfor acoustic and elastic wave problems.part I:Semidiscreteerror estimates[R].TICAM report,2001,0:1-12.
[44]Qiu J,Shu C W.Hermite WENO schemes and their applicationas limiters for Runge-Kutta discontinuous Galerkin methodII:two dimensional case[J].Comput.Fluids,2005,34:642-663.
[45]Luo H,Baum J D,Lhner R.A Hermite WENO-based limiterfor discontinuous Galerkin method on unstructured grids[J].J.Comput.Phys.,2007,225:686-713.
[46]Hesthaven J S,Warburton T.High-order nodal methods onunstructured grids,I.Time-domain solution of Maxwell’sequations[J].J.Comput.Phys.,2002,181:186-221.
[47]Cockburn B,Li F,Shu C W.Local divergence-free discontinuousGalerkin methods for the Maxwell equations[J].J.Comput.Phys.,2004,194:588-610.
[48]Fezoui L,Lanteri S,Lohrengel S,Piperno S.Convergenceand stability of a discontinuous Galerkin time-domain methodfor the 3Dheterogeneous Maxwell equations on unstructuredmeshes[J].M2AN,Math.Model Numer.Anal.,2005,39:1149-1176.
[49]Grote M,Schneebeli A,Schtzau D.Discontinuous Galerkinfinite element method for the wave equation[J].SIAM J.Num.Anal.,2006,44:2408-2431.
[50]Monk P,Richter G R.A discontinuous Galerkin method forlinear symmetric hyperbolic systems in inhomogeneous media[J].J.Sci.Comput.,2005,22:443-477.
[51]Chung E T,Engquist B.Optimal discontinuous Galerkinmethods for wave propagation[J].SIAM J.Sci.Comput.2006,44:2131-2158.
[52]Etienne V,Chaljub E,Virieux J,Glinsky N.An hp-adaptivediscontinuous Galerkin finite-element method for 3-D elasticwave modelling[J].Geophys.J.Int.,2010,183:941-962
[53]Komatitsch D,Martin R.An unsplit convolutional perfectlymatched layer improved at grazing incidence for the seismicwave equation[J].Geophysics,2007,72(5):SM155-SM167.
[54]Peraire J,Persson P O.The compact discontinuous Galerkin(CDG)method for elliptic problems[J].SIAM J.Sci.Comput.,2008,30(4):1806-1824.
[55]Aagaard B T,Hall J F,Heaton T H.Characterization ofnearsource ground motion with earthquake simulations[J].Earthq.Spectra.,2001,17:177-207.
[56]Akcelik V,et al.High resolution forward and inverse earthquakemodeling on terascale computers,in SC’03:Proceedings ofthe 2003 ACM/IEEE Conference on Supercomputing[C].IEEE Computer Society,Washington,DC,USA,2003,52.
[57]Ichimura T,Hori M,Kuwamoto H.Earthquake motionsimulation with multi-scale finite element analysis on hybridgrid[J].Bull.seism.Soc.Am.,2007,97(4):1133-1143.
[58]Chin-Joe-Kong M J S,Mulder W A,Van Veldhuizen M.Higherorder triangular and tetrahedral finite elements withmass lumping for solving the wave equation[J].J.Eng.Math.,1999,35:405-426.
[59]Chaljub E,Komatitsch D,Vilotte J P,Capdeville Y,ValetteB,Festa G.Spectral-element analysis in Seismology[J].Advances in Geophysics,2007,48:365-419.
[60]Komatitsch D,Tromp J.Spectral-element simulations of globalseismic wave propagation-II.3-D models,oceans,rotation,and self-gravitation[J].Geophys.J.Int.,2002,150:303-318.
[61]Komatitsch D,Labarta J,Michéa D.A simulation of seismicwave propagation at high resolution in the inner core of theEarth on 2166 processors of MareNostrum[C],LectureNotes Comput.Sci.,2008,5336:364-377.
[62]Mercerat E D,Vilotte J P,Sanchez-Sesma F J.Triangularspectral element simulation of two-dimensional elastic wavepropagation using unstructured triangular grids[J].Geophys.J.Int.,2006,166:679-698.
[63]Pasquetti R,Rapetti F.Spectral element methods onunstructured meshes:comparisons and recent advances[J].J.Scient.Comput.,2006,27:377-387.
[64]Hughes T J R.The finite element method-linear static anddynamic finite element analysis[M].Prentice Hall,1987.
[65]Igel H.Wave propagation in three-dimensional spherical sectionsby the Chebyshev spectral method[M].Geophys.J.Int.,1996,136:559-566.
[66]Thompson L L,Pinsky P.A space-time finite element methodfor structural acoustics in infinite domains.Part 1:formulation,stability and convergence[J].Comp.Meth.Appl.Mech.Eng.,1996,132:195-227.
[67]Dumbser M,Schwartzkopff T,Munz C D.Arbitrary highorder finite volume schemes for linear wave propagation,toappear in Notes on Numerical Fluid Mechanics andMultidisciplinary Design(NNFM)[C],Springer,2005.
[68]Toro E F,Millington A C,Nejad L A.Towards very highorder Godunov schemes,in Godunov methods:Theory andapplications[M].Kluwer Academic/Plenum Publishers,Oxford,2001:907-940.
[69]Toro E F,Titarev V A.Solution of the generalized Riemannproblem for advection-reaction equations[J].Proc.Roy.Soc.London,2002:271-281.
[70]Schwartzkopff T,Munz C D,Toro E F.ADER:A high-Order approach for linear hyperbolic systems in 2D[J].J.Sci.Comput.,2002,17:231-240.
[71]Schwartzkopff T,Dumbser M,Munz C D.Fast high orderADER schemes for linear hyperbolic equations[J].J.Comput.Phys.,2004,197:532-539.
[72]Dumbser M,Munz C D.ADER discontinuous Galerkin schemesfor Aeroacoustics[J].Comptes Rendus-Mecanique,2005,333:683-687.
[73]Dumbser M,Munz C D.Arbitrary high order discontinuousGalerkin schemes,in numerical methods for hyperbolic andkinetic problems,eds.Cordier S,Goudon T,Gutnic M,Sonnendrucker E,IRMA series in mathematics andtheoretical physics[J].EMS Publishing House,2005:295-333.
[74]Butcher J C.The numerical analysis of ordinary differentialequations:Runge-Kutta and general linear methods[M].Wiley,1987.
[75]Dumbser M,Kser M.An arbitrary high order discontinuousGalerkin method for elastic waves on unstructured meshes II:The three-dimensional isotropic case[J].Geophys.J.Int.2006,167:319-336.
[76]Kser M,Dumbser M,de la Puente J,Igel H.An arbitraryhigh order discontinuous Galerkin method for elastic waves onunstructured meshes III:Viscoelastic attenuation[J].Geophys.J.Int.,2007,168:224-242.
[77]de la Puente J,Kser M,Dumbser M,Igel H.An arbitraryhigh order discontinuous Galerkin method for elastic waves onunstructured meshes IV:Anisotropy[J].Geophys.J.Int.,2007,169:1210-1228.
[78]de la Puente J,Dumbser M,Kser M,Igel H.DiscontinuousGalerkin methods for wave propagation in poroelastic media[J].Geophysics,2008,73:77-97.
[79]Dumbser M,Kser M,Toro E F.An arbitrary high orderdiscontinuous Galerkin method for elastic waves onunstructured meshes V:Local Time Stepping and p-Adaptivity[J].Geophys.J.Int.,2007,171(3):695-717.
[80]Castro C E,Kser M,Brietzke G B.Seismic waves inheterogeneous material:subcell resolution of thediscontinuous Galerkin method[J].Geophys.J.Int.2010,182:250-264.
[81]de la Puente J,Ampuero J-P,Kser M.Dynamic rupturemodeling on unstructured meshes using a discontinuousGalerkin method,J.Geophys.Res.,2009,114,B10302.
[82]Pelties C,de la Puente J,Ampuero J-P,Brietzke G B,et al.Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructuredtetrahedral meshes,J.Geophys.Res.,2012,117,B02309.
[83]Delcourte S,Fezoui L,Glinsky-Olivier N.A high-orderdiscontinuous Galerkin method for the seismic wavepropagation[J].ESAIM:Proceedings,2009,27:70-89.
[84]Piperno S,Remaki M,Fezoui L.A non-diffusive finite volumescheme for the 3DMaxwell equations on unstructured meshes[J].SIAM J.Numer.Anal.,2002,39:2089-2108.
[85]Benjemaa M,Glinsky-Olivier N,Cruz-Atienza V M,VirieuxJ,et al.Dynamic non-planar crack rupture by a finite volumemethod[J].Geophys.J.Int.,2007,171:271-285.
[86]Kser M,Dumbser M.A highly accurate discontinuousGalerkin method for complex interfaces between solids andmoving fluids[J].Geophysics,2008,73(3):23-35.
[87]Moczo P,Kristek J,Galis M,Chaljub E,Etienne V.3-Dfinite-difference,finite-element,discontinuous-Galerkin andspectral-element schemes analysed for their accuracy withrespect to P-wave to S-wave speed ratio[J].Geophys.J.Int.2011,187:1645-1667.
[88]Bernacki M,Fezoui L,Lanteri S,Piperno S.Parallel discontinuousGalerkin unstructured mesh solvers for the calculation ofthree-dimensional wave propagation problems[J].AppliedMathematical Modeling,2006,30:744-763.
[89]Kser M,Hermann V,de la Puente J.Quantitative accuracyanalysis of the discontinuous Galerkin method for seismicwave propagation[J].Geophys.J.Int.,2008,173(2):990-999.
[90]Warburton T.Rocky mountain mathematics consortium summer2008:Parallel numerical methods for partial differentialequations[OL].2008,(http://www.caam.rice.edu/timwar/RMMC/MIDG.html).
[91]Hu F,Hussaini M,Rasitarinera P.An analysis of thediscontinuous Galerkin method for wave propagationproblems[J].J.Comput.Phys.,1999,151(2):921-946.
[92]Stanescu D,Kopriva D,Hussaini M.Dispersion analysis fordiscontinuous spectral element methods[J].J.Scient.Comput.,2000,15(2):149-171.
[93]Ainsworth M.Dispersive and dissipative behaviour of high orderdiscontinuous Galerkin finite element methods[J].J.Comput.Phys.,2004,198(1):106-130.
[94]Ainsworth M,Monk P,Muniz W.Dispersive and dissipativeproperties of discontinuous Galerkin finite element methodsfor the second-order wave equation[J].J.Scient.Comput.,2006,27(1):5-40.
[95]De Basabe,J D,Sen M K.Grid dispersion and stabilitycriteria of some common finite-element methods for acousticand elastic wave equations[J].Geophysics,2007,72:T81-T95.
[96]De Basabe J D,Sen M,Wheeler M.The interior penaltydiscontinuous Galerkin method for elastic wave propagation:grid dispersion[J].Geophys.J.Int.,2008,175:83-93.
[97]龙桂华,李小凡,江东辉.基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案[J].地球物理学报,2010,53(12):2964-2971.Long G H,Li X F,Jiang D H.Accelerating seismic modelingwith staggered-grid Fourier Pseudo-spectral differentiationmatrix operator method on graphics processing unit[J].Chinese J.Geophys.(in Chinese),2010,53(12):2964-2971.
[98]刘红伟,刘洪,李博,等.起伏地表叠前逆时偏移理论及GPU加速技术[J].地球物理学报,2011,54(7):1883-1892.Liu H W,Liu H,Li B,et al.Pre-stack reverse timemigration for rugged topography and GPU accelerationtechnology[J].Chinese J.Geophys.(in Chinese),2011,54(7):1883-1892.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心