4.4 Article

Ensemble-Based Electrical Resistivity Tomography with Data and Model Space Compression

Journal

PURE AND APPLIED GEOPHYSICS
Volume 178, Issue 5, Pages 1781-1803

Publisher

SPRINGER BASEL AG
DOI: 10.1007/s00024-021-02730-1

Keywords

Electrical resistivity tomography; Bayesian inversion

Funding

  1. Universita di Pisa within the CRUI-CARE Agreement
  2. Ministero dell'Ambiente e della Tutela del Territorio e del Mare, Italy, project DILEMMA-Imaging, Modeling, Monitoring and Design of Earthen Levees

Ask authors/readers for more resources

The proposed alternative ERT method utilizes ensemble smoother with multiple data assimilation to provide multiple realizations of resistivity values in the subsurface and reduce ill-conditioning of the inversion process. It allows for accurate uncertainty appraisal and can be applied to nonlinear inverse problems.
Inversion of electrical resistivity tomography (ERT) data is an ill-posed problem that is usually solved through deterministic gradient-based methods. These methods guarantee a fast convergence but hinder accurate assessments of model uncertainties. On the contrary, Markov Chain Monte Carlo (MCMC) algorithms can be employed for accurate uncertainty appraisals, but they remain a formidable computational task due to the many forward model evaluations needed to converge. We present an alternative approach to ERT that not only provides a best-fitting resistivity model but also gives an estimate of the uncertainties affecting the inverse solution. More specifically, the implemented method aims to provide multiple realizations of the resistivity values in the subsurface by iteratively updating an initial ensemble of models based on the difference between the predicted and measured apparent resistivity pseudosections. The initial ensemble is generated using a geostatistical method under the assumption of log-Gaussian distributed resistivity values and a Gaussian variogram model. A finite-element code constitutes the forward operator that maps the resistivity values onto the associated apparent resistivity pseudosection. The optimization procedure is driven by the ensemble smoother with multiple data assimilation, an iterative ensemble-based algorithm that performs a Bayesian updating step at each iteration. The main advantages of the proposed approach are that it can be applied to nonlinear inverse problems, while also providing an ensemble of models from which the uncertainty on the recovered solution can be inferred. The ill-conditioning of the inversion procedure is decreased through a discrete cosine transform reparameterization of both data and model spaces. The implemented method is first validated on synthetic data and then applied to field data. We also compare the proposed method with a deterministic least-square inversion, and with an MCMC algorithm. We show that the ensemble-based inversion estimates resistivity models and associated uncertainties comparable to those yielded by a much more computationally intensive MCMC sampling.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.4
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available