TY - JOUR
T1 - Inversion of Multiconfiguration Complex EMI Data with Minimum Gradient Support Regularization
T2 - A Case Study
AU - Deidda, Gian Piero
AU - Díaz de Alba, Patricia
AU - Rodriguez, Giuseppe
AU - Vignoli, Giulio
N1 - Funding Information:
The authors wish to thank the reviewers for their comments, which led to improvements in the presentation. This research was supported in part by the Fondazione di Sardegna 2017 research project “Algorithms for Approximation with Applications [Acube]”, the INdAM-GNCS research project “Metodi numerici per problemi mal posti”, the INdAM-GNCS research project “Discretizzazione di misure, approssimazione di operatori integrali ed applicazioni”, the Regione Autonoma della Sardegna research project “Algorithms and Models for Imaging Science [AMIS]” (RASSR57257, intervento finanziato con risorse FSC 2014-2020 - Patto per lo Sviluppo della Regione Sardegna), the Visiting Scientist Programme 2015/2016 (University of Cagliari), and an RAS/FBS Grant (Grant No. F71/17000190002).
Funding Information:
The authors wish to thank the reviewers for their comments, which led to improvements in the presentation. This research was supported in part by the Fondazione di Sardegna 2017 research project ?Algorithms for Approximation with Applications [Acube]?, the INdAM-GNCS research project ?Metodi numerici per problemi mal posti?, the INdAM-GNCS research project ?Discretizzazione di misure, approssimazione di operatori integrali ed applicazioni?, the Regione Autonoma della Sardegna research project ?Algorithms and Models for Imaging Science [AMIS]? (RASSR57257, intervento finanziato con risorse FSC 2014-2020 - Patto per lo Sviluppo della Regione Sardegna), the Visiting Scientist Programme 2015/2016 (University of Cagliari), and an RAS/FBS Grant (Grant No. F71/17000190002).
Publisher Copyright:
© 2020, International Association for Mathematical Geosciences.
PY - 2020/10/1
Y1 - 2020/10/1
N2 - Frequency-domain electromagnetic instruments allow the collection of data in different configurations, that is, varying the intercoil spacing, the frequency, and the height above the ground. Their handy size makes these tools very practical for near-surface characterization in many fields of applications, for example, precision agriculture, pollution assessments, and shallow geological investigations. To this end, the inversion of either the real (in-phase) or the imaginary (quadrature) component of the signal has already been studied. Furthermore, in many situations, a regularization scheme retrieving smooth solutions is blindly applied, without taking into account the prior available knowledge. The present work discusses an algorithm for the inversion of the complex signal in its entirety, as well as a regularization method that promotes the sparsity of the reconstructed electrical conductivity distribution. This regularization strategy incorporates a minimum gradient support stabilizer into a truncated generalized singular value decomposition scheme. The results of the implementation of this sparsity-enhancing regularization at each step of a damped Gauss–Newton inversion algorithm (based on a nonlinear forward model) are compared with the solutions obtained via a standard smooth stabilizer. An approach for estimating the depth of investigation, that is, the maximum depth that can be investigated by a chosen instrument configuration in a particular experimental setting, is also discussed. The effectiveness and limitations of the whole inversion algorithm are demonstrated on synthetic and real data sets.
AB - Frequency-domain electromagnetic instruments allow the collection of data in different configurations, that is, varying the intercoil spacing, the frequency, and the height above the ground. Their handy size makes these tools very practical for near-surface characterization in many fields of applications, for example, precision agriculture, pollution assessments, and shallow geological investigations. To this end, the inversion of either the real (in-phase) or the imaginary (quadrature) component of the signal has already been studied. Furthermore, in many situations, a regularization scheme retrieving smooth solutions is blindly applied, without taking into account the prior available knowledge. The present work discusses an algorithm for the inversion of the complex signal in its entirety, as well as a regularization method that promotes the sparsity of the reconstructed electrical conductivity distribution. This regularization strategy incorporates a minimum gradient support stabilizer into a truncated generalized singular value decomposition scheme. The results of the implementation of this sparsity-enhancing regularization at each step of a damped Gauss–Newton inversion algorithm (based on a nonlinear forward model) are compared with the solutions obtained via a standard smooth stabilizer. An approach for estimating the depth of investigation, that is, the maximum depth that can be investigated by a chosen instrument configuration in a particular experimental setting, is also discussed. The effectiveness and limitations of the whole inversion algorithm are demonstrated on synthetic and real data sets.
KW - Depth of investigation (DOI)
KW - Electromagnetic induction
KW - Minimum gradient support stabilizer
KW - Nonlinear inverse problems
KW - Truncated generalized singular value decomposition (TGSVD)
UR - http://www.scopus.com/inward/record.url?scp=85079400334&partnerID=8YFLogxK
U2 - 10.1007/s11004-020-09855-4
DO - 10.1007/s11004-020-09855-4
M3 - Article
AN - SCOPUS:85079400334
SN - 1874-8961
VL - 52
SP - 945
EP - 970
JO - Mathematical Geosciences
JF - Mathematical Geosciences
IS - 7
ER -