Numerical Modeling of 2D Elastic Wave Propagation in Porous VTI Medium Using Implicit Rotated Staggered Grid FD Method

Summary We combine the rotated staggered grid finite difference method with the implicit staggered grid finite difference method together to perform numerical modeling of the porous elastic anisotropic medium. The combined method requires nearly the same amount of computation and occupies nearly the same amount of memory as those of the explicit staggered-grid method under the same order but can obtain higher accuracy when additional cost of visiting arrays is ignored. A high order explicit FD can be replaced by some lower order implicit FD so that lots of computational cost may be saved while maintaining the precision. Finally, we show two examples and a comparison of the precision between the rotated staggered grid finite difference and the combined method to demonstrate advantages of the combined method.