An implicit spatial and high-order temporal finite difference scheme for 2D acoustic modelling
Enjiang Wang 1 2 Yang Liu 1 2 3
+ Author Affiliations
- Author Affiliations
1 State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China.
2 CNPC Key Laboratory of Geophysical Prospecting, China University of Petroleum, Beijing 102249, China.
3 Corresponding author. Email: wliuyang@vip.sina.com
Exploration Geophysics 49(2) 187-201 https://doi.org/10.1071/EG16094
Submitted: 27 July 2016 Accepted: 1 December 2016 Published: 16 January 2017
Abstract
The finite difference (FD) method exhibits great superiority over other numerical methods due to its easy implementation and small computational requirement. We propose an effective FD method, characterised by implicit spatial and high-order temporal schemes, to reduce both the temporal and spatial dispersions simultaneously. For the temporal derivative, apart from the conventional second-order FD approximation, a special rhombus FD scheme is included to reach high-order accuracy in time. Compared with the Lax-Wendroff FD scheme, this scheme can achieve nearly the same temporal accuracy but requires less floating-point operation times and thus less computational cost when the same operator length is adopted. For the spatial derivatives, we adopt the implicit FD scheme to improve the spatial accuracy. Apart from the existing Taylor series expansion-based FD coefficients, we derive the least square optimisation based implicit spatial FD coefficients. Dispersion analysis and modelling examples demonstrate that, our proposed method can effectively decrease both the temporal and spatial dispersions, thus can provide more accurate wavefields.
Key words: acoustic wave equation, high-order finite difference, implicit scheme, least squares, spatial dispersion, Taylor-series expansion, temporal dispersion.
References
Chang, S., and Liu, Y., 2013, A truncated implicit high-order finite-difference scheme combined with boundary conditions: Applied Geophysics, 10, 53–62| A truncated implicit high-order finite-difference scheme combined with boundary conditions:Crossref | GoogleScholarGoogle Scholar |
Chen, J., 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., 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, An implicit finite-difference operator for the Helmholtz equation: Geophysics, 77, T97–T107
| An implicit finite-difference operator for the Helmholtz equation:Crossref | GoogleScholarGoogle Scholar |
Dablain, M. A., 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 |
De Hoop, A. T., 1960, A modification of Cagniard’s method for solving seismic pulse problems: Applied Scientific Research Section B, 8, 349–356
| A modification of Cagniard’s method for solving seismic pulse problems:Crossref | GoogleScholarGoogle Scholar |
Du, X., Fowler, P. J., and Fletcher, R. P., 2014, Recursive integral time-extrapolation methods for waves: A comparative review: Geophysics, 79, T9–T26
| Recursive integral time-extrapolation methods for waves: A comparative review:Crossref | GoogleScholarGoogle Scholar |
Fang, G., Fomel, S., Du, Q., and Hu, J., 2014, Lowrank seismic-wave extrapolation on a staggered grid: Geophysics, 79, T157–T168
| Lowrank seismic-wave extrapolation on a staggered grid:Crossref | GoogleScholarGoogle Scholar |
Finkelstein, B., and Kastner, R., 2007, Finite difference time domain dispersion reduction schemes: Journal of Computational Physics, 221, 422–438
| Finite difference time domain dispersion reduction schemes:Crossref | GoogleScholarGoogle Scholar |
Fomel, S., Ying, L., and Song, X., 2013, Seismic wave extrapolation using lowrank symbol approximation: Geophysical Prospecting, 61, 526–536
| Seismic wave extrapolation using lowrank symbol approximation:Crossref | GoogleScholarGoogle Scholar |
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 |
Kosloff, D. D., and Baysal, E., 1982, Forward modeling by a Fourier method: Geophysics, 47, 1402–1412
| Forward modeling by a Fourier method:Crossref | GoogleScholarGoogle Scholar |
Kosloff, D., Filho, A. Q., Tessmer, E, and Behle, A, 1989, Numerical solution of the acoustic and elastic wave equations by a new rapid expansion method: Geophysical Prospecting, 37, 383–394
| Numerical solution of the acoustic and elastic wave equations by a new rapid expansion method:Crossref | GoogleScholarGoogle Scholar |
Lax, P., and Wendroff, B., 1960, Systems of conservation laws: Communications on Pure and Applied Mathematics, 13, 217–237
| Systems of conservation laws:Crossref | GoogleScholarGoogle Scholar |
Lele, S. K., 1992, Compact finite difference schemes with spectral-like resolution: Journal of Computational Physics, 103, 16–42
| Compact finite difference schemes with spectral-like resolution:Crossref | GoogleScholarGoogle Scholar |
Levander, A. R., 1988, Fourth-order finite-difference P-SV seismograms: Geophysics, 53, 1425–1436
| Fourth-order finite-difference P-SV seismograms: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., 2014, Optimal staggered-grid finite-difference schemes based on least-squares for wave equation modelling: Geophysical Journal International, 197, 1033–1047
| Optimal staggered-grid finite-difference schemes based on least-squares for wave equation modelling:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2009a, 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 | 1:CAS:528:DC%2BD1MXht1GnsLvL&md5=7ce9f8d8c31ba8af95d6592ae7b10951CAS |
Liu, Y., and Sen, M. K., 2009b, A practical implicit finite-difference method: examples from seismic modelling: Journal of Geophysics and Engineering, 6, 231–249
| A practical implicit finite-difference method: examples from seismic modelling:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXivFeitbk%3D&md5=36d3f56dc8893848e6acea06640f3881CAS |
Liu, Y., and Sen, M. K., 2009c, An implicit staggered-grid finite-difference method for seismic modelling: Geophysical Journal International, 179, 459–474
| An implicit staggered-grid finite-difference method for seismic modelling:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2010, A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation: Geophysics, 75, A1–A6
| A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2013, Time–space domain dispersion-relation-based finite-difference method with arbitrary even-order accuracy for the 2D acoustic wave equation: Journal of Computational Physics, 232, 327–345
| Time–space domain dispersion-relation-based finite-difference method with arbitrary even-order accuracy for the 2D acoustic wave equation:Crossref | GoogleScholarGoogle Scholar |
Liu, H., Dai, N., Niu, F., and Wu, W., 2014, An explicit time evolution method for acoustic wave propagation: Geophysics, 79, T117–T124
| An explicit time evolution method for acoustic wave propagation:Crossref | GoogleScholarGoogle Scholar |
Long, G., Zhao, Y., and Zou, J., 2013, A temporal fourth-order scheme for the first-order acoustic wave equations: Geophysical Journal International, 194, 1473–1485
| A temporal fourth-order scheme for the first-order acoustic wave equations:Crossref | GoogleScholarGoogle Scholar |
McMechan, G. A., 1983, Migration by extrapolation of time-dependent boundary values: Geophysical Prospecting, 31, 413–420
| Migration by extrapolation of time-dependent boundary values:Crossref | GoogleScholarGoogle Scholar |
Pestana, R. C., and Stoffa, P. L., 2010, Time evolution of the wave equation using rapid expansion method: Geophysics, 75, T121–T131
| Time evolution of the wave equation using rapid expansion method:Crossref | GoogleScholarGoogle Scholar |
Shan, G., 2009, Optimized implicit finite-difference and Fourier finite-difference migration for VTI media: Geophysics, 74, WCA189–WCA197
| Optimized implicit finite-difference and Fourier finite-difference migration for VTI media:Crossref | GoogleScholarGoogle Scholar |
Song, X., and Fomel, S., 2011, Fourier finite-difference wave propagation: Geophysics, 76, T123–T129
| Fourier finite-difference wave propagation:Crossref | GoogleScholarGoogle Scholar |
Song, X., Fomel, S., and Ying, L., 2013, Lowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation: Geophysical Journal International, 193, 960–969
| Lowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation:Crossref | GoogleScholarGoogle Scholar |
Sun, W., Zhou, B., and Fu, L., 2015, A staggered-grid convolutional differentiator for elastic wave modeling: Journal of Computational Physics, 301, 59–76
Sun, J., Fomel, S., and Ying, L., 2016, Low-rank one-step wave extrapolation for reverse time migration: Geophysics, 81, S39–S54
| Low-rank one-step wave extrapolation for reverse time migration:Crossref | GoogleScholarGoogle Scholar |
Tan, S., and Huang, L., 2014a, 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 |
Tan, S., and Huang, L., 2014b, 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 |
Tarantola, A., 1984, Inversion of seismic reflection data in the acoustic approximation: Geophysics, 49, 1259–1266
| Inversion of seismic reflection data in the acoustic approximation: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., and Operto, S., 2009, An overview of full-waveform inversion in exploration geophysics: Geophysics, 74, WCC1–WCC26
| An overview of full-waveform inversion in exploration geophysics:Crossref | GoogleScholarGoogle Scholar |
Virieux, J., Calandra, H., and Plessix, R. É., 2011, A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging: Geophysical Prospecting, 59, 794–813
| A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging: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 |
Wang, E., Liu, Y., and Sen, M. K., 2016, Effective finite-difference modelling methods with 2-D acoustic wave equation using a combination of cross and rhombus stencils: Geophysical Journal International, 206, 1933–1958
| Effective finite-difference modelling methods with 2-D acoustic wave equation using a combination of cross and rhombus stencils: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 |
Zhang, Y., and Zhang, G., 2009, One-step extrapolation method for reverse time migration: Geophysics, 74, A29–A33
| One-step extrapolation method for reverse time migration:Crossref | GoogleScholarGoogle Scholar |
Zhou, B., and Greenhalgh, S. A., 1992, Seismic scalar wave equation modeling by a convolutional differentiator: Bulletin of the Seismological Society of America, 82, 289–303
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 |
Zhou, B., Greenhalgh, S. A., 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 |