Our previous paper [1] described and analyzed the performance and some of the advantages obtained when a reservoir simulator is upgraded with: the possibility of flow between any given grid blocks, especially from one to several blocks;a multiple time step calculation. These improvements were appreciated by modelling engineers. The upgraded simulator was used for the simulation of actual reservoirs with complex geometry (such as faulted reservoirs), and for its possibility of having a composite discretization grid (i.e. a grid with several nested levels of local sub-gridding). This practical need justified further extensive study of the flow behaviour near a sub-gridding boundary in a composite grid. Composite grids possess some features which do not exist in regular classical grids. In a regular grid system, when two blocks are neighbours in a certain direction, their grid centers are always drawn up in that direction. In a composite grid, this is not always the case for neighbouring blocks belonging to different levels of local sub-gridding. As a consequence, in a composite grid, approximation of material exchange between two blocks by the classical two point finite difference scheme may be less accurate than in a regular grid. Tests have shown that a significant grid effect may be induced by computing flows between neighbouring blocks belonging to different sub-gridding levels with this scheme. To eliminate these effects, a more sophisticated numerical scheme, involving more than two points, must be used for flow computation. This paper reports an analysis of this phenomenon. It also gives a description of some schemes especially designed for this purpose. These schemes may be viewed as being the classical two point scheme with one of the two points being a grid point but not a grid block center. Reservoir data are then approximated at this new grid point by interpolating data values at the grid block centers. More or less sophisticated interpolations are considered. It also describes some numerical examples performed to compare these schemes when they are used for the simulation of multiphase flows with a composite grid. One of the schemes is selected for its accuracy, its reliability, and its simplicity. On numerical examples, it is shown to nearly reproduce, with a composite grid, results obtained with a fine regular grid having a much larger number of blocks. Moreover, this scheme is implemented on a black oil simulator using a Sequential Solution Method and designed for the simulation of three dimensional three phase flows in large reservoirs with composite grids and several thousands of blocks.

Our previous paper [I] described and analyzed the performance and some of the advantages obtained when a reservoir simulator is upgraded with :