1887
Volume 60, Issue 5
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

In this paper, we propose a nearly‐analytic central difference method, which is an improved version of the central difference method. The new method is fourth‐order accurate with respect to both space and time but uses only three grid points in spatial directions. The stability criteria and numerical dispersion for the new scheme are analysed in detail. We also apply the nearly‐analytic central difference method to 1D and 2D cases to compute synthetic seismograms. For comparison, the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method are used to model acoustic wavefields. Numerical results indicate that the nearly‐analytic central difference method can be used to solve large‐scale problems because it effectively suppresses numerical dispersion caused by discretizing the scalar wave equation when too coarse grids are used. Meanwhile, numerical results show that the minimum sampling rate of the nearly‐analytic central difference method is about 2.5 points per minimal wavelength for eliminating numerical dispersion, resulting that the nearly‐analytic central difference method can save greatly both computational costs and storage space as contrasted to other high‐order finite‐difference methods such as the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2011.01033.x
2012-01-17
2024-04-25
Loading full text...

Full text loading...

References

  1. AkiK. and RichardsP.G.1980. Quantitative Seismology: Theory and Methods . W.H. Freeman & Co., San Francisco , CA .
    [Google Scholar]
  2. AlfordR.M., KellyK.R. and BooreD.M.1974. Accuracy of finite‐difference modelling of the acoustic wave equation. Geophysics39, 834–842.
    [Google Scholar]
  3. AppeloD. and PeterssonN.A.2009. A stable finite‐difference method for the elastic wave equation on complex geometries with free surfaces. Communications in Computational Physics5, 84–107.
    [Google Scholar]
  4. BlanchJ.O. and RobertssonA.1997. A modified Lax‐Wendroff correction for wave propagation in media described by Zener elements. Geophysical Journal International131, 381–386.
    [Google Scholar]
  5. DablainM.A.1986. The application of high‐order differencing to scalar wave equation. Geophysics51, 54–66.
    [Google Scholar]
  6. DuS.T.1997. Seismic Waves Dynamics. Petroleum University Press, Shandong , China (in Chinese).
    [Google Scholar]
  7. FornbergB.1990. High‐order finite‐differences and pseudo‐spectral method on staggered grids. SIAM Journal on Numerical Analysis27, 904–918.
    [Google Scholar]
  8. GellerR.J. and TakeuchiN.1998. Optimally accurate second‐order time‐domain finite‐difference scheme for the elastic equation of motion: One‐dimensional case. Geophysical Journal International135, 48–62.
    [Google Scholar]
  9. HigdonR.L.1991. Absorbing boundary conditions for elastic waves. Geophysics56, 231–241.
    [Google Scholar]
  10. IgelH., MoraP. and RiolletB.1995. Anisotropic wave propagation through finite‐difference grids. Geophysics60, 1203–1216.
    [Google Scholar]
  11. KellyK., WardR., TreitelS. and AlfordR.1976. Synthetic seismograms: A finite‐difference approach. Geophysics41, 2–27.
    [Google Scholar]
  12. KolsloffD. and BaysalE.1982. Forward modelling by a Fourier method. Geophysics56, 231–241.
    [Google Scholar]
  13. LaxP.D. and WendroffB.1964. Difference schemes for hyperbolic equations with high order of accuracy. Communications on Pure and Applied Mathematics17, 381–398.
    [Google Scholar]
  14. LeleS.K.1992. Compact finite‐difference schemes with spectral‐like resolution. Journal of Computational Physics103, 16–42.
    [Google Scholar]
  15. LisitsaV. and VishnevskiyD.2010. Lebedev scheme for the numerical simulation of wave propagation in 3D anisotropic elasticity. Geophysical Prospecting58, 619–635.
    [Google Scholar]
  16. LiuE., DobsonA., PanD.M. and YangD.H.2008. The matrix formulation of boundary integral modelling of elastic wave propagation in 2D multi‐layered media with irregular interfaces. Journal of Computational Acoustics16, 381–396.
    [Google Scholar]
  17. MarfurtK.J.1984. Accuracy of finite‐difference and finite‐element modelling of the scalar and elastic wave equations. Geophysics49, 533–549.
    [Google Scholar]
  18. MattssonK. and NordströmJ.2006. High order finite‐difference methods for wave propagation in discontinuous media. Journal of Computational Physics220, 249–269.
    [Google Scholar]
  19. MizutaniH., GellerR.J. and TakeuchiN.2000. Comparison of accuracy and efficiency of time‐domain schemes for calculating synthetic seismograms. Physics of the Earth and Planetary Interiors119, 75–97.
    [Google Scholar]
  20. MoczoP., KristekJ. and HaladaL.2000. 3D fourth‐order staggered‐grid finite‐difference schemes: Stability and grid dispersion. Bulletin of the Seismological Society of America90, 587–603.
    [Google Scholar]
  21. MoczoP., KristekJ., VavrycukV., ArchuletaR.J. and HaladaL.2002. 3D heterogeneous staggered‐grid finite‐difference modelling of seismic motion with volume harmonic and arithmetic averaging of elastic modulus and densities. Bulletin of the Seismological Society of America92, 3042–3066.
    [Google Scholar]
  22. MoczoP., RobertssonJ.O.A. and EisnerL.2007. The finite‐difference time‐domain method for modelling of seismic wave propagation. Advances in Geophysics48, 421–516.
    [Google Scholar]
  23. RobertssonJ.O.A., BlanchJ.O. and SymesW.W.1994. Viscoelastic finite‐difference modelling. Geophysics59, 1444–1456.
    [Google Scholar]
  24. SaengerE.H. and BohlenT.2004. Finite‐difference modelling of viscoelastic and anisotropic wave propagation using the rotated staggered grid. Geophysics69, 583–591.
    [Google Scholar]
  25. SaengerE.H., GoldN. and ShapiroS.A.2000. Modelling the propagation of elastic waves using a modified finite‐difference grid. Wave Motion31, 77–92.
    [Google Scholar]
  26. SantosaF. and PaoY.H.1986. Accuracy of a Lax‐Wendroff scheme for the wave equation. Journal of the Acoustical Society of America80, 1429–1437.
    [Google Scholar]
  27. SeiA. and SymesW.1994. Dispersion analysis of numerical wave propagation and its computational consequences. Journal of Scientific Computing10, 1–27.
    [Google Scholar]
  28. TakeuchiN. and GellerR.J.2000. Optimally accurate second order time‐domain finite‐difference scheme for computing synthetic seismograms in 2‐D and 3‐D media. Physics of the Earth and Planetary Interiors119, 99–131.
    [Google Scholar]
  29. VersteegR.J. and GrauG.1991. The Marmousi experience. Proceedings of the EAGE workshop on Practical Aspects of Seismic Data Inversion (Copenhagen, 1990), European Assocication of Exploration Geophysicists, Zeist .
  30. VichnevetskyR.1979. Stability charts in the numerical approximation of partial differential equations: A review. Mathematics and Computers in Simulation21, 170–177.
    [Google Scholar]
  31. VirieuxJ.1986. P‐SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method. Geophysics51, 889–901.
    [Google Scholar]
  32. YangD.H., LiuE. and ZhangZ.J.2008. Evaluation of the u‐W finite element method in anisotropic porous media. Journal of Seismic Exploration17, 273–299.
    [Google Scholar]
  33. YangD.H., LiuE., ZhangZ.J. and TengJ.2002. Finite‐difference modelling in two‐dimensional anisotropic media using a flux‐corrected transport technique. Geophysical Journal International148, 320–328.
    [Google Scholar]
  34. YangD.H., PengJ.M., LuM. and TerlakyT.2006. Optimal nearly‐analytic discrete approximation to the scalar wave equation. Bulletin of the Seismological Society of America96, 1114–1130.
    [Google Scholar]
  35. YangD.H., SongG.J., HuaB.L. and CalandraH.2010. Simulation of acoustic wave‐fields in heterogeneous media: A robust method for automatic suppression of numerical dispersion. Geophysics75, T99–T110.
    [Google Scholar]
  36. YangD.H., TengJ.W., ZhangZ.J. and LiuE.2003a. A nearly‐analytic discrete method for acoustic and elastic wave equations in anisotropic media. Bulletin of the Seismological Society of America93, 882–890.
    [Google Scholar]
  37. YangD.H., WangL. and DengX.Y.2010. An explicit split‐step algorithm of the implicit Adams method for solving 2‐D acoustic and elastic wave equations. Geophysical Journal International180, 291–310.
    [Google Scholar]
  38. YangD.H., WangS.Q., ZhangZ.J. and TengJ.W.2003b. n‐times absorbing boundary conditions for compact finite‐difference modelling of acoustic and elastic wave propagation in the 2‐D TI Medium. Bulletin of the Seismological Society of America93, 2389–2401.
    [Google Scholar]
  39. ZengY.Q. and LiuQ.H.2001. A staggered‐grid finite‐difference method with perfectly matched layers for poroelastic wave equations. Journal of the Acoustical Society of America109, 2571–2580.
    [Google Scholar]
  40. ZhangZ.J., WangG.J. and HarrisJ M.1999. Multi‐component wave‐field simulation in viscous extensively dilatancy anisotropic media. Physics of Earth and Planetary Interiors114, 25–38.
    [Google Scholar]
  41. ZhengH.S., ZhangZ.J. and LiuE.2006. Non‐linear seismic wave propagation in anisotropic media using the flux‐corrected transport technique. Geophysical Journal International165, 943–956.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2011.01033.x
Loading
/content/journals/10.1111/j.1365-2478.2011.01033.x
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Numerical dispersion; Wave equation

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error