An optimized 2D/3D finite-difference seismic wave propagator using rotated staggered grid for complex elastic anisotropic structures

The synergy of computing power and physical simulation has enabled deeper insights into geological processes and properties. In geophysics, seismic anisotropy is one such crucial yet often oversimplified property, where wave propagation velocity varies with direction. To simulate the seismic wave propagation in complex anisotropic media, the rotated staggered grid (RSG) scheme was introduced decades ago. However, publicly available software for this purpose has been scarce. To address this gap, we present a newly implemented wave solver, integrated with the open-source finite-difference package Devito, that supports the simulation of seismic wave propagation in both 2D and 3D complex anisotropic media at variable spatial orders. Our implementation includes strategies to mitigate checkerboard artifacts and optimizations to reduce the number of derivative operations, thereby enhancing performance and efficiency. This wave solver aims to assist the geophysics community in more accurately modeling seismic wave propagation in intricate materials, ultimately improving our understanding of geological processes, physical properties of earth material, and subsurface structures.