Efficient modeling of wave propagation in a vertical transversely isotropic attenuative medium based on fractional Laplacian

To efficiently simulate wave propagation in a vertical transversely isotropic (VTI) attenuative medium, we have developed a viscoelastic VTI wave equation based on fractional Laplacian operators under the assumption of weak attenuation (Q≫1), where the frequency-independent Q model is used to mathematically represent seismic attenuation. These operators that are nonlocal in space can be efficiently computed using the Fourier pseudospectral method. We evaluated the accuracy of numerical solutions in a homogeneous transversely isotropic medium by comparing with theoretical predictions and numerical solutions by an existing viscoelastic-anisotropic wave equation based on fractional time derivatives. To accurately handle heterogeneous Q, we select several Q values to compute their corresponding fractional Laplacians in the wavenumber domain and interpolate other fractional Laplacians in space. We determined its feasibility by modeling wave propagation in a 2D heterogeneous attenuative VTI medium. We concluded that the new wave equation is able to improve the efficiency of wave simulation in viscoelastic-VTI media by several orders and still maintain accuracy.