Non-linear optimization for seismic traveltime tomography

This paper presents a non-linear algorithmic approach for seismic traveltime. It is based on large-scale optimization using non-linear least-squares and trust-region methods. These methods provide a natural way to stabilize algorithms based on Newton's iteration for non-linear minimization. They also correspond to an alternative (and often more efficient) view of the Levenberg-Marquardt method. Numerical experience on synthetic data and on real borehole-to-borehole problems are presented. In particular, results produced by the new algorithm are compared with those of Ivansson (1985) for the Kra˙kema˙la experiment

[1]  Stephen J. Wright,et al.  A new non-linear least squares algorithm for the seismic inversion problem , 1986, Geophysical Journal International.

[2]  Jonathan M. Lees,et al.  Tomographic inversion for three‐dimensional velocity structure at Mount St. Helens using earthquake data , 1989 .

[3]  Ray Buland,et al.  The mechanics of locating earthquakes , 1976, Bulletin of the Seismological Society of America.

[4]  Sven Ivansson,et al.  Tomographic velocity estimation in the presence of low velocity zones , 1984 .

[5]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[6]  Manfred Koch,et al.  A numerical study on the determination of the 3-D structure of the lithosphere by linear and non-linear inversion of teleseismic travel times , 1985 .

[7]  S. Treitel,et al.  Fast l p solution of large, sparse, linear systems: application to seismic travel time tomography , 1988 .

[8]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[9]  Guust Nolet,et al.  Solving or resolving inadequate and noisy tomographic systems , 1985 .

[10]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[11]  Keiiti Aki,et al.  Determination of three‐dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes: 1. A homogeneous initial model , 1976 .

[12]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[13]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[14]  R. P. Bording,et al.  Applications of seismic travel-time tomography , 1987 .

[15]  C. H. Chapman,et al.  Crosshole seismic tomography , 1989 .

[16]  Trond Steihaug.  Local and superlinear convergence for truncated iterated projections methods , 1983, Math. Program..

[17]  Guust Nolet,et al.  Seismic tomography : with applications in global seismology and exploration geophysics , 1987 .

[18]  Keiiti Aki,et al.  Determination of the three‐dimensional seismic structure of the lithosphere , 1977 .

[19]  R. Crosson,et al.  Crustal structure modeling of earthquake data: 1. Simultaneous least squares estimation of hypocenter and velocity parameters , 1976 .

[20]  T. Steihaug.  The Conjugate Gradient Method and Trust Regions in Large Scale Optimization , 1983 .

[21]  Harold Jeffreys,et al.  An Alternative to the Rejection of Observations , 1932 .

[22]  J. Virieux.  Fast and accurate ray tracing by Hamiltonian perturbation , 1991 .

[23]  T. J. Moser.  The shortest path method for seismic ray tracing in complicated media , 1992 .