Two-dimensional Time-Domain Finite-Difference Modeling for Viscoelastic Seismic Wave Propagation

Na Fan,Lian-Feng Zhao,Xiao-Bi Xie,Zengxi Ge,Zhen-Xing Yao
DOI: https://doi.org/10.1093/gji/ggw228
IF: 3.352
2016-01-01
Geophysical Journal International
Abstract: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-stress wave 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 pre-calculated 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 the MacCormack algorithm. The new method 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.
What problem does this paper attempt to address?