Abstract
Bayesian modelling and computational inference by Markov chain Monte Carlo (MCMC) is a principled framework for large-scale uncertainty quantification, though is limited in practice by computational cost when implemented in the simplest form that requires simulating an accurate computer model at each iteration of the MCMC. The delayed acceptance Metropolis–Hastings MCMC leverages a reduced model for the forward map to lower the compute cost per iteration, though necessarily reduces statistical efficiency that can, without care, lead to no reduction in the computational cost of computing estimates to a desired accuracy. Randomizing the reduced model for the forward map can dramatically improve computational efficiency, by maintaining the low cost per iteration but also avoiding appreciable loss of statistical efficiency. Randomized maps are constructed by a posteriori adaptive tuning of a randomized and locally-corrected deterministic reduced model. Equivalently, the approximated posterior distribution may be viewed as induced by a modified likelihood function for use with the reduced map, with parameters tuned to optimize the quality of the approximation to the correct posterior distribution. Conditions for adaptive MCMC algorithms allow practical approximations and algorithms that have guaranteed ergodicity for the target distribution. Good statistical and computational efficiencies are demonstrated in examples of calibration of large-scale numerical models of geothermal reservoirs and electrical capacitance tomography.
Original language | English |
---|---|
Article number | 26 |
Journal | GEM - International Journal on Geomathematics |
Volume | 11 |
Issue number | 1 |
DOIs | |
Publication status | Published - 1 Dec 2020 |
Keywords
- Adaptive MCMC
- Capacitance tomography
- Delayed acceptance
- Geothermal reservoir
- Inverse problem
- Markov chain Monte Carlo (MCMC)
- Reduced model
ASJC Scopus subject areas
- General Earth and Planetary Sciences
- Modelling and Simulation