Population MCMC methods for history matching and uncertainty quantification

Lina Mahgoub Yahya Mohamed, Ben Calderhead, Maurizio Filippone, Mike Christie, Mark Girolami

    Research output: Contribution to journalArticlepeer-review

    14 Downloads (Pure)


    This paper presents the application of a population Markov Chain Monte Carlo (MCMC) technique to generate history-matched models. The technique has been developed and successfully adopted in challenging domains such as computational biology but has not yet seen application in reservoir modelling. In population MCMC, multiple Markov chains are run on a set of response surfaces that form a bridge from the prior to posterior. These response surfaces are formed from the product of the prior with the likelihood raised to a varying power less than one. The chains exchange positions, with the probability of a swap being governed by a standard Metropolis accept/reject step, which allows for large steps to be taken with high probability. We show results of Population MCMC on the IC Fault Model-a simple three-parameter model that is known to have a highly irregular misfit surface and hence be difficult to match. Our results show that population MCMC is able to generate samples from the complex, multi-modal posterior probability distribution of the IC Fault model very effectively. By comparison, previous results from stochastic sampling algorithms often focus on only part of the region of high posterior probability depending on algorithm settings and starting points.

    Original languageEnglish
    Pages (from-to)423-436
    Number of pages14
    JournalComputational Geosciences
    Issue number2
    Publication statusPublished - Mar 2012


    • History matching
    • MCMC methods
    • Population MCMC
    • Uncertainty quantification


    Dive into the research topics of 'Population MCMC methods for history matching and uncertainty quantification'. Together they form a unique fingerprint.

    Cite this