Wave front evolution in strongly heterogeneous layered media using the fast marching method

The fast marching method (FMM) is a grid based numerical scheme for tracking the evolution of monotonically advancing interfaces via finite-difference solution of the eikonal equation. Like many other grid based techniques, FMM is only capable of finding the first-arriving phase in continuous media; however, it distinguishes itself by combining both unconditional stability and rapid computation, making it a truly practical scheme for velocity fields of arbitrary complexity. The aim of this paper is to investigate the potential of FMM for finding later arriving phases in layered media. In particular, we focus on reflections from smooth subhorizontal interfaces that separate regions of continuous velocity variation. The method we adopt for calculating reflected phases involves two stages: the first stage initializes FMM at the source and tracks the incident wave front to all points on the reflector surface; the second stage tracks the reflected wave front by reinitializing FMM from the interface point with minimum traveltime. Layer velocities are described by a regular grid of velocity nodes and layer boundaries are described by a set of interface nodes that may be irregularly distributed. A triangulation routine is used to locally suture interface nodes to neighbouring velocity nodes in order to facilitate the tracking of wave fronts to and from the reflector. A number of synthetic tests are carried out to assess the accuracy, speed and robustness of the new scheme. These include comparisons with analytic solutions and with solutions obtained from a shooting method of ray tracing. The convergence of traveltimes as grid spacing is reduced is also examined. Results from these tests indicate that wave fronts can be accurately tracked with minimal computational effort, even in the presence of complex velocity fields and layer boundaries with high curvature. Incident wave fronts containing gradient discontinuities or shocks also pose no difficulty. Further development of the wave front reinitialization scheme should allow other later arrivals such as multiples to be successfully located.