Ensemble Subsurface Modeling Using Grid Computing Technology

Ensemble Kalman filter (EnKF) uses a randomized ensemble of subsurface models for error and uncertainty estimation. However, the complexity of geological models and the requirement of a large number of simulation runs make routine applications extremely difficult due to expensive computation cost. Grid computing technologies provide a cost-efficient way to combine geographically distributed computing resources to solve large-scale data and computation intensive problems. Hence, we design and implement a grid-enabled EnKF solution to ill-posed model inversion problems for subsurface modeling. It has been integrated into the ResGrid, a problem solving environment aimed at managing distributed computing resources and conducting subsurface-related modeling studies. Two synthetic cases in reservoir studies indicate that the enhanced ResGrid efficiently performs EnKF inversions to obtain accurate, uncertainty-ware predictions on reservoir production. This grid-enabled EnKF solution is also being applied for data assimilation of large-scale groundwater hydrology nonlinear models. The ResGrid with EnKF solution is open-source and available for downloading.

[1]  Dean S. Oliver,et al.  History Matching of the PUNQ-S3 Reservoir Model Using the Ensemble Kalman Filter , 2005 .

[2]  Dean S. Oliver,et al.  Three-dimensional reservoir description from multiwell pressure data and prior information , 1997 .

[3]  Ian T. Foster,et al.  The Anatomy of the Grid: Enabling Scalable Virtual Organizations , 2001, Int. J. High Perform. Comput. Appl..

[4]  Frank T.-C. Tsai,et al.  Geophysical parameterization and parameter structure identification using natural neighbors in groundwater inverse problems , 2005 .

[5]  Xin Li,et al.  ResGrid: A Grid-aware Toolkit for Reservoir Uncertainty Analysis , 2006, Sixth IEEE International Symposium on Cluster Computing and the Grid (CCGRID'06).

[6]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[7]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[8]  Roland N. Horne,et al.  A Procedure to Integrate Well Test Data, Reservoir Performance History and 4-D Seismic Information into a Reservoir Description , 1997 .

[9]  A. Reynolds,et al.  Optimization Algorithms for Automatic History Matching of Production Data , 2002 .

[10]  Dean S. Oliver,et al.  Sensitivity Coefficients for Three-Phase Flow History Matching , 2001 .

[11]  Daniel S. Katz,et al.  A comparison of two methods for building astronomical image mosaics on a grid , 2005, 2005 International Conference on Parallel Processing Workshops (ICPPW'05).

[12]  Peter S. Maybeck,et al.  Stochastic Models, Estimation And Control , 2012 .

[13]  R. Olea Geostatistics for Natural Resources Evaluation By Pierre Goovaerts, Oxford University Press, Applied Geostatistics Series, 1997, 483 p., hardcover, $65 (U.S.), ISBN 0-19-511538-4 , 1999 .

[14]  Yan Chen,et al.  Dynamic Reservoir Data Assimilation With an Efficient, Dimension-Reduced Kalman Filter , 2007 .