The conventional reason for upscaling in reservoir simulation is the computational limit of the simulator. However, we argue that, from a system-theoretical point of view, a more fundamental reason is that there is only a limited amount of information (output) that can be observed from production data, while there is also a limited amount of control (input) that can be exercised by adjusting the well parameters; in other words, the input-output behavior is usually of much lower dynamical order than the number of grid blocks in the model. Therefore, we propose an upscaling approach to find a coarse model that optimally describes the input-output behavior of a reservoir system. In this control-relevant method, the coarse-scale model parameters are calculated as the solution of an optimization problem that minimizes the distance between the input-output behavior of the fine- and coarse-scale models. This distance is measured with the aid of the Hankel- or energy norms, where Hankel singular values are used as a measure of the combined controllability and observability of the system. The method is particularly attractive to scale up simulation models in flooding optimization and/or history matching studies for a given configuration of wells. An advantage of our upscaling method is that it intrinsically relies most heavily on those parameter values that directly influence the input-output behavior. It is a global method, in the sense that it relies on the system properties of the entire reservoir. It does not, however, require any forward simulation, neither of the full nor of the upscaled model. it also does not depends on a particular control strategy, but instead uses the dynamical system equations directly. A drawback, however, is its dependency on well positions, which implies that it should be (partially) repeated when those positions are changed. We tested the method on several examples and for nearly all cases obtained coarse scale models with a superior input-output behavior compared to common upscaling algorithms.
Introduction
Simulation of flow in porous media is an important tool to predict and optimize reservoir performance. As a precursor to the simulation, reservoir data are collected from different sources with various temporal and spatial scales, which are usually integrated into detailed geological models with a large number of cells (107 to 108). The large size of geo-models makes it often impossible to undertake reasonable multiphase flow simulations, sensitivity analysis, uncertainty assessment of a large number of realizations, computer-assisted history matching or life-cycle optimization. In order to perform these processes within a practical time frame and/or a reasonable cost, an upscaling procedure is required to transfer the flow and transport processes from the detailed fine-scale model to a lower-order, coarser representation. The approximated coarse-scale model (typically up to 106 cells) should reduce the complexity, while retaining the main properties of the reservoir system.
Upscaling background
Various upscaling techniques have been developed in both hydrology and petroleum engineering. These methods vary from simple averaging methods on Cartesian cells to sophisticated flow-based techniques on unstructured grids. Extensive reviews on upscaling were written by e.g. Wen et al. (1996), Renard and Marsily (1997) and Durlofsky (2005). Upscaling methods are classified based on the type of the parameters that are scaled up. In single-phase upscaling, we only consider the pressure (flow) equation and calculate equivalent permeability and porosity values, whereas in two-phase upscaling both pressure and transport equations are considered. Therefore, in addition to the single-phase parameters, equivalent relative permeability and capillary pressure curves are calculated. In principle, upscaling includes both parameters and equations. However, in most single-phase upscaling procedures, coarse-scale equations are chosen to have the same mathematical form as the fine-scale equations.