Efficient History Matching Using a Multiscale Technique
- Sigurd I. Aanonsen (U. of Bergen)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- February 2008
- Document Type
- Journal Paper
- 154 - 164
- 2008. Society of Petroleum Engineers
- 5.1.5 Geologic Modeling, 7.6.2 Data Integration, 5.5.8 History Matching, 5.5.3 Scaling Methods, 5.5 Reservoir Simulation, 4.1.2 Separation and Treating, 5.1 Reservoir Characterisation, 7.2.2 Risk Management Systems, 5.2.1 Phase Behavior and PVT Measurements, 4.3.4 Scale, 5.6.9 Production Forecasting
- 2 in the last 30 days
- 747 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
It is demonstrated that a method for multiscale history matching can be used to improve efficiency and/or quality of the solution when achieving a fine-scale match as compared to history-matching directly on the fine scale. Starting from a given fine-scale realization, coarser models are generated using a global upscaling technique in which the coarse models are "history matched?? with respect to the solution at the fine scale. Conditioning to dynamic data is performed by history matching a coarse model, and this model is then successively refined using a combination of downscaling and history matching until a model that matches dynamic data is obtained at the finest scale. Bias in predicted data because of upscaling errors may be taken into account. The advantage of this procedure is that the large-scale corrections are obtained using fast models which—combined with proper downscaling procedures—provide a better initial model for the final adjustment on the fine scale. Coarse-scale history matching also provides a regularization of the fine-scale match, making the process less dependent on a correct prior model. With the proposed methodology, a series of models with varying degrees of complexity—all being consistent with both static and dynamic data—may be generated without additional cost. Effects of including a priori information and different initial downscaling techniques, such as sampling or block kriging with sequential Gaussian simulation (BKSGS), are investigated using two synthetic reservoir models.
In Aanonsen and Eydinov (2006), a multiscale method was proposed for more-effective history matching of petrophysical properties. Starting from a given fine-scale realization, coarser models are generated by use of a global upscaling technique in which the coarse models are "history matched?? with respect to the solution at the fine scale. Conditioning to dynamic data is accomplished by history matching a coarse model, and this model is then successively refined using a combination of downscaling and history matching until a model that matches dynamic data is obtained at the finest scale. It was demonstrated through some simple examples that this procedure may reduce the computational effort and/or improve the quality of the match as compared to history matching directly on the fine scale. In this paper, the technique is taken further, by taking into account bias in the coarse-scale solution. Also, effects of including a priori information and different initial downscaling techniques, such as sampling and BKSGS, are investigated.
In traditional multiscale estimation, a series of estimations is performed in which the resolution of the zonation is increased for each step in the sequence for examples, see Liu (1993); Chavent and Bissel (1998); and Yoon et al. (1999). This approach reduces both the computational effort and the overparameterization problem. The multiscale estimation was improved by Ben Ameur et al. (2002); Grimstad et al. (2003); and Grimstad and Mannseth (2004), who introduced a method for adaptive multiscale estimation in which the resolution is increased only in some regions of the reservoir of each stage in the sequence. Other types of multiscale approaches include updating a fine model on the basis of simulations performed on a coarse model (Mezghani and Roggero 2001).
Common for all of the above-mentioned multiscale techniques is that the computational grid is kept the same while the parameter resolution is changed. In the method used here, the forward problem is also computed on grids of varying resolution. The main challenge with the method is that it requires an efficient way for tuning the properties in every grid cell in the model. We have used commercial software in which the number of parameters (which, with this method, is equal to the number of active grid cells) is limited to 500. However, with software based on the adjoint method (Li et al. 2003; Gao and Reynolds 2006), this procedure may be applied to much larger cases—up to real field-model size.
|File Size||5 MB||Number of Pages||11|
Aanonsen, S.I. and Eydinov, D. 2006. A Multiscale Method forDistributed Parameter Estimation With Application to Reservoir HistoryMatching. Computational Geosciences 10 (1): 97-117. DOI:10.1007/s10596-005-9012-4.
Behrens, R.A., MacLeod, M.K., Tran, T.T., and Alimi, A.C. 1998. Incorporating Seismic Attribute Mapsin 3D Reservoir Models. SPEREE 1 (2): 122-126. SPE-36499-PA.DOI: 10. 2118/36499-PA.
Ben Ameur, H., Chavent, G., and Jaffré, J. 2002. Refinement and CoarseningIndicators for Adaptive Parameterization: Application to the Estimation ofHydraulic Transmissivities. Inverse Problems 18 (3): 775-794.DOI: 10.1088/0266-5611/18/3/317.
Chavent, G. and Bissel, R. 1998. Indicator for the Refinement ofParameterization. In Inverse Problems in Engineering Mechanics(Proc., International Symposium on Inverse Problems in EngineeringMechanics, Nagano, Japan). M. Tanaka and G.S. Dulikravich, eds. Amsterdam:Elsevier Science B.V., 309-314.
Chen, Y., Durlofsky, L.J., and Wen, X.H. 2004. Robust Coarse Scale Modelingof Flow and Transport in Heterogeneous Reservoirs. Proc., EuropeanConference on the Mathematics of Oil Recovery, Cannes, France, 30 August-2September.
Christie, M.A., Subbey, S., and Sambridge, M. 2002. Prediction UnderUncertainty in Reservoir Modelling. Proc., European Conference on theMathematics of Oil Recovery, Freiberg, Germany, 3-6 September.
Floris, F.J.T., Bush, M.D., Cuypers, M., Roggero, F., and Syversveen, A.-R.2001. Methods for Quantifying the Uncertainty of Production Forecasts: AComparative Study. Petroleum Geoscience 7: S87-S96.
Gao, G. and Reynolds, A.C. 2006. An Improved Implementation of theLBFGS Algorithm for Automatic History Matching. SPEJ 11 (1):5-17. SPE 90058-PA. DOI: 10.2118/90058-PA.
Grimstad, A.-A. and Mannseth, T. 2004. Comparison of Methods for Downscalingof Coarse Scale Permeability Estimates. Proc., European Conference onthe Mathematics of Oil Recovery, Cannes, France, 30 August-2 September.
Grimstad, A.-A., Mannseth, T., Nævdal, G., and Urkedal, H. 2003. Adaptive MultiscalePermeability Estimation. Computational Geosciences 7 (1):1-25. DOI: 10.1023/A:1022417923824.
Holden, L. and Nielsen, B.F. 2000. Global Upscaling ofPermeability in Heterogeneous Reservoirs; The Output Least Squares (OLS)Method. Transport in Porous Media 40 (2): 115-143. DOI:10.1023/A:1006657515753.
Isaaks, E.H. and Srivastava, R.M. 1989. An Introduction to AppliedGeostatistics. New York City: Oxford University Press.
Kitanidis, P.K. 1995. Quasi-Linear Geostatistical Theoryfor Inversing. Water Resources Research 31 (10): 2411-2420.DOI: 10.1029/95WR01945.
Li, R., Reynolds, A.C., and Oliver, D.S. 2003. History Matching of Three-Phase FlowProduction Data. SPEJ 8 (4): 328-340. SPE-87336-PA. DOI:10.2118/87336-PA.
Liu, J. ,1993. A MultiresolutionMethod for Distributed Parameter Estimation. SIAM J. on ScientificComputing 4 (2): 389-405. DOI: 10.1137/0914024.
Mezghani, M. and Roggero, F. 2001. Combining Gradual Deformation andUpscaling Techniques for Direct Conditioning of Fine Scale Reservoir Models toDynamic Data. Paper SPE 71334 presented at the SPE Annual TechnicalConference and Exhibition, New Orleans, 30 September-3 October. DOI:10.2118/71334-MS.
Oliver, D.S., He, N., and Reynolds, A.C. 1996. Conditioning PermeabilityFields to Pressure Data. Proc., European Conference on the Mathematicsof Oil Recovery, Leoben, Austria, 3-6 September.
Omre, H. and Lødøen, O.P. 2004. Improved Production Forecasts andHistory Matching Using Approximate Fluid-Flow Simulators. SPEJ9 (3): 339-351. SPE-74691-PA. DOI: 10.2118/74691-PA.
Pickup, G.E., Monfared, H., Zhang, P., and Christie, M.A. 2004. A New Way ofLooking at Upscaling. Proc., European Conference on the Mathematics ofOil Recovery, Cannes, France, 30 August-2 September.
Sun N.-Z. 1994. Inverse Problems in Groundwater Modeling. Dordrecht,The Netherlands: Kluwer Academic Publishers.
Yoon, S., Malallah, A.H., Datta-Gupta, A., Vasco, D.W., and Behrens, R.A.1999. A Multiscale Approach toProduction Data Integration Using Streamline Models. Paper SPE 56653presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6October. DOI: 10.2118/56653-MS.