Accuracy and efficiency of stencils for the eikonal equation in earth modelling

Motivated by the needs for creating fast and accurate models of complex geological scenarios, accuracy and efficiency of three stencils for the isotropic eikonal equation on rectangular grids are evaluated using a fast marching implementation. The stencils are derived by direct modelling of the wave front, resulting in new and valuable insight in terms of improved upwind and causality conditions. After introducing a method for generalising first-order upwind stencils to higher order, a new second-order diagonal stencil is presented. Similarly to the multistencil fast marching approach, the diagonal stencil makes use of nodes in the diagonal directions, whereas the traditional Godunov stencil uses solely edge-connected neighbours. The diagonal stencil uses nodes close to each other, reaching upwind, to get a more accurate estimate of the angle of incidence of the arriving wave front. Although the stencils are evaluated in a fast marching setting, they can be adapted to other efficient eikonal solvers. All first- and second-order stencils are evaluated in a range of tests. The first test case models a folded structure from the Zagros fold belt in Iran. The other test cases are constructed to investigate specific properties of the examined stencils. The numerical investigation considers convergence rates and CPU times for non-constant and constant speed first-arrival computations. In conclusion, the diagonal stencil is the most efficient and accurate of the three alternatives.

[1]  S. A. Petersen,et al.  Earth Recursion, an Important Component in Shared Earth Model Builders , 2008 .

[2]  Ron Kimmel,et al.  Optimal Algorithm for Shape from Shading and Path Planning , 2001, Journal of Mathematical Imaging and Vision.

[3]  Jun-Wei Huang,et al.  Joint transmission and reflection traveltime tomography using the fast sweeping method and the adjoint-state technique , 2012 .

[4]  Stanley Bak,et al.  Some Improvements for the Fast Sweeping Method , 2010, SIAM J. Sci. Comput..

[5]  S. A. Petersen,et al.  Optimization Strategy for Shared Earth Modeling , 2004 .

[6]  J. Mallet.  Space–Time Mathematical Framework for Sedimentary Geology , 2004 .

[7]  Knut-Andreas Lie,et al.  Fast computation of arrival times in heterogeneous media , 2006 .

[8]  Qingfen Lin,et al.  Enhancement, Extraction, and Visualization of 3D Volume Data , 2001 .

[9]  Stefan Turek,et al.  THE EIKONAL EQUATION: NUMERICAL EFFICIENCY VS. ALGORITHMIC COMPLEXITY ON QUADRILATERAL GRIDS , 2005 .

[10]  T. Gillberg.  A Semi-Ordered Fast Iterative Method (SOFI) for monotone front propagation in simulations of geological folding , 2011 .

[11]  W. Jeong.  UUCS-07-010 A Fast Iterative Method for a Class of Hamilton-Jacobi Equations on Parallel Systems , 2007 .

[12]  J. Sethian,et al.  3-D traveltime computation using the fast marching method , 1999 .

[13]  Jean-François Rainaud,et al.  Knowledge-driven applications for geological modeling , 2005 .

[14]  Aly A. Farag,et al.  MultiStencils Fast Marching Methods: A Highly Accurate Solution to the Eikonal Equation on Cartesian Domains , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[15]  S. A. Petersen,et al.  Earth Modelling Using Distance-Fields Derived by Fast Marching , 2007 .

[16]  Are Magnus Bruaset,et al.  Turning Rocks into Knowledge , 2010, Simula Research Laboratory.

[17]  Alexander M. Bronstein,et al.  Parallel algorithms for approximation of distance maps on parametric surfaces , 2008, TOGS.

[18]  Qingfen Lin,et al.  A Modified Fast Marching Method , 2003, SCIA.

[19]  J. Tsitsiklis,et al.  Efficient algorithms for globally optimal trajectories , 1994, Proceedings of 1994 33rd IEEE Conference on Decision and Control.

[20]  Anthony J. Yezzi,et al.  Fully Isotropic Fast Marching Methods on Cartesian Grids , 2010, ECCV.

[21]  J. Sethian,et al.  Ordered upwind methods for static Hamilton–Jacobi equations , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[23]  Hongkai Zhao,et al.  A fast sweeping method for Eikonal equations , 2004, Math. Comput..

[24]  Chi-Wang Shu,et al.  Uniformly Accurate Discontinuous Galerkin Fast Sweeping Methods for Eikonal Equations , 2011, SIAM J. Sci. Comput..

[25]  Mohamed-Jalal Fadili,et al.  Dual Front Evolution Model and Its Application in Medical Imaging , 2004, MICCAI.

[26]  Ø. Hjelle,et al.  A Hamilton–Jacobi Framework for Modeling Folds in Structural Geology , 2011 .

[27]  Ross T. Whitaker,et al.  A Fast Iterative Method for Eikonal Equations , 2008, SIAM J. Sci. Comput..

[28]  B. Alaei,et al.  Geological modelling and finite difference forward realization of a regional section from the Zagros fold-and-thrust belt , 2007, Petroleum Geoscience.

[29]  Tariq Alkhalifah,et al.  Implementing the fast marching eikonal solver: spherical versus Cartesian coordinates , 2001 .

[30]  James A. Sethian,et al.  3-D imaging using higher order fast marching traveltimes , 2002 .

[31]  Christopher M. Kuster,et al.  Computational Study of Fast Methods for the Eikonal Equation , 2005, SIAM J. Sci. Comput..

[32]  Hongkai Zhao,et al.  High Order Fast Sweeping Methods for Static Hamilton–Jacobi Equations , 2006, J. Sci. Comput..

[33]  Benjamin J. Kadlec,et al.  Leveraging graphics processing units (GPUs) for real-time seismic interpretation , 2010 .

[34]  Song-Gang Xu,et al.  A Fast Sweeping Method for Computing Geodesics on Triangular Manifolds , 2010, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[35]  Sergey Fomel,et al.  A second-order fast marching eikonal solver , 2000 .

[36]  S. Petersen.  Compound Modelling - a Geological Approach to the Construction of Shared Earth Models , 1999 .

[37]  Gilles Lambaré,et al.  Two-dimensional multivalued traveltime and amplitude maps by uniform sampling of a ray field , 1996 .

[38]  Seongjai Kim,et al.  An O(N) Level Set Method for Eikonal Equations , 2000, SIAM J. Sci. Comput..

[39]  Bobby Bodenheimer,et al.  Synthesis and evaluation of linear motion transitions , 2008, TOGS.

[40]  M. Sambridge,et al.  Seismic ray tracing and wavefront tracking in laterally heterogeneous media , 2008 .

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

[42]  Gerard T. Schuster,et al.  Finite‐difference solution of the eikonal equation along expanding wavefronts , 1992 .

[43]  Seongjai Kim.  3-D eikonal solvers: First-arrival traveltimes , 2002 .

[44]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[45]  M. Sambridge,et al.  Multiple reflection and transmission phases in complex layered media using a multistage fast marching method , 2004 .