Numerical simulation of elastic wave propagation in inhomogeneous media

A generalized Haskell matrix method is proposed for modeling elastic wave propagation in inhomogeneous media based on mixed variable elastic wave equations and the exponential of matrices. The method presented can be used to accurately simulate all kinds of waves propagating, solving effectively problems of artificial reflection and incident wave excitation. This method is especially suited to simulate elastic wave propagation in laterally homogeneous and laterally inhomogeneous mixed layer models. The numerical results demonstrate the good accuracy of the method.