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

• Main Authors: Chao Ma, Yan Gao, Cheng Lu
• Format: Article
• Language: English
• Published: Yandy Scientific Press 2019-12-01
• Series: Advances in Geo-Energy Research
• Subjects: staggered grid; frequency-domain; finite-difference scheme; impedance matrix; heterogeneous medium
• Online Access: https://www.yandy-ager.com/index.php/ager/article/view/175

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.