Publications des scientifiques de l'IRD

Legchenko Anatoli, Comte J. C., Ofterdinger U., Vouillamoz Jean-Michel, Lawson F. M. A., Walsh J. (2017). Joint use of singular value decomposition and Monte-Carlo simulation for estimating uncertainty in surface NMR inversion. Journal of Applied Geophysics, 144, p. 28-36. ISSN 0926-9851.

Titre du document
Joint use of singular value decomposition and Monte-Carlo simulation for estimating uncertainty in surface NMR inversion
Année de publication
2017
Type de document
Article référencé dans le Web of Science WOS:000408781800003
Auteurs
Legchenko Anatoli, Comte J. C., Ofterdinger U., Vouillamoz Jean-Michel, Lawson F. M. A., Walsh J.
Source
Journal of Applied Geophysics, 2017, 144, p. 28-36 ISSN 0926-9851
We propose a simple and robust approach for investigating uncertainty in the results of inversion in geophysics. We apply this approach to inversion of Surface Nuclear Magnetic Resonance (SNMR) data, which is also known as Magnetic Resonance Sounding (MRS). Solution of this inverse problem is known to be non-unique. We inverse MRS data using the well-known Tikhonov regularization method, which provides an optimal solution as a trade-off between the stability and accuracy. Then, we perturb this model by random values and compute the fitting error for the perturbed models. The magnitude of these perturbations is limited by the uncertainty estimated with the singular value decomposition (SVD) and taking into account experimental errors. We use 106 perturbed models and show that the large majority of these models, which have all the water content within the variations given by the SVD estimate, do not fit data with an acceptable accuracy. Thus, we may limit the solution space by only the equivalent inverse models that fit data with the accuracy close to that of the initial inverse model. For representing inversion results, we use three equivalent solutions instead of the only one: the "best" solution given by the regularization or other inversion technic and the extreme variations of this solution corresponding to the equivalent models with the minimum and the maximum volume of water. For demonstrating our approach, we use synthetic data sets and experimental data acquired in the framework of investigation of a hard rock aquifer in the Ireland (County Donegal).
Plan de classement
Sciences fondamentales / Techniques d'analyse et de recherche [020] ; Hydrologie [062] ; Géologie et formations superficielles [064]
Description Géographique
IRLANDE
Localisation
Fonds IRD [F B010070971]
Identifiant IRD
fdi:010070971
Contact