Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion

Astronomy and Astrophysics – Astronomy

Scientific paper

Rate now

  [ 0.00 ] – not rated yet Voters 0   Comments 0

Details

118

Diffraction, Finite-Difference Methods, Inversion, Numerical Techniques, Seismic Velocities, Wave Equation

Scientific paper

By specifying a discrete matrix formulation for the frequency-space modelling problem for linear partial differential equations (`FDM' methods), it is possible to derive a matrix formalism for standard iterative non-linear inverse methods, such as the gradient (steepest descent) method, the Gauss-Newton method and the full Newton method. We obtain expressions for each of these methods directly from the discrete FDM method, and we refer to this approach as frequency-domain inversion (FDI). The FDI methods are based on simple notions of matrix algebra, but are nevertheless very general. The FDI methods only require that the original partial differential equations can be expressed as a discrete boundary-value problem (that is as a matrix problem). Simple algebraic manipulation of the FDI expressions allows us to compute the gradient of the misfit function using only three forward modelling steps (one to compute the residuals, one to backpropagate the residuals, and a final computation to compute a step length). This result is exactly analogous to earlier backpropagation methods derived using methods of functional analysis for continuous problems. Following from the simplicity of this result, we give FDI expressions for the approximate Hessian matrix used in the Gauss-Newton method, and the full Hessian matrix used in the full Newton method. In a new development, we show that the additional term in the exact Hessian, ignored in the Gauss-Newton method, can be efficiently computed using a backpropagation approach similar to that used to compute the gradient vector. The additional term in the Hessian predicts the degradation of linearized inversions due to the presence of first-order multiples (such as free-surface multiples in seismic data). Another interpretation is that this term predicts changes in the gradient vector due to second-order non-linear effects. In a numerical test, the Gauss-Newton and full Newton methods prove effective in helping to solve the difficult non-linear problem of extracting a smooth background velocity model from surface seismic-reflection data.

No associations

LandOfFree

Say what you really think

Search LandOfFree.com for scientists and scientific papers. Rate them and share your experience with other people.

Rating

Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion will most certainly appreciate the feedback.

Rate now

     

Profile ID: LFWR-SCP-O-1369134

  Search
All data on this website is collected from public sources. Our data reflects the most accurate information available at the time of publication.