Least squares staggered-grid finite-difference for elastic wave modelling
Lei Yang 1 2 Hongyong Yan 1 3 Hong Liu 11 Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, 100029, China.
2 University of Chinese Academy of Sciences, Beijing, 100049, China.
3 Corresponding author. Email: yanhongyong@163.com
Exploration Geophysics 45(4) 255-260 https://doi.org/10.1071/EG13087
Submitted: 17 October 2013 Accepted: 14 January 2014 Published: 24 February 2014
Abstract
Staggered-grid finite-difference (SFD) methods have been used widely in seismic wave numerical modelling and migration. The conventional way to calculate the high-order SFD coefficients on spatial derivatives is the Taylor-series expansion method, which generally leads to great accuracy on just a small frequency zone. In this paper, we first derive the SFD coefficients of arbitrary even-order accuracy for the first-order spatial derivatives by the dispersion relation and the least-squares method, which can satisfy the specific numerical solution accuracy of the derivative on a wide frequency zone. Then we use the SFD coefficients based on the least-squares to solve the first-order spatial derivatives and analyse the accuracy of the numerical solution. Finally, we perform elastic wave numerical modelling with the least-squares staggered-grid finite-difference (LSSFD) method. Meanwhile, the numerical dispersion, the modelling accuracy and the computing costs of the new method are compared with that of the Taylor-series expansion staggered-grid finite-difference (TESFD) method. The numerical dispersion analysis and elastic wavefield modelling results demonstrate that the LSSFD method can efficiently suppress the numerical dispersion and has greater modelling accuracy than the conventional TESFD method under the same discretisations and without extra computing costs.
Key words: dispersion, elastic wave, least squares, numerical modelling, staggered-grid finite-difference (SFD).
References
Bérenger, J. P., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of Computational Physics, 114, 185–200| A perfectly matched layer for the absorption of electromagnetic waves:Crossref | GoogleScholarGoogle Scholar |
Chu, C., and Stoffa, P. L., 2012a, Implicit finite-difference simulations of seismic wave propagation: Geophysics, 77, T57–T67
| Implicit finite-difference simulations of seismic wave propagation:Crossref | GoogleScholarGoogle Scholar |
Chu, C., and Stoffa, P. L., 2012b, 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 |
Dong, L. G., Ma, Z. T., Cao, J. Z., Wang, H. Z., Gao, J. H., Lei, B., and Xu, S. Y., 2000a, A staggered-grid high-order difference method of one-order elastic wave equation: Chinese Journal of Geophysics, 43, 411–419
Dong, L. G., Ma, Z. T., and Cao, J. Z., 2000b, A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation: Chinese Journal of Geophysics, 43, 856–864
Etgen, J. T., 2007, A tutorial on optimizing time domain finite-difference schemes: ‘Beyond Holberg’: Stanford Exploration Project Report, 129, 33–43.
Graves, R. W., 1996, Simulating seismic wave propagation in 3D elastic media using staggered-grid finite difference: Bulletin of the Seismological Society of America, 86, 1091–1106
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., Riollet, B., and Mora, P., 1992, Accuracy of staggered 3-D finite-difference grids for anisotropic wave propagation: 62nd Annual International Meeting, SEG, Expanded Abstracts, 1244–1246.
Kim, J. W., and Lee, D. J., 1996, Optimized compact finite difference schemes with maximum resolution: AIAA Journal, 34, 887–893
| Optimized compact finite difference schemes with maximum resolution:Crossref | GoogleScholarGoogle Scholar |
Kindelan, M., Kamel, A., and Sguazzero, P., 1990, On the construction and efficiency of staggered numerical differentiators for the wave equation: Geophysics, 55, 107–110
| On the construction and efficiency of staggered numerical differentiators for the wave equation:Crossref | GoogleScholarGoogle Scholar |
Lee, C., and Seo, Y., 2002, A new compact spectral scheme for turbulence simulations: Journal of Computational Physics, 183, 438–469
| A new compact spectral scheme for turbulence simulations:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., 2013, Globally optimal finite-difference schemes based on least squares: Geophysics, 78, T113–T132
| Globally optimal finite-difference schemes based on least squares:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2009a, A practical implicit finite-difference method: examples from seismic modeling: Journal of Geophysics and Engineering, 6, 231–249
| A practical implicit finite-difference method: examples from seismic modeling:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXivFeitbk%3D&md5=ba1dd9481927d9f8b7f4c928eea65ef5CAS |
Liu, Y., and Sen, M. K., 2009b, An implicit staggered-grid finite-difference method for seismic modeling: Geophysical Journal International, 179, 459–474
| An implicit staggered-grid finite-difference method for seismic modeling:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2011, Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes: Bulletin of the Seismological Society of America, 101, 141–159
| Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes:Crossref | GoogleScholarGoogle Scholar |
Madariaga, R., 1976, Dynamics of an expanding circular fault: Bulletin of the Seismological Society of America, 66, 639–666
Nakata, N., Tsuji, T., and Matsuoka, T., 2011, Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit: Exploration Geophysics, 42, 98–104
| Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3MXisVKrtL4%3D&md5=5b847e5a7f11d813dcd9698790c77587CAS |
Pei, Z., 2004, Numerical modeling using staggered-grid high order finite difference of elastic wave equation on arbitrary relief surface: Oil Geophysical Prospecting, 39, 629–634
Tam, C. K. W., and Webb, J. C., 1993, Dispersion-relation-preserving finite difference schemes for computational acoustics: Journal of Computational Physics, 107, 262–281
| Dispersion-relation-preserving finite difference schemes for computational acoustics:Crossref | GoogleScholarGoogle Scholar |
Tessmer, E., 2011, Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration: Geophysics, 76, S177–S185
| Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration:Crossref | GoogleScholarGoogle Scholar |
Virieux, J., 1984, SH-wave propagation in heterogeneous media: velocity stress finite-difference method: Geophysics, 49, 1933–1942
| SH-wave propagation in heterogeneous media: velocity stress finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Virieux, J., 1986, P-SV wave propagation in heterogeneous media: velocity stress finite difference method: Geophysics, 51, 889–901
| P-SV wave propagation in heterogeneous media: velocity stress finite difference method:Crossref | GoogleScholarGoogle Scholar |
Yan, H., and Liu, Y., 2012, High-order finite-difference numerical modeling of wave propagation in viscoelastic TTI media using rotated staggered grid: Chinese Journal of Geophysics, 55, 252–265
| High-order finite-difference numerical modeling of wave propagation in viscoelastic TTI media using rotated staggered grid:Crossref | GoogleScholarGoogle Scholar |
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., and Yao, Z., 2013, 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, H., and Zhang, G., 2011, Prefactored optimized compact finite-difference schemes for second spatial derivatives: Geophysics, 76, WB87–WB95
| Prefactored optimized compact finite-difference schemes for second spatial derivatives:Crossref | GoogleScholarGoogle Scholar |