Search Documents

By Harold L. Overton, Leonard B. Lipson

The first paper in this series1 outlined practical methods for applying the theory of steady-state flow of an ideal Bingham plastic liquid through a circular pipe and axially through a stationary concentric annulus to the enginekring analysis of friction loss problems. In the present paper the properties of another useful rheological model, the pseudoplastic generalized Newtonian liquid, are considered. The terminology applied to the rheological models is derived from the following classification of rheological models. The term "generalized Newtonian" is a class designation applied to models which do not have yield points but which exhibit a dependence on shear rate. The term "pseudoplastic" applies to a sub-group of models within this class which exhibit a decrease in "viscosity" with increasing shear rate. The term "power model" refers to a mathematical model proposed for describing the behavior of the pseudoplastic liquid. This model is of interest since certain of the salt saturated, oil emulsion, and "low solids" drilling fluids, inverted emulsion- and oil-base drilling fluids, aqueous gels, gelled crudes and blocking agents employed in hydraulic fracturing operations are of this type. The present treatment considers the equations describing steady-state Poiseuille flow through a circular pipe and a stationary concentric annul us, and Couette flow between concentric rotating cylinders for a particular pseudoplastic model. As in the previous paper' it is demonstrated how these equations can be applied to engineering calculations of friction losses. Graphical procedures for recovering the various flow curves and cxample calculations are included. A method for simplifying turbulent-flow correlations is also proposed. PROPERTIES OF PSEUDOPLASTIC GENERALIZED NEWTONIAN LIQUID In contrast to the Bingham plastic the generalized Newtonian liquid is yield point-free. However, in the flow region the relation between stress and rate of shear is not linear as it is with ordinary Newtonian fluids. If the stress increases with shear rate at a less than linear rate the system is called pseudoplastic. While there is no yield point in this system, the Bingham advocates' coined the term pseudoplastic to distinguish this kind of behavior from Bingham plastic behavior. For the pseudoplastic a plot of flow rate vs friction loss or rpm vs torque, from pipe flow or rotational viscometer measurements, respectively, generally has three qualities: (1) it can be clearly extrapolated toward the origin; (2) there is very little slope at the origin. This accounts for the high apparent viscosities at low mixing and pumping rates; and (3) the flow curve possesses a continuous curvature which gradually diminishes as the shear increases. Rheological measurements on a pseudoplastic liquid can be easily misinterpreted as indicating Bingham plastic behavior. A typical case in point is illustrated in Fig. 1. These data were obtained from measurements on an oil-emulsion mud using a commercially available rotational-type viscometer (Model 35 Fann V-G meter). Applying the "two-point"* method3 gives a plastic viscosity of 17 cp and a yield point of 24 lb/100 ft2; at the same time weighting material settles very easily in this system. This illustrates that it is not possible to distinguish between the pseudoplastic and the Bingham plastic from only two data points. It is possible to determine qualitatively which model is more appropriate by a combination of dynamic and static measurements. Experience indicates that the simplest criterion for the Bingham plastic liquid is that the calculated dynamic yield point should be of the same magnitude as the measured gel strengths. For the oil-emulsion drilling fluid illustrated in Fig. 1, a value of 4 lb/100 ft2 was obtained for gel strengths measured

By Kenneth F. Whinery, John M. Campbell

Basically, the Buckley-Leverett theory involves two systems which are similar in nature but are differentiated by time. These systems may be described by the fractional flow and frontal advance equations which essentially characterize the mechanics of oil movement while being expelled from the reservoir. The fractional flow equation originally developed by Leverettl may be expressed in the more usable form The development of this equation is based on Darcy's law describing fluid flow through porous media and applies to the flow at only one point (it is a point function). For simplification, the fractional flow equation is written in the above form because the capillary forces always increase the fractional flow of the displacing phase regardless of the direction of flow or the displacing phase. If, for simplicity, the effects due to gravity and capillary pressure differences are neglected, the fraction of the displacing fluid, f, at any point in the flowing stream is related to the properties of the system by 1 Thus, it is seen that in the absence of capillary and gravitational effects, fd for a given sand and fluids varies only with saturation and pressure. The magnitude of the viscosity ratio, —a has an effective range (range of about 30) in the system where gas is displacing oil and a much smaller range for water displacing oil. In order to make the fractional flow equation more versatile, it is necessary to connect the fractional flow at a given point and saturation with time. This problem was approached by Buckley and Leverett" who developed the frontal advance equation, Eq. 3 states that the rate of advance of a plane that has a fixed saturation, Sd, is proportional to the change in composition of the flow stream caused by a small change in the saturation of the displacing fluid. It is, essentially, a transformation of a material balance equation representing the net rate of accumulation of the displacing phase within a homogeneous sand block. This accumulation is proportional to the difference between the rate at which the displacing fluid enters the sand and that at which it leaves. Eq. 3 describes the velocity with which a plane of constant displacing phase saturation advances through a porous system. Buckley and Leverett," Babson,' Kern," Welge, and others have adequately discussed the basic mechanism and application of the fractional flow and frontal advance equations. APPLICATIONS OF THE FRONTAL ADVANCE THEORY TO PETROLEUM RECOVERY Two general applications of the Buckley-Leverett frontal advance theory involve the system in which the oil is being displaced by an expanding gas cap overlying the oil zone and that in which the oil is being displaced by water. Displacement by Gas in the Presence of an Immobile Water Saturation A system in which gas is the displacing phase may be thought of as having two forces effecting the displacement process. These forces are the gravitational force and that force exerted by the displacing gas. The gravitational effects control the displacing efficiency of the gas. The gravitational effect will be less at higher rates of flow, thereby reducing the effectiveness of the displacement of the oil by the gas. The more efficient displacements occur at flow rates which are less than the gravity free fall rate. Capillary forces can be neglected without materially changing the magnitude of the gas saturation. The Mile Six pool is used herein to illustrate the calculating procedures in evaluating gas drive-gravity drainage field perfomlance. These calculations represent the determination of the gas-oil contact when the distribution of the hydrocarbon pore volume is considered. Two methods

By B. H. Caudle, A. B. Dyes

In a recent publication' it was shown that wells with a free surface in a homogeneous gravity-drainage reservoir have a hyperbolic decline with index n '. This paper reports efforts to extend this theory to practical field cases. As a first step in making the extension, a large number of lease production curves from gravity-drainage fields in the stripper stage were fitted with the general hyperbolic decline equation, q/q, = 1/(1 + y) y and n in each case being determined by statistical methods. It was found that this equation is reliable for decline curve extrapolation for periods up to at least 25 years. The resrilts obtained were disturbing in the sense that values for the index, n, were often smaller than one, which fact leads to the prediction that a well would produce all infinite amount of oil in infinite time. A theoretical investigation based on the results of the previous paper provides a suitable explanation for these low values of n. It suggests that these are obtained if production occurs from two or more layers of diflerent permeabilitv and thicktress or two lnyers having different skin effects. It is predicted that as the more permeable layers are exhausted, the index of decline, n, will rise from a value smaller than 2, to n = 2. A practical method is given for fitting the general hyperbolic decline Eq. 6 to production decline data. When this method is used, the index of decline observed in gravity-drainage reservoirs in the stripper stage when wing only data over early years 01 decline usually agrees quite closely with that determined using data over 25 years or more. This indicates that the hyperbolic decline equation is reliable for predicting behavior of such reservoirs. A method is also suggested for taking into account any changes in the value of the index INTRODUCTION As a result of both experimental work with models and theoretical deductions, it was recently concluded' that wells having a free surface in a homogeneous reservoir producing by pure gravity drainage should have a hyperbolic decline with n = 2. The present report details efforts to extend this work on ideal homogeneous reservoirs to actual field cases. It should be remembered that both the previous work and the present extension are based on the assumption that the reservoir is at very low pressure, and hence, nearly all the potential available for driving fluid to the well is the gravity head. Attempts to apply conclusions of this study to reservoirs which have considerable pressures supplementing the gravity drive may lead to difficulty. WELLS WITH NO FREE SURFACE In a dipping gravity-drainage reservoir such as that shown in Fig. 1, some of the down-dip wells have no free surface. The previous study' showed that the up-dip wells, which have a free surface, should follow Eq. 1. The following approximate analysis, of a type which was justified experimentally for up-dip wells in the work just mentioned, leads to results for the down-dip wells.

By Stephen G. Dardaganian

Pore volume compresibilitieu measured in the laborutory on core .samples are reported for typical reservoir .sarrrl.rtorres at reservoir pressure. These cornpressibilities ore different for each reseriloir sample and cannot be correlated to porosity. The compressihilities are a1so functions of presslire. In recent years several authors have pointed out the importance of including pore volume compressibility of the reservoir rock in certain reservoir engineering calculations. Hall' showed that in calculating the hydrocarbon volume of an undersaturated reservoir from the production per unit change in reservoir pressure, the neglect of rock compressibility could, in the extreme case of low porosity rock, lead to results in error by a factor of two. Hawkinsl' and Hobson and Mrosovsky showed that both rock and interstitial water compressibilities must be included to achieve satisfactory accuracy in this calculation. Although rock compressibility is an important reservoir property in these calculations, very few data are available. Hall' made compressibility measurements on a few samples in a small pressure range. No other direct pore volume compressibility measurements have been reported. Geertsma' discussed the elasticity of reservoir rock but did not present any new data. In this paper, rock compressibilities suitable for reservoir calculations are reported for a variety of sandstones typical of U. S. reservoirs. The measurements were made through a range of pressures which is also typical of reservoirs. EXPERIMENTAL PROCEDURE Plugs 1 in. in diameter and 2- to 3-in. long were diamond drilled, 14 from cores from seven oil-bearing or potentially oil-bearing sandstone, one from a quarried sandstone, and another from a shallow, oil-free, pure orthoquartzite. The plugs were extracted in toluene, vacuum dried, and then weighed. The plugs were then saturated with kerosene and reweighed. This gave pore volume. They were then again extracted in toluene. vacuum dried, and jacketed in 0.004-in. thick copper foil. A short length of -in high pressure tubing led from the jacket. The jacketed sample was then saturated with kerosene and placed in the hydraulic pressure cell as shown in Fig. 1. Movement of the mercury slug showed the change in pore volume of the sample under a given pressure change. Temperature was about 72" F for all measurements. Before taking any measurements, the external pressure on the sample was raised to 10,000 psig and kept there for several hours. This formed the thin copper jacket around the plug. 7-he pressure was then removed and the sample left at least overnight, during which time it recovered its original pore volume. An external pressure of 12,000 psig was then applied and held constant. The change in pore volume was noted as the internal pressure was raised in 1,000 psig increments to 10,000 psig. The internal pressure was then reduced to zero, again noting pore volume changes,

By William Hurst

Muskat's depletion performance equation is here derived considering the expansion behavior of the reservoir hydrocarbon system and a simple fractional-flow equation. This nietkod of derivation leads logically to the two extensions that follow. The first of these is concerned with gravity segregation in a depletion-drive reservoir. The second is concerned with including an empirically determined tern? for water influx in the performance equation. The more general equation for gravity segregation when there is a primary gas cap and empirically-determined water influx is stated for completeness. These equations have been found useful in reservoir performance calculations in Eastern Venezuela. A discussion on the methods of solving these equations follows, and considers firstly the effect of taking finite intervals in the numerical integration, and sec-ondly, methods of incorporating the time functions involved in segregation in with the expansion behavior. The paper concludes with a brief general discussion on further extensions to the depletion performance equation. INTRODUCTION The two fundamental sources of energy by which oil is produced from a reservoir result from pressure depletion inside the boundaries of the reservoir and fluid encroachment across the boundaries of the reservoir. The wells in either case form low pressure outlets through which oil and gas may be produced by the expansive force of the reservoir fluids and associated encroaching fluids. When the reservoir pressure is higher than the bubble-point pressure of the oil, so that there is no free gas in the reservoir, these expansive forces are the only ones available for the production of oil. However, when the reservoir pressure is less than the bubble-point pressure of the oil, free gas is vaporized as the pressure falls. With both oil and free-gas phases present, the additional forces of gravity and capillarity may operate on the gas-oil system, as they have previously operated on the oil-interstitial water system. Gravity tends to segregate the free gas from the oil due to their density difference. Capillarity opposes and eventually balances gravity as the more extreme free gas and oil saturations are reached, preventing the independent move- ment of free gas until it is above a certain saturation, and the independent movement of oil when it is below a certain saturation. The type of depletion performance equation chosen for predicting the future performance of a reservoir depends on the amount of past history available. When the reservoir is somewhere past the halfway mark in depletion, some form of decline curve is often used. With less past history, material balance equations which incorporate empirical factors based on the past performance are often used. When, however, the amount of past history is small, the Muskat depletion performance equation will usually be used. The distinguishing feature of this type of equation is that empirical factors based on the over-all or macroscopic reservoir behavior are almost or entirely absent. Each parameter affecting the reservoir performance is ascribed an independent set of values based on measurements made on laboratory samples; that is, incorporating microscopic empirical factors. In establishing Muskat-type depletion performance equations, it is necessary to consider the reservoir as consisting of a number of associated blocks, in each of which the saturations and pressures may be considered uniform, and in each of which all substances have uniform pressure-volume characteristics. Thus, a primary gas cap can usually be considered as one block and an aquifer as another. Gravity segregation may be negligible for practical purposes when the rock and oil properties are adverse and/or the dip or thickness of the reservoir is too small. In this case the whole oil leg may be considered as one block, except in very large reservoirs. In very large reservoirs the fluid and rock properties may vary enough, particularly in the dip direction, for it to be necessary to divide the oil leg into a number of blocks, in each of which the relevant quantities may be considered uniform. When gravity segregation of the oil and free gas is not negligible, it is necessary to consider the space occupied by the initial oil leg as divided into two blocks, a secondary gas cap and an effective oil leg, in each of which saturations may be considered to be uniform. The total volume of these two blocks is thus constant, but the secondary gas cap grows continuously at the expense of the oil leg. Muskat1 derived depletion performance equations for the basic case of an oil reservoir with closed boundaries and without segregation, and for the case of an oil

By T. G. Richardson, F. M. Perkins

In the recent paper of Richardson and Perkins entitled "A Laboratory Investigation of the Effect of Rate on Recovery of Oil by Water Flooding,"' the authors found very little appzrent effect on oil recovery of having free gas saturation present prior to water flooding in unconsolidated sands. We agree that the effect in this particular system was small but do not feel that the data presented should be regarded as disproving work of others which has shown an appreciable effect. The authors do not make such a claim, but neither do they advance an explanation of the difference. In our experience the following factors must be carefully considered when studying the relation between residual oil saturation after water flooding and initial free gas saturation: (1) the magnitude of the inherent residual oil saturation, i.e., the residual oil saturation with no free gas; (2) the increase in oil-water viscosity ratio with pressure lowering if the gas saturation is obtained by solution gas drive; and (3) the oil shrinkage which occurs with pressure lowering and gas evolution. The system studied by Richardson and Perkins gave residual oil saturations as low as 20 per cent without free gas present, and one should perhaps not expect much further lowering. As an illustration of the magnitude of the effect when the residual oil saturation is inherently higher we submit the accompanying figure. These data were obtained on an unconsolidated sand under elevated pressure in a system similar to the one described by Richardson and Perkins. Gas saturation was established by gas drive before water flooding. The viscosity ratio was 600:1 under the conditions of the displacement. Here the gas was extremely effective in reducing residual oil saturation. This fact was also indicated by Kyte, et al.2 Their plots of residual oil saturation vs initial gas saturation show trends of greater free gas effect with higher oil viscosities. The pressure-dependent effects were not discussed by Richardson and Perkins. They should be separated from the free gas effect if the objective of the study is a correlation of residual oil saturation vs initial free gas saturation. The increase in oil viscosity when gas is released from solution tends to increase residual oil saturation. This effect is rather small in the viscosity range of kerosene. A Buckley-Leverett calculation, assuming an increase in viscosity ratio from 1.2 to 1.8, indicates an increase in residual oil saturation of about 1 per cent of pore volume. Thus, one would expect this system, if operated at constant viscosity ratio, to show a slightly greater free gas effect. As Dyes- as shown, oil shrinkage can offset the beneficial effect of free gas on waterflood oil recovery, if the gas saturation is established by depletion. Since

Jan 1, 1958

By Charles R. Stewart, William W. Owens

Reservoir performance predictions based on laboratory core test data assume that fluid flow is laminar for the laboratory test. A study has been made to determine the validity of this assumption for laboratory tests on various types of porosity found in producing limestone formation. Data are presented which show that turbulence and slippage can occur during laboratory tests on hetero-geneous-type porosity limestones, thus causing serious errors in measured single-phase permeabilities and two-phase relative permeability characteristics. In single-phase flow tests it is possible to eliminate turbulence and correct for slippage or to eliminate both factors by controlling test conditions. It is not always possible to control test conditions and thereby eliminate turbulence and slippage in two-phase .flow tests. A correction method is presented which can be used to calculate the true two-phase laminar flow relative permeability characteristic even though furbulence and slippage exist. .INTRODUCTION It is customary to make use of Darcy's law and modifications of this law, together with laboratory data on formation core samples to predict the performance of producing reservoirs. Such predictions are based on an assumption that fluid flow is in the laminar or streamline region for the laboratory test. It was the purpose of this inves- tigation to determine the extent to which turbulent flow may occur in laboratory fluid flow tests on hetero-geneous porosity limestones. Considering that turbulent flow conditions might exist in some laboratory fluid flow tests, additional emphasis was placed on the development of a method to correct for turbulence when laminar flow conditions could not be attained. FLUID FLOW CONCEPTS FOR POROUS MEDIA The Influence of Pore Geometry on Fluid Flow One of the more important factors influencing fluid flow in porous media is the geometry of the pore space which includes such characteristics of the pores as size, shape, distribution, roughness, uniformity, etc. In general, oil- and gas-producing formations can be divided into two broad types on the basis of pore geometry. One has been called sandstone-type porosity media, which is characterized by a small range in pore size, uniformity in shape of the pores, smooth pore surfaces and a regular and uniform distribution of pores. The other type has been called heterogeneous porosity media and is usually limited to the dolomites and limestones. This type is characterized by a wide variation in the size, shape, and distribution of the pores and rough, irregular pore surfaces. It is therefore apparent that conditions are much more favorable for turbulent flow* in heterogeneous-type porosity media than in sandstone-type porosity media. Interrelationship Between Turbulence and Gas Slippage In studying the problem of turbulent flow in laboratory tests on porous media, it is necessary to be aware of the interrelationship between slippage and turbulence for gas flow. As a result of slippage or the Klinken-berg effecta, apparent perrneabilities to gas are greater than the true value because there is no stationary layer of gas in contact with the walls of the flow channels. Gas slippage decreases as the mean free path of the gas molecules decreases. Since the mean free path of any gas decreases with increasing density, increases in static pressure result in lower apparent gas permeabilities. However, a reduction in gas permeability can also be due to turbulence. Therefore, in studying only turbulent flow in porous media, it is necessary to hold gas density, and slippage, constant or to reduce slippage to a negligible value by operating at high static pressures. Presentation of Laminar and Turbulent Flow Data A graphical relationship between permeability and a pseudo-Reynolds number, N,, will be used to show the two types of fluid flow, i.e., laminar and turbulent. The usual graphical method for such a description has been the use of friction factor-Reynolds number charts4. On such a logarithmic diagram, the laminar region appears as a straight line having a slope of 45 degrees. As the friction factor decreases and the Reynolds number increases, the turbulent region is reached and appears as a deviation from the 45-degree slope line. However, in petroleum engineering literature resistance of por-

By J. C. Melrose, W. R. Foster, J. G. Savins, E. R. Parish

Many differences can be imagined between gas-oil flow in which the gas is supplied at the face of the core and gas-oil flow in which the flowing gas was originally dissolved in the oil. If capillary pressure characteristics and flow requirements control gas saturation distribution, the gas would be expected to be located at preferred sites within the porous medium as determined by pore sizes. On the other hand, during solution gas drive the gas first appears as bubbles through a nuclea-tion process. Nothing in self-nucleation theory specifies at which sites the first bubbles should be formed. In all probability they will be randomly distributed throughout the porous medium. Furthermore, it is not at all certain that even at low rates of production the gas will redistribute itself after nucleation to the channels normally occupied by gas in simple gas flow. Stewart, et at, have shown that at least for some limestone samples, oil recoveries could not be predicted for all rates of production using any one set of relative gas and oil permeabilities. An important factor in controlling recoveries during solution gas drive was the rate of bubble formation, higher rates giving higher recoveries. Stewart, et al, attributed the increase in recovery to a better distribution of the gas phase in heterogeneous limestone samples than is obtained by simple external gas drive. Differences in recovery from these causes were not reported for sandstone cores. In the experiments to be reported here, oil recovery, pressure and producing GOR history were measured during solution gas drive for a 5-ft sandstone core. The results were compared with predictions from the Muskat method for computing solution gas-drive behavior using external gas-drive relative permeability. The effects of changing the rate of production and oil viscosity were studied. At high laboratory rates of average pressure decline, two observations were made which would not have been predicted by Muskat's depletion theory: (1) oil recovery increased with increasing rate of production for a given viscosity oil, and (2) oil recovery increased with increasing oil viscosity for a given high rate of production. Both of these observations are explained as consequences of diffusion control of gas saturations superimposed on the normal gas-oil flow requirements, Fur- thermore, discontinuous gas phase flow appears to be significant during solution gas drive. The laboratory tests were performed at rates of average pressure decline many times greater than the maximum possible rate of average pressure decline in an actual oil field. It is, therefore, not possible to draw any direct conclusions regarding the effect of rate on recovery for the solution gas-drive mechanism under actual field conditions. However, at the lower laboratory rates, recoveries were nearly independent of rate and could be predicted by the Muskat method, using external gas-drive relative permeability data. These results suggest that at normal oilfield rates the effect of rate on recovery for the solution gas-drive mechanism is negligible. EXPERIMENTAL PROCEDURES The core material for the pressure depletion studies was Bandera sandstone from an outcrop in the Mid-Continent. This sandstone was selected because of its low permeability (about 10 md), which would permit the development of substantial pressure gradients in the corn at moderate flow rates. The core was 5-ft long and 2-in. in diameter. Its properties are listed in Table 1. Relative gas-to-oil permeability ratios were measured by an external gas-drive method2. The results are shown for a short 2-in. core and for the 5-ft core in Fig. 8. Oils used in the pressure depletion experiments were kerosene and a highly refined white oil (standard white oil No. 3) with gas-free viscosities of 1.8 and 25 cp, respectively. The gas was a naturally occurring methane from Gough field, Inglewood, Calif. The oil viscosities, gas solubilities and formation volume factors are plotted as functions of pressure at 75°F in Figs. 1 and 2. Methane viscosities and compressibilities were obtained from the literature"'. A core mounting was required which could withstand up to 2,500 psi internal pressure. This was obtained by first encasing the core completely in a plastic resin (Scotch Cast, manufactured by Minnesota Mining & Manufacturing Co.) The plastic covered core was then inserted into a steel pipe equipped with screw caps so that the plastic coating could be pressured from

By J. G. Richardson, J. W. Graham

The purpose of this paper is to present the results of theoretical and experimental studies of water imbibition. The imbibition processes are involved in recovery of oil from stratified and fractured-matrix formations in natural water drives and water flooding. An understanding of the role of inhibition in implementing the recovery of oil from such formations is deemed essential to proper control of these reservoirs to achieve maximum recovery. The theoretical studies involved development of the differential equations which describe the spontaneous imbibition of water by an oil-saturated rock. The dependence of the rate of water intake by the rock on the permeability, interfacial tension, contact angles, fluid viscosities and fluid saturatiorls is discussed. A few experiments were performed using core samples to determine the effects of core length and presence of a free gas suturation. The role of water imbibition in recovery of oil from a fractured-matrix reservoir by water flooding was investigated by use of a laboratory model. This model was scaled to represent one element of a frac-tured-matrix formation. Water floods were made at various rates with several fracture widths. Interpretations were made of the behavior expected in a system containing many matrix blocks. The presence of a free gas sntu.ration was found to reduce the rate of water imbibition. In the reservoir prototype of the fractured-matrix model, water imbibition rather than direct displacement by water was the dominant mechanism in the recovery of oil at low rates. INTRODUCTION Imbibition may be defined as the spontaneous taking up of a liquid by a porous solid. The spontaneous process of imbibition occurs when the fuid-filled solid is immersed or brought in contact with another fluid which preferentially wets the solid. In the process of wetting and flowing into the solid, the imbibing fluid displaces the non-wetting resident fluid. Common examples of this phenomenon are dry bricks soaking up water and expelling air, a blotter soaking up ink and expelling air and reservoir rock soaking up water and expelling oil. As increasingly better lithological descriptions have been made of the characteristics of petroleum-bearing formations, it has become obvious that imbibition phenomena which were once considered laboratory curiosities are of practical importance. For instance, in reservoirs composed of water-wet sand strata of different permeability in intimate contact, the tendency of water to channel through the more permeable stratum is offset by the tendency for water to imbibe into the tight sand and expel oil into the coarse sand. Also, in fractured-matrix formations the tendency of water to channel through the fractures is offset by water-wet matrix blocks. As some imbibition of the water into the of the largest fields in the world are fractured-matrix reservoirs, it has become increasingly important to understand all the factors involved in the imbibition process. Examples of fractured-matrix reservoirs are the Spraberry field in West Texas which produces from a fractured sandstone', the giant Kirkuk field in Iran', the Dukhan field in Qatar, Persian Gulf2, and the Masjid-I-Sula-main and the Haft-Kel fields in Southwestern Iran, which produce from fissured limestone3. Research into recovery of oil from fractured-matrix formations was stimulated by the rapid decline of oil productivity of wells in the Spraberry formation. One result of this research was the water imbibition process developed by the Atlantic Refining Co.4 Another idea was that much of the Spraberry oil could be recovered by conventional water-flooding procedures5. Subsequently, pilot floods were conducted in this field to test the feasibility of these ideas. It was felt that an understanding of the role played by imbibition processes in displacement of oil from a fractured-matrix reservoir could not be obtained from field data alone because of the many complicating factors and uncertainties involved. Therefore, theoretical and laboratory studies were undertaken to provide this understanding. Study of the equations which describe the linear, countercurrent imbibition process provided an insight into the role of various factors in the process, such as the permeability of rock and inter-facial tension. In addition to the theoretical studies, imbibition experiments were conducted with core samples to determine the effect on the rate of imbibition of such variables as core length and free gas saturation. The principal experimental studies were conducted by water flooding a scaled model of an clement of a frac-tu red-matrix reservoir to evaluate

By O. C. Holbrook, George G. Bernard

A new theoretical treatment has been obtained for the behavior of pattern waterflood injection wells when closed in. Two cases are treated: Case I where oil and water are assumed to have the same properties, and Case 2 where they arc different. In applying the method, one plots log (p — p,) vs closed-in time, where p is well-bore pressure at any tims and p, is static pressure. The value of p. is determined by trial and error as that value which makes the plot linear at large time. A value for the permeability-thickness product can be determined from the intercept of this linear part, and a value of the skin factor from the injection pressure at time of closing in. Application of the method to data from water floods at three fields seems to give reasonable results. For the case of unit mobility ratio, it is proved that this new method should give the same value for permeability-thickness product as the conventional pressure build-up method. In addition, the new method gives correct values for static pressure, whereas the conventional method does not, often indicating negative static pressures. The new method may be used in cases where the surface pressure persists after closing in as well as in cases where it does not. INTRODUCTION It is of considerable interest and importance to be able to determine the characteristics of the reservoir in an area surrounding a water injection well. Thus, if we can determine early in the life of an injection well that there is a considerable "skin effect", remedial measures can be started before a full-scale pattern flood begins. Similarly, if it can be shown that a gradual buildup of skin effect is occurring with time, measures to free the water of plugging material can be taken. Determination of static pressure in the water-injection well may show that the water is entering a thief zone and not the desired reservoir. Finally, determination of the permeability of the sand around the injection well will allow estimation of the future relation between injection pressure and rate. It should be possible to determine average reservoir permeability, skin effect and static pressure from pressure fall-off data. However, at the time we began work on this subject, it was thought that no adequate theory on which to base such determinations' was available. According to the conventional method which considers the reservoir to be filled with one fluid of small compressibility (see Van Everdingen, Joers2, and Nowak2), shut-in pressure is plotted vs log where is injection time, and At is closed-in time. The physical significance of injection time, may well be questioned in this case, since in a reservoir completely filled with a single fluid (as required by this theory) and with input and output rates equal, the pressure behavior after an initial transient is independent of t,. Attempts by our Tulsa area to use this theory led to negative values of static pressure in most cases. Because of these limitations of the method discussed above, it was decided to attempt to develop a new theory of pressure behavior in water injection wells, one which would apply when there is a gas saturation, as is so often the case in water floods. In the following treatment the assumptions and basic equations are given first, then the method of application of the equations. A complete example is given to clarify details of application. All difficult mathematics has been placed in the appendices so that the reader can follow the text without difficulty. However, if he wishes only to apply the results without knowing the basis for them, he can learn how to do this from reading only the sections entitled "Plotting of Experimental Results" and "Example." ASSUMPTIONS AND BASIC EQUATIONS Statement of Problem It will be assumed that a horizontal layer of constant thickness contains in its pore system a mixture of oil, gas and water. While water is being injected into this pore- system through a well at constant rate, an oil bank is built up, gas being expelled from the space taken by the oil as shown in Fig. 1. The saturations within each

By C. S. Matthews, H. C. Lefkovits

The performance of depletion-type reservoirs at the stage in which the gas is at such low pressure that gravity is essentially the sole driving force has been investigated both theoretically and by use of scaled models. A simple theory is developed which completely explains the scaled model experiments. It also agrees with some field data. When the range of production rates to be considered is less than 25, the theory predicts that the production rate decline curve for homogeneous reservoirs in which there is a free surface will be of the hyperbolic type with n = 2: The theory gives a physical significance to the parameter y. INTRODUCTION This investigation was undertaken to determine a method for prediction of performance of reservoirs where gravity is the sole driving agent. Depletion-type reservoirs fall into this category after the pressure has fallen to some low value. Some examples of this type are found at Healdton ('Oklahoma), Oklahoma City (Oklahoma), Midway-Sunset (California), and East Coalinga (California). Generally the production rate per unit thickness is quite low in such reservoirs, with the conse- quence that this stage is sometimes called the stripper stage. However, per well production rates may be of the order of 50 B/D because of the considerable sand thickness open in some of the reservoirs. In the past, predictions of behavior of such reservoirs have been based on empirical extrapolation of the rate of production to the economic limit. Difficulties arise when new wells are drilled in the field or when production is artificially curtailed, since in such cases the decline rate of individual wells or leases changes and no clear trend can be detected. It was hoped that the present study would lead to improved methods for evaluation of such reservoirs. The study was initiated by building a scaled model for a particular field and by running a number of experiments on this model. Later a simple theory was developed which explained the laboratory results completely. The theory led, as shown below, to a hyperbolic-type decline equation which had been used previously' in a completely empirical manner. THEORY The General Case The problem to be treated here is the following: We are given a homogeneous cylindrical porous reservoir with exterior radius r, wellbore radius r., and height h,, containing, from the bottom up to a height hit oil, some gas, and interstitial water (see Fig. 1 A). The portion of the reservoir between K, and h, is filled with gas at low pressure, with residual oil, and with interstitial water. If at a time t = 0, the level of the oil in the wellbore is lowered to a height h,, what is the production rate from the reservoir at a time to The configuration within the reservoir at some instant

Jan 1, 1957

By Denton R. Wieland, Harvey T. Kennedy

The frequency of bubble formation is measured for oils from the East Texas field and the Slaughter field in cores from these fields. East Texas oil was also tested in a Slaughter core. The data checked previous results in that a definite supersatrtration may exist without the formation of bubb1es in periods of observation totaling 158 hours. Measured frequencies varied from zero at supersaturations below 14 psi to 3.1 bubbles/sec/ cu ft of rock at 40 psi supersaturation INTRODUCTION The formation of bubbles is essential in a solution gas-drive reservoir in the displacement of oil. In formations where the basic pore structure consists of fine pores or fractures connected to larger pores or crevices, bubble formation in the fine pores may be the most effective means by which the oil may be displaced, even if secondary methods are applied. The most efficient oil displacement would occur where at least one bubble would form in every pore. As the pressure decreased, the bubbles would expand, due to the gas diffusing into the bubbles, and displace the oil toward the lower pressure area around the wellbore. Kennedy and Olson' have shown that the total number of bubbles formed in a reservoir is primarily dependent upon the rate of pressure drop, the rate of diffusion of gas through oil, and surface area of the rock. The surface tension of the liquid and other characteristics of the rock may also exert an important effect, which would give a different frequency to different rocks and oils and various combinations of different rocks and oils. The purpose of this paper is to present frequency measurements on two rocks and two oils. A hydrocarbon liquid is at equilibrium with gas when it contains the maximum volume of gas in solution that it can have at the prevailing temperature and pressure. If the temperature is maintained constant and the pressure is reduced below the equilibrium pressure, the liquid is supersaturated to the extent of the difference in the two pressures until the gas comes out of solution as a free phase. As shown by previous work, and confirmed in the present investigation, the number of bubbles formed per second increases with an increase in supersaturation. Kennedy and Olson showed that a tenfold increase in the rate of pressure decline resulted in a tenfold increase in the number of bubbles formed per cubic foot of reservoir rock, e.g., for decline rates of 0.1, 1, and 10 psi/day, the number of bubbles formed per cubic foot would be 40, 400, and 4,000, respectively. For another hypothetical system, Stewart, Hunt, and Berry' estimated that for pressure declines of 0.1, 1, and 10 psi/day, the number of bubbles formed would be approximately 10, 100, and 1,000 cu ft of reservoir rock. Their value of bubbles formed is lower than given by Kennedy and Olson; however, the order of magnitude agreed satisfactorily. The rate of diffusion will also affect the number of bubbles formed. If the diffusion coefficient is high, the gas in solution in the oil surrounding a single bubble will tend to diffuse into the bubble rather than forming other bubbles. On the other hand, the formation of more bubbles will tend to occur if the diffusion coefficient is very low. In his thesis dated Aug., 1953, Wood" reports on the frequency of bubble formation in a system consisting of oil and rock from the Rangely field. In the present work, the apparatus of Wood was redesigned to allow greater precision, especially at low supersaturations, where bubble growth is slower, and detection is more difficult.

Jan 1, 1958

By J. R. Kyte, L. A. Rapoport, R. J. Stanclift, S. C. Stephan

Experimental studies covering a wide range of core materials and fluid properties have been conducted to determine the mechanism of oil displacement by water in a partially gas-saturated porous medium. In all instances, the presence of a gas phase was found to have a beneficial effect in reducing residual oil saturations. The practical significance of this benefit is discussed, and a simplified procedure is outlined for evaluating the effects of free gas on water flooding by means of short core tests. INTRODUCTION In addition to oil and water, reservoirs subjected to water flooding frequently contain, also, a gas phase. Common engineering procedures account for the presence of this "free gas" only from the viewpoint of volumetric balance, implying that the only role of the gas consists of providing "fill up" space. It is usually visualized that during the initial stages of the water invasion, the oil, moving ahead of the water, displaces part of the gas and that subsequently the remaining portion of the gas phase is totally compressed and dissolved in the advancing oil bank. Thus, consideration of a two-phase, water-oil flooding mechanism supplies an adequate basis for predictions of oil recovery if the pressure build-up caused by the flood is sufficiently great, so as to reduce the free gas saturation to a negligible value in any portion of the reservoir by the time that portion is reached by the advancing flood water. In many in- stances, however, waterflooding operations are carried out when the reservoir pressure is still relatively high SO that the pressure build-up associated with the water flood may not result in complete dissipation of the gas phase. Under such circumstances it is necessary to ascertain whether or not the presence of free gas has an effect on waterflood behavior and, if so, to account for any such effect in the evaluation of field operations. The difficulties encountered in attempting to arrive at a comprehensive picture of the displacement mechanism of oil by water in the presence of free gas stem from the fact that most of the published data on this subject were obtained under conditions that were either not clearly representative of those existing in a reservoir or not sufficiently controlled to permit evaluation of the effects specifically caused by the gas phase. Some of the most precise data found in the literature were obtained in tests using a "static" gas phase created prior to each water flood by circulating oil through partly gas-saturated cores until a "trapped," residual gas saturation was established throughout the system.' Other investigations were conducted under more representative conditions where the gas phase was free to move during the water flood: but in most cases the results of these tests were obscured by auxiliary phenomena which arise from the compressibility and solubility of the gas and preclude quantitative conclusions.2,3,4 Prior to water flooding, a reservoir usually contains oil and gas, distributed in a manner in which both are free to move. Thus, the principal objective of the present studies was to determine the type of displacement mechanism occurring when water invades a porous medium containing oil and a "mobile" gas phase. Further-

Jan 1, 1957

By H. J. Ramey, I. S. Bousaid

Experimental results on the oxidation reaction kinetics in the forward combustion oil recovery process are presented. A total of 48 runs were made wherein a stationary thin layer of coked, unconsoli-dated sand was burned isothermally in a combustion cell. Individual runs were made at various temperature levels to permit determination of the effect of temperature upon the reaction. An expression was obtained for the burning rate of carbon as a function of carbon concentration, combustion temperature and oxygen partial pressure. The carbon burning rate for two types of crude oil indicated a first order reaction with respect to both carbon concentration and oxygen partial pressure. The effect of combustion temperature on the reaction rate constant matched the Arrhenius equation. The activation energy was similar for the two crude oils examined. The activation energy decreased for a porous media containing clay. The rate of oxidation of crude oil at reservoir temperature was found to be significant. Other significant findings included information on hydrogen-carbon content of fuel residues, fuel reactivity and the products of combustion. INTRODUCTION The production of crude oil by underground combustion has been studied in the laboratory by many investigators.1-5 Results of laboratory and field experiments have been reported in the literature describing the forward combustion process. But as yet, no qualitative or quantitative study of the kinetics of fuel combustion involved in this process has been reported. The fuel concentration and the rate at which fuel is burned at the front are important factors governing the air requirement in a forward combustion operation. Although the fuel is essentially unrecoverable crude, the air required to burn the fuel is an important economic factor in this process. Because fuel is burned, the heat transport associated with forward combustion is a key and unique feature of this new oil recovery method. Many investiga-tors5-9 have presented information on the heat transmission and fluid mechanics involved in forward combustion. Berry and Parrish10 demonstrated the utility of considering reaction kinetics in reverse burning. From differential thermal analysis, Tadema11 presented a qualitative discussion of the nature of reactions between oil and oxygen in combustion oil recovery. Although little quantitative work has been done on the reaction kinetics involved in forward combustion oil recovery, an extensive literature does exist on combustion of carbons and oils, and carbonaceous residues from cracking catalyst pellets. Dart, et al., 12 studied the combustion rate for oxidation of carbonaceous residues on clay catalyst pellets, and found the reaction to be second-order with respect to carbon concentration, and first-order with respect to oxygen partial pressure for carbon concentrations less than 2 weight percent of the catalyst weight. The reaction appeared to be first-order with respect to carbon concentration for concentrations greater than 2 percent. Metcalfe13 noted that other workers had found that aging of the fuel during the combustion process was responsible for changing coke properties, and it accounted for the apparent second-order carbon concentration effect found by Dart, et al. It appears that burning of residues from cracking pellets is first-order with respect to both carbon concentration and oxygen partial pressure. Dart, et al., also observed that hydrogen in the hydrocarbon residue appeared to react faster than the carbon. Lewis, et al.,14 studied oxidation of charcoal, coke and graphite in a fluidized bed. Gas velocities were high enough to partially lift and circulate the carbon particles. Their results indicated first-order reaction dependency with respect to both carbon concentration and oxygen partial pressure. One

Jan 1, 1969

By R. O. Leach, O. W. Wagner

Because of unfavorable wetting conditions much residual oil is left when a porous material is Pushed by water. Methods suggested to change reservoir wetting to improve oil displncernrnt efficiency are generally expensitlr. The present 1aborator.y study was undertaken to gain an under.standinx of the factors which determine reservoir wettability, arid to find out if oil displacement efficiency might be improved by a wettahility change accomplished at low cost in on oil reservoir. Contact angle measurements were made on mineral surfaces using sevc.r~zl sets of reservoir oil and water samp1es. Results of the contact angle studies suggest that reservoir wetta-hility may he primarily determined by natural surface-active substances present in the reservoir fluids. The effect of changing sa1inity and pH of the water phase was studied. The re.suits suggest that gross changes in preferential wettability might be acc~o~npli.shed by injection of water containing simple chernicnls to alter pH or salinity in the reservoir. Such treatment could he much less expensive than injection of commercial surface-active agents. Waterflood tests have also been made using synthetic cores and oil and water having wening characteristics similar to those of reservoir fluids. Cores initially oil-wet were flooded in such a way that they were made prefermtial1y water-wct by the advancing flood water. This reversal in preferential wettability achieved greater oil displacement efficiency than when either oil-wet or water-wet conditions were maintained throughout the flood. For the systems studied, the higher the oil viscosity the greater the percentage improvement obtained over conventional waterflood recovery. This suggests that a flooding process making use of wettability-reversal may extend the oil viscosity range over which water flooding is attractive. Because a precise adjustment of reservoir wettability does not seem to be required, and because altering the pH or salinity in some reservoirs may be inexpensive, it appears that a waterflooding process employing wet-[ability-reversal could find .succesful field application. I NTRODUCTION The efficiency with which water will displace oil from a porous material is related to the nature of the capillary forces present. These in turn are controlled by the preferential wetting of the solid by the two fluids. Because of unfavorable wetting conditions, 30 per cent or more of the original oil in place may remain unrecovered in that portion of a reservoir flushed by water. This paper is concerned with the possibility of improving waterflood oil displacement efficiency by alterations in the wettability of the porous material. A laboratory study was made to gain a better understanding of the factors which control reservoir wettability, and to determine if the oil displacement efficiency could be improved by some inexpensive means of manipulating wettability of the porous medium. Contact angle measurements were made with several natural and synthetic oil, water and solid systems (1) to obtain a better understanding of how to duplicate reservoir wettability in the laboratory, and (2) to discover possible means for changing preferential wettability of natural reservoir systems. Flooding tests were also made in synthetic systems to determine if oil displacement efficiency could be improved by those wettability manipulations suggested by the contact angle measurements. Based on these studies a possible method for improving waterflood oil displacement efficiency is presented. This method involves causing an originally oil-wet porous material to become preferentially water-wet during the course of a water flood. The purpose of this paper is to present results of the laboratory studies. THEORY Rock surfaces in some oil reservoirs are believed to be covered with a firmly attached bituminous or other organic coating. Such surfaces would be preferentially oil-wet in the presence of both oil and water, regardless of composition of reservoir fluids. Other reservoirs are believed to contain rock surfaces not permanently coated with such materials, and which would be preferentially wet by water in the presence of water and oil free from surface-active substances. However, when the reservoir fluids do contain certain natural surface-active materials in sufficient quantity, rock surfaces acquire a degree of preferential oil wettability caused by adsorption of these natural surface-active materials on the solid. The equilibrium amount of these materials adsorbed per unit surface area is believed to depend upon their concentration in the bulk liquid phases.

By L. W. Holm

This study shows that in the presence of foam, gas and liquid flow separately through porous media representative of reservoir rock. These results were obtained by using tracer techniques to measure the flow of the gas and liquid comprising the foam. Foam does not flow through the porous medium as a body even when the liquid and gas are combined outside the system and injected as foam. Instead, the liquid and gas forming the foam separate as the foam films break and then re-form in the porous system. Liquid moves through the porous medium via the film network of the bubbles and gas moves progressively through the system by breaking and re-forming bubbles throughout the length of the flow path. The flow rates of the gas and liquid are a function of the number and strength of the films in the porous medium. There is no free flow of gas; i.e., no continuous gas phase. On the basis of these results, foam can be expected to improve a water flood or gas drive by decreasing the permeability of the reseruoir rock to a displacing liquid or gas. This improves the mobility ratio and thus the conformance of the flood. INTRODUCTION Foam is formed when gas and a solution of a surface active agent are injected into a porous medium either simultaneously or intermittently. During the past few years, several papers have been published on the subject of foam flow in Dorous media. Foam has been used successfully in the removal of capillary water blocks from producing formations.''2 The use of foam in gas storage reservoirs to reduce gas leaks and to increase storage capacity has been considered in recent years. Foam has also been investigated as an oil displacing agent,3-5 and as an agent to improve the mobility ratio in a waterflood.6-10 However, the mechanism by which the gas and liquid phases comprising the foam flow through a porous medium has nor been described adequately. Normally, when two immiscible phases (gas and liquid) flow concurrently through a porous medium, each phase follows separate paths or channels. At given saturations of the two phases, a certain number of channels are available to each phase, and as saturations change, the number and configuration of the channels available for each phase also change. The effective permeability of each phase is a function of the saturation of that phase only, and the flow of each phase can be described by Darcy's law. When foam is present, the effective permeability of the porous medium to each phase is greatly reduced compaced with permeabilities measured in the absence of foam.6,11-13 Based upon the observed flow of surfactant solutions and gas in capillaries, it has been concluded that the gas and liquid may flow separately or they may flow combined as foam. At least four mechanisms have been postulated to explain how fluids flow with foam present: 1. A large portion of the gas is trapped in the porous medium and a small fraction flows as free gas, following Darcy's law.9 2. The foam structure moves as a body; the rate of gas flow is the same as the rate of liquid flow. 6,13 3. Gas flows as a discontinuous phase by breaking and re-forming films. Liquid flows as a free phase.12 4. A portion of the liquid and gas move as a foam body while excess surfactant solution moves as a free phase.14 It also has been suggested that different flow mechanisms exist for high quality (dry) foams made from dilute surfactant solutions and for foams made from more concentrated solutions. Studies conducted on the flow of foam through capillaries have shown that a plug-type flow occurs and that foam flows as a body.6,15 Fried also has described the flow of foam through porous media as

Jan 1, 1969

By W. B. Hickman, J. J. Pickell, B. F. Swanson

Many physical properties of the porous media-immiscible liquid system are dependent upon the distribution of fluids within the pores; this in turn, is primarily a function of pore structure, liquid-liquiri interfacial tcnsion and liquid-solid wetting conditions. The cnpillary pressure hysteresis process provides a means of investigating the influence of pore structure upon fluid distribution for consistent sur/acc conditions. Invcstigntiot2s indicate that residual non-wetting-phase saturations /ollozuir,g the imbibition process (i.e.,wetting phase displacing nor2-wctting phase) are dependent upon both pore structure and initial non-raetting phase saturation and suggest that residual fluid is distributed as Discontinuous ,globules, one to a few pore sizes in dimension, thrnugh the entire range 01 pore sizes originally occupied. It clppears that air-mercury capillary pressure dala adequntely rellect the distribution of fluiris in a water-oil system when strong wetting condition preliczil. An oil-air counter-current imbibition technique has also been found to provide a rapid rneans of obtaining residual-initial saturation data. In a majority of cnses, rcsidual saturations detf-rniinecl from the oil-air or air-mercury process reasonably approximate residual oil saturation following water drive of a strongly water-wet medium. INTRODUCTION A reliable estimate of recoverable reserves depends not only on the amount of original oil-in-place but also on pore geometry and distribution of fluids within the pores. A critical parameter determining the recovery from a reservoir under waterflood, for example, is the amount and distribution of residual oil within the various rock types present. The purpose of this paper is to investigate the mechanism of capillary trapping and assess its importance in laboratory measureqents of residual oil saturation. The degree of wettability of a reservoir rock is recognized as an important factor in waterflood or imbibition experiments. In this paper, however, only the water-wet case has been considered. Considerable experimental evidence1 suggests that for water-wet rocks, capillary forces predominate in tile distribution of fluids and that viscous forces in the range normallv of interest in the reservoir have a minimum influence on residual oil saturation. It follows that if the ultimate recovery is controlled by pore geometry, a unique residual non-wetting phase saturation should exist for a given set of initial conditions. Two laboratory procedures.found to be extremely useful in the study of pore structure and degree of fluid interconnection at various saturations are decribed. Although air-nercury capillary injection curves have been used2 previously to characterize the drainage case, the withdrawal or imbibition case can provide valuable supplementary data. The air-mercury process, however, has several disadvantages; it is difficult to run in a sufficiently accurate manner, mercury does not always act as a strongly non-wetting liquid and in the air-mercury process the sample is rendered unsuitable for future analyses. An alternative process is described in which air is the non-wetting phase and naptha, hentane, octane or toluene is the wetting phase. INTERFACIAL TENSION AND CAPILLARY PRESSURE Interfacial tension between immiscible fluids is due to the difference in attraction of like molecules as compared with their attraction to molecules of the neighbouring fluid. This net attraction results in a tension at the interface. To extend the interface; thus, interfacial tension s can also be thought of as free surface energy. Interfacial tension is normally expressed as dynes/cm, and interfacial energy is measured in ergs/cm2 hence, both have dimensions ml A 2 and are numerically equal.

Jan 1, 1967

By P. Datta, L. L. Handy

For a wetting phase displacing a nonwetting phase from a porous medium the distribution of the residual fluid may depend on displacement conditions. Although this subject has been debated in the literature, only a few experiments have been cited to support the various conclusions. Experimental results presented in this paper show that fluid distributions are dependent on imbibition pmcedures. Results agree qualitatively with predictions from the pore doublet model. If the rate of water imbibition is restricted, the nonwetting phase is trapped preferentially in .the larger pores as expected. But if the rate of water imbibition is unrestricted, trapping occurs somewhat more in the smaller pores. These conclusions were deduced primarily from relative permeability measurements. INTRODUCTION Relative permeabilities are known to depend on the saturation history of the porous medium. For either continuously increasing or continuously decreasing wetting phase saturations, however, they have normally been assumed to be single-valued functions of saturation. Severat investigators have compared relative permeabilities measured by different methods and have reported acceptable agreement.1-5 Studies of simple models of the displacement process suggest that the distribution of residual fluids can be influenced by the displacement method. If this is so, relative permeabilities also should depend on the method of displacement. These predictions have not been supported by experimental data. The object of this paper is to investigate experimentally some of the predictions of these model studies. The particular question of importance is whether steady-state and unsteady-state methods should be expected to give the same values of relative permeability. Much of the reported data showing agreement between steady-state and unsteady-state methods have been obtained on un consolidated sand. Johnson, Bossler and Naumann, however, compared steady-state and unsteady-state results for Weiler sandstone and found no significant differences.5 Levine made a detailed study of pressure and saturation distributions for a laboratory waterflood in a consolidated system but made no attempt to calculate performance from steady-state relative permeability data Bail and Marsden in a somewhat similar set of experiments did attempt a comparison of observations with predictions but the results were inconclusive.7 Excluding gravity, two forces affect the distribution of wetting and nonwetting phases during an unsteady-state, immiscible displacement: viscous and capillary forces. If relative permeabilities are indeed the same when measured by steady- or unsteady-state methods, the fluid distributions at the same fluid saturations must be similar for both methods. Furthermore, the implication is that the relation between capillary and viscous forces is the same for both methods insofar as the effect on microscopic fluid distributions is concerned. Unsteady-state displacements are normally run at high rates to maximize the ratio of the viscous to capillary forces. The objective is to reduce the effect of capillary forces on macroscopic fluid distributions. For floods in which the phase which wets the porous medium displaces the nonwetting phase, capillary forces compete with viscous forces in determining fluid distributions.8 The residual nonwetting phase is discontinuous. Competitive aspects of viscous and capillary forces and the resulting effect on water-oil distribution in porous media have been illustrated in an elementary way using a pore doublet model. This model and predictions from it have been discussed at considerable length by Rose and Witherspoon, Rose and Cleary and by Moore and Slobod.9-13 Rose and Witherspoon show that so long as water invades the model at a rate equal to or greater than the free imbibition rate, the water will move through the larger capillary of a doublet first and trap oil in the smaller capillary. On the other hand, if a negative pressure is applied to the water

Jan 1, 1967

By C. Dirksen

This study investigates whether oxygen consumption during segregated forward combustion may be affected by natural convection. Linear theory indicates that thermal instability occurs in a horizontal porous medium when the modified Rayleigh number N 'Ra exceeds a critical value of about 40. Most experimental results, including those described here, indicate this value to be about an order of magnitude smaller. N ', is derived from "differential similarity" and the proper time scale factor is also obtained. The critical value of N'ra for a horizontal, water-saturated porous medium at atmospheric pressure subjected to a uertical temperature gradient was found to be about 3.0. The ambiguity in the value of N'R, arising when fluid properties cannot be assumed constant is indicated. Natural convection was observed for oblique systems below N 'Ra , while simultaneous horizontal flow did not affect N Racrit. The expected range of values for the various parameters pertaining to segregated forward combustion is indicated and the corresponding values of N 'Ra calculated. Only under very favorable conditions, such as high permeability, high pressure and low top temperature, can the critical value of N 'R, be exceeded. Thickness of the convection layer should not become too large, since the time required to obtain significant mixing is proportional to the square of the thickness. For an oblique combustion frant the conditions for significant mixing may be less strenuous. LNTRODUCTION Laboratory1,2 and field314 experiments have shown that a markedly upgraded oil can be produced by in situ forward combustion. Under certain conditions injected air may move through a high permeability, gas-saturated channel over the oil-saturated layer, while combustion takes place, exclusively, at the horizontal interface (Fig. 1). Such a channel may be a natural gas cap, or it may result from a forward combustion process in which injected air and combustion gases fingered through the top of the formation, causing a premature breakthrough at the production well. In this particular case of segregated forward combustion, the air stream moves parallel to the combustion front. Molecular diffusion and transverse mechanical dispersion will be insufficient to cause all the oxygen in the air stream to reach the combustion front. When the oxygen in the air stream is not completely consumed in the combustion process, it will reach the production well and create a serious fire hazard. In addition, the reduced efficiency may make the process economically unprofitable while also serious corrosion problems may be experienced in the production we11.3 This undesirable situation should, therefore, be prevented. However, in a situation as depicted in Fig. 1, an unstable temperature gradient exists in the gas phase. The combustion is taking place at the lower boundary at a nearly constant temperature that may be anywhere from 600 to 1,200F, depending on the conditions. The upper boundary of the gas phase is overlying, impermeable rock; its temperature might be around 100F, but may increase to 300F or higher due to conduction and convection as the combustion process continues. This represents a steep temperature gradient and it is conceivable that under those conditions free or natural convection will take place in the gas

Jan 1, 1967

By H. Y. Jennings

An extensive field and laboratory experimental program was carried out to compare the waterflood behavior of carefully preserved soft and unconsolidated cores with measurements on the same cores after extraction. Results obtained from using idealized consolidated and unconsolidated porous media in which wettability could be carefully controlled were contrasted with the preserved core data. The controlled tests on idealized porous media investigated the effect of wettability, flood rate, core length, core permeability and consolidation on the displacement of high viscosity oils. It was concluded from these studies that waterfloods are more favorable when carried out with crude oil in preserved soft and unconsolidated cores than with the same cores after they have been extracted and re-saturated. Waterfloods are usually more favorable when carried out with crude oil in extracted .soft and unconsolidated cores than with refined oil of the .same viscosity in the same cores. The less favorable behavior of extracted soft and unconsolidated cores compared to preserved cores is due to alteration of the core by extraction. Preserved cores saturated with native water and oil should be used for laboratory displacement experiments because they more accurately reflect true reservoir behavior INTRODUCTION The demand for low gravity crude oil created by refinery modernization has focused attention on increasing the production of this viscous crude oil. Billions of barrels are in place in fields that are depleted, or nearly depleted, by primary production mechanisms. Since low gravity reservoirs are relatively recent, geologically, the solid matrix material is usually soft and unconsolidated sand. Such formations are also characterized by a high clay content. Evaluation of sophisticated oil recovery processes with the associated high capital investments has increased the demand for special core analysis tests on material from these soft and unconsolidated sand reservoirs. Data in this paper have resulted from an extensive field and laboratory experimental program. The initial objective vras to provide soft and unconsolidated sand cores for laboratory measurements with as little alteration as possible from their reservoir condition. A secondary objective was to compare the waterflood behavior of these carefully preserved cores with measurements on the same cores after they had been extracted and resaturated. When the comparison showed markedly different behavior the final objective was to attempt to explain the difference by making measurements on idealized porous media free of clay in which initial wettability could be carefully controlled. EXPERIMENTS MATERIALS Core Material The preserved cores used in this study were obtained with as little alteration as possible from their reservoir conditions. Special techniques were developed to satisfy this objective. The cores were cut with native crude whenever the crude had the necessary properties to satisfy the minimum drilling fluid requirements. A pure hydrocarbon chromatographic tracer was added to the crude to provide a simple, safe and inexpensive method to distinguish between oil filtrate and formation oil. Details of the tracer technique used in this study and the procedures used to handle and process the soft, unconsolidated cores have been published.' The core samples were then preserved and packaged at the well site with a dip-applied strippable plastic, or by using a rubber sleeve from a rubber-sleeve core barrel. The final step was to insure that the carefully taken and preserved samples were not altered by processing in the laboratory. Some soft and poorly consolidated sands were carefully shaped in the lab and encased in a protective mounting without disturbing their three-dimensional integrity. Many samples were so poorly consolidated that they literally flowed from the package when the seal was broken. A technique was developed to obtain cylindrical plugs from these samples by using liquid nitrogen as the drilling fluid in a conventional core-drill, drill-press assembly. The idealized porous media consisted of sintered rods of Alundum RA 139, outcrop sandstone identified as Al-hambra sandstone and packs of No. 130 Nevada sand. This sand is 95 per cent 140-200 mesh. (Physical properties of these porous media with those representative of preserved cores are in Table 1.) The natural cores and most of the idealized porous media were 1 1/2-in. in diameter and 3-in. long. The samples used to study the effect of core length were 1 1/2-in. in diameter and 12-in. long. Liquids Water and oil produced from the formations being studied were protected from the atmosphere and collected in carefully cleaned glass or plastic-lined containers. The production was free of chemical additives such as emulsion breakers and corrosion inhibitors. The crude oil was

Jan 1, 1967