1887
  1. Home
  2. A-Z Publications
  3. Geophysical Prospecting
  4. Volume 65, Issue S1
  5. Article
Volume 65, Issue S1
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Classical least‐squares techniques (Moore–Penrose pseudoinverse) are covariance based and are therefore unsuitable for the solution of very large‐scale linear systems in geophysical inversion due to the need of diagonalisation. In this paper, we present a methodology to perform the geophysical inversion of large‐scale linear systems via the discrete wavelet transform. The methodology consists of compressing the linear system matrix using the interesting properties of covariance‐free orthogonal transformations, to design an approximation of the Moore–Penrose pseudoinverse. We show the application of the discrete wavelet transform pseudoinverse to well‐conditioned and ill‐conditioned linear systems. We applied the methodology to a general‐purpose linear problem where the system matrix has been generated using geostatistical simulation techniques and also to a synthetic 2D gravimetric problem with two different geological set‐ups, in the noise‐free and noisy cases. In both cases, the discrete wavelet transform pseudoinverse can be applied to the original linear system and also to the linear systems of normal equations and minimum norm. The results are compared with those obtained via the Moore–Penrose and the discrete cosine transform pseudoinverses. The discrete wavelet transform and the discrete cosine transform pseudoinverses provide similar results and outperform the Moore–Penrose pseudoinverse, mainly in the presence of noise. In the case of well‐conditioned linear systems, this methodology is more efficient when applied to the least‐squares system and minimum norm system due to their higher condition number that allows for a more efficient compression of the system matrix. Also, in the case of ill‐conditioned systems with very high underdetermined character, the application of the discrete cosine transform to the minimum norm solution provides very good results. Both solutions might differ on their regularity, depending on the wavelet family that is adopted. These methods have a general character and can be applied to solve any linear inverse problem arising in technology, particularly in geophysics, and also to non‐linear inversion by linearisation of the forward operator.

© 2017 European Association of Geoscientists & Engineers © 2017 European Association of Geoscientists & Engineers
Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12595
2017-12-26
2024-04-25

  • From This Site

    /content/journals/10.1111/1365-2478.12595
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. AbediM., GholamiA., NorouziG.H. and FathianpourN.2013. Fast inversion of magnetic data using Lanczos bidiagonalization method. Journal of Applied Geophysics9, 126–137.
     [Google Scholar]
  2. AbediM., SiahkoohiH.R., GholamiA. and NorouziG.H.2015. 3D inversion of magnetic data through wavelet based regularization method. International Journal of Mining and Geo‐Engineering49(1), 1–18.
     [Google Scholar]
  3. AsterR.C., Borchers, B. and ThurberC.H.2005. Parameter Estimation and Inverse Problems. Elsevier Academic Press.
     [Google Scholar]
  4. ChamoliA.2009. Wavelet analysis of geophysical time series. E‐Journal Earth Science India 2 (IV), 258–275.
     [Google Scholar]
  5. ChiaoL.Y. and LiangW.T.2003. Multiresolution parameterization for geophysical inverse problems. Geophysics68, 199–209.
     [Google Scholar]
  6. CoifmanR.R., MeyerY. and WickerhauserV.1992. Wavelet analysis and signal processing. In: Wavelets and Their Applications (eds M.B.Ruskaiet al.), pp. 153–178. Boston, MA: Jones and Bartlett.
     [Google Scholar]
  7. CoifmanR.R. and DonohoD.L.1995. Translation‐Invariant De‐noising. Technical Report 475. Stanford, CA: Stanford University.
     [Google Scholar]
  8. DaubechiesI.1992. Ten lectures on wavelets, SIAM lectures. CBMS‐NSF Regional Conference Series in Applied Mathematics.
     [Google Scholar]
  9. De VoreR.A., JawerthB. and LucierB.1992. Image compression through wavelet transform coding. IEEE Transactions on Information Theory38, 719–746.
     [Google Scholar]
  10. DoiT., HayanoS. and SaitoY.1997. Wavelet solution of inverse source problems. IEEE Transactions on Magnetics33, 1935–1938.
     [Google Scholar]
  11. Fernández‐MuñizZ., Fernández‐MartínezJ.L., SrinivasanS. and MukerjiT.2015. Comparative analysis of the solution of linear continuous inverse problems using different basis expansions. Journal of Applied Geophysics113, 92–102.
     [Google Scholar]
  12. Fernández‐MartínezJ.L., Fernández‐MuñizZ., PalleroJ.L.G. and Bonvalot, S.2017. Linear geophysical inversion via the DCT pseudoinverse. Application to potential fields. Geophysical Prospecting.
     [Google Scholar]
  13. GholamiA. and SiahkoohiH.R.2009. A two‐step wavelet‐based regularization for linear inversion of geophysical data. Geophysical Prospecting57, 847–862.
     [Google Scholar]
  14. HansenP.C., PereyraV. and SchererG.2013. Least Squares Data Fitting with Applications. Johns Hopkins University Press.
     [Google Scholar]
  15. KaneJ.A.2003. Wavelets domain inversion and joint deconvolution/interpolation of geophysical data. PhD thesis, Massachusetts Institute of Technology, USA.
  16. KlannE.2005. Regularization of linear ill‐posed problems in two steps: combination of data smoothing and reconstruction methods. PhD thesis, University of Bremen, Germany.
  17. KumarP. and Foufoula‐GeorgiouE.1997. Wavelet analysis for geophysical applications. Review of Geophysics35(4), 385–412.
     [Google Scholar]
  18. LangM., GuoH., OdegartJ.E., BurrusC.S. and WellR.O.1995. Nonlinear processing of a shift invariant DWT for noise reduction. Proceedings of SPIE Mathematical Imaging: Wavelet Applications for Dual Use, Orlando, FL.
  19. LiY. and OldenburgD.W.2003. Fast inversion of large‐scale magnetic data using wavelet transforms and a logarithmic barrier method. Geophysical Journal International152, 251–265.
     [Google Scholar]
  20. MallatS.1989. A theory for multiresolution signal decomposition: the wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence11(7), 674–693.
     [Google Scholar]
  21. MarteletG., SailhacP., MoreauF. and DiamentM.2001. Characterization of geological boundaries using 1‐D wavelet transform on gravity data: theory and application to the Himalayas. Geophysics66, 1116–1129.
     [Google Scholar]
  22. MoreauF., GibertD., HolschneiderM. and SaraccoG.1997. Wavelet analysis of potential fields. Inverse Problem13, 165–178.
     [Google Scholar]
  23. MoreauF., GibertD., HolschneiderM. and SaraccoG.1999. Identification of sources of potential fields with the continuous wavelet transform: basic theory. Journal of Geophysical Research104, 5003–5013.
     [Google Scholar]
  24. OuadfeulS., AliouaneL. and EladjS.2010. Multiscale analysis of geomagnetic data using the continuous wavelet transform: a case study from Hoggar (Algeria). SEG, Expanded Abstracts, 1222–1225.
  25. PolikarR.1999. The story of wavelets. IMACS/IEEE CSCC'99 Proceedings, 5481–5486.
  26. TamboliS.S. and UdupiV.R.2012. Compression methods using wavelet transform. International Journal of Computer Engineering & Technology (IJCET)3(1), 314–321.
     [Google Scholar]
  27. TangY.Y., YangL.H., LinJ. and MaH.2000. Wavelet Theory and Its Applications to Pattern Recognition. Series in Machine Perception Artificial Intelligence, Vol. 36.
  28. Van FleetP.J.2008. Discrete Wavelet Transformations: An Elementary Approach with Applications. John Wiley & Sons, Inc.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12595
Loading
Linear inversion via the discrete wavelet transform pseudoinverse
Geophysical Prospecting 65, 131 (2017); https://doi.org/10.1111/1365-2478.12595
/content/journals/10.1111/1365-2478.12595
/content/journals/10.1111/1365-2478.12595
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Discrete wavelet transform pseudoinverse; Geophysical inversion; Inverse problem

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