The Calculation Of The Seismic Wave-field With High-order Staggered-grid Finite-difference Scheme

Finite-difference method is an extensively applied method in seismic wave-field computing, because of its characteristics of the quick calculation speed and high accuracy and wide application. In this paper I make a brief summary of the dominating computing methods and research results, and compare the merits and defects of the methods. Based on this, I choose high-order staggered-grid finite-difference scheme to do the calculation of the seismic wave-field.In this paper I derive completely the any order accurate time and space formulation of 2D staggered-grid finite-difference scheme based on elastic equations and Taylor formula. The scheme has a high accuracy and can be adopted model any media with arbitrarily distributed material properties. I also analyze its stability, dispersion properties and boundary condition. I model the elastic wave for both isotropic and anisotropic media based on elastic equation, the results can reflect the characters of dynamics and kinematics of waves in seismic wave-field. According to the practical status, I make seismic modeling and reverse-time migration based on the acoustic equation.Numerical tests for both seismic modeling and reverse-time migration are performed on various theoretical models and a model based on a practical structure. The test results demonstrate that the high-order staggered-grid finite-difference scheme can be efficient and sufficiently accurate. It may also be an effective and robust tool in wide areas involving the numerical computation of partial differential equations.