Journal of Canadian Petroleum Technology, Vol.54, No.6, 424-441, 2015
Initial Sampling of Ensemble for Steam-Assisted-Gravity-Drainage-Reservoir History Matching
History matching of a steam-assisted-gravity-drainage (SAGD) reservoir requires a large ensemble size for proper uncertainty assessment, which ultimately results in high computational cost. Therefore, it is necessary to reduce the number of realizations for SAGD-reservoir simulation purposes. In this paper, a novel sampling method (based on the probability-distance-minimization method) to generate an initial ensemble of reduced size is discussed. This method considers multiple static measurements and geological properties and uses Kantorovich distance to quantify the probability distance between the original ensemble and the reduced ensemble, which is later optimized by use of the mixed-integer linear-optimization (MILP) technique. To show the effectiveness of the method, we have shown history matching of an SAGD reservoir using the smaller size initial ensemble derived from the proposed method and compared with the original ensemble. For history matching, the ensemble Kalman filter (EnKF) has been used because of its ability to assimilate data for large-scale nonlinear systems. Results are compared with several other methods, such as importance sampling, kernel K-means clustering, and sampling by use of orthogonal ensemble members. The robustness and usefulness of each method for generating an improved initial ensemble of reduced size are analyzed on the basis of two criteria: (1) Does the smaller ensemble retain the same statistical distribution characteristics as the original ensemble, and (2) does the smaller ensemble improve the performance of history matching? In general, we conclude that the improved, smaller initial ensemble created by use of the proposed method retains the best statistical characteristics of the original ensemble. Also, it provides better performance compared with other ranking methods in sampling and history matching using EnKF. Finally, the proposed method can reduce the computing cost significantly without compromising uncertainty in the forecast model, which allows for real-time updating at smaller time intervals.