Comparison of numerical dispersion in acoustic finite-difference algorithms
Wen-Quan Liang 1 Yan-Fei Wang 1 2 Chang-Chun Yang 1
+ Author Affiliations
- Author Affiliations
1 Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China.
2 Corresponding author. Email: yfwang@mail.iggcas.ac.cn
Exploration Geophysics 46(2) 206-212 https://doi.org/10.1071/EG13072
Submitted: 11 October 2013 Accepted: 25 February 2014 Published: 16 April 2014
Abstract
Numerical simulation of the acoustic wave equation is widely used to synthesise seismograms theoretically, and is also the basis of the reverse time migration. With some stability conditions, grid dispersion often exists because of the discretisation of the time and the spatial derivatives in the wave equation. How to suppress the grid dispersion is therefore a key problem for finite-difference approaches. Different methods are proposed to address the problem. The commonly used methods are the high order Taylor expansion methods and the optimised methods. In this paper, we compare the performance of these methods in the space and time–space domains. We demonstrate by dispersion analysis and numerical simulation that a linear method without iteration performs comparably to the optimised methods, but with reduced computational effort.
Key words: acoustic wave-equation modelling, finite-difference scheme, optimised methods, Taylor expansion method.
References
Alford, R., Kelly, K., and Boore, D., 1974, Accuracy of finite-difference modeling of the acoustic wave equation: Geophysics, 39, 834–842| Accuracy of finite-difference modeling of the acoustic wave equation:Crossref | GoogleScholarGoogle Scholar |
Baysal, E., Kosloff, D. D., and Sherwood, J. W. C., 1983, Reverse time migration: Geophysics, 48, 1514–1524
| Reverse time migration:Crossref | GoogleScholarGoogle Scholar |
Billette, F., and Brandsberg-Dahl, S., 2005, The 2004 BP velocity benchmark, in Proceedings of the 67th EAGE Conference and Exhibition, B035.
Chen, J. B., 2007, High-order time discretizations in seismic modeling: Geophysics, 72, SM115–SM122
| High-order time discretizations in seismic modeling:Crossref | GoogleScholarGoogle Scholar |
Chen, J. B., 2011, A stability formula for Lax-Wendroff methods with fourth-order in time and general-order in space for the scalar wave equation: Geophysics, 76, T37–T42
| A stability formula for Lax-Wendroff methods with fourth-order in time and general-order in space for the scalar wave equation:Crossref | GoogleScholarGoogle Scholar |
Chu, C., and Stoffa, P. L., 2012, Determination of finite-difference weights using scaled binomial windows: Geophysics, 77, W17–W26
| Determination of finite-difference weights using scaled binomial windows:Crossref | GoogleScholarGoogle Scholar |
Dablain, M., 1986, The application of high‐order differencing to the scalar wave equation: Geophysics, 51, 54–66
| The application of high‐order differencing to the scalar wave equation:Crossref | GoogleScholarGoogle Scholar |
Daudt, C., Braile, L., Nowack, R., and Chiang, C., 1989, A comparison of finite-difference and Fourier method calculations of synthetic seismograms: Bulletin of the Seismological Society of America, 79, 1210–1230
Etgen, J. T., 2007, A tutorial on optimizing time domain finite-difference scheme: ‘Beyond Holberg’: Stanford Exploration Project Report, 129, 33–43.
Fornberg, B., 1987, The pseudospectral method: comparisons with finite differences for the elastic wave equation: Geophysics, 52, 483–501
| The pseudospectral method: comparisons with finite differences for the elastic wave equation:Crossref | GoogleScholarGoogle Scholar |
Fornberg, B., 1988, The pseudospectral method: accurate representation of interfaces in elastic wave calculations: Geophysics, 53, 625–637
| The pseudospectral method: accurate representation of interfaces in elastic wave calculations:Crossref | GoogleScholarGoogle Scholar |
Fornberg, B., 1990, High-order finite differences and the pseudospectral method on staggered grids: SIAM Journal on Numerical Analysis, 27, 904–918
| High-order finite differences and the pseudospectral method on staggered grids:Crossref | GoogleScholarGoogle Scholar |
Holberg, O., 1987, Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena: Geophysical Prospecting, 35, 629–655
| Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena:Crossref | GoogleScholarGoogle Scholar |
Igel, H., Käser, M., and Stupazzini, M., 2009, Simulation of seismic wave propagation in media with complex geometries, in W. H. K. Lee, ed., Encyclopedia of complexity and system science: Springer Verlag, 7891–7914.
Liang, W. Q., Yang, C. C., Wang, Y. F., and Liu, H. W., 2013, Acoustic wave equation modeling with new time–space domain finite difference operators: Chinese Journal of Geophysics, 56, 3497–3506
Liang, W. Q., Wang, Y. F., and Yang, C. C., 2014, Determining the finite difference weights for the acoustic wave equation by a new dispersion-relationship-preserving method: Geophysical Prospecting, ,
Liu, Y., and Sen, M. K., 2009, A new time–space domain high-order finite-difference method for the acoustic wave equation: Journal of Computational Physics, 228, 8779–8806
| A new time–space domain high-order finite-difference method for the acoustic wave equation:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2010, Acoustic VTI modeling with a time–space domain dispersion-relation-based finite-difference scheme: Geophysics, 74, T67–T73
Liu, Y., and Sen, M. K., 2011, Finite-difference modeling with adaptive variable-length spatial operators: Geophysics, 76, T79–T89
| Finite-difference modeling with adaptive variable-length spatial operators:Crossref | GoogleScholarGoogle Scholar |
Liu, H. W., Li, B., Liu, H., Tong, X. L., Liu, Q., Wang, X. W., and Liu, W. Q., 2012, The issues of prestack reverse time migration and solutions with graphic processing unit implementation: Geophysical Prospecting, 35, 629–655
Marfurt, K. J., 1984, Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations: Geophysics, 49, 533–549
| Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations:Crossref | GoogleScholarGoogle Scholar |
Moczo, P., Kristek, J, and Halada, L, 2000, 3D 4th-order staggered-grid finite-difference schemes: stability and grid dispersion: Bulletin of the Seismological Society of America, 90, 587–603
| 3D 4th-order staggered-grid finite-difference schemes: stability and grid dispersion:Crossref | GoogleScholarGoogle Scholar |
Mora, P., 1986, Elastic finite difference with convolutional operators: Stanford Exploration Project Report, 48, 277–289.
Robertsson, J., Blanch, J., and Symes, W., 1994, Viscoelastic finite‐difference modeling: Geophysics, 59, 1444–1456
| Viscoelastic finite‐difference modeling:Crossref | GoogleScholarGoogle Scholar |
Wang, Y. F., Stepanova, I. E., Titarenko, V. N., and Yagola, A. G., 2011, Inverse problems in geophysics and solution methods: Higher Education Press (Beijing).
Yan, H., Liu, Y., and Zhang, H., 2013, Prestack reverse-time migration with a time–space domain adaptive high-order staggered-grid finite-difference method: Exploration Geophysics, 44, 77–86
| Prestack reverse-time migration with a time–space domain adaptive high-order staggered-grid finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Zhang, J. H., and Yao, Z. X., 2013a, Optimized explicit finite-difference schemes for spatial derivatives using maximum norm: Journal of Computational Physics, 250, 511–526
| Optimized explicit finite-difference schemes for spatial derivatives using maximum norm:Crossref | GoogleScholarGoogle Scholar |
Zhang, J. H., and Yao, Z. X., 2013b, Optimized finite-difference operator for broadband seismic wave modeling: Geophysics, 78, A13–A18
| Optimized finite-difference operator for broadband seismic wave modeling:Crossref | GoogleScholarGoogle Scholar |
Zhou, B., and Greenhalgh, S., 1992, Seismic scalar wave equation modeling by a convolutional differentiator: Bulletin of the Seismological Society of America, 82, 289–303
Zhou, B., Greenhalgh, S., and Zhe, J., 1993, Numerical seismogram computations for inhomogeneous media using a short, variable length convolutional differentiator: Geophysical Prospecting, 41, 751–766
| Numerical seismogram computations for inhomogeneous media using a short, variable length convolutional differentiator:Crossref | GoogleScholarGoogle Scholar |