arrow
Volume 33, Issue 5
Adaptive Ensemble Kalman Inversion with Statistical Linearization

Yanyan Wang, Qian Li & Liang Yan

Commun. Comput. Phys., 33 (2023), pp. 1357-1380.

Published online: 2023-06

Export citation
  • Abstract

The ensemble Kalman inversion (EKI), inspired by the well-known ensemble Kalman filter, is a derivative-free and parallelizable method for solving inverse problems. The method is appealing for applications in a variety of fields due to its low computational cost and simple implementation. In this paper, we propose an adaptive ensemble Kalman inversion with statistical linearization (AEKI-SL) method for solving inverse problems from a hierarchical Bayesian perspective. Specifically, by adaptively updating the unknown with an EKI and updating the hyper-parameter in the prior model, the method can improve the accuracy of the solutions to the inverse problem. To avoid semi-convergence, we employ Morozov’s discrepancy principle as a stopping criterion. Furthermore, we extend the method to simultaneous estimation of noise levels in order to reduce the randomness of artificially ensemble noise levels. The convergence of the hyper-parameter in prior model is investigated theoretically. Numerical experiments show that our proposed methods outperform the traditional EKI and EKI with statistical linearization (EKI-SL) methods.

  • AMS Subject Headings

62F15, 65C35, 65N21, 65N75

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-33-1357, author = {Wang , YanyanLi , Qian and Yan , Liang}, title = {Adaptive Ensemble Kalman Inversion with Statistical Linearization}, journal = {Communications in Computational Physics}, year = {2023}, volume = {33}, number = {5}, pages = {1357--1380}, abstract = {

The ensemble Kalman inversion (EKI), inspired by the well-known ensemble Kalman filter, is a derivative-free and parallelizable method for solving inverse problems. The method is appealing for applications in a variety of fields due to its low computational cost and simple implementation. In this paper, we propose an adaptive ensemble Kalman inversion with statistical linearization (AEKI-SL) method for solving inverse problems from a hierarchical Bayesian perspective. Specifically, by adaptively updating the unknown with an EKI and updating the hyper-parameter in the prior model, the method can improve the accuracy of the solutions to the inverse problem. To avoid semi-convergence, we employ Morozov’s discrepancy principle as a stopping criterion. Furthermore, we extend the method to simultaneous estimation of noise levels in order to reduce the randomness of artificially ensemble noise levels. The convergence of the hyper-parameter in prior model is investigated theoretically. Numerical experiments show that our proposed methods outperform the traditional EKI and EKI with statistical linearization (EKI-SL) methods.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2023-0012}, url = {http://global-sci.org/intro/article_detail/cicp/21764.html} }
TY - JOUR T1 - Adaptive Ensemble Kalman Inversion with Statistical Linearization AU - Wang , Yanyan AU - Li , Qian AU - Yan , Liang JO - Communications in Computational Physics VL - 5 SP - 1357 EP - 1380 PY - 2023 DA - 2023/06 SN - 33 DO - http://doi.org/10.4208/cicp.OA-2023-0012 UR - https://global-sci.org/intro/article_detail/cicp/21764.html KW - Ensemble Kalman inversion, statistical linearization, adaptive, Bayesian inverse problem. AB -

The ensemble Kalman inversion (EKI), inspired by the well-known ensemble Kalman filter, is a derivative-free and parallelizable method for solving inverse problems. The method is appealing for applications in a variety of fields due to its low computational cost and simple implementation. In this paper, we propose an adaptive ensemble Kalman inversion with statistical linearization (AEKI-SL) method for solving inverse problems from a hierarchical Bayesian perspective. Specifically, by adaptively updating the unknown with an EKI and updating the hyper-parameter in the prior model, the method can improve the accuracy of the solutions to the inverse problem. To avoid semi-convergence, we employ Morozov’s discrepancy principle as a stopping criterion. Furthermore, we extend the method to simultaneous estimation of noise levels in order to reduce the randomness of artificially ensemble noise levels. The convergence of the hyper-parameter in prior model is investigated theoretically. Numerical experiments show that our proposed methods outperform the traditional EKI and EKI with statistical linearization (EKI-SL) methods.

Yanyan Wang, Qian Li & Liang Yan. (2023). Adaptive Ensemble Kalman Inversion with Statistical Linearization. Communications in Computational Physics. 33 (5). 1357-1380. doi:10.4208/cicp.OA-2023-0012
Copy to clipboard
The citation has been copied to your clipboard