Two-dimensional time-domain finite-difference modeling for viscoelastic seismic wave propagation

Real Earth media are not perfectly elastic. Instead, they attenuate propagating mechanical waves. This anelastic phenomenon in wave propagation can be modeled by a viscoelastic mechanical model consisting of several standard linear solids. Using this viscoelastic model, we approximate a constant Q over a frequency band of interest. We use a four-element viscoelastic model with a trade-off between accuracy and computational costs to incorporate Q into 2-D time-domain first-order velocity-stresswave equations. To improve the computational efficiency, we limit the Q in the model to a list of discrete values between 2 and 1000. The related stress and strain relaxation times that characterize the viscoelastic model are precalculated and stored in a database for use by the finite-difference calculation. A viscoelastic finite-difference scheme that is second order in time and fourth order in space is developed based on theMacCormack algorithm. The newmethod is validated by comparing the numerical result with analytical solutions that are calculated using the generalized reflection/transmission coefficient method. The synthetic seismograms exhibit greater than 95 per cent consistency in a two-layer viscoelastic model. The dispersion generated from the simulation is consistent with the Kolsky-Futterman dispersion relationship.