Abstract
In this paper, we focus on parameter estimation for an elliptic inverse problem. We consider a 2D steady-state single-phase Darcy flow model, where permeability and boundary conditions are uncertain. Permeability is parameterized by the Karhunen-Loeve expansion and thus assumed to be Gaussian distributed. We employ two ensemble-based data assimilation methods: ensemble Kalman filter and ensemble transform particle filter. The formal one approximates mean and variance of a Gaussian probability function by means of an ensemble. The latter one transforms ensemble members to approximate any posterior probability function. Ensemble Kalman filter considered here is employed with regularization and localization—R(L)EnKF. Ensemble transform particle filter is also employed with a form of regularization called tempering and localization—T(L)ETPF. Regularization is required for highly non-linear problems, where prior is updated to posterior via a sequence of intermediate probability measures. Localization is required for small ensemble sizes to remove spurious correlations. We have shown that REnKF outperforms TETPF. We have shown that localization improves estimations of both REnKF and TETPF. In numerical experiments when uncertainty is only in permeability, TLETPF outperforms RLEnKF. When uncertainty is both in permeability and in boundary conditions, TLETPF outperforms RLEnKF only for a large ensemble size 1000. Furthermore, when uncertainty is both in permeability and in boundary conditions but we do not account for error in boundary conditions in data assimilation, RLEnKF outperforms TLETPF.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Doucet, A., de Freitas, N., Gordon, N.: Sequential Monte-Carlo Methods in Practice. Springer, New York (2001)
Emerick, A.A., Reynolds, A.C.: Ensemble smoother with multiple data assimilation. Comput. Geosci. 55, 3–15 (2013)
Iglesias, M.A.: Iterative regularization for ensemble data assimilation in reservoir models. Comput. Geosci. 19(1), 177–212 (2015)
Iglesias, M.A.: A regularizing iterative ensemble Kalman method for PDE-constrained inverse problems. Inverse Probl. 32, 025002 (2016)
Acevedo, W., de Wiljes, J., Reich, S.: Second-order accurate ensemble transform particle filters. SIAM J. Sci. Comput. 39(5), A1834–A1850 (2017)
Daniela, C., Dunlop, M., Somersalo, E., Stuart A.: Iterative updating of model error for Bayesian inversion. Inverse Probl. 34(2), 025008 (2018)
Oliver, D.: Metropolized randomized maximum likelihood for improved sampling from multimodal distributions. SIAM/ASA J. Uncertain. Quant. 5(1), 259–277 (2017)
Evensen, G.: Analysis of iterative ensemble smoothers for solving inverse problems. Comput. Geosci. 22(3), 885–908 (2018)
Gaspari, G., Cohn, S.E.: Construction of correlation functions in two and three dimensions. Q. J. Roy. Meteorol. Soc. 125(554), 723–757 (1999)
Penny, S.G.: Miyoshi T. A local particle filter for high-dimensional geophysical systems. Nonlin. Processes Geophys. 23, 391–405 (2016)
Carrera, J., Neuman, S.P.: Estimation of aquifer parameters under transient and steady state conditions: III. Application to synthetic and field data. Water Resour. Res., 22 (1986)
Kaipio, J., Somersalo, E.: Statistical and Computational Inverse Problems. Springer, New York (2004)
Poterjoy, J.A.: Localized particle filter for high-dimensional nonlinear systems. Mon. Wea. Rev. 144(1), 59–76 (2016)
Cotter, S.L., Roberts, G.O., Stuart, A.M., White, D.: MCMC methods for functions: modifying old algorithms to make them faster. Stat. Sci. 28, 424–446 (2013)
Houtekamer, P.L., Lefaivre, L., Derome, J., Ritchie, H., Mitchell, H.L.: A system simulation approach to ensemble prediction. Mon. Wea. Rev. 124, 1225–1242 (1996)
Houtekamer, P.L., Mitchell H. L.: A sequential ensemble Kalman filter for atmospheric data assimilation. Mon. Wea. Rev. 129, 123–137 (2001)
Luo, X., Stordal, A.S., Lorentzen, R.J., Naevdal, G.: Iterative ensemble smoother as an approximate solution to a regularized minimum-average-cost problem: Theory and applications. Soc. Petrol. Eng. J. 20(5), 962–982, 10 (2015)
Luo, X.: Ensemble-based kernel learning for a class of data assimilation problems with imperfect forward simulators. PLOS ONE 14(7), 1–40, 07 (2019)
Hamill, T.M., Whitaker, J.S.: Distance-dependent filtering of background error covariance estimates in an ensemble Kalman. Filter. Mon. Wea Rev. 129, 2776–2790 (2001)
Hanke, M.: A regularizing Levenberg-Marquardt scheme, with applications to inverse groundwater filtration problems. Inverse Probl. 13, 79–95 (1997)
Neal, R.M.: Annealed importance sampling. Stat. Comput. 11(2), 125–139 (2001)
Ma, X., Bi, L.: A robust adaptive iterative ensemble smoother scheme for practical history matching applications. Comput. Geosci. 23(3), 415–442 (2019)
Kantas, N., Beskos, A., Jasra, A.: Sequential Monte Carlo for inverse problems: a case study for the Navier Stokes equation. SIAM/ASA J. Uncertain. Quant. 2, 464–489 (2014)
Pele, O., Werman, M.: Fast and robust earth mover’s distances. In: IEEE 12th international conference on computer vision, pp. 460–467. IEEE (2009)
Sakov, P., Haussaire, J.-M., Bocquet, M.: An iterative ensemble Kalman filter in the presence of additive model error. Q. J. Roy. Meteorol. Soc. 144(713), 1297–1309 (2018)
Dubinkina, S.: Relevance of conservative numerical schemes for an ensemble Kalman filter. Q. J. Roy. Meteorol. Soc. 144(711), 468–477 (2018)
Reich, S., Cotter, C.: Probabilistic Forecasting and Bayesian Data Assimilation. Cambridge University Press, Cambridge (2015)
Ruchi, S., Dubinkina, S.: Application of ensemble transform data assimilation methods for parameter estimation in reservoir modeling Nonlin. Process. Geophys. 25, 731–746 (2018)
Ruchi, S., Dubinkina, S., Iglesias, M.A.: Transform-based particle filtering for elliptic Bayesian inverse problems (2018)
Acknowledgements
We would like to acknowledge Marco Iglesias (U. of Nottingham) for fruitful discussions about sources of model error and for providing a code to model the Darcy flow and an MCMC code.
Funding
This work is part of the research programme Shell-NWO/FOM Computational Sciences for Energy Research (CSER) with project number 14CSER007 which is partly financed by the Netherlands Organization for Scientific Research (NWO).
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Dubinkina, S., Ruchi, S. Comparison of regularized ensemble Kalman filter and tempered ensemble transform particle filter for an elliptic inverse problem with uncertain boundary conditions. Comput Geosci 24, 149–160 (2020). https://doi.org/10.1007/s10596-019-09904-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10596-019-09904-w