An effective numerical method to remove the edge effects for seismic wave modeling in arbitrary viscoelastic anisotropic media

The generalized stiffness reduction method (GSRM) was previously proposed to absorb the artificial reflections in seismic wave modeling in arbitrary elastic anisotropic media. However, the traditional PML technique and the original GSRM lose their effectiveness of removing the artificial edge effects in viscoelastic anisotropic media. To overcome this issue, we propose a new version of the GSRM to attenuate the viscoelastic waves arriving the edges of the grid. The frequency-domain wave modeling is carried out in VTI media and shows the results used the PML and two GSRMs. Comparing the numerical solutions with the analytic solutions, we show perfect matchings of the wavefields by using the new version of the GSRM comparing to the use of the PML and the original GSRM, which exhibit serious deviations from the true solution. Moreover, the new version of the GSRM maintains the simple numerical implementation of the original GRSM.