Search Documents

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

Immiscible displacement tests were performed in a consolidated sandstone core over the interfacial tension range from less thdn 0.01 to 5 dynes/cm to better define how interfacial tension (IFTJ reduction can lead to increased oil recovery. Data obtained were displacement efficiency at breakthbrough vs IFT for both drainage and imbibition conditions. These tests simulate water flooding under oil-wet and water-wet conditions, respectively. Results of the study have shown that displacement efficiency under both oil-wet and water-wet conditions can be markedly improved by a sufficient reduction in IFT. In the particular porous media used and for the low pressure gradients employed, the IFT must be reduced to a value less than about 0.07 dynes/cm to achieve increased recovery at the time of breakthrough of the injected phase. Below 0.07 dynes /cm, further small reductions in IFT result in large increases in displacement efficiency. Observed increases in recovery were obtained at pressure gradients which are well below those which can exist in the interwell area of a reservoir under water flood. The effect of pressure gradient on recovery is discussed. INTRODUCTION A residual oil saturation remains in rock which has been water flooded because, under usual reservoir conditions, the driving force which can be generated is inadequate to expel oil trapped by capillary forces. Since these capillary forces can be reduced by reducing the IFT a frequently studied method of increasing oil recovery has been the use of surfactants to reduce the water-oil interfacial tension. Mungan l observed improved recovery in both water-wet and oil-wet systems at 1.1 dyne /cm IFT, finding that the amount of improvement was greater in oil-wet systems. Moore and Blum,2 working with visual flow cell micro-models, con- cluded that recovery cannot be improved in water-wet systems by injecting surfactant solutions. They calculated that, for a pressure gadient of 1 psi/ft in their model, the IFT must be reduced to 0.03 dynes/cm to release oil trapped under water-wet conditions. Berkeley et al.3 indicated that, at representative field flooding rates, the IFT must be reduced to 0.001 dyne/cm to improve recovery, Thus, there is a wide variation of opinion about the IFT levels needed to improve recovery under field conditions. These previous investigators all have used as their experimental method the addition of surfactants to the injected water. The use of surfactant solutions to reduce IFT creates two experimental problems: (1) the loss of surfactant through adsorption on reservoir rock can obscure the true IFT value which exists at the displacement front, and (2) at very low values of IFT emulsifica-tion of oil and water commonly occurs. It is difficult, therefore, to determine whether the increased recovery is caused by IFT reduction as such, or is instead caused by emulsification. Also the properties of the available surfactants limit the IFT range which can be studied. In the present study the purpose is to better define how interfacial tension reduction can lead to increased oil recovery. A matter of great interest is the amount of recovery improvement potentially achievable in this way. The study was made using very low pressure gradients which was well within the range achievable by water flooding in the interwell region of petroleum reservoirs. A unique experimental approach was chosen to avoid adsorption and emulsification problems, and to allow convenient control of IFT. The fluid phases used were the equilibrium vapor and liquid phases of the methane-n-pentane system. The interfacial tension level was varied by changing the equilibrium pressure of the methane-pentane system over the range from 1,200 psia to near the critical pressure (2,420 psia). All tests were performed at a controlled temperature of 100F. The IFT-vs-pressure relationship for the methane-pentane system was based on the data of Stegemeier and Hough,6 and on new data obtained in this study. With this experimental approach it has been possible to study displacement at lower values of IFT than have been previously investigated.

Jan 1, 1967

By D. E. Menzie, D. L. Dauben

This paper discusses the physical parameters involved in the slow flow of high molecular weight polymer solutions in porous media. The interacting effects of polymer properties and porous media properties on flow performance are considered. Experiments were conducted with the polyethylene oxides, with molecular weights ranging from 200,000 to over 5,000,000. Frontal advance velocities ranged from I to 30 ft/day. The porous matrix consisted of a flow cell packed with glass beads. Polymer solutions were characterized by viscosity and normal stress measurements. Under certain conditions, unexpectedly high flow resistance was observed. This behavior was observed to be a function of flow rate, pore size, polymer molecular weight and concentration. The polymer solutions exhibit "dilatant" flow behavior in porour media in contrast with the pseudo-plastic behavior in simple flow systems. A theoretical ex-planation of such behavior is presented. INTRODUCTION The use of polymers in the injected water of a waterflood increases oil displacement efficiency by reducing the mobility (k,/pu) of the driving phase. A reduced driving phase mobility results in improvements in the areal sweep efficiency and in the vertical coverage in stratified reservoirs. This mobility reduction may be achieved by a permeability reduction, a viscosity increase or by a combination of the two. Early attempts at increasing the injected water viscosity were not successful because of the poor economics involved. The use of such materials as glycerin, sugar or glycols to increase water viscosity was not economically feasible. Attempts in using certain naturally occurring polymers were not too successful because of the high polymer losses to the rock. A high molecular weight, partially hydrolyzed polyacry-lamide was introduced 3 years ago as a waterflood additive.',' Initial work by Pye' indicated that the presence of these polymers in dilute concentrations decreases the water mobility 5 to 20 times more than would be expected from measurements of the solution viscosity. Such an effect would be of obvious economic value since only a small polymer concentration would be required to accomplish a large reduction in water mobility. This research was designed to study the basic flow mech anisms of polymer solutions in porous media. The interacting effects of polymer properties and porous media properties on flow performance are considered. This study indicates some of the conditions under which these interactions can lead to significant mobility effects. is valid only for Newtonian fluids. For polymer solutions and other non-Newtonian fluids the equation must be modified to consider that viscosity is a variable quantity. Modification of the Darcy equation to include non-Newtonian effects has been the subject of several recent investigations."" The modifications generally adopt a rheological model. such as an Ellis or Dower law model. to Dorous media by defining some characteristic channel radius. Most of these studies showed the porous media flow behavior to be predictable from viscometric data. Investigations involving the flow of high molecular weight polyacrylamide solutions through cores have generally encountered the high flow resistances reported previously by Pye. This high flow resistance has generally been attributed to an in-depth permeability reduction, as evidenced by a reduced permeability to water which has displaced polymer solution from the porous media. Marshall and Metzner" report high flow resistances, which they attribute to viscoelastic effects. In this paper, the effects of polymer molecular weight, pore size, flow rate and concentration are considered. The polymers used are the polyethylene oxides known as Poly-ox. This class of polymers was used because of its availability in a wide range of molecular weights and because of its known ability to propagate well through porous media. THEORY PHYSICAL PROPERTIES OF POLYMER SOLUTIONS The polymer molecule dissolves in water by means of hydrogen bonding, but retains some of its own structural identity while in solution. Nonionic polymers, such as polyethylene oxide, are generally considered to have a random coiling configration. This type of molecule has the ability to "sequester" or hold a large volume of solvent within its coils in a manner similar to that of a sponge. The random coil is easily deformed and under an applied stress

By R. R. Goddard

By use of the frequency response method with a radioactive tracer, it was possible to determine fluid dispersion and distribution in a natural consolidated and an unconsolidated medium. Measurements were made in a linear flow system at oleic saturations of 69 per cent in the consolidated medium, and 100 per cent in both media. Dispersion and distribution were obtained by measuring the amplitude attenuation and the phase velocity of sinusoidal waves with a dual monitor apparatus. The gamma ray emissions permitted in situ measurements at any distance along the porous samples. One result of importance was that the effective diffusivity increased as the wave length increased. As a consequence, a dispersion coefficient appropriate for the injection of large slugs might exceed the value measured by use of small slugs. Since flow models based solely on fluid velocity and an effective diffusivity coefficient imply that the diffusivity should be independent of frequency, such representations were not adequate for the data of this study. A comparison was made with a capacitance model of porous media with deadend PVs, but even this model was not completely adequate. By using attenuation and phase velocity data, fluid dispersion can be predicted without postulating a differential equation satisfied by the tracer concentration, thereby eliminating the need of a complicated model to represent dispersion. INTRODUCTION The flow of similar miscible fluids through a porous medium can be fairly adequately described by two parameters: the average fluid velocity and the effective diffusivity.1-3 It has been pointed out recently, however, that significant discrepancies exist between this representation and the experimental data.4-7 An improved agreement can be obtained by introducing additional parameters based on the concept of dead-end pores. The purpose of the present investigation was to find out whether the frequency response method could be used to measure the relevant parameters. The method was used in the following form. A stream of fluid was flowed at a constant rate through a sample of porous material and the concentration of a radioactive tracer in the fluid was varied sinusoidally at a fixed frequency. The effects of flow through a porous medium are a decrease in the amplitude of the concentration wave and an increase in the velocity of the peaks of the waves above the average velocity. Attenuation and phase veIocity of the waves were measured as a function of frequency and fluid velocity. The simple two-parameter model implies that the diffusivity should be independent of frequency. Data reported in this paper show that the diffusivity decreases as the frequency increases. Hence, as shown also by many others, the two-parameter model is not completely adequate. Coats and Smith 5 used two additional parameters in their model: the volume of the dead-end pores and the rate of mass transfer between dead-end pores and the flowing stream. Their capacitance model of a porous medium containing some stagnant fluid, to which transfer occurs by molecular diffusion, did not explain the dispersion results of either the present study or of theirs. Instead, the capacitance effect can be better described as the result of extreme velocity variations within the pores of the medium, with transfer between the velocity zones by convection. EXPERIMENTAL PREPARATION OF POROUS MEDIA Data reported here were obtained with a natural Berea sandstone, and with two packings of 20/25 mesh sand. Pertinent information about these media is listed in Table 1. The sandstone core was first sealed with a non-

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 R. W. Parsons, P. R. Chaney

Oil recovery by the imbibition mechanism can be important in fractured carbonate reservoirs with a bottom water drive. Laboratory experiments were performed on water-wet carbonate rocks to model this process. The results are interpreted in light of the applicable scaling laws. The imbibition production behavior of a pillar of rock subjected to a slowly rising water table can be synthesized from data on total immersion experiments on small rock samples. If the rate of water table rise is slow enough and if the physical properties are such that the imbibition zone is of a small vertical extent, then only the final oil recovery on the small samples need be known. Applying laboratory imbibition data to a real reservoir should be done with caution because of possible wettability problems and the unknown behavior of two-phase flow in fractures. INTRODUCTION The spontaneous taking up of a preferential wetting fluid into a porous medium with the simultaneous expulsion of the contained fluid is termed imbibition. One of the many physical imbibition processes is the production of oil by the taking up of water into a water-wet reservoir rock. In certain types of heterogeneous reservoirs, e.g., fractured, it is realized that imbibition alone may be the dominant mechanism. Theoretical studies of imbibition oil recovery have followed several lines of attack. The behavior of certain mathematically tractable and conceptually simple models have been analyzed. Birks1 considered a bundle of vertical capillaries model, while Perotti et al.2 and Kelemen3 used a model of a series of vertically oriented fractures. If the oil production behavior from a block of reservoir rock is known, then the cumulative behavior of a vertical stack of these blocks subjected to a rising water table can be determined by a method outlined by Aronofsky et al. 4 Knowing the capillary pressure and relative permeability functions should uniquely set the imbibition behavior for a given sample. Solutions of the partial differential equations describing the process, however, present some formidable problems. Digital computer techniques have been applied to displacement processes in which the capillary pressure has been included.5"8 More recently, similar techniques have been used by Blair 9 to calculate countercurrent imbibition behavior. The fact that imbibition behavior is described by the usual fluid flow equations (Darcy's Law and continuity equation) implies that the derived scaling laws should apply equally to this situation. Use of model experiments to study imbibition has been employed by Graham and Richardson 10 and Mattax and Kyte,11 who actually tested the scaling law applicability. This study is concerned with scaled model imbibition experiments on water-wet carbonate rocks. Of particular interest is the slowly rising water table experiment. This simulates a highly fractured reservoir with the water table rising uniformly through the vertical fracture system. Experiments were performed on different sized samples and with different rates of rise of the water-oil interface surrounding the rock pillar. Total immersion imbibition experiments on smaller rock samples were performed to test the basic linear scaling laws and to see if these simpler tests could be used to synthesize the rising water table experimental results. THEORY IMBIBITION SCALING LAWS Scaling laws for the flow of two incompressible immiscible fluids in a porous medium have been presented by Rapoport.12 These are derived directly from Darcy's Law for the individual phases and the continuity equation. Results from a laboratory model experiment will duplicate those from some reservoir prototype if a constant proportion is maintained among the three forces acting on the fluids — the capillary pressure gradient, the gravitational gradient and the flowing pressure gradient. Operation of the model (i.e., fluid injection or withdrawal) must be conducted in accordance with specific equations. If the

Jan 1, 1967

By N. Mungan

A study was made of the effects of wettability and interfacial tension on the immiscible displacement of a Iiquid by another Iiquid from porous media. The influence of viscosity ratio was also investigated. Porous media used were polytetrafluoroethylene (TFE) cores prepared by compressing TFE powder under different pressures. lt is shown that displacement of a wetting by a nonwetting Iiquid is always less efficient than the displacement of a nonwetting by a wetting fluid, all other things being equal. In the former case, the recovery efficiency can be increased substantially by either reducing the interfacial tension or increasing the viscosity of the displacing fluid. A qualitative discussion is given on the implications of this work to the recovery of crude oil by water flooding. INTRODUCTION The high cost of oil exploration and new recovery schemes makes it imperative that waterflooding be conducted under conditions favoring most efficient oil recovery. To improve oil recovery by waterflooding, it is essential that the role played by interfacial forces in the entrapment of residual oil be studied and understood. Interfacial phenomena in natural rock, connate water and crude oil systems are very complicated because of the complexity of the natural liquids found in petroleum reservoirs, because of our inability to adequately describe the geometrical structure of the porous media and because of a lack of understanding of physical and chemical interactions between the liquids and surface of the pores. The problem becomes further complicated when one tries to elucidate the role of interfacial phenomena in fluid flow. Numerous studies of the displacement of oil by water under different interfacial tension or wet-tability conditions have been made.lm3 These studies have been performed in silica, alundum or sandstone systems using water and paraffin oil and also some surface active material to control the interfacial tension s and the contact angle 0. Unfortunately, the high energies of various interfaces involved favor adsorption and orientation of the surface active material at the interfaces. Also the surface active material concentration at the interfaces exceeds that in the bulk of the liquid phases. Such surface excess may cause the surfactant distribution, the contact angle and the interfacial tension to differ from their measured static equilibrium values and makes interpretation of the displacement experiments difficult. Furthermore, as changes in s also lead to changes in cos 0, the role played individually by one of these parameters in the displacement becomes obscured by the effect of the other. To circumvent these difficulties, a low surface energy solid and true solutions or pure liquids should be used. Use of a low surface energy solid minimizes adsorption and orientation effects at the solid-liquid interfaces. By controlling s and cos 0 through use of selected pairs of pure liquids or true solutions rather than by surfactants, the adsorption effects at liquid-liquid interfaces are eliminated. In the present study TFE cores were used as the porous media. Liquids used were water, sucrose solutions, paraffin oils and benzyl, n-butyl and isobutyl alcohols. The interfacial tension was varied from 40 to 1.1 dynes/cm by suitably choosing the liquid pair. A surface active material was added to the water-oil system only in the case where interfacial tension of 0.5 dyne/ cm was desired. No precise changes of cos 0 were attempted. However, either the displaced or the displacing liquid could be made the one which preferentially wets the TFE surface. Using sucrose solutions and blends of paraffin oils proved to be a convenient way of changing the viscosity ratio between the displaced and displacing liquids. The present investigation examines the effect of interfacial tension, wettability and viscosity ratio on the immiscible liquid-liquid displacement from porous media. Observed data permit a

Jan 1, 1967

By W. B. Gogarty

With the use of polymer solutions in secondary recovery operations, the need has developed to understand the mobility control mechanism. This study investigated mobility control by considering both permeability and rheological effects. Experiments used a high molecular weight, partially hydrolyzed polyacrylamide polymer. Flow studies took place in reservoir and Berea cores having zero oil saturation. Effective size of the polymer flow unit was inferred from Nuclepore filter tests. Clay studies indicated the particle size capable of decreasing the core permeability. Flushed permeabilities measured the approximate core permeabilities with flowing polymer solutions. These permeabilities were considerably lower than original values. With mobility data and the flushed permeability, maximum effective viscosities were determined for polymer solution flow in a core. Effective viscosities showed that rheological properties play an important part in mobility control with polymer solutions. The study showed that permeabilities decrease and stabilize with polymer flow. At the lower permeabilities, high shear rates exist in the cores. Because of the pseudoplastic character of the polymer solution, the high shear rates caused low effective viscosities. This condition pointed to the inefficient use of the potentially high viscosity of the polymer solution at low shear rates. INTRODUCTION In the oil industry, a great deal of interest is being shown in the use of polymer solutions for secondary recovery and a number of polymer floods are being performed in the United States.1 Some of these floods have become commercial while others have been reported as failures. A number of floods are still in progress and remain to be evaluated. With the advent of polymer flooding, the need developed to understand the mobility control mechanism in porous media. Rheologically, the polymer solutions behave as pseudoplastic fluids. Investigators have studied the rheology for this type of non-Newtonian fluid in porous media?-6 Results have also been reported on the ability of polymer solutions to decrease mobility.7-10 Some investigators indicated that polymer solutions show a different rheological behavior in cores than on the viscometer.7,8 This study is representative of an initiatory investigation on mobility control with polymer solutions. Some effects are considered on a qualitative basis. Mobility control is considered from the standpoint of both permeability and rheology and, to the extent possible, permeability and rheological effects are considered separately. Results show that mobility control occurs from a reduced permeability caused by polymer retention and an effective viscosity determined by the average shear rate in the core. In this respect, rheological behavior is as important in mobility control as permeability. PROCEDURES A high molecular weight, partially hydrolyzed, polyacrylamide polymer was used in the investigation. Solutions of various concentrations were made with this polymer in different types of water. Some solutions used distilled water and others used tap water. The tap water contained approximately 500 ppm dissolved solids. As used in this paper, concentrations in ppm refer to parts by weight solute to one million parts by weight of solvent. In the tap water, calcium and magnesium accounted for from 60 to 80 ppm of the solids, the remainder being mostly sodium chloride. All solutions contained a bactericide, and anti-oxidants and antiprecipitants were added to the tap water. Capillary and concentric-cylinder* viscometers determined the rheological properties of solutions. The capillary viscometer is described elsewhere.6 Berea and reservoir sandstone were used as core materials for this investigation. Cores were 1 in. in diameter by approximately 2 to 3 in. long. Prior to mounting, the Berea cores were heat-treated for 3 hours at 850F. This treatment decreased the water-sensitivity of the clays in the cores. The reservoir cores were cleaned with solvent and

By A. R. Gregory

A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy. S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King, 1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. 5- and P-wave velocity measurements on Bandera sandstone and Solenholen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2 and by Gregory 3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fate sand by Gregory, 3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King 1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8 A new shear wave velocity laboratory apparatus has been developed for testing rock samples under elevated pressures. This apparatus employs liquid-solid interfaces as mode converters. The technique for generating pure shear waves at liquid-solid interfaces was developed from studies of critical angle methods.6 Two identical transducers are placed in contact with opposite ends of the rock sample, one transducer acting as an ultrasonic energy emitter or generator and the other acting as an ultrasonic energy receiver or motor. The P-wave energy emitted from the piezoelectric crystal in the generator transducer is converted into pure S-wave energy at the Drakeol-aluminum

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 I. H. Silberberg, R. L. Reed, O. K. Kimbler

interfacial films have frequently been observed at interfaces between certain crude oils and water. Several investigators have postulated that the presence of these films should influence the efficiency of oil recovery in water drive or waterflood operations. They may also influence the stability of emulsions which are sometimes a problem in petroleum production, and may be a factor in the fonation of paraffin deposits in oil well tubing and flow lines. This paper presents a technique with which a modified Langmuir film balance may be used to study the compressibility and collapse pressure of these natural interfacial films. Experimental data are presented for several nude oil-water systems. Data developed are used to infer the phase state of the film as a function of such variables as rate of reduction of interfacial area, ionic composition of the subtrate and pH of the subtrate. A film of known physical characteristics is shown to have a sinnificant effect on oil recovery from an unconsolidated sand pack. Possible applications of these results to petroleum production are discussed. INTRODUCTION The use of water to displace petroleum from reservoir rocks is of major importance both as a primary and a secondary recovery process. As water invades the rock, oil is completely displaced from some pores and left as a discontinuous phase in other pores. The manner in which water moves from pore to pore is strongly influenced by capillary forces. In view of the complexity of reservoir fluid systems, there can be little doubt that complicated interactions take place at both the liquid-solid and oil-water interfaces. One of the more interesting, and least understood, of the phenomena which take place at the oil-water interface is the formation of interfacial films. These films are believed to result from the adsorption of high molecular weight polar molecules at the interface.l.2 Presence of such molecules may cause a striking alteration in interfacial tension. When the oil-water interfacial area of certain crudes is rapidly reduced, a thin region (film) about the interface assumes the appearance of a solid membrane, and striations, wrinkles and gross distortions may occur. If such a film is solid, it should greatly alter the interfacial tension normally assumed to exist between the oil and water phases. If the membrane is continuous, a solid phase would separate the oil and water. Interfacial films between crude oil and water were observed in 1949 by Bartell and Niederhausers who commented upon the apparent rigidity of the films and their possible importance in the petroleum industry. Morrell and Egloff4 had earlier attributed the extreme stability of emulsions of sea water in fuel oil to very stable asphaltic films. Numerous investigators have observed rigid films in the course of crude oil-water interfacial tension determinations by the pendent drop method. Several investigators5,6,2 have separated inter-facially active materials from crudes and attempted to characterize them chemically. Reisberg and Doscher,2 using Ventura crude, showed the interfacial tension against water (as measured by the pendent drop method) to be affected by aging, contraction and expansion of the interface, and the pH of the water. These investigators attributed the adhesion of oil to a water-wetted surface and the distortions of flow paths in glass capillaries to the presence of rigid films. Dodd7 has studied the interfacial viscosity of adsorbed films and found them to be non-Newtonian in behavior. Craighead and Harvey8 reported a series of displacements in tubes packed with 60 mesh glass beads. They interpreted the results as indicating an effect of stearic acid films on waterflood recovery and imply that natural films may produce similar results.

Jan 1, 1967

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By J. F. Evers, S. Sadiq, F. W. Preston, G. W. Swift

Proper scaling of unsteady-state flow in laboratory gas reservoir models can be achieved if the permeability/fractional porosity ratio is less than 0.5 md. Porous media meeting this specification can be obtained from mixtures of sand and Portland cement provided (1) the water/cement ratio is in the range of 0.35 to 0.40, (2) the sand/cement ratio is in the range of 1/1 to 2/1, (3) the sand grain size is in the range of 48-100 to 18-28 mesh, (4) the mortar curing time is 7 days or more and (5) the cured mortar is dried under vacuum to reduce the temperature level employed during the drying process. The areal homogeneity of the mortars tested is reasonably good, although significant trends in vertical permeability were noted. The Klinkenberg and non-Darcy coefficients of mortar are such that proper scaling of these effects in the laboratory model could be expected. INTRODUCTION A long term deliverability prediction for a gas well is generally based on data determined from flow tests that are of short duration relative to the producing life of the well. These data, together with an estimate of gas reserves in place, form the basis for establishing the parameters of the mathematical model used to predict performance at future times. Results obtained from a mathematical model represent only an interpretation of the physical process being modeled. Experimental verification of the production process from discovery to abandonment is needed where production occurs under carefully controlled conditions. Ideally, one might consider such an experiment to be conducted in the field. Even if the cost of such an experiment could be justified, the time required would make it impractical to consider. However, if the time scaling problem could be handled properly, an experimental model of the gas reservoir production process could provide the necessary data to study the mechanisms of the production process in a reasonable time period. An experimental model has been designed and constructed using the data of this study as a basis. The model reservoir system consists of one producing well located at the center of a right-circular cylinder of porous material. Ten to 15 observation wells are located at strategic offset positions from the producing well. These wells, and the producing well, are fitted with pressure transducers so that the pressure-time history can be recorded at various radial and azimuthal positions in the drainage area during the entire production life of the system. Time is scaled in the model so that the time required to produce gas from discovery pressure to abandonment pressure is about 2 hours. This compares with a 20-year producing life in the model's field counterpart and gives a time scale factor of approximately 1:100,000. DESIGN CRITERIA To establish design criteria for the experimental model, the reservoir is considered to be a finite right-circular cylinder of radius rr, and height h with a producing well of radius rw located at the center and completed through the producing formation. This reservoir consists of a porous medium which has a uniform permeability k and gas porosity and the gas contained in the reservoir has density p, viscosity µ and molecular weight M. Initially, the pressure at all points within the reservoir is uniform at po. and the reservoir temperature is T. Since the reservoir production process has been demonstrated to occur in two distinct transient modes for either Darcy or non-Darcy flow,l it is

By W. B. Gogarty

With flow of non-Newtonian fluids in porous media, effective viscosities are needed for use in the Darcy equation. These viscosities depend on the rock parameters and flowing conditions. In this investigation, rheological and flow results were correlated to develop an expression for calculating effective viscosities. Surfactant-stabilized dispersions of water in hydrocarbon were used with consolidated sandstone cores of different permeabilities. The average shear rate in a core was related to the permeability and porosity of the core and the frontal velocity through the core. By using the correlation obtained for determining effective viscosity, experimental and theoretical values were compared for each fluid system by determining average and maximum errors. For the fluid system showing the largest average error, the correlation fit 38 experimental points with an average and maximum percentage error of 4.3 and 14.8, respectively. For the fluid system showing the least error, average and maximum percentage errors of 1.8 and 4.3, respectively, were calculated using 29 experimental points. INTRODUCTION Use of non-Newtonian fluids in the petroleum industry is not new. Fluids of this type have been used for many years as fracturing agents and drilling muds. Recently, attention has focused on the use of high molecular weight polymer solutions for secondary recovery. Results have been reported on studies relating to the flow behavior of polymer solutions in porous media.1-4 Each study gives the viscosity in the Darcy equation as a function of the rheological properties of the fluid, the characteristics of the porous medium and the pressure gradient. The functional forms of the reported expressions differ somewhat in each study. The first three studies used unconsolidated porous media while the fourth used a single consolidated core. Because of the limited information available on the flow of non-Newtonian fluids in porous media, further investigation was undertaken in this area. We selected non-Newtonian fluids that did not contain high molecular weight polymers in order to minimize core plugging. The over-all objective was to gain a better understanding of the flow behavior of non-Newtonian fluids in consolidated porous media. Specifically, our goal was to determine the relationship between the average shear rate in a core and the frontal velocity. Because of this goal, it was necessary to determine the effect of permeability on the shear rate-frontal velocity relationship. We wanted to replace the viscosity in the Darcy equation with an expression relating to the average shear rate in the core and the rheological properties of the non-Newtonian fluid. THEORY For a time-independent, non-Newtonian fluid, the apparent viscosity is a function of the shear rate. This relationship is given by µ = F(v')...........(1) With Eq. 1, a single value of shear rate determines the apparent viscosity at that shear rate. With flow in porous media, a spectrum of shear rates occurs at any average throughput. An average shear rate v' can be used to represent the spectrum of shear rates that exists within the core. By using v' and the rheological data given by Eq. 1, an effective viscosity is determined for the flowing fluid. If v' is the correct average, µe becomes the appropriate viscosity for use in the Darcy equation. Under these conditions, flow rate-pressure drop relationships can be determined for non-Newtonian fluids. HOW the correct average shear rate might be related to the frontal velocity through the core and the characteristics of the porous media will now be considered. Later, the experimental data are presented that indicate the validity of these

By N. Mungan, F. W. Smith, J. L. Thompson

Adsorption of polymers and transport, rheology and oil recovery efficiency of their solutions were studied in the laboratory to evaluate the use of polymers in waterflood-ing. While a tenfold mobility reduction was obtained with polymer concentrations as low as 0.05 per cent by weight, the mobility reduction depended on the type of polymer, molecular weight, salinity and pH of water, crude oil and capillary properties of the porous media. Choice of a suitable polymer and a workable concentration will have to be tailored for each application. Little reduction in the residual oil saturation can be expected from polymer flooding. Improvement in the volumetric sweep efficiency is possible hut the extent of the improvement can best be evaluated by properly designed field testing. Some aspects of the field use of polymer floods are discussed. INTRODUCTION Waterflooding is a simple, inexpensive secondary recovery method and is being used widely. Innumerable laboratory studies have been made to unravel the fundamentals of the displacement of oil by water and to find the ways of most efficient oil recovery. These studies and a great many field case histories have revealed that the prime cause of poor oil recovery is the inefficient and incomplete sweep of reservoir volume by the injected water. Sweep efficiency is affected by many factors of which the mobility ratio is an important one. Mobility ratio M is defined here as the ratio of water to oil mobilities: M = (k»/y,r)/(k,JJJJ........(1) In Eq. I, the permeabilities are the effective permeabilities and depend on fluid saturations and, hence, change during the different depletion stages in a flood. A wide practice is to use the effective water permeability at residual oil saturation and the effective oil permeability at interstitial water saturation in Eq. 1. If the mobility ratio is greater than one, the mobility ratio is unfavorable and water, being more mobile than oil, would finger through the oil zone resulting in poor oil recovery efficiency. If the mobility ratio is favorable (one or less) the displacement of oil by water occurs more or less in a pistonlike fashion. In some waterfloods. the mobility ratio is unfavorable and any additives by which the mobility of water can be decreased would favor more efficient oil recovery. The thing to bear in mind, however, is whether or not the improvement in oil recovery is sufficient to more than pay out the cost of the additives needed. For example, materials like sugars, alcohols and glycerine reduce water mobility by thickening the water, but the cost of material requirement precludes any field application. For an additive to be useful in water-flooding, it must bring about a large reduction in water mobility at low concentrations; it must be adsorbed only negligibly; and it must not completely plug up the formation. Some synthetic organic polymers have shown promise of meeting these requirements and have been used in the field.'-W owever, no in-depth studies of the rheological, adsorption and oil displacement characteristics of polymer solutions have been reported. The present work is a study of these properties. EXPERIMENTAL In this work, concentrations are given on a weight per volume basis; 0.5 per cent concentration means 0.5 gm of polymer is dissolved in enough water (or NaCl solution) to make 100 ml. A bactericide, usually 0.1 to 0.2 per cent by volume of 38 per cent formaldehyde solution, was used in the polymer solutions. The NaCl solution was 30,000 ppm. Some properties of the polymers studied are given in Table I. Physical properties of all cores used are in Table 2. Flow behavior of polymer solutions was studied by three consecutive flow tests in cores. First, water (or brine) was injected at constant rate of about 1 ft/D to obtain the water mobility. Then, filtered polymer solution (prepared in water or brine) was flowed through the core. Since the rate was constant, increase in the pressure drop across the core reflects decrease in the mobility. Finally, the core was flooded with water (or brine) to study recovery of mobility. The Alundum cores which were used in

Jan 1, 1967

By G. W. Thomas

Jan 1, 1967

By D. L. Kelley

High-viscosity water injection has been proposed for use in reservoirs containing high-viscosity crude oils. Previous publications have largely ignored the possible effects of the connate water on the proposed process. This paper describes experimental work which indicates that the connate water will be forced ahead of the injected water to form a bank of low-viscosity water. This decreases the oil recovery which would be expected if such a bank were not formed. These effects are shown for a range of fluid mobilities and connate-water saturations for a five-spot injection system. In general, oil recoveries using viscous water are significantly greater than for untreated water even though they are less than would be expected if no connate water bank were formed. INTRODUCTION The effect of mobility ratio on the oil recovery of wa-terfloods has been known for many years. Muskat first pointed out that the fluid mobilities (k/µ) in the oil and water regions would affect the performance of the water-flood, and he estimated the general effect of these variables.' Since this early work, studies of the effect of mobility ratio on secondary recovery have been reported where mathematical,' potentiometric3 and scaled flow models' were used. These studies have shown that a reduction in the mobility ratio between the oil and the displacing fluid would cause additional oil recovery when water-flooding reservoirs containing viscous crude oils. Studies reported by Pye- nd Sandiford 8 have indicated that chemicals to increase injection water viscosity are now available and can be used to reduce the over-all mobility ratio of a waterflood. Where mobility ratios are controlled by the injection of viscous fluids, the connate water of the reservoir can play an important part in the displacement of the reservoir oil. The purpose of this study was to determine the effect of the connate-water saturation in waterfloods where viscous waters are used for injection. DISPLACEMENT OF THE CONNATE WATER Russell, Morgan and Muskat7 were the first to recognize the mobility of connate waters in waterflooding. They conducted waterfloods on oil-saturated cores containing 20 and 35 per cent irreducible water saturations, and found that from 80 to 90 per cent of the "irreducible" water was produced after only one pore volume of water was injected. However, their experiments were conducted at rates of flow significantly higher than those ordinarily occurring in waterfloods. Also, the cores were only from 4.0 to 8.5 cm long. Brown 4 studied a 100-cm linear sand pack which had been prepared to contain connate water and oil. He used 140- and 1.8-cp oils with injection water of essentially the same viscosity as the connate water. He found that all of the connate water was displaced by the injection water in both cases. However, the injection volumes required for complete displacement of the connate water were considerably higher in the case of the more viscous oil. To verify the results of the foregoing experiment, a 10-ft-long linear model was constructed by packing 250-300 mesh sand in a 1/2-in. diameter nylon tube. The model was evacuated, saturated with a brine of 1-cp viscosity, and flooded with a 41-cp mineral oil to the irreducible water saturation of 10.9 per cent. The model was then waterflooded by the injection of a water solution which had an apparent viscosity of 42.6 cp. The solution consisted of 0.5 per cent methylcellulose in distilled water. The viscosities of the oil and connate water were measured with an Ostwald viscosimeter. The viscosity of the polymer solution was calculated by Darcy's law using pressures measured during actual flow conditions. The ratio of the mobility in the oil region to the mobility in the inject ion-water region was approximately 0.32. The mobility ratio of the oil region to the connate-water bank was approximately 14. The mobility ratio between the connate-water bank and the injection water region was 0.024. Approximately 84.5 per cent of the recoverable oil was produced before water breakthrough. Immediately following breakthrough, oil and connate water were produced at an increasing water-oil ratio until the viscous injection water broke through. At viscous-water breakthrough, 96 per cent of the original connate water had been produced. After breakthrough of the viscous water, there was no additional production of connate water or oil. The near-

Jan 1, 1967

By S. S. Marsden, S. A. Khan

Externally generated foam was injected continuously into short porous media. Both flow rate and pressure drop were measured. Liquid saturation was determined by electrical conductivity. Foam yuality I, expressed as the ratio of gas volume to total volume, was varied from 0.70 to 0.96. As measured with a modified fann VG Meter, apparent viscosity of this foam µa decreases with increasing shear rate but usually falls within the range of 50 to 500 cp. At a given shear rate, µa increases almost linearly with I-. When measured with a Bendix Ultraviscoson, kinematic pa is independent of V but absolute pa increases with T from about 3 to 8 cp. The effective permeability -apparent viscosity ratio ke/µa decreases almost linearly with V for porous media of high permeability, but the rate of decrease becomes less for tighter ones. The relative permeability-apparent viscosity ratio kr/µa us V data does not fall on a single line. The kr/µa ratio increases with liquid saturation in the porous medium and with surfactant concentration. Estimates of pa for foam in porous media vary from 30 to 100 cp. INTRODUCTION Although research on the development of a foam-drive, oil recovery process has been going on for almost a decade, most of the significant publications have appeared within the last several years. This illustrates well the rate at which interest in this process is accelerating. Bond and Holbrook 1 were the first to describe the use of foam to improve oil recovery in their patent of 1958. They proposed that an aqueous foaming agent slug be injected into the formation and that this be followed by gas to produce a foam in situ. Fried 2 studied the injection of foam into porous media which had already been subjected to conventional gas or water drives and found that gas could be used to drive a foam bank which would, in turn, displace additional oil in the form of an oil bank. He attributed the increased oil recovery to the high effective viscosity of foam flowing in porous media. His microscopic observations showed the importance of foam generation and regeneration within the porous medium. By injecting both air and aqueous surfactant solution, Bernard3 generated foams within the porous medium in which oil displacement was being studied. In a separate empirical test, he also measured the dynamic foaming characteristics of the same surfactants in water and/or oil. With some exceptions and for the seyen surfactants studied, there seems to be a qualitative relationship between the efficiency of liquid rlisplacrment and the dynamic foaming testtsed. This relationship was not consistent enough to eliminate the necessity of actual foam flood tests in porous media for surfactant selection. In a study basic to gas storage in aquifers, Bennett4 described the displacement of brine by foam in consolidated porous media. Among other things, he stated that the ability of a surfactant solution to foam is more important than the stability of its foam. The presence of a foam bank between the displacing air and the displaced brine improved both breakthrough and ultimate recovery. In a continuation of this work Kolb5 attributed the great reduction in surfactant solution production rate as displacement by air progressed to a decrease in relative permeability to gas. These several effects reported by both Bennett and Kolb can all be attributed to the high apparent viscosity of foam which was obviously flowing in the porous media. Bernard and Holm6 observed that gas permeability in the presence of foam was much less than in the absence of foam when both were measured at the same gas saturations. This reduction in permeability was proportionately greater for more permeable sands and sandstones than for the less permeable. Viscosities of the surfactant solution and air were used for their calculations rather than apparent foam viscosities. The foam stability increased with decreasing permeability, which was probably due to the natural regenerative properties of porous media, and the permeability reduction increased

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

By D. L. Archer, W. W. Owens

Conventional waterflooding often is uneconomic in highly fractured reservoirs because of the gross bypassing of the reservoir oil by injected water. Imbibition and pressure pulse flooding have been used in at least one fractured reservoir in an attempt to achieve better oil recovery performance. This paper presents the results of laboratory flow tests conducted on large cores to evaluate the possible applicability of these methods (particularly pressure pulse flooding) to different types of reservoir systems. Test data were obtained on both water-wet and oil-wet systems, and on systems having two widely different levels of compressibility and flow capacity. Results indicate that pressure pulse recoveries from fractured reservoirs will likely not exceed 5 to 10 per cent pore space with maximum response achieved during the first pulse cycle. Improved recovery by this method is possible from both oil-wet and water-wet reservoirs. Comparable saturation distributions during imbibition and pressure pulse production suggest that an initial pressure pulse cycle to speed production response would not interfere with subsequent imbibition flooding in water-wet reservoirs. INTRODUCTION There is a steady increase in the number of fluid injection projects being initiated each year to improve oil and gas recovery over that obtained by primary production mechanisms. Experience gained in different geographical areas and with different recovery methods is teaching that reservoir anatomy is one of the more important factors controlling the success or failure of such projects. Fractured formations appear to be the rule rather than the exception, especially in carbonate reservoirs. Conventional methods of fluid injection, whether it is waterflooding, gas injection or miscible flooding, have limited applicability to highly fractured reservoirs because of the severe bypassing of reservoir fluids. The result is early breakthrough of the injected fluid and rapidly increasing ratios of injected to in-place fluid with an undesirable effect on return on investment. Thus, innovations to conventional methods must be developed if recovery from fractured reservoirs is to be optimized. One new method proposed for application to highly fractured reservoirs is pressure pulse waterflooding. Although pressure pulsing with gas has been suggested' and rested in at least one reservoir; the Sohio Oil Co. is generally credited with the development and testing of pressure pulsing in conjunction with waterflooding. Publications" of Sohio's waterflood operations in the Spra-berry Trend indicate that the idea for pressure pulse flooding evolved from a critical analysis of the disappointing early performance of their flood initiated in a portion of the trend in April, 1961. This flood was planned to take advtantage of imbibiltion and high pressure gradients across the reservoir matrix blocks (unfractured Mocks of the reservoir rock which are visualized to be generally surrounded 'by the fracture system) evaluated in pilot tests by the Atlantic Richfield Co.5 and Humble Oil & Refining Co.' However, the early recovery and pressure performance of this flood indicalted that high pressure gradients induced between the fracture system and the centers of the matrix blocks were forcing water into the periphery of the blocks and temporarily interfering with the counter-current flow of oil due to imbibition. Fill-up in the blocks was occurring with re-solution of the free gas phase as pressures climbed above the original reservoir bubble point. It appeared that cessation of water injection to permit the capillary forces to become dominant and that expansion of the rock and its contained fluids during pressure reduction might aid in expulsion of oil from the rock matrix into the fractures. Subsequent production performance. aflter cessation of injection, proved this hypotheslis correct. Thus, pressure pulse waterflooding was born. Pressure pulse flooding appears to have several advan'tages over imbibition type flooding in highly fractured reservoirs. First, all wells may be used for water injection which should hasten fill-up and achieve a more rapid increase in reservoir pressure. Second, because of increased injection pressures and rates, flooding gradients will force water into the reservoir matrix more rapidly than would be achieved by imbibition alone. Third, during the pressure depletion cycle of the flood, the compressibility of the system, which has been increased by resolution of the free gas phase, provides energy for the displacement of oil at a rate greater than would correspond to countercurrent flow during imbibition. Fourth, during depletion all wells may be put on production, thus contributing to higher oil withdrawal rates on a reservoir-wide basis. One possible disadvantage of pressure pulse flooding as compared to imbibition floodling, however, is that the outer periphery of each matrix block is flooded to a residual or near residual oil saturation during the Water lnjection stage. This zone of reduced permeability to oil m'ay offer greater restriction to the production of

Jan 1, 1967

By J. Lohrenz, R. J. Francis, G. C. Clark

Compositional material balance methods have been extended to solution gas-gas cap-water drive reservoirs. Partial or complete pressure maintenance as well as reservoir repressuring by gas and/or water injection may be considered. Increrments of reservoir withdrawals and injections equivalent to a given time increment are used as the basic iterative step to allow the extension. An example field study shows how the reservoir volumes. extent of the water drive, and relative permeability ratio function may be estimated. Knowing these, the methods can be used to predict reservoir behavior. INTRODUCTION The compositional material balance is similar to conventional material balances in many aspects. The assumption of a pressure-homogenous reservoir with no transients is common to all conventional material balances1-1 although numerical methods are different. The compositional material balance also contains this assumption. Inclusion of phase equilibria directly in reservoir calculations is the key difference and improvement in the compositional material balance. A continuous account of the compositions of the reservoir fluids and effluents is determined. Calculations of fluid phases and properties can be based on these compositions. Excepting fluid phases and properties determined in this manner, the compositional material balance uses procedures exactly analogous to conventional methods. Brinkman and Weinaug5 and Jacoby and Berry6 have given compositional material balances for simple solution gas drive reservoirs; subsequently Jacoby and Berry' extended their procedure to include pressure maintenance by gas injection at constant reservoir pressure. This paper extends compositional material balance techniques to a more general combination drive reservoir with various types of pressure maintenance operations. The procedure for this extended compositional material balance and an example of its field application are presented. PROCEDURE The numerical calculations of this compositional mate- rial balance are performed on a high-speed digital computer. The most comprehensive case the present compositional material balance program can handle is a solution gas-gas cap-water drive reservoir with (1) water production, (2) water influx, (3) separate or commingled production from the oil and gas zones, (4) any surface processing arrangement including LPG and gasoline production, and (5) gas injection into the gas and/or oil zones." This description of the procedure applies to this most comprehensive case. Appropriate deletions in the procedure can be made for more simple mechanisms. ASSUMPTIONS In common with conventional material balances, the compositional material balance makes the following assumptions: (1) the total reservoir is at a uniform pressure at a given time; (2) the volumetric ratio of gas to liquid withdrawal at reservoir conditions is directly proportional to their mobilities; and (3) instead of assuming a pressure-dependent total gas-oil equilibrium, the cornpbsitional material balance makes a sianificantly different and more rigorous assumption as follows. Thermodynamic equilibrium exists for each component in the liquid-gas mixtures in the reservoir and in surface processing. These three assumptions are fundamental in the compositional material balance procedure. Other assumptions involve arbitrary simplifications of the reservoir mechanisms. This compositional material balance uses the following simplifying assumptions: (1) no bulk mixing occurs between the gas and oil zone; (2) total mixing occurs between any injected gas and the contents of the injection zone; (3) no in situ hydrocarbon liquid is produced from the gas zone; (4) no reservoir gas is retained behind any net bottom-water influx; (5) the residual saturation of hydrocarbon liquid at reservoir conditions behind any net bottom-water influx is constant; and (6) fluids retained behind any net bottom-water influx exert no further influence on the reservoir. These six assumptions are not fundamental. Modifications of the present procedure can be made for different assumptions of the reservoir mechanisms. THE TIME INCREMENT Previous compositional material balances have used a pressure increment for simple solution gas drive reservoirs"" and an increment of injected gas for pressure maintenance.' These increments are suitable for these simple cases. They are, however, inadequate for reservoirs with time dependent drive mechanisms including water drive. A time increment removes this inadequacy and is