Open Access Articles

Contact Us

Ausasia Science and Technology Press, Australia

Press address: Unit 3/3 Peace St, Box Hill South, Melbourne, Victoria 3128, Australia

Press email:

Press phone: 0421672696

Journal office email:

Our Visitors

269908 visitors Look details

Numerical modeling of elastic wave in frequency-domain by using staggered grid fourth-order finite-difference scheme

Chao Ma, Yan Gao, Cheng Lu

(Published: 2019-12-01)

Corresponding Author and Email:Cheng Lu,; ORCID:

Citation:Ma, C., Gao, Y., Lu, C. Numerical modeling of elastic wave in frequency-domain by using staggered grid fourth-order finite-difference scheme. Advances in Geo-Energy Research, 2019, 3(4): 410-423, doi: 10.26804/ager.2019.04.08.

Article Type:Original article


Simulation of elastic wave propagation is an important method for oil and gas exploration. Accuracy and efficiency of elastic wave simulation in complex geological environment are always the focus issue. In order to improve the accuracy and efficiency in numerical modeling, a staggered grid fourth-order finite-difference scheme of modeling elastic wave in frequency-domain is developed, which can provide stable numerical solution with fewer number of grid points per wavelength. The method is implemented on first-order velocity-stress equation and a parsimonious spatial staggered-grid with fourth-order approximation of the first-order derivative operator. Numerical tests show that the accuracy of the fourth-order staggered-grid stencil is superior to that of the mixed-grid and other conventional finite difference stencils, especially in terms of shear-wave phase velocity. Measures of mass averaging acceleration and optimization of finite difference coefficients are taken to improve the accuracy of numerical results. Meanwhile, the numerical accuracy of the finite difference scheme can be further improved by enlarging the mass averaging area at the price of expanding the bandwidth of the impedance matrix that results in the reduction of the number of grid points to 3 per shear wavelength and computer storage requirement in simulation of practical models. In our scheme, the phase velocities of compressional and shear wave are insensitive to Poisson's ratio does not occur conventional finite difference scheme in most cases, and also the elastic wave modeling can degenerate to acoustic case automatically when the medium is pure fluid or gas. Furthermore, the staggered grid scheme developed in this study is suitable for wave propagation modeling in media with coupling fluid-solid interfaces that are not resolved for previous finite difference method.

Keywords:Staggered grid, frequency-domain, finite-difference scheme, impedance matrix, heterogeneous medium.


Download this Manuscript:PDF

Copyright ©2016. All Rights Reserved