Use of Ensemble Kalman Filter to 3-Dimensional Reservoir Characterization during Waterflooding

Author(s):  
Kwangwon Park ◽  
Jonggeun Choe
SPE Journal ◽  
2006 ◽  
Vol 11 (02) ◽  
pp. 181-192 ◽  
Author(s):  
Erlend H. Vefring ◽  
Gerhard H. Nygaard ◽  
Rolf J. Lorentzen ◽  
Geir Naevdal ◽  
Kjell K. Fjelde

Summary Two methods for characterizing reservoir pore pressure and reservoir permeability during UBD while applying active tests are presented and evaluated. Both methods utilize a fast, dynamic well fluid-flow model that is extended with a transient reservoir model. Active testing of the well is applied by varying the bottomhole pressure in the well during the drilling operations. The first method uses the Levenberg-Marquardt optimization algorithm to estimate the reservoir parameters by minimizing the difference between measurements from the drilling process and the corresponding model states. The method is applied after the drilling process is finished, using all the recorded measurements. The second method is the ensemble Kalman filter, which simulates the drilling process using the dynamic model while drilling is performed, and updates the model states and parameters each time new measurements are available. Measurements are used that usually are available while drilling are used, such as pump rates, pump pressure, bottomhole pressure, and outlet rates. The methods are applied to different cases, and the results indicate that active tests might improve the estimation results. The results also show that both estimation methods give useful results, and that the ensemble Kalman filter calculates these results during the UB operation. Introduction During UBD, the well pressure is kept below the reservoir pore pressure, and reservoir fluids flow into the well. The flow rate from the reservoir depends on the pressure difference between the reservoir pore pressure and the well pressure, in addition to other reservoir parameters, such as permeability and porosity. The viscosity and compressibility of the reservoir fluids also influence the influx rate. The influx of reservoir fluids causes variations in the annulus section of the well, because of changes in well fluid composition and well fluid-flow rate. By measuring some of the fluid-flow parameters of the well, such as pressures changes and rate changes, the reservoir parameters causing the influx might be identified. This is the principal idea that also is the basis for well testing and transient reservoir analysis. Identification of the reservoir properties close to the well gives important information for planning the well-completion design. If highly productive zones can be located, then the use of smart completion can be better utilized. Reservoir characterization during UBD has received attention from several research groups in recent years. Kardolus and van Kruijsdijk (1997) developed a transient reservoir model based on the boundary-element method. This model was compared with a transient analytical reservoir model. One of their findings was that the transient analytical reservoir model could be used for evaluation of the parameters in the reservoir. In a following study, van Kruijsdijk and Cox (1999) presented a method for identifying the permeability in a horizontal reservoir based on measurements of the reservoir inflow. The flow effects caused by the reservoir boundaries were included in the flow calculations.


SPE Journal ◽  
2010 ◽  
Vol 15 (04) ◽  
pp. 1062-1076 ◽  
Author(s):  
A.. Seiler ◽  
S.I.. I. Aanonsen ◽  
G.. Evensen ◽  
J.C.. C. Rivenæs

Summary Although typically large uncertainties are associated with reservoir structure, the reservoir geometry is usually fixed to a single interpretation in history-matching workflows, and focus is on the estimation of geological properties such as facies location, porosity, and permeability fields. Structural uncertainties can have significant effects on the bulk reservoir volume, well planning, and predictions of future production. In this paper, we consider an integrated reservoir-characterization workflow for structural-uncertainty assessment and continuous updating of the structural reservoir model by assimilation of production data. We address some of the challenges linked to structural-surface updating with the ensemble Kalman filter (EnKF). An ensemble of reservoir models, expressing explicitly the uncertainty resulting from seismic interpretation and time-to-depth conversion, is created. The top and bottom reservoir-horizon uncertainties are considered as a parameter for assisted history matching and are updated by sequential assimilation of production data using the EnKF. To avoid modifications in the grid architecture and thus to ensure a fixed dimension of the state vector, an elastic-grid approach is proposed. The geometry of a base-case simulation grid is deformed to match the realizations of the top and bottom reservoir horizons. The method is applied to a synthetic example, and promising results are obtained. The result is an ensemble of history-matched structural models with reduced and quantified uncertainty. The updated ensemble of structures provides a more reliable characterization of the reservoir architecture and a better estimate of the field oil in place.


SPE Journal ◽  
2007 ◽  
Vol 12 (03) ◽  
pp. 382-391 ◽  
Author(s):  
Mohammad Zafari ◽  
Albert Coburn Reynolds

Summary Recently, the ensemble Kalman Filter (EnKF) has gained popularity in atmospheric science for the assimilation of data and the assessment of uncertainty in forecasts for complex, large-scale problems. A handful of papers have discussed reservoir characterization applications of the EnKF, which can easily and quickly be coupled with any reservoir simulator. Neither adjoint code nor specific knowledge of simulator numerics is required for implementation of the EnKF. Moreover, data are assimilated (matched) as they become available; a suite of plausible reservoir models (the ensemble, set of ensemble members or suite or realizations) is continuously updated to honor data without rematching data assimilated previously. Because of these features, the method is far more efficient for history matching dynamic data than automatic history matching based on optimization algorithms. Moreover, the set of realizations provides a way to evaluate the uncertainty in reservoir description and performance predictions. Here we establish a firm theoretical relation between randomized maximum likelihood and the ensemble Kalman filter. Although we have previously generated reservoir characterization examples where the method worked well, here we also provide examples where the performance of EnKF does not provide a reliable characterization of uncertainty. Introduction Our main interest is in characterizing the uncertainty in reservoir description and reservoir performance predictions in order to optimize reservoir management. To do so, we wish to generate a suite of plausible reservoir models (realizations) that are consistent with all information and data. If the set of models is obtained by correctly sampling the pdf, then the set of models give a characterization of the uncertainty in the reservoir model. Thus, by predicting future reservoir performance with each of the realizations, and calculating statistics on the set of outcomes, one can evaluate the uncertainty in reservoir performance predictions.


2015 ◽  
Vol 137 (4) ◽  
Author(s):  
Amit Panwar ◽  
Japan J. Trivedi ◽  
Siavash Nejadi

Distributed temperature sensing (DTS), an optical fiber down-hole monitoring technique, provides a continuous and permanent well temperature profile. In steam assisted gravity drainage (SAGD) reservoirs, the DTS plays an important role to provide depth-and-time continuous temperature measurement for steam management and production optimization. These temperature observations provide useful information for reservoir characterization and shale detection in SAGD reservoirs. However, use of these massive data for automated SAGD reservoir characterization has not been investigated. The ensemble Kalman filter (EnKF), a parameter estimation approach using these real-time temperature observations, provides a highly attractive algorithm for automatic history matching and quantitative reservoir characterization. Due to its complex geological nature, the shale barrier exhibits as a different facies in sandstone reservoirs. In such reservoirs, due to non-Gaussian distributions, the traditional EnKF underestimates the uncertainty and fails to obtain a good production data match. We implemented discrete cosine transform (DCT) to parameterize the facies labels with EnKF. Furthermore, to capture geologically meaningful and realistic facies distribution in conjunction with matching observed data, we included fiber-optic sensor temperature data. Several case studies with different facies distribution and well configurations were conducted. In order to investigate the effect of temperature observations on SAGD reservoir characterization, the number of DTS observations and their locations were varied for each study. The qualities of the history-matched models were assessed by comparing the facies maps, facies distribution, and the root mean square error (RMSE) of the predicted data mismatch. Use of temperature data in conjunction with production data demonstrated significant improvement in facies detection and reduced uncertainty for SAGD reservoirs. The RMSE of the predicted data is also improved. The results indicate that the assimilation of DTS data from nearby steam chamber location has a significant potential in significant reduction of uncertainty in steam chamber propagation and production forecast.


Sign in / Sign up

Export Citation Format

Share Document