Abstract: Elastic wave reverse time migration with Gaussian beams is a flexible and adaptable imaging tool for seismic data, and it can image subsurface structures more accurately.It adopts elastic wave dynamic Gaussian beams to characterize the Green function, then constructs forward and backward continuation wavefields, and finally calculates the imaging results through cross correlation.As the earth is essentially anisotropic, isotropic processing will result in artificial imaging.To solve the problem of rapid divergence of Gaussian beams with increasing propagation distance, this paper proposes a method based on the elastic wave reverse time migration with Gaussian beams.An operator of dynamic focused beams is derived for elastic waves to constrain seismic energy within the effective range.An anisotropic ray-tracing algorithm based on phase velocity is used to correct the influence of anisotropy on seismic imaging.This algorithm has high computational efficiency and can mitigate the uncertainties of obtaining weakly anisotropic parameters of subsurface media.The test results for a complex VTI structural model and a TTI fault model show that the proposed method can constrain seismic energy within the effective range to improve the imaging quality of deep complex structures.