Three-dimensional travel-time computation using the Fast marching method

We present a fast algorithm for solving the eikonal equation in three dimensions, based on the Fast Marching Method (FMM). The algorithm is of order O(N log N), where N is the total number of grid points in the computational domain. The algorithm can be used in any orthogonal coordinate system, and globally constructs the solution to the eikonal equation for each point in the coordinate domain. The method is unconditionally stable, and constructs solutions consistent with the exact solution for arbitrarily large gradient jumps in velocity. In addition, the method resolves any overturning propagation wavefronts. We begin with the mathematical foundation for solving the eikonal equation using the FMM, and follow with the numerical details. We show examples of traveltime propagation through the SEG/EAGE Salt Model, and the use of these first arrival traveltimes to image 3D prestack data.

[1]  J. Vidale.  Finite‐difference calculation of traveltimes in three dimensions , 1990 .

[2]  W. Schneider.  Robust and efficient upwind finite-difference traveltime calculations in three dimensions , 1995 .

[3]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[4]  J. Sethian.  AN ANALYSIS OF FLAME PROPAGATION , 1982 .

[5]  Håvar Gjøystdal,et al.  Traveltime and amplitude estimation using wavefront construction , 1993 .

[6]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[7]  James A. Sethian,et al.  Theory, algorithms, and applications of level set methods for propagating interfaces , 1996, Acta Numerica.

[8]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[9]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[10]  E. Rouy,et al.  A viscosity solutions approach to shape-from-shading , 1992 .

[11]  William W. Symes,et al.  Upwind finite-difference calculation of traveltimes , 1991 .

[12]  D. Nichols.  Maximum energy traveltimes calculated in the seismic frequency band , 1996 .

[13]  J. Sethian.  Numerical algorithms for propagating interfaces: Hamilton-Jacobi equations and conservation laws , 1990 .

[14]  P. Podvin,et al.  Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools , 1991 .

[15]  J. Vidale.  Finite-difference calculation of travel times , 1988 .

[16]  J. Sethian.  Curvature and the evolution of fronts , 1985 .

[17]  J. Sethian,et al.  An O(N log N) algorithm for shape modeling. , 1996, Proceedings of the National Academy of Sciences of the United States of America.