Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

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 - https://doi.org/10.1071/EG16094
Submitted: 27 July 2016  Accepted: 1 December 2016   Published online: 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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

Kosloff, D. D., and Baysal, E., 1982, Forward modeling by a Fourier method: Geophysics, 47, 1402–1412
Forward modeling by a Fourier method:CrossRef |

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 |

Lax, P., and Wendroff, B., 1960, Systems of conservation laws: Communications on Pure and Applied Mathematics, 13, 217–237
Systems of conservation laws:CrossRef |

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 |

Levander, A. R., 1988, Fourth-order finite-difference P-SV seismograms: Geophysics, 53, 1425–1436
Fourth-order finite-difference P-SV seismograms:CrossRef |

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 |

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 |

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 | 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 | 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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

Song, X., and Fomel, S., 2011, Fourier finite-difference wave propagation: Geophysics, 76, T123–T129
Fourier finite-difference wave propagation:CrossRef |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |

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 |



Export Citation