Wave propagation in multi-layered viscoelastic media with dipping interfaces

Fast and efficient seismic modeling methods are usually needed for seismic inverse problems. In this work, we propose a flexible algorithm for simulating seismic wave propagation in dipping layered viscoelastic media, which has taken the frequency dispersion into account. Since the vector wave equations are independent on the coordinate system, we study the wave propagation in the constant and local coordinate systems alternately. Firstly, for single dipping interface, the reflection and transmission coefficients of plane wave and the expressions of secondary waves are derived by coordinate transforming. Then the reflection and transmission coefficients of multi-layered media with dipping interfaces are obtained by recursive approach. Lastly, the wavefield for point source is synthesized with a fast integral method. A novel integral path is chose through comparison. Analysis shows this method is stable for all frequencies. Synthetic records with our method are consistent with the results of finite difference method, and our method is faster.