Adaptive regularisation for ensemble Kalman inversion

Abstract

We propose a new regularisation strategy for the classical ensemble Kalman inversion (EKI) framework. The strategy consists of: (i) an adaptive choice for the regularisation parameter in the update formula in EKI, and (ii) criteria for the early stopping of the scheme. In contrast to existing approaches, our parameter choice does not rely on additional tuning parameters which often have severe effects on the efficiency of EKI. We motivate our approach using the interpretation of EKI as a Gaussian approximation in the Bayesian tempering setting for inverse problems. We show that our parameter choice controls the symmetrised Kulback-Leibler divergence between consecutive tempering measures. We further motivate our choice using a heuristic statistical discrepancy principle. We test our framework using electrical impedance tomography with the complete electrode model. Parameterisations of the unknown conductivity are employed which enable us to characterise both smooth or a discontinuous (piecewise-constant) fields. We show numerically that the proposed regularisation of EKI can produce efficient, robust and accurate estimates, even for the discontinuous case which tends to require larger ensembles and more iterations to converge. We compare the proposed technique with a standard method of choice and demonstrate that the proposed method is a viable choice to address computational efficiency of EKI in practical/operational settings.

References in corpus (11)

• The Ensemble Kalman Filter for Inverse Problems
• Quantifying and Reducing Model-Form Uncertainties in Reynolds-Averaged Navier-Stokes Simulations: A Data-Driven, Physics-Based Bayesian Approach
• Ensemble Kalman Inversion: A Derivative-Free Technique For Machine Learning Tasks
• On the Convergence of the Ensemble Kalman Filter
• Parameterizations for Ensemble Kalman Inversion
• Convergence of the Square Root Ensemble Kalman Filter in the Large Ensemble Limit
• Well Posedness and Convergence Analysis of the Ensemble Kalman Inversion
• Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems
• Bayesian inversion in resin transfer molding
• Never look back - A modified EnKF method and its application to the training of neural networks without back propagation
• Accelerating Parallel Tempering: Quantile Tempering Algorithm (QuanTA)