Adaptive sequential Monte Carlo for posterior inference and model selection among complex geological priors

Details

Serval ID
serval:BIB_5E4F452BF1DF
Type
Article: article from journal or magazin.
Collection
Publications
Institution
Title
Adaptive sequential Monte Carlo for posterior inference and model selection among complex geological priors
Journal
Geophysical Journal International
Author(s)
Amaya M, Linde N, Laloy E
ISSN
0956-540X
1365-246X
Publication state
Published
Issued date
19/05/2021
Volume
226
Number
2
Pages
1220-1238
Language
english
Abstract
Bayesian model selection enables comparison and ranking of conceptual subsurface models described by spatial prior models, according to the support provided by available geophysical data. Deep generative neural networks can efficiently encode such complex spatial priors, thereby, allowing for a strong model dimensionality reduction that comes at the price of enhanced non-linearity. In this setting, we explore a recent adaptive sequential Monte Carlo (ASMC) approach that builds on annealed importance sampling (AIS); a method that provides both the posterior probability density function (PDF) and the evidence (a central quantity for Bayesian model selection) through a particle approximation. Both techniques are well suited to parallel computation and rely on importance sampling over a sequence of intermediate distributions, linking the prior and the posterior PDF. Each subsequent distribution is approximated by updating the particle weights and states, compared with the previous approximation, using a small pre-defined number of Markov chain Monte Carlo (MCMC) proposal steps. Compared with AIS, the ASMC method adaptively tunes the tempering between neighboring distributions and performs resampling of particles when the variance of the particle weights becomes too large. We evaluate ASMC using two different conceptual models and associated synthetic cross-hole ground penetrating radar tomography data. For the most challenging test case, we find that the ASMC method is faster and more reliable in locating the posterior PDF than state-of-the-art adaptive MCMC. The evidence estimates are found to be robust with respect to the choice of ASMC algorithmic variables and much less sensitive to the model proposal type than MCMC. The variance of the evidence estimates are best estimated by replication of ASMC runs, while approximations based on single runs provide comparable estimates when using a sufficient number of proposal steps in approximating each intermediate distribution.
Keywords
Geochemistry and Petrology, Geophysics
Funding(s)
Swiss National Science Foundation / 184574
Create date
03/08/2021 21:42
Last modification date
04/08/2021 6:37
Usage data