Multiple reflection and transmission phases in complex layered media using a multistage fast marching method

Traditional grid-based eikonal schemes for computing traveltimes are usually confined to obtaining first arrivals only. However, later arrivals can be numerous and of greater amplitude, making them a potentially valuable resource for practical applications such as seismic imaging. The aim of this paper is to introduce a grid-based method for tracking multivalued wavefronts composed of any number of reflection and refraction branches in layered media. A finite-difference eikonal solver known as the fast marching method (FMM) is used to propagate wavefronts from one interface to the next. By treating each layer that the wavefront enters as a separate computational domain, one obtains a refracted branch by reinitializing FMM in the adjacent layer and a reflected branch by reinitializing FMM in the incident layer. To improve accuracy, a local grid refinement scheme is used in the vicinity of the source where wavefront curvature is high. Several examples are presented which demonstrate the viability of the new method in highly complex layered media. Even in the presence of velocity variations as large as 8:1 and interfaces of high curvature, wavefronts composed of many reflection and transmission events are tracked rapidly and accurately. This is because the scheme retains the two desirable properties of a single-stage FMM: computational speed and stability. Local grid refinement about the source also can increase accuracy by an order of magnitude with little increase in computational cost.

[1]  S. Chapman.  Notes on isomagnetic charts: I—Their singular points and contour-lines , 1940 .

[2]  Klaus H. Jacob,et al.  Three‐dimensional seismic ray tracing in a laterally heterogeneous spherical Earth , 1970 .

[3]  Victor Pereyra,et al.  Solving two-point seismic-ray tracing problems in a heterogeneous medium. Part 1. A general adaptive finite difference method , 1980 .

[4]  Kaoru Yamaguchi,et al.  A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE , 1986 .

[5]  B. Barsky,et al.  An Introduction to Splines for Use in Computer Graphics and Geometric Modeling , 1987 .

[6]  V. Červený.  Ray tracing algorithms in three-dimensional laterally varying layered structures , 1987 .

[7]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[8]  Raul Madariaga,et al.  Non-Linear Reflection Tomography , 1988 .

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

[10]  William A. Prothero,et al.  A fast, two-point, three-dimensional raytracing algorithm using a simple step search method , 1988 .

[11]  Malcolm Sambridge,et al.  Non-linear arrival time inversion: constraining velocity anomalies by seeking smooth models in 3-D , 1990 .

[12]  Robert L. Nowack,et al.  Inversion for Crustal Structure Using Reflections From the PASSCAL Ouachita Experiment , 1990 .

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

[14]  T. Moser.  Shortest path calculation of seismic rays , 1991 .

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

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

[17]  Robert B. Smith,et al.  Seismic traveltime inversion for 2-D crustal velocity structure , 1992 .

[18]  Jonathan M. Lees,et al.  Shortest path ray tracing with sparse graphs , 1993 .

[19]  Christopher Juhlin,et al.  3-D interpretation of reflected arrival times by finite-difference techniques , 1994 .

[20]  John A. Hole,et al.  3-D finite-difference reflection travel times , 1995 .

[21]  Petr Bulant,et al.  Two-point ray tracing in 3-D , 1996 .

[22]  Leigh House,et al.  Minimum traveltime calculation in 3-D graph theory , 1996 .

[23]  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.

[24]  George A. McMechan,et al.  3-D two-point ray tracing for heterogeneous, weakly transversely isotropic media , 1996, 1902.01969.

[25]  Dimitri Bevc,et al.  Imaging complex structure with semirecursive Kirchhoff migration , 1997 .

[26]  R. Kubota,et al.  Crustal Structure in the Southern Kanto-Tokai Region Derived from Tomographic Method for Seismic Explosion Survey , 1997 .

[27]  William W. Symes,et al.  A slowness matching finite difference method for traveltimes beyond transmission caustics , 1998 .

[28]  J. Benamou.  Direct computation of multivalued phase space solutions for Hamilton-Jacobi equations , 1999 .

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

[30]  Richard I. Cook,et al.  3-D traveltime computation using second‐order ENO scheme , 1999 .

[31]  Ettrich,et al.  Comparison of six different methods for calculating traveltimes , 1999 .

[32]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[33]  Kazuki Koketsu,et al.  Finite difference traveltime calculation for head waves travelling along an irregular interface , 2000 .

[34]  N. Rawlinson,et al.  Inversion of seismic refraction and wide-angle reflection traveltimes for three-dimensional layered crustal structure , 2001 .

[35]  David L. Chopp,et al.  Some Improvements of the Fast Marching Method , 2001, SIAM J. Sci. Comput..

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

[37]  C. Baag,et al.  Rapid and Accurate Two-Point Ray Tracing Based on a Quadratic Equation of Takeoff Angle in Layered Media with Constant or Linearly Varying Velocity Functions , 2002 .

[38]  J. Sethian,et al.  Fast-phase space computation of multiple arrivals , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Stanley Osher,et al.  Geometric Optics in a Phase-Space-Based Level Set and Eulerian Framework , 2002 .

[40]  Björn Engquist,et al.  High-Frequency Wave Propagation by the Segment Projection Method , 2002 .

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

[42]  Jianliang Qian,et al.  An adaptive finite-difference method for traveltimes and amplitudes , 2002 .

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

[44]  M. Sambridge,et al.  Wave front evolution in strongly heterogeneous layered media using the fast marching method , 2004 .