An alternate algorithm of numerical simulation in anisotropic media

A possible problem of anisotropic medium at the internal interface for finite-difference method is investigated in this paper. Firstly, we show a finite-difference scheme for the 2-dimensional wave equations with 15 elastic parameters. Secondly, we mainly analyse the stress-continuity conditions for the finite-difference scheme at discontinuities in the 2-D transversely isotropic medium and suggest an alternate simulation technique. Finally, the seismic wavefield of theoretical model in transversely isotropic medium is simulated. The numerical computation shows that both computational results and theoretical analyse for the model are the same and that the algorithm is feasible and satisfactory. This study is supported by the National Natural Science Foundation of China.