NDLI: Iterative regularization for ensemble data assimilation in reservoir models

Please wait, while we are loading the content...

Iterative regularization for ensemble data assimilation in reservoir models

Content Provider	Springer Nature Link
Author	Iglesias, Marco A.
Copyright Year	2014
Abstract	We propose the application of iterative regularization for the development of ensemble methods for solving Bayesian inverse problems. In concrete, we construct (i) a variational iterative regularizing ensemble Levenberg-Marquardt method (IR-enLM) and (ii) a derivative-free iterative ensemble Kalman smoother (IR-ES). The aim of these methods is to provide a robust ensemble approximation of the Bayesian posterior. The proposed methods are based on fundamental ideas from iterative regularization methods that have been widely used for the solution of deterministic inverse problems (Katltenbacher et al. de Gruyter, Berlin 2008). In this work, we are interested in the application of the proposed ensemble methods for the solution of Bayesian inverse problems that arise in reservoir modeling applications. The proposed ensemble methods use key aspects of the regularizing Levenberg-Marquardt scheme developed by Hanke (Inverse Problems 13,79–95 1997) and that we recently applied for history matching in Iglesias (Comput. Geosci. 1–21 2013). Unlike most existing methods where the stopping criteria and regularization parameters are typically selected heuristically, in the proposed ensemble methods, the discrepancy principle is applied for (i) the selection of the regularization parameters and (ii) the early termination of the scheme. The discrepancy principle is key for the theory of iterative regularization, and the purpose of the present work is to apply this principle for the development of ensemble methods defined as iterative updates of solutions to linear ill-posed inverse problems. The regularizing and convergence properties of iterative regularization methods for deterministic inverse problems have long been established. However, the approximation properties of the proposed ensemble methods in the context of Bayesian inverse problems is an open problem. In the case where the forward operator is linear and the prior is Gaussian, we show that the tunable parameters of the proposed IR-enLM and IR-ES can be chosen so that the resulting schemes coincide with the standard randomized maximum likelihood (RML) and the ensemble smoother (ES), respectively. Therefore, the proposed methods sample from the posterior in the linear-Gaussian case. Similar to RML and ES methods, in the nonlinear case, one may not conclude that the proposed methods produce samples from the posterior. The present work provides a numerical investigation of the performance of the proposed ensemble methods at capturing the posterior. In particular, we aim at understanding the role of the tunable parameters that arise from the application of iterative regularization techniques. The numerical framework for our investigations consists of using a state-of-the art Markov chain Monte Carlo (MCMC) method for resolving the Bayesian posterior from synthetic experiments. The resolved posterior via MCMC then provides a gold standard against to which compare the proposed IR-enLM and IR-ES. Our numerical experiments show clear indication that the regularizing properties of the regularization methods applied for the computation of each ensemble have significant impact of the approximation properties of the proposed ensemble methods at capturing the Bayesian posterior. Furthermore, we provide a comparison of the proposed regularizing methods with respect to some unregularized methods that have been typically used in the literature. Our numerical experiments showcase the advantage of using iterative regularization for obtaining more robust and stable approximation of the posterior than unregularized methods.
Starting Page	177
Ending Page	212
Page Count	36
File Format	PDF
ISSN	14200597
Journal	Computational Geosciences
Volume Number	19
Issue Number	1
e-ISSN	15731499
Language	English
Publisher	Springer International Publishing
Publisher Date	2014-12-09
Publisher Place	Cham
Access Restriction	One Nation One Subscription (ONOS)
Subject Keyword	Bayesian data assimilation Inverse problems Reservoir models Earth Sciences Geotechnical Engineering & Applied Earth Sciences Hydrogeology Mathematical Modeling and Industrial Mathematics Soil Science & Conservation
Content Type	Text
Resource Type	Article
Subject	Computers in Earth Sciences Computational Theory and Mathematics Computer Science Applications Computational Mathematics

Sl.	Authority	Responsibilities	Communication Details
1	Ministry of Education (GoI), Department of Higher Education	Sanctioning Authority	https://www.education.gov.in/ict-initiatives
2	Indian Institute of Technology Kharagpur	Host Institute of the Project: The host institute of the project is responsible for providing infrastructure support and hosting the project	https://www.iitkgp.ac.in
3	National Digital Library of India Office, Indian Institute of Technology Kharagpur	The administrative and infrastructural headquarters of the project	Dr. B. Sutradhar bsutra@ndl.gov.in
4	Project PI / Joint PI	Principal Investigator and Joint Principal Investigators of the project	Dr. B. Sutradhar bsutra@ndl.gov.in Prof. Saswat Chakrabarti will be added soon
5	Website/Portal (Helpdesk)	Queries regarding NDLI and its services	support@ndl.gov.in
6	Contents and Copyright Issues	Queries related to content curation and copyright issues	content@ndl.gov.in
7	National Digital Library of India Club (NDLI Club)	Queries related to NDLI Club formation, support, user awareness program, seminar/symposium, collaboration, social media, promotion, and outreach	clubsupport@ndl.gov.in
8	Digital Preservation Centre (DPC)	Assistance with digitizing and archiving copyright-free printed books	dpc@ndl.gov.in
9	IDR Setup or Support	Queries related to establishment and support of Institutional Digital Repository (IDR) and IDR workshops	idr@ndl.gov.in

Evaluation of Gaussian approximations for data assimilation in reservoir models

Iterative Bayesian inversion with Gaussian mixtures: finite sample implementation and large sample asymptotics

An iterative ensemble Kalman filter for reservoir engineering applications

The regularizing Levenberg–Marquardt scheme for history matching of petroleum reservoirs

Multimodal ensemble Kalman filtering using Gaussian mixture models

Assessment of ordered sequential data assimilation

Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data

An iterative version of the adaptive Gaussian mixture filter

Multiscale ensemble filtering for reservoir engineering applications

Iterative regularization for ensemble data assimilation in reservoir models

Similar Documents

Evaluation of Gaussian approximations for data assimilation in reservoir models

Iterative Bayesian inversion with Gaussian mixtures: finite sample implementation and large sample asymptotics

An iterative ensemble Kalman filter for reservoir engineering applications

The regularizing Levenberg–Marquardt scheme for history matching of petroleum reservoirs

Multimodal ensemble Kalman filtering using Gaussian mixture models

Assessment of ordered sequential data assimilation

Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data

An iterative version of the adaptive Gaussian mixture filter

Multiscale ensemble filtering for reservoir engineering applications

Iterative regularization for ensemble data assimilation in reservoir models