Two efficient algorithms for iterative linearized inversion of seismic waveform data

We present two equivalent algorithms for iterative linearized waveform inversion for 3-D Earth structure with respect to an arbitrary 3-D starting model; one is a matrix formulation, and the second is a wavefield formulation. Both algorithms require the computation of accurate synthetic seismograms, but neither requires that any particular method be used to compute the synthetics. The matrix formulation is equivalent to our previously published algorithm (Hara, Tsuboi & Geller 1991), but requires less than 10 per cent of the CPU time of the previous algorithm. The wavefield algorithm is equivalent to that of Tarantola (1986) and Mora (1987), but appears to be substantially more efficient.

[1]  John H. Woodhouse,et al.  Mapping the upper mantle: Three‐dimensional modeling of earth structure by inversion of seismic waveforms , 1984 .

[2]  R. Pratt Inverse theory applied to multisource cross-hole tomography, Part2 : Elastic wave-equation method , 1990 .

[3]  R. Geller,et al.  Partial derivatives of synthetic seismograms for a laterally heterogeneous Earth model , 1987 .

[4]  G. Nolet Partitioned waveform inversion and two‐dimensional structure under the network of autonomously recording seismographs , 1990 .

[5]  A. Tarantola Inversion of seismic reflection data in the acoustic approximation , 1984 .

[6]  Peter Mora,et al.  Inversion = Migration + Tomography , 1988, Shell Conference.

[7]  A. Tarantola A strategy for nonlinear elastic inversion of seismic reflection data , 1986 .

[8]  R. Pratt,et al.  INVERSE THEORY APPLIED TO MULTI‐SOURCE CROSS‐HOLE TOMOGRAPHY.: PART 1: ACOUSTIC WAVE‐EQUATION METHOD1 , 1990 .

[9]  R. Geller,et al.  Inversion for laterally heterogeneous upper mantle S-wave velocity structure using iterative waveform inversion , 1993 .

[10]  Z. Alterman,et al.  Rotational Multiplets in the Spectrum of the Earth , 1961 .

[11]  Gene H. Golub,et al.  Matrix computations , 1983 .

[12]  R. Geller,et al.  Variational free oscillation computations for three laterally heterogeneous Earth models , 1987 .

[13]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[14]  Robert J. Geller,et al.  Computation of synthetic seismograms and their partial derivatives for heterogeneous media with arbitrary natural boundary conditions using the Direct Solution Method , 1994 .