Bayesian tomography with prior-knowledge-based parametrization and surrogate modeling

TL;DR

This paper introduces a Bayesian tomography method that leverages prior knowledge-based parametrization and surrogate modeling to improve inversion speed and accuracy in subsurface imaging, demonstrated with GPR travel-time data.

Contribution

The paper presents a novel Bayesian tomography framework combining prior-knowledge parametrization with surrogate models for efficient inversion.

Findings

Reduced-dimensional parametrization improves inversion constraints.

Surrogate models enable reliable MCMC with fewer data sets.

Uncertainty quantification is effectively achieved through component reintroduction.

Abstract

We present a Bayesian tomography framework operating with prior-knowledge-based parametrization that is accelerated by surrogate models. Standard high-fidelity forward solvers solve wave equations with natural spatial parametrizations based on fine discretization. Similar parametrizations, typically involving tens of thousand of variables, are usually employed to parameterize the subsurface in tomography applications. When the data do not allow to resolve details at such finely parameterized scales, it is often beneficial to instead rely on a prior-knowledge-based parametrization defined on a lower dimension domain (or manifold). Due to the increased identifiability in the reduced domain, the concomitant inversion is better constrained and generally faster. We illustrate the potential of a prior-knowledge-based approach by considering ground penetrating radar (GPR) travel-time tomography in a crosshole configuration. An effective parametrization of the input (i.e., the permittivity distributions) and compression in the output (i.e., the travel-time gathers) spaces are achieved via data-driven principal component decomposition based on random realizations of the prior Gaussian-process model with a truncation determined by the performances of the standard solver on the full and reduced model domains. To accelerate the inversion process, we employ a high-fidelity polynomial chaos expansion (PCE) surrogate model. We show that a few hundreds design data sets is sufficient to provide reliable Markov chain Monte Carlo inversion. Appropriate uncertainty quantification is achieved by reintroducing the truncated higher-order principle components in the original model space after inversion on the manifold and by adapting a likelihood function that accounts for the fact that the truncated higher-order components are not completely located in the null-space.