Although upscaling has been extensively investigated, a quantitative relationship between the multi-million block geological model and the upscaled flow model is likely to be lost during the history match process, primarily due to parameter modification on an upscaled model. In this paper, we propose a new history match procedure combined with an innovative faster upscaling technique.
In the proposed procedure, after each history run, parameter modification for the next run was done on the detailed geological model in close collaboration with geologist. This loop was made possible with a new faster upscaling technique, distinct against the existing method, by adopting empirical correlations of computing upscaled absolute and relative permeabilities. A capillary-limit upscaling technique was applied for the case study, since a target reservoir is capillary dominant.
In the case study, history match was conducted on a sector of an oil-wet carbonate reservoir under waterflooding for 20 years. In each modification, the revised geological model of three million blocks was upscaled to a flow model of 60 thousand blocks in a process independent manner with a scale-up factor of 50. The CPU speed-up factor was approximately 200: 10 days for the detailed model with massive parallel computing and an hour for the upscaled model with a personal computer.
Even with such significant CPU time reduction, it was confirmed that the upscaled model still maintained the simulation quality of the detailed geological model.
Introduction
Recently, the size of a geological model has been becoming huge, e.g., multi-million-gridblock model, in order to describe the heterogeneity as much as possible. Direct transfer of such a huge geological model (detailed geological model) to a flow model would result in a too long CPU time in a flow simulation and be impossible even for nowadays computer capacity in terms of the required memory.
In general, the objective of upscaling is to reduce the number of gridblocks of a flow model so that a flow simulation can be conducted in a reasonable range of a computation time. The requirement for upscaling is that the important characteristics of heterogeneity and multiphase flow are to be retained in the resultant upscaled flow model made up of the fewer number of gridblocks than the detailed geological model. In other words, upscaling is one of such techniques that assign representative properties to each simulation block of a coarse grid model given that a finely gridded model is available.
In the literature, there are two major upscaling approaches: dynamic approaches and effective properties approaches (Ekrann and Dale 1992). Dynamic approaches adopt pseudo- functions based on the simulation results of an entire reservoir model or representative portion of a fine grid model. Thus, these approaches give accurate results limited to an imposed flow condition from which pseudofunctions are derived. However, the derivation is too time-consuming to upscale very large geological models. In addition, the pseudofunctions have many other deficiencies as criticized by Barker and Thiebeau (1997).
On the other hand, effective properties approaches average flow properties from the information on sub-grid scale heterogeneity alone in which the local flow regime assumed to be extreme: viscous, gravity, or capillary dominated. These approaches, which work on the isolated coarse grid block, enable us to compute effective properties much faster than dynamic approaches and give accurate results if the appropriate assumption is made for the actual flow field.