A new staggered grid finite difference scheme optimised in the space domain for the first order acoustic wave equation
Wenquan Liang 1 Xiu Wu 1 Yanfei Wang 2 5 Jingjie Cao 3 Chaofan Wu 1 Baoqing He 41 College of Resource Engineering, Longyan University, Longyan 364000, China.
2 Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China.
3 Hebei GEO University, Shijiazhuang, Hebei 050031, China.
4 Bureau of Geophysical Prospecting, China National Petroleum Corporation, Zhuozhou, Hebei 072750, China.
5 Corresponding author. Email: yfwang@mail.iggcas.ac.cn
Exploration Geophysics 49(6) 898-905 https://doi.org/10.1071/EG17088
Submitted: 2 July 2017 Accepted: 30 January 2018 Published: 27 March 2018
Abstract
Staggered grid finite difference (FD) methods are widely used to synthesise seismograms theoretically, and are also the basis of reverse time migration and full waveform inversion. Grid dispersion is one of the key problems for FD methods. It is desirable to have a FD scheme which can accelerate wave equation simulation while still preserving high accuracy. In this paper, we propose a totally new staggered grid FD scheme which uses different staggered grid FD operators for different first order spatial derivatives in the first order acoustic wave equation. We determine the FD coefficient in the space domain with the least-squares method. The dispersion analysis and numerical simulation demonstrated the effectiveness of the proposed method.
Key words: acoustic wave equation, finite difference scheme, finite difference time domain, numerical dispersion, staggered grid.
References
Bohlen, T., and Wittkamp, F., 2016, Three-dimensional viscoelastic time-domain finite-difference seismic modeling using the staggered Adams–Bashforth time integrator: Geophysical Journal International, 204, 1781–1788| Three-dimensional viscoelastic time-domain finite-difference seismic modeling using the staggered Adams–Bashforth time integrator: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 |
Di Bartolo, L., Dors, C., and Mansur, W. J., 2012, A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation: Geophysics, 77, T187–T199
| A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation:Crossref | GoogleScholarGoogle Scholar |
Etemadsaeed, L., Moczo, P., Kristek, J., Ansari, A., and Kristekova, M., 2016, A no-cost improved velocity–stress staggered-grid finite-difference scheme for modelling seismic wave propagation: Geophysical Journal International, 207, 481–511
| A no-cost improved velocity–stress staggered-grid finite-difference scheme for modelling seismic wave propagation:Crossref | GoogleScholarGoogle Scholar |
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, 840–850
| Acoustic wave equation modeling with new time-space domain finite difference operators: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 |
Ren, Z., and Liu, Y., 2015, Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-difference schemes: Geophysics, 80, T17–T40
| Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-difference schemes:Crossref | GoogleScholarGoogle Scholar |
Ren, Z., Liu, Y., and Sen, M. K., 2017, Least-squares reverse time migration in elastic media: Geophysical Journal International, 208, 1103–1125
| Least-squares reverse time migration in elastic media:Crossref | GoogleScholarGoogle Scholar |
Robertsson, J. O., Blanch, J. O., and Symes, W. W., 1994, Viscoelastic finite-difference modeling: Geophysics, 59, 1444–1456
| Viscoelastic finite-difference modeling:Crossref | GoogleScholarGoogle Scholar |
Tan, S., and Huang, L., 2014a, A staggered-grid finite-difference scheme optimized in the time-space domain for modeling scalar-wave propagation in geophysical problems: Journal of Computational Physics, 276, 613–634
| A staggered-grid finite-difference scheme optimized in the time-space domain for modeling scalar-wave propagation in geophysical problems:Crossref | GoogleScholarGoogle Scholar |
Tan, S., and Huang, L., 2014b, An efficient finite-difference method with high-order accuracy in both time and space domains for modelling scalar-wave propagation: Geophysical Journal International, 197, 1250–1267
| An efficient finite-difference method with high-order accuracy in both time and space domains for modelling scalar-wave propagation: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 |
Wang, Y., Liang, W., Nashed, Z., Li, X., Liang, G., and Yang, C., 2014, Seismic modeling by optimizing regularized staggered-grid finite-difference operators using a time-space-domain dispersion-relationship-preserving method: Geophysics, 79, T277–T285
| Seismic modeling by optimizing regularized staggered-grid finite-difference operators using a time-space-domain dispersion-relationship-preserving method:Crossref | GoogleScholarGoogle Scholar |
Wu, Z., and Alkhalifah, T., 2014, The optimized expansion based low-rank method for wavefield extrapolation: Geophysics, 79, T51–T60
| The optimized expansion based low-rank method for wavefield extrapolation:Crossref | GoogleScholarGoogle Scholar |
Yan, H., and Liu, Y., 2013, Acoustic VTI modeling and pre-stack reverse-time migration based on the time-space domain staggered-grid finite-difference method: Journal of Applied Geophysics, 90, 41–52
| Acoustic VTI modeling and pre-stack reverse-time migration based on the time-space domain staggered-grid finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Zhang, J. H., and Yao, Z. X., 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 |