Discrete cosine transform for parameter space reduction in Bayesian electrical resistivity tomography

A. Vinciguerra,M. Aleardi,A. Hojat,M. H. Loke,E. Stucchi

DOI: https://doi.org/10.1111/1365-2478.13148

2021-10-01

Geophysical Prospecting

Abstract:Electrical Resistivity Tomography (ERT) is a non-linear and ill-posed geophysical inverse problem that is usually solved through gradient-descent methods. This strategy is computationally fast and easy to implement but impedes accurate uncertainty appraisals. We present a probabilistic approach to 2D ERT in which a Markov Chain Monte Carlo (MCMC) algorithm is used to numerically evaluate the posterior probability density function that fully quantifies the uncertainty affecting the recovered solution. The main drawback of MCMC approaches is related to the considerable number of sampled models needed to achieve accurate posterior assessments in high-dimensional parameter spaces. Therefore, to reduce the computational burden of the inversion process, we employ the Differential Evolution Markov Chain (DEMC), a hybrid method between non-linear optimisation and MCMC sampling, which exploits multiple and interactive chains to speed up the probabilistic sampling. Moreover, the Discrete Cosine Transform (DCT) reparameterization is employed to reduce the dimensionality of the parameter space removing the high-frequency components of the resistivity model which are not sensitive to data. In this framework, the unknown parameters become the series of coefficients associated with the retained DCT basis functions. First, synthetic data inversions are used to validate the proposed method and to demonstrate the benefits provided by the DCT compression. To this end, we compare the outcomes of the implemented inversion approach with those provided by a DEMC algorithm running in the full, un-reduced model space. Then, we apply the proposed approach to invert field data acquired along a river embankment. The results yielded by the implemented approach are also benchmarked against a standard local inversion algorithm. The proposed Bayesian ERT inversion provides posterior mean models in agreement with the predictions achieved by the gradient-based inversion, but the implemented approach also provides model uncertainties, which can be used for penetration depth and resolution limit identification.This article is protected by copyright. All rights reserved

geochemistry & geophysics