A Hierarchical Multilevel Markov Chain Monte Carlo Algorithm with Applications to Uncertainty Quantification in Subsurface Flow

Abstract

In this paper we address the problem of the prohibitively large computational cost of existing Markov chain Monte Carlo methods for large--scale applications with high dimensional parameter spaces, e.g. in uncertainty quantification in porous media flow. We propose a new multilevel Metropolis-Hastings algorithm, and give an abstract, problem dependent theorem on the cost of the new multilevel estimator based on a set of simple, verifiable assumptions. For a typical model problem in subsurface flow, we then provide a detailed analysis of these assumptions and show significant gains over the standard Metropolis-Hastings estimator. Numerical experiments confirm the analysis and demonstrate the effectiveness of the method with consistent reductions of more than an order of magnitude in the cost of the multilevel estimator over the standard Metropolis-Hastings algorithm for tolerances $\epsilon <10^{-2}$.

References in corpus (1)

• Explicit error bounds for Markov chain Monte Carlo

Cited by in corpus (7)

• Multilevel ensemble Kalman filtering
• Importance sampling correction versus standard averages of reversible MCMCs in terms of the asymptotic variance
• Multilevel Sequential Monte Carlo Samplers for Normalizing Constants
• Multilevel ensemble Kalman filtering for spatially extended models
• A multilevel Monte Carlo method for a class of McKean-Vlasov processes
• Optimal fidelity multi-level Monte Carlo for quantification of uncertainty in simulations of cloud cavitation collapse
• Scalable iterative methods for sampling from massive Gaussian random vectors