To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.
Traditional AVO forward modeling only considers the impact of reflection coefficients at the interface on seismic wave field amplitude and ignores various propagation effects. Introducing wave propagation effects including geometric spreading, transmission loss, attenuation into seismic wave propagation, multi-wave amplitude-preserved AVO forward modeling for horizontally layered media based on ray theory is proposed in this paper. We derived the multi-wave geometric spreading correction formulas for horizontally layered media in order to describe the geometric spreading effect of multi-wave propagation. Introducing the complex traveltime directly, we built the relationship between complex traveltime and quality factor without the help of complex velocity to describe the attenuation of viscoelastic media. Multi-wave transmission coefficients, obtained by solving the Zoeppritz equations directly, is used to describe the transmission loss. Numerical results show that the effects of geometric spreading, attenuation, and transmission loss on multi-wave amplitude varies with offset and multi-wave amplitude-preserved AVO forward modeling should consider the reconstructive effect of wave propagation on reflection amplitude.
