Multifidelity adaptive sequential Monte Carlo for geophysical inversion
Détails
Télécharger: ggae040 (5).pdf (2156.73 [Ko])
Etat: Public
Version: Final published version
Licence: CC BY 4.0
Etat: Public
Version: Final published version
Licence: CC BY 4.0
ID Serval
serval:BIB_1BB11C809A87
Type
Article: article d'un périodique ou d'un magazine.
Collection
Publications
Institution
Titre
Multifidelity adaptive sequential Monte Carlo for geophysical inversion
Périodique
Geophysical Journal International
ISSN
0956-540X
1365-246X
1365-246X
Statut éditorial
Publié
Date de publication
07/03/2024
Peer-reviewed
Oui
Volume
237
Numéro
2
Pages
788-804
Langue
anglais
Résumé
In the context of Bayesian inversion, we consider sequential Monte Carlo (SMC) methods that provide an approximation of the posterior probability density function and the evidence (marginal likelihood). These particle approaches build a sequence of importance sampling steps between gradually tempered distributions evolving from the prior to the posterior PDF. To automate the definition of the tempering schedule, adaptive SMC (ASMC) allows tuning the temperature increments on-the-go. One general challenge in Bayesian inversions is the computational burden associated with expensive, high-fidelity forward solvers. Lower-fidelity surrogate models are interesting in this context as they can emulate the response of expensive forward solvers at a fraction of their cost. We consider surrogate modelling within ASMC and introduce first an approach involving surrogate modelling only, in which either prior samples are used to train the surrogate, or the surrogate model is retrained by updating the training set during the inversion. In our implementation, we rely on polynomial chaos expansions for surrogate modelling, principal component analysis for model parametrization and a ground-penetrating radar cross-hole tomography problem with either an eikonal or finite-difference time-domain solver as high-fidelity solver. We find that the method based on retraining the surrogate during the inversion outperforms the results obtained when only considering prior samples. We then introduce a computationally more expensive multifidelity approach including a transition to the high-fidelity forward solver at the end of the surrogate-based ASMC run leading to even more accurate results. Both methods result in speed-ups that are larger than one order of magnitude compared to standard high-fidelity ASMC inversion.
Web of science
Open Access
Oui
Financement(s)
Fonds national suisse
Création de la notice
01/11/2024 12:07
Dernière modification de la notice
01/11/2024 14:05