Inversion of electromagnetic induction data using a novel wavelet-based and scale-dependent regularization term
- 6 May 2021
- journal article
- research article
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 226 (3) , 1715-1729
Abstract
SUMMARY: The inversion of electromagnetic induction data to a conductivity profile is an ill-posed problem. Regularization improves the stability of the inversion and a smoothing constraint is typically used. However, the conductivity profiles are not always expected to be smooth. Here, we develop a new inversion scheme in which we transform the model to the wavelet space and impose a sparsity constraint. This sparsity constrained inversion scheme will minimize an objective function with a least-squares data misfit and a sparsity measure of the model in the wavelet domain. A model transform to the wavelet domain allows to investigate the temporal resolution (periodicities at different frequencies) and spatial resolution (location of the peaks) characteristics of the model, and penalizing small-scale coefficients effectively reduces the complexity of the model. The novel scale-dependent regularization term can be used to favour either blocky or smooth structures, as well as high-amplitude models in globally smooth structures in the inversion. Depending on the expected conductivity profile, a suitable wavelet basis function can be chosen. The scheme supports multiple types of regularization with the same algorithm and is thus flexible. Finally, we apply this new scheme on a frequency domain electromagnetic sounding data set, but the scheme could equally apply to any other 1-D geophysical method.Keywords
All Related Versions
Funding Information
- FWO
- Fonds De La Recherche Scientifique - FNRS
This publication has 49 references indexed in Scilit:
- Electrical Conductivity Depth Modelling with a Multireceiver EMI Sensor for Prospecting Archaeological FeaturesArchaeological Prospection, 2011
- Realistic Subsurface Anomaly Discrimination Using Electromagnetic Induction and an SVM ClassifierEURASIP Journal on Advances in Signal Processing, 2010
- Evaluating the multiple coil configurations of the EM38DD and DUALEM‐21S sensors to detect archaeological anomaliesArchaeological Prospection, 2009
- For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solutionCommunications on Pure and Applied Mathematics, 2006
- Decoding by Linear ProgrammingIEEE Transactions on Information Theory, 2005
- An iterative thresholding algorithm for linear inverse problems with a sparsity constraintCommunications on Pure and Applied Mathematics, 2004
- Environmental geophysics at Scott Base: ground penetrating radar and electromagnetic induction as tools for mapping contaminated ground at Antarctic research basesCold Regions Science and Technology, 2003
- Line search algorithms with guaranteed sufficient decreaseACM Transactions on Mathematical Software, 1994
- A theory for multiresolution signal decomposition: the wavelet representationIEEE Transactions on Pattern Analysis and Machine Intelligence, 1989
- THE MAGNETIC DIPOLE OVER THE HORIZONTALLY STRATIFIED EARTHCanadian Journal of Physics, 1951