Core Analysis - The Effect of Permeability Stratification in Complete Water-Drive Systems

A theory is presented for calculating the performance history of complete water-drive systems producing from idealized stratified formations. The general equations are applied to systems where the permeability stratification is either of the exponential or linear type. Calculations were carried through for different degrees of permeability stratification, but with special emphasis on the effect of the mobility ratio between the produced oil and the invading water on the resultant performance. These results are also expressed graphically as curves for the initial water breakthrough recovery, for the different degrees of stratification, as a function of the mobility ratio, and of the composition of the produced fluid stream as a function of the cumulative oil recovery. For several typical cases the latter has also been plotted as a function of the cumulative oil and water throughflow. The general result is that when the mobility of the oil is lower than that of the invading water the channelling tendency resulting from the permeability stratification becomes aggravated as the higher permeability zones become flooded out. Situations of this type would obtain when producing low gravity or highly viscous oils. Conversely, if the mobility of the oil is high compared to that of the invading water, the flooding of the high permeability zones will lead to a retarding and choking effect, and the gross bypassing phenomena will be partially suppressed. These conditions would correspond to those of flooding high gravity or low viscosity oils. A discussion is given of the various basic assumptions made in the analysis, including that of ignoring the stripping phase of the production history as implied by relative permeability concepts. INTRODUCTION The physical ultimate recoveries from oil reservoirs are basically determined and limited by the physical oil displacement processes associated with the reservoir producing mechanism. In practice, however, the economic ultimate recoveries are further limited by the mobility of the reservoir fluids and the uniformity and continuity of the producing formation. In fact, it is the differential depletion between the component parts of the composite reservoir which ultimately determines the total recovery at the time of field abandonment. While this observation applies to both the solution gas drive and gravity drainage mechanisms, in which use is made only of the energy contained within the original oil-bearing reservoir, it is of even more paramount importance under operations wherein the energy associated with extraneous fluids provides the ultimate oil expulsion mechanism. Whether the invading fluid is the water from an edgewater drive, water injected for pressure maintenance, gas injected for pressure maintenance, or gas returned to the formation in a cycling program, it is often the continuity and uniformity of the producing section which will control the economic efficiency of the operations. The importance of the problem of reservoir non-uniformity does not, of course, lessen its complexity or the difficulties of its solution. In fact, these are inherently such that the concept of a "general" solution is virtually meaningless. About all that can be reasonably hoped for is the analysis of specific and well-defined types of non-uniformity which may give some degree of approximation to actual reservoir conditions. Since variations in the nature of the reservoir which depend only on the position along the streamlines will not lead directly to major differential depletion development within the reservoir, the types of non-uniformity considered thus far have involved stratification assumptions. That is, the producing section has been replaced by a multi-layer "sandwich," each uniform areally, and differing from the others only in its basic physical constants as to thickness, porosity and permeability. The fluid motion in the composite system is thus approximated by a parallel superposition of the independent fluid movements in the individual strata. For the specific application to cycling operations the theory of the effect of permeability stratification has been developed for both discontinuous' and continuous types of permeability stratification. Among the latter, treatments have been given of systems in which permeability distributions are governed by exponential, linear' or probability3 functions. In all these studies complete dynamical equivalence was assumed between the injected dry gas and the displaced wet gas. The overall effective permeability of each stratum was therefore considered as constant and independent of the degree of invasion of the injected fluid. In the case of the displacement of oil by water, the assumption of dynamical equivalence between the water and oil will be strictly valid only by accident. Even if the oil viscosity should be the same as that of the water, the effective permeability to the oil in the presence of the connate water will in general be quite different from that of the water behind the water-oil interface flowing past the trapped residual oil. As a result the effective permeability for the stratum as a whole and rate of water invasion will change with time as the intrusion continues. The differential fluid motion in the individual strata will thus also vary with time. Qualitatively, it is easy to predict the resultant effects. If the permeability to viscosity ratio of the invading fluid exceeds that of the fluid displaced, the stratification and bypassing effect of the perme-