Numerical Simulation of Water-Oil Flow in Naturally Fractured Reservoirs
Authors: Kazemi H. , Merrill L.S. , Porterfield K.L. , Zeman P.R. 
Published: December 2007

Abstract
A three-dimensional, multiple-well, numerical simulator for simulating single- or two-phase flow of water and oil is developed for fractured reservoirs. The simulator equations are two-phase flow extensions of the single-phase flow equations derived by Warren and Root. The simulator accounts for relative fluid mobilities, gravity force, imbibition, and variation in reservoir properties. The simulator handles uniformly and nonuniformly properties. The simulator handles uniformly and nonuniformly distributed fractures and for no fractures at all. The simulator can be used to simulate the water-oil displacement process and in the transient testing of fractured reservoirs. The simulator was used on the conceptual models of two naturally fractured reservoirs: a quadrant of a five-spot reservoir and a live-well dipping reservoir with water drive. These results show the significance of imbibition in recovering oil from the reservoir rock in reservoirs with an interconnected fracture network.

Introduction
Numerical reservoir simulators are being used extensively to simulate multiphase, multicomponent flow in "single-porosity" petroleum reservoirs. Such simulators generally cannot be used to petroleum reservoirs. Such simulators generally cannot be used to study flow behavior in the naturally fractured reservoirs that are usually classified as double-porosity systems. In the latter, one porosity is associated with the matrix blocks and the other porosity is associated with the matrix blocks and the other represents that of the fractures and vugs. If fractures provide the main path for fluid flow from the reservoir, then usually the oil from the matrix blocks flows into the fracture space, and the fractures carry the oil to the wellbore. When water comes in contact with the oil zone, water may imbibe into the matrix blocks to displace oil. Combinations of large flow rates, low matrix permeability, and weak imbibition may result in water fingering permeability, and weak imbibition may result in water fingering through the fractures into the wellbore. Once fingering of water occurs, the water-oil ratio may increase to a large value. None of the published theoretical work on multiphase flow in naturally fractured systems has been applied directly to the simulation of a reservoir as a whole. Usually, only a segment of the reservoir was simulated, and the results were extrapolated to the entire reservoir. To simulate a reservoir as a whole, we have developed a mathematical formulation of the flow problem that has been programmed as a three-dimensional, compressible, water-oil reservoir simulator. The simulator equations are two-phase flow extensions of the single-phase flow equations derived by Warren and Root. The theory is based on the assumption of double porosity at each point in a manner that the fractures form a continuum filled by the noncontinuous matrix blocks. In other words, the fractures are the boundaries of the matrix blocks. The flow equations are solved by a finite difference method. A typical finite-difference grid cell usually contains one or several matrix blocks. In this case, all the matrix blocks within the finite-difference grid cell have the same pressure and saturation. Gravity segregation within individual matrix blocks is not calculated, but the over-all gravity segregation from one grid cell to another is accounted for. In many practical problems, this approximation is acceptable. In some situations, a matrix block encloses several finite-difference grid cells. In this case, the gravity segregation within the matrix block is calculated. To include heterogeneity, a redefinition of local porosities and permeabilities provides a method for simulating situations where part of the reservoir is fractured and where part is not fractured. The above description points to the complexity of the situations that one encounters. Therefore, the judicious choice of the number of finite difference grid cells with respect to the number of matrix blocks becomes a critical engineering decision. Later sections will provide insight to alleviate such decisions.
SPEJ
P. 317