The geometric optics approximation to high frequency anisotropic wave propagation reduces the anisotropic wave equation to a static Hamilton–Jacobi equation. This equation is known as the anisotropic eikonal equation and has three different coupled wave modes as solutions. We introduce here a level set-based Eulerian approach that captures all three of these wave propagations. In particular, our method is able to accurately reproduce the quasi-transverse, or quasi-S, waves with cusps, which form a class of multi-valued solutions. The level set formulation we use is borrowed from one for moving curves in three spatial dimensions, with the velocity fields for evolution following from the method of characteristics on the anisotropic eikonal equation. We present here our derivation of the algorithm and numerical results to illustrate its accuracy in different cases of anisotropic wave propagations related to seismic imaging.
Read full abstract