1887
Volume 56 Number 1
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Linearized inversion methods such as Gauss‐Newton and multiple re‐weighted least‐squares are iterative processes in which an update in the current model is computed as a function of data misfit and the gradient of data with respect to model parameters. The main advantage of those methods is their ability to refine the model parameters although they have a high computational cost for seismic inversion. In the Gauss‐Newton method a system of equations, corresponding to the sensitivity matrix, is solved in the least‐squares sense at each iteration, while in the multiple re‐weighted least‐squares method many systems are solved using the same sensitivity matrix. The sensitivity matrix arising from these methods is usually not sparse, thus limiting the use of standard preconditioners in the solution of the linearized systems. For reduction of the computational cost of the linearized inversion methods, we propose the use of preconditioners based on a partial orthogonalization of the columns of the sensitivity matrix. The new approach collapses a band of co‐diagonals of the normal equations matrix into the main diagonal, being equivalent to computing the least‐squares solution starting from a partial solution of the linear system. The preconditioning is driven by a bandwidth which can be interpreted as the distance for which the correlation between model parameters is relevant. To illustrate the benefit of the proposed approach to the reduction of the computational cost of the inversion we apply the multiple re‐weighted least‐squares method to the 2D acoustic seismic waveform inversion problem. We verify the reduction in the number of iterations in the conjugate'gradient algorithm as the bandwidth of the preconditioners increases. This effect reduces the total computational cost of inversion as well.

Loading

Article metrics loading...

/content/journals/10.1111/j.1365-2478.2007.00662.x
2007-11-15
2024-03-28
Loading full text...

Full text loading...

References

  1. AlfordR., KellyK. and BooreD.1974. Accuracy of finite‐difference modeling of the acoustic wave equation. Geophysics39, 834–842.
    [Google Scholar]
  2. BenziM. and TumaM.2003. A robust preconditioner with low memory requirements for large sparse least squares problems. SIAM Journal on Scientific and Statistical Computing25, 499–512.
    [Google Scholar]
  3. GallivanK.A., PlemmonsR.J. and SamehA.H.1990. Parallel algorithms for dense linear algebra computations. SIAM Review32, 54–135.
    [Google Scholar]
  4. HeathM.T.1984. Numerical methods for large sparse linear least squares problems. SIAM Journal on Scientific and Statistical Computing5, 497–513.
    [Google Scholar]
  5. KennettB.L.N. and SambridgeM.S.1992. Earthquake location ‐ genetic algorithms for teleseisms. Physics of the Earth and Planetary Interiors75, 103–110.
    [Google Scholar]
  6. KosloffD. and BaysalE.1982. Forward modelling by a Fourier method. Geophysics47, 1402–1412.
    [Google Scholar]
  7. MenkeW.1989. Geophysical Data Analysis: Discrete Inverse Theory. International Geophysics Series, V . 45. Academic Press. ISBN 0124909313.
    [Google Scholar]
  8. PorsaniM.J., NiwasS. and FerreiraN.R.2001. Robust inversion of vertical electrical sounding data using a multiple re‐weighted least‐squares method. Geophysical Prospecting49, 255–264.
    [Google Scholar]
  9. PorsaniM.J. and OliveiraS.P.2005. Modified incomplete Cholesky preconditioners for 2D resistivity modelling. 67th EAGE meeting, Madrid , Spain , Expanded Abstract.
  10. PrattR.G., ShinC. and HicksG.J.1998. Gauss‐Newton and full Newton methods in frequency domain seismic waveform inversion. Geophysical Journal International133, 341–362.
    [Google Scholar]
  11. RavautC., OpertoS., ImprotaL., VirieuxJ., HerreroA. and Dell'AversanaP.2004. Multiscale imaging of complex structures from multifold wide‐aperture seismic data by frequency‐domain full‐waveform tomography: Application to a thrust belt. Geophysical Journal International159, 1032–1056.
    [Google Scholar]
  12. RothmanD.H.1985. Nonlinear inversion, statistical mechanics and residual statics estimation. Geophysics50, 2784–2796.
    [Google Scholar]
  13. SenM.K. and StoffaP.L.1995. Global Optimization Methods in Geophysical Inversion . Elsevier Science. ISBN 0444817670.
    [Google Scholar]
  14. ScalesJ.A. and GersztenkornA.1988. Robust methods in inverse theory. Inverse Problems. 4, 1071–1091.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/j.1365-2478.2007.00662.x
Loading
/content/journals/10.1111/j.1365-2478.2007.00662.x
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error