Research - Calculation of Initial Fluid Distribution in Oil Reservoirs (TP 2405, Petr. Tech., July 1948)

It is pointed out that the application of capillary pressure curves obtained by drainage or desaturation processes to the calculation of the fluid distribution in interphase transition zones involves a number of difficulties; namely: (I) the development of very low nonwetting phase saturations appears to be in contradiction with the lack of mobility of such distributions indicated by permeability-saturation curves, (2) dispersed nonwetting phases are thermodynamically unstable, and (3) discontinuous phases should not be subject to hydrostatic equilibrium requirements. While these difficulties could be obviated by assuming that the capillary pressure drainage curve has an initial horizontal segment, they are automatically circumvented by application of the imbibition capillary pressure curve to the lower part of the water-oil transition zone. These generally show zero displacement pressure at only partial wetting-phase saturations. The countercurrent upward flow of oil into the main oil-saturated pay and downward drainage of water also suggests that wetting phase imbibition processes will control the saturation distribution immediately above the water-saturated section. Similar considerations, with a generalized interpretation of the apparent wetting-phase behavior of the oil and gas phases, provides a basis for constructing the curve for fluid distribution in the oil-gas transition zone. In the transition zones so derived, the oil begins with a nonvanishing saturation at the water-oil contact and terminates with a similar saturation at the top of the gas-oil contact. The gas-oil transition zone begins with an equivalent nonvanishing gas saturation. Introduction It has been a common assumption for some years that the nature of the fluid dis- tribution in virgin reservoirs, and in particular that in the transition zones between the oil and water and between the oil and gas sections, can be computed by a simple application of capillary pressure data. The results of such calculations appear to have been first reported by Leverett.' The latter, however, indicated only the numerical values of the parameters used in the computations, without explicitly describing the procedure. Leverett's illustrative calculated transition zones are reproduced in Fig I. Although no critical study of the apparently obvious method of calculation has been published, one cannot proceed very far in the computation without encountering some rather fundamental questions not yet answered in the literature. This paper is not written under the pretense that the whole transition-zone problem has been fully and satisfactorily solved. Its purpose, rather, is to discuss suggestions for its treatment, much of it admittedly hypothetical, and emphasize the nature of the uncertainties arising therein; for the literature, as now available, gives little indication that there are any problems still outstanding. The Transition-zone Equation The basic equation, which presumably gives directly the fluid distribution without recourse to any further consideration, is: Ap = (r1 - r2)gh [1] where r1, y2 are the densities of the two contiguous fluids, h is the height above the "capillary-free" zone of complete saturation with fluid I, g the acceleration of gravity, and Ap is the capillary pressure across the average interface at the height