摘要
Inthispaper,weproposeanearlyanalyticexponentialtimedifference(NETD)methodforsolvingthe2Dacousticandelasticwaveequations.Inthismethod,weusethenearlyanalyticdiscreteoperatortoapproximatethehigh-orderspatialdifferentialoperatorsandtransformtheseismicwaveequationsintosemi-discreteordinarydifferentialequations(ODEs).Then,theconvertedODEsystemissolvedbytheexponentialtimedifference(ETD)method.WeinvestigatethepropertiesofNETDindetail,includingthestabilityconditionfor1-Dand2-Dcases,thetheoreticalandrelativeerrors,thenumericaldispersionrelationforthe2-Dacousticcase,andthecomputationalefficiency.Inordertofurthervalidatethemethod,weapplyittosimulatingacoustic/elasticwavepropagationinmultilayermodelswhichhavestrongcontrastsandcomplexheterogeneousmedia,e.g.,theSEGmodelandtheMarmousimodel.Fromourtheoreticalanalysesandnumericalresults,theNETDcansuppressnumericaldispersioneffectivelybyusingthedisplacementandgradienttoapproximatethehigh-orderspatialderivatives.Inaddition,becauseNETDisbasedonthestructureoftheLiegroupmethodwhichpreservesthequantitativepropertiesofdifferentialequations,itcanachievemoreaccurateresultsthantheclassicalmethods.
出版日期
2014年01月11日(中国期刊网平台首次上网日期,不代表论文的发表时间)