Reconstructing the structure of the soil using non-invasive techniques is a very relevant problem in many scientific fields, like geophysics and archaeology. This can be done, for instance, with the aid of frequency domain electromagnetic (FDEM) induction devices. Inverting FDEM data is a very challenging inverse problem, as the problem is extremely ill-posed, i.e. sensible to the presence of noise in the measured data, and non-linear. Regularization methods substitute the original ill-posed problem with a well-posed one whose solution is an accurate approximation of the desired one. In this paper we develop a regularization method to invert FDEM data. We propose to determine the electrical conductivity of the ground by solving a variational problem. The minimized functional is made up by the sum of two term: the data fitting term ensures that the recovered solution fits the measured data, while the regularization term enforces sparsity on the Laplacian of the solution. The trade-off between the two terms is determined by the regularization parameter. This is achieved by minimizing an ℓ2 − ℓq functional with 0 < q ⩽ 2. Since the functional we wish to minimize is non-convex, we show that the variational problem admits a solution. Moreover, we prove that, if the regularization parameter is tuned accordingly to the amount of noise present in the data, this model induces a regularization method. Some selected numerical examples on synthetic and real data show the good performances of our proposal.

A variational non-linear constrained model for the inversion of FDEM data

Díaz de Alba, Patricia
2022-01-01

Abstract

Reconstructing the structure of the soil using non-invasive techniques is a very relevant problem in many scientific fields, like geophysics and archaeology. This can be done, for instance, with the aid of frequency domain electromagnetic (FDEM) induction devices. Inverting FDEM data is a very challenging inverse problem, as the problem is extremely ill-posed, i.e. sensible to the presence of noise in the measured data, and non-linear. Regularization methods substitute the original ill-posed problem with a well-posed one whose solution is an accurate approximation of the desired one. In this paper we develop a regularization method to invert FDEM data. We propose to determine the electrical conductivity of the ground by solving a variational problem. The minimized functional is made up by the sum of two term: the data fitting term ensures that the recovered solution fits the measured data, while the regularization term enforces sparsity on the Laplacian of the solution. The trade-off between the two terms is determined by the regularization parameter. This is achieved by minimizing an ℓ2 − ℓq functional with 0 < q ⩽ 2. Since the functional we wish to minimize is non-convex, we show that the variational problem admits a solution. Moreover, we prove that, if the regularization parameter is tuned accordingly to the amount of noise present in the data, this model induces a regularization method. Some selected numerical examples on synthetic and real data show the good performances of our proposal.
2022
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11386/4775513
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 4
  • ???jsp.display-item.citation.isi??? 3
social impact