Transform-based particle filtering for elliptic Bayesian inverse problems

Abstract

We introduce optimal transport based resampling in adaptive SMC. We consider elliptic inverse problems of inferring hydraulic conductivity from pressure measurements. We consider two parametrizations of hydraulic conductivity: by Gaussian random field, and by a set of scalar (non-)Gaussian distributed parameters and Gaussian random fields. We show that for scalar parameters optimal transport based SMC performs comparably to monomial based SMC but for Gaussian high-dimensional random fields optimal transport based SMC outperforms monomial based SMC. When comparing to ensemble Kalman inversion with mutation (EKI), we observe that for Gaussian random fields, optimal transport based SMC gives comparable or worse performance than EKI depending on the complexity of the parametrization. For non-Gaussian distributed parameters optimal transport based SMC outperforms EKI.

References in corpus (6)

• A regularizing iterative ensemble Kalman method for PDE-constrained inverse problems
• Parameterizations for Ensemble Kalman Inversion
• Well-Posed Bayesian Geometric Inverse Problems Arising in Subsurface Flow
• Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems
• Solving Large-Scale PDE-constrained Bayesian Inverse Problems with Riemann Manifold Hamiltonian Monte Carlo
• Long-time stability and accuracy of the ensemble Kalman-Bucy filter for fully observed processes and small measurement noise