Elastic reverse time migration based on nested triangular mesh in 2D isotropic media

Abstract Reverse time migration (RTM) based on two-way wave equation is an important imaging tool, which has been widely used in exploration geophysics. It offers distinct advantages in handling complex geological structures. Moreover, elastic reverse time migration (ERTM) can produce more physically meaningful images for subsurface and has received increasing attention in recent years. However, rugged topography presents challenges in the reconstruction of subsurface structures for ERTM. To overcome this problem, we introduce the grid method for topography ERTM, which is an unstructured mesh based algorithm. However, the conventional gird method for RTM is implementing by discretizing model into small meshes, leading to high computational cost. To address this problem, we introduce the grid method based on nested unstructured meshes to implement ERTM with complex surface topography. Low-order small meshes, due to their fine scale, can accurately depict complex surface topography. However, their generation process is relatively time-consuming. In contrast, high-order large meshes contain many similar small meshes internally, which can be processed in parallel, resulting in higher computational efficiency. Therefore, we propose the following strategy: for the receiver wavefield, several layers on surface discreted by the small meshes are used to place receivers and remaining layers discreted by the large meshes to improve computational efficiency. For the source wavefield, only the large meshes are used because the spacing of shot is relatively large. In this way, we ensure the accurate description of complex surface topography for the implementation of ERTM with high computational efficiency. Finally, the effectiveness of our method is validated through experiments on three models with complex topography.