Efficient and accurate computation of seismic traveltimes and amplitudes

We describe two practicable approaches for an efficient computation of seismic traveltimes and amplitudes. The first approach is based on a combined finite-difference solution of the eikonal equation and the transport equation (the ‘FD approach’). These equations are formulated as hyperbolic conservation laws; the eikonal equation is solved numerically by a third-order ENO–Godunov scheme for the traveltimes whereas the transport equation is solved by a first-order upwind scheme for the amplitudes. The schemes are implemented in 2D using polar coordinates. The results are first-arrival traveltimes and the corresponding amplitudes. The second approach uses ray tracing (the ‘ray approach’) and employs a wavefront construction (WFC) method to calculate the traveltimes. Geometrical spreading factors are then computed from these traveltimes via the ray propagator without the need for dynamic ray tracing or numerical differentiation. With this procedure it is also possible to obtain multivalued traveltimes and the corresponding geometrical spreading factors. Both methods are compared using the Marmousi model. The results show that the FD eikonal traveltimes are highly accurate and perfectly match the WFC traveltimes. The resulting FD amplitudes are smooth and consistent with the geometrical spreading factors obtained from the ray approach. Hence, both approaches can be used for fast and reliable computation of seismic first-arrival traveltimes and amplitudes in complex models. In addition, the capabilities of the ray approach for computing traveltimes and spreading factors of later arrivals are demonstrated with the help of the Shell benchmark model.

[1]  P. Hubral,et al.  PARABOLIC AND HYPERBOLIC PARAXIAL TWO‐POINT TRAVELTIMES IN 3D MEDIA1 , 1993 .

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

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

[4]  Traveltime Based True-amplitude Migration , 2000 .

[5]  Determination of geometrical spreading from traveltimes , 2003 .

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

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

[8]  Moshe Reshef,et al.  Migration of common‐shot gathers , 1986 .

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

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

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

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

[13]  M. Kvasnička,et al.  3-D network ray tracing , 1994 .

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

[15]  Martin Tygel,et al.  3-D true‐amplitude finite‐offset migration , 1993 .

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

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

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

[19]  C. Vanelle,et al.  Traveltime Based True Amplitude Migration of PS Converted Wave , 2001 .

[20]  W. Symes,et al.  Finite‐difference quasi‐P traveltimes for anisotropic media , 2001 .

[21]  John E. Vidale,et al.  Rapid calculation of seismic amplitudes , 1990 .

[22]  Håvar Gjøystdal,et al.  Estimation of multivalued arrivals in 3D models using wavefront construction—Part I , 1996 .

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

[24]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

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

[26]  Determining the Ray Propagator From Traveltimes , 1998 .

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

[28]  Jean Brac,et al.  Can we image complex structures with first‐arrival traveltime? , 1993 .