A multiscale strategy for Bayesian inference using transport maps

Abstract

In many inverse problems, model parameters cannot be precisely determined from observational data. Bayesian inference provides a mechanism for capturing the resulting parameter uncertainty, but typically at a high computational cost. This work introduces a multiscale decomposition that exploits conditional independence across scales, when present in certain classes of inverse problems, to decouple Bayesian inference into two stages: (1) a computationally tractable coarse-scale inference problem; and (2) a mapping of the low-dimensional coarse-scale posterior distribution into the original high-dimensional parameter space. This decomposition relies on a characterization of the non-Gaussian joint distribution of coarse- and fine-scale quantities via optimal transport maps. We demonstrate our approach on a sequence of inverse problems arising in subsurface flow, using the multiscale finite element method to discretize the steady state pressure equation. We compare the multiscale strategy with full-dimensional Markov chain Monte Carlo on a problem of moderate dimension (100 parameters) and then use it to infer a conductivity field described by over 10,000 parameters.

References in corpus (1)

• Likelihood-informed dimension reduction for nonlinear inverse problems

Cited by in corpus (6)

• An introduction to sampling via measure transport
• Principal component analysis and sparse polynomial chaos expansions for global sensitivity analysis and model calibration: application to urban drainage simulation
• Residual-based error correction for neural operator accelerated infinite-dimensional Bayesian inverse problems
• Some models are useful, but how do we know which ones? Towards a unified Bayesian model taxonomy
• Recurrent Localization Networks applied to the Lippmann-Schwinger Equation
• Bayesian nonparametric generative modeling of large multivariate non-Gaussian spatial fields