Finite difference traveltime calculation for anisotropic media

SUMMARY Seismic traveltimes for anisotropic elastic media are determined by numerical solution of the eikonal equation using a 2-D finite difference approach. Most previous examples of this technique have been restricted to the geometrical-optics form of the eikonal equation associated with isotropic media. Here, a general anisotropic form is used, with particular solutions for media with transversely isotropic (TI) symmetry. The algorithm tracks arrivals corresponding to a particular wave type (i.e. 4P, 4s) using an expanding-wavefront scheme on a hexagonal mesh. In a homogeneous example, the method is demonstrated to be accurate for 4P waves, but prone to large errors for 4s waves in the vicinity of wavefront cuspidal edges and along shear wave singularities. Two other examples that incorporate both anisotropy and inhomogeneity demonstrate the potential usefulness of this algorithm for calculating 4P traveltimes in more complex models.

[1]  Vlastislav Červený,et al.  Ray method in seismology , 1977 .

[2]  F. K. Levin Seismic velocities in transversely isotropic media , 1979 .

[3]  A. Balch,et al.  A dynamic programming approach to first arrival traveltime computation in media with arbitrarily distributed velocities , 1992 .

[4]  V. Červený Seismic Rays and Ray Intensities in Inhomogeneous Anisotropic Media , 1972 .

[5]  J. Hudson Wave speeds and attenuation of elastic waves in material containing cracks , 1981 .

[6]  Dirk Gajewski,et al.  Computation of high-frequency seismic wavefields in 3-D laterally inhomogeneous anisotropic media , 1987 .

[7]  Samuel Fox,et al.  41st ANNUAL MEETING , 1964 .

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

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

[10]  Daniel H. Rothman,et al.  Cellular‐automaton fluids: A model for flow in porous media , 1988 .

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

[12]  Hideki Sailo,et al.  Traveltimes And Raypaths of First Arrival Seismic Waves: Computation Method Based On Huygens’ Principle , 1989 .

[13]  Herbert F. Wang,et al.  Ultrasonic velocities in Cretaceous shales from the Williston basin , 1981 .

[14]  C. Thomson,et al.  A comment on the form of the geometrical spreading equations, with some numerical examples of seismic ray tracing in inhomogeneous, anisotropic media , 1989 .

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

[16]  S. Crampin,et al.  Seismic anisotropy - the state of the art , 1984 .

[17]  Joe Dellinger Anisotropic Finite-Difference Traveltimes , 1991 .

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

[19]  M. Musgrave,et al.  Crystal Acoustics: Introduction to the Study of Elastic Waves and Vibrations in Crystals , 1970 .

[20]  F. Hron,et al.  Reflection and transmission coefficients for transversely isotropic media , 1977, Bulletin of the Seismological Society of America.

[21]  W. Beydoun,et al.  Elastic Ray‐Born L2‐Migration/Inversion , 1989 .