Finite-difference frequency-domain modeling of viscoacoustic wave propagation in 2D tilted transversely isotropic (TTI) media

A 2D finite-difference, frequency-domain method was developed for modeling viscoacoustic seismic waves in transversely isotropic media with a tilted symmetry axis. The medium is parameterized by the P-wave velocity on the symmetry axis, the density, the attenuation factor, Thomsen's anisotropic parameters delta and epsilon, and the tilt angle. The finite-difference discretization relies on a parsimonious mixed-grid approach that designs accurate yet spatially compact stencils. The system of linear equations resulting from discretizing the time-harmonic wave equation is solved with a parallel direct solver that computes monochromatic wavefields efficiently for many sources. Dispersion analysis shows that four grid points per P-wavelength provide sufficiently accurate solutions in homogeneous media. The absorbing boundary conditions are perfectly matched layers (PMLs). The kinematic and dynamic accuracy of the method wasassessed with several synthetic examples which illustrate the propagation of S-waves excited at the source or at seismic discontinuities when epsilon