Lithological tomography with the correlated pseudo-marginal method

Friedli, L.; Linde, N.; Ginsbourger, D.; Doucet, A. (2022). Lithological tomography with the correlated pseudo-marginal method. Geophysical journal international, 228(2), pp. 839-856. Oxford University Press 10.1093/gji/ggab381

Text ggab381.pdf - Published Version Restricted to registered users only Available under License Publisher holds Copyright. Author holds Copyright Download (3MB) | Request a copy

We consider lithological tomography in which the posterior distribution of (hydro)geological parameters of interest is inferred from geophysical data by treating the intermediate geophysical properties as latent variables. In such a latent variable model, one needs to estimate the intractable likelihood of the (hydro)geological parameters given the geophysical data. The pseudo-marginal (PM) method is an adaptation of the Metropolis–Hastings algorithm in which an unbiased approximation of this likelihood is obtained by Monte Carlo averaging over samples from, in this setting, the noisy petrophysical relationship linking (hydro)geological and geophysical properties. To make the method practical in data-rich geophysical settings with low noise levels, we demonstrate that the Monte Carlo sampling must rely on importance sampling distributions that well approximate the posterior distribution of petrophysical scatter around the sampled (hydro)geological parameter field. To achieve a suitable acceptance rate, we rely both on (1) the correlated PM (CPM) method, which correlates the samples used in the proposed and current states of the Markov chain and (2) a model proposal scheme that preserves the prior distribution. As a synthetic test example, we infer porosity fields using crosshole ground-penetrating radar (GPR) first-arrival traveltimes. We use a (50 × 50)-dimensional pixel-based parametrization of the multi-Gaussian porosity field with known statistical parameters, resulting in a parameter space of high dimension. We demonstrate that the CPM method with our proposed importance sampling and prior-preserving proposal scheme outperforms current state-of-the-art methods in both linear and non-linear settings by greatly enhancing the posterior exploration.

Interest & Impact

Downloads

1 since deposited on 11 Apr 2022
1 in the past 12 months

Citations

5 Citations in Web of Science ®
6 Citations in Scopus

Search

in Google Scholar™

Services

Actions (login required)

Edit item

Actions (login required)

Edit item