This paper presents the results of laboratory measurements of relative permeabilities to oil and gas on small core samples of reservoir rock by five methods, and describes the influences of such factors as boundary effect, hysteresis, and rate upon these measurements. The five methods used were the "Penn State," the "single core dynamic," the "gas drive," the "stationary liquid," and the "Hassler" techniques.In those methods in which the results are subject to error because of the boundary effect, the error may be minimized by the use of high rates of flow. In order to avoid complexities introduced by hysteresis, it is necessary to approach each saturation unidirectionally. Observed deviations of relative permeabilities with rate can be explained as a manifestation of the boundary effect, and disappear as the boundary effect vanishes.The results indicate that all five methods yield essentially the same relative permeabilities to gas. Of the four methods applicable to the determination of relative permeability to oil, three, the Penn State, single core dynamic, and gas drive, gave relative permeabilities to oil which were in close agreement. The Hassler method gave relative permeabilities to oil which were consistently lower than the results obtained by the other methods.

This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. In the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. Introduction The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradient s, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942. Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers. The influence of capillary forces in laboratory water floods has also been discussed by several authors. While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured- matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been published on the experimental techniques involved. Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. SPEJ P. 195ˆ

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. Introduction Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of particle differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow pressure and saturation. Thus, solution of the flow problems by difference equations involves solving problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems or problems of gas percolation. Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the net, level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. SPEJ P. 417

Published in Petroleum Transactions, AIME, Vol. 213, 1958, pages 96–102.Paper presented at 32nd Annual Fall Meeting of SPE in Dallas, Texas, Oct. 6–9, 1957. ABSTRACT The calculation of the behavior of an oil reservoir during a water flood has long been an important problem to reservoir engineers. Buckley and Leverett derived the differential equation which describes the displacement of oil from a linear porous medium by an immiscible fluid, but this equation could not be solved by the methods of classic mathematics. Consequently, in order to integrate the equation over the length of the reservoir, they neglected the effects of capillary pressure. In the present paper, a numerical method has been developed for determining the behavior of a linear water flood with the inclusion of capillary pressure. The differential equation which was derived for the case of incompressible fluids is second order and non-linear. This differential equation was approximated by an implicit form of difference equation which is second order correct in both time and distance. An electronic digital computer was used to perform the numerical solution of the difference equation. INTRODUCTION The problem of calculating the flow and distribution of fluids in an oil reservoir subjected to a water flood has long challenged the reservoir engineer. The ability to solve this problem would provide a valuable tool for the design and study of field waterflooding programs. One of the first contributions in this field was made by Buckley and Leverett, who developed a method of calculating waterflood performance in a linear reservoir. Their technique was limited by the practical necessity of excluding quantitative consideration of capillary pressure. It is the purpose of this paper to describe a method for calculating the behavior of a linear water flood with capillary pressure considered. This method, although limited to the linear case, should serve as a step toward the solution of the two- or three-dimensional waterflooding problem which would better describe actual reservoirs.

The shape and position of the gas-oil transition zone during downdip displacement of oil by gas has been calculated using flow equations which include the effects of gravity, relative permeability, capillary pressure and compressibility of the fluids. The calculations treat the problem in two space dimensions, and results are compared with data from a laboratory model tilted at 30 degrees and 60 degrees from the horizontal on displacements near and above the maximum rate at which gravity segregation prevents channeling of the gas along the top of the stratum. The good agreement between calculated and experimental results demonstrates the validity of the technique as well as that of the flow equations. Introduction Knowledge of the fluid distribution and movement in and oil reservoirs important in producing operations and estimation of reserves. The history of the oil industry has included steady progress in improving the accuracy of calculations which provide the required knowledge. The earliest method of calculating reservoir performance consisted of material-balance equations based on the assumption that all properties were uniform throughout a reservoir. For many reservoirs such a simple formulation is still the most useful. However, when large pressure and saturation gradients exist in a reservoir, the assumption of uniform values throughout may lead to significant error. To reduce these errors, Buckley and Leverett introduced a displacement equation which considers pressure and saturation gradients. Methods available at that time permitted solutions to the Buckley-Leverett equation in one space dimension; these solutions have been very useful in solving many problems related to the production of oil. However, the one-dimensional methods are not adequate for systems in which saturations vary in directions other than the direction of flow. An example of such a system is the case of gas displacing oil down a dipping stratum in which the gas-oil contact becomes significantly tilted. Of course, the Buckley-Leverett displacement method cannot predict the tilt of the gas-oil contact. Recent improvements of the one-dimensional Buckley-Leverett method achieve some success in predicting the tilt of the gas-oil contact at sufficiently low flow rates. However, at rates high enough that the viscous pressure gradient nearly equals or exceeds the gravity gradient, even these improved one-dimensional methods incorrectly predict the shape and velocity of the contact. Further progress in estimating such fluid movements in a reservoir appears to require consideration of the problem in more than one space dimension. The recent two-dimensional method of Douglas, Peaceman and Rachford appears adaptable to calculate changes with time of the saturation distribution in a vertical cross-section of a reservoir. The movement of saturation contours should represent the moving fluid contacts and include the effects of crossflow due to gravity, as well as variations in the rock and fluid properties. The nonlinear nature of the equations used in the method has prevented proof of the validity of the solutions. Douglas, Peaceman and Rachford made some comparisons with experiment but did not include cases in which gravity was important nor cases involving displacement by the nonwetting phase. Forthesereasons, atestof the two-dimensional method for a case in which these factors are included would be very desirable. The test selected was a comparison of calculated results with those from a carefully controlled laboratory experiment on a model with measured physical properties. The model selected was one in which gas displaced oil down a tilted, rectangular sand pack. The model can be thought of as representing a vertical cross-section taken parallel to the dip of a reservoir. The displacement thus simulates gas displacing oil downdip that might result from gas-cap expansion or gas injection. SPEJ P. 19^