摘要
在地震波正演模拟中,常采用有限差分法(finite difference method,FDM)、伪谱法(pseudo-spectral method,PSM)等方法。为研究2种方法的特性,对二维声波方程数值模拟中的有限差分法和伪谱法在计算效率、计算精度、内存占用方面做了详细的对比。发现2种方法均可达到高精度,有限差分法在计算效率上占优,伪谱法占用的内存更小。
To study the characters of finite difference method and pseudo-spectral method in numerical simulation of seismic waves,the two methods are compared in the areas of computational efficiency,computational accuracy and memory usage.It is found that both methods can achieve high accuracy.Finite difference method is good at computational efficiency while pseudo-spectral method takes up less memory.
引文
[1]MITCHELL A R,GRIFFITHS D F.The finite difference method in partial differential equations[M].Hoboken:John Wiley,1980.
[2]ALTERMAN Z,KARAL F C.Propagation of elastic waves in layered media by finite difference methods[J].Bulletin of the Seismological Society of America,1968,58(1):367-398.
[3]SUN G,TRUEMAN C W.Analysis and numerical experiments on the numerical dispersion of two-dimensional ADI-FDTD[J].IEEE Antennas&Wireless Propagation Letters,2003,2(1):78-81.
[4]李锐坚.VSP高阶有限差分正演模型[J].西安石油大学学报(自然科学版),1996(2):38-41.LI Ruijian.High order finite forward modeling for VSPof acoustic wave equation[J].Journal of Xi’an Petroleum Institute,1996(2):38-41.(in Chinese)
[5]ALFORD R M,KELLY K R,BOORE D M.Accuracy of finite-difference modeling of the acoustic wave equation[J].Geophysics,2012,39(6):834-842.
[6]DAN D K,BAYSAL E.Forward modeling by a Fourier method[J].Geophysics,1982,47(10):1402-1412.
[7]FORNBERG B.The pseudospectral method:comparisons with finite differences for the elastic wave equation[J].Geophysics,1986,52(4):483.
[8]LIU Q H.The PSTD algorithm:a time-domain method requiring only two cells per wavelength[J].Microwave&Optical Technology Letters,2015,15(3):158-165.
[9]ENGQUIST B,RUNBORG O.Computational high frequency wave propagation[J].Acta Numerica,2003,12(12):181-266.
[10]孙卫涛.弹性波动方程的有限差分数值方法[M].北京:清华大学出版社,2009.SUN Weitao.Finite difference numerical method for elastic wave equations[M].Beijing:Tsinghua University Press,2009.(in Chinese)