Search Documents

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 N. Mungan

Laboratory imbibition and displacement experiments were performed using crude oil and cores drilled with water and preserved under anaerobic conditions. The purpose of these tests was to determine reservoir rock wettability and to find out if more oil could be recovered by use of NaOH solution than by conventional waterflooding. The preserved cores were found to be oil-wet. Contrary to work in the literature, these cores changed to water-wet upon contact with air. After exposure to air for a week, the cores yielded more oil by waterflooding than when preserved under exclusion of air. At reservoir temperature of 160F, flooding the preserved cores with 0.5N NaOH solution recovered more oil than an ordinary wa-terflood, and additional oil when following a waterflood. When the caustic solution was used from the beginning, all the extra oil was obtained before breakthrough; when the caustic followed a conventional waterflood, the extra oil was produced in the form of an oil bank ahead of the injected caustic. The increase in oil recovery resulted from wettability reversal. Also, use of caustic reduced the volume of injection required to flood out the cores. At room temperature, however, the caustic solution did not reverse the wettability and gave no additional oil recovery. Cores which had become water-wet by air exposure or caustic flooding were restored to their original oil-wet state when saturated with crude oil and allowed to equcilibrate at reservoir temperature for two weeks. Therefore, in the absence of preserved cores, it may be possible to restore weathered cores to their original wettability for use in laboratory floods. INTRODUCTION Waterflooding has been in use since 1865, and is by far the simplest of secondary recovery methods. Unfortunately, most waterfloods are inefficient in recovering oil, often leaving half or more of the original oil in place un-recovered. The low oil recovery generally results from low sweep efficiency and low displacement efficiency. Consequently, to increase oil recovery by waterflooding, sweep and displacement efficiencies should be improved. Sweep efficiency is primarily affected by reservoir heterogeneities and mobility ratio, while displacement efficiency is affected by the capillary forces between fluids and rock surfaces. For petroleum reservoirs, the capillary forces are expressed in terms of interfacial tension and wettability. If oil recovery is to be improved significantly in water- flooding, the capillary forces holding the oil in the raervoir porous matrix must be reduced or eliminated. One way to reduce capillary forces is to inject commercial surfactants ahead of the injection water into the reservoir. Laboratory tests of this method have shown no promise of an economical process yet, and no increase in oil recovery was obtained in the field trials which have been reported. Work is continuing in many companies to find surface-active agents which, in workable concentrations, can yield substantial added oil recovery. Another way to change capillary forces operating in petroleum reservoirs is by changing the pH of the injected water. Wagner et al.' showed that change in the pH sometimes activates the surface-active materials natural to some crudes and brings about gross wettability change. Since pH alteration can be obtained with cheap chemicals, such as hydrochloric acid or sodium hydroxide, the process shows promise of being economical in a field application. Pan American Oil Corp. reported oil recovery by use of caustic solution from a flooded-out reservoir.' Their test, conducted at a small additional cost, yielded results which were so sufficiently favorable and encouraging that the wettability reversal flood was expanded to portions of the field not previously flooded.13 It is important to bear in mind that changes in the pH of the water not only can reverse wettability but also can lower the interfacial tension between water and crude oil. Reisberg and Doscher4 have studied the pH dependency of the interfacial tension of Venture crude using sodium hydroxide solutions of various concentrations. Their data show that the interfacial tension was lowered from 23.0 to 0.02 dynes/cm by increasing the NaOH concentration from 0.005 to 0.5 per cent by weight. Thus, the use of NaOH may lead to additional oil recovery due to both wettability reversal and lowering of interfacial tension. Whether alteration of pH results in wettability reversal from oil-wet to water-wet and increases oil recovery depends on wetting properties of the reservoir rock and the crude. This necessitates delicate laboratory experiments, with suitable core and fluid samples from a field. Although many investigators have studied wettability reversal floods in the laboratory,1,2,5,6 these studies have been carried out with synthetic porous media, refined laboratory fluids and surface-active chemicals to simulate the process. The study presented in this paper is the first time that wettability reversal by pH alteration has been accomolished in laboratory core floods using carefully preserved natural cores, live crude and with experiments performed at reservoir pressure and temperature.

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 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 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. 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 E. J. Koval

Practical miscible displacement processes will be characterized by lingering of the solvent into the oil. The lingering process is brought on by viscosity differences, and can be accentuated by channeling and longitudinal dispersion. The effects of these factors on the efficiency of unstable completely miscible displacements are combined in what is called the K-factor method. This method, analogous to the Buckley-Leverett method, predicts recovery and solvent cut as a function of pore volumes of solvent injected. Experimental data are included and show excellent agreement with theory for a wide variety of sandstone cores and viscosity ratios. INTRODUCTION Theoretical considerations,1 laboratory experiments,2 and pilot tests3 lead to the conclusion that miscible displacements in the field will be unstable In an unstable miscible displacement, the solvent fingers through the oil. This fingering leads to early breakthrough of the solvent and an extended period during which both oil and solvent are produced. For such a system, there appear to be four principal factors which bring about or accentuate the effects of instability: longitudinal dispersion (including geometrical effects), channeling, viscosity differences and gravity differences. Other factors, such as diffusion and flooding rate, can also influence the effects of instability, but at the flooding rates considered in this report, circa 15 ft/D, they are unimportant. Longitudinal dispersion can be thought of as a spreading of the solvent front caused by the presence of microscopic inhomogeneities. Channeling of the solvent occurs when a porous medium has macroscopic inhomogeneities; i.e., gross permeability variations. Viscosity differences lead to fingering of the less viscous solvent. This difference in viscosity accelerates the growth of fingers along paths previously developed because of permeability variations. Gravity differences lead to overriding by the usually less dense solvent. Al- though gravity effects are generally small at a flooding rate of 15 ft/D, they would, nevertheless, unnecessarily complicate the interpretation of flooding experiments. For this reason, all the experiments reported herein were done with matched density fluids. Fingering and the resultant poor areal sweep were recognized early as the dominant influences on the efficiency and the economics of miscible displacement processes. Much research effort has been spent on ways to minimize fingering and increase areal sweep such as the use of graded viscosity slugs4 or water slugs. Some researchers attempted to work out ways to prevent fingering completely; i.e., to achieve a stable displacement through gravity or latitudinal dispersion stabilization. l Others did not attempt to control fingering but obtained an economic process by merely recycling the solvent and sweeping pattern by pattern.6 During this period, all aspects of fingering came under close scrutiny. Some researchers reported on how fingering looks2 and how it is affected by viscosity ratio,7 geometry,7 and slug size.2,* Peaceman and Rachfordg suggested a mathematical approach to the prediction of unstable miscible displacements in relatively homogeneous sand packs, but their work cannot be extended conveniently to heterogeneous systems. Hence for a heterogeneous system, no method is presently available for predicting solvent cut and recovery as functions of pore volumes of solvent injected. The purpose of this investigation was to attempt to fill in the gap in our knowledge concerning the prediction of performance of unstable miscible dis placements. Necessarily, the system selected for study was a relatively simple one. The restrictions placed on the system were: (1) The system was linear; (2) The solvent was miscible in all proportions with the oil in place; (3) The solvent was continuously injected into the porous medium; (4) Gravitational effects were eliminated by matching densities and (5) All the flood rates were high and constant at 15 ft/D to avoid any small rate effect and to minimize any diffusion effects. To simplify and to indicate that both longitudinal dispersion and channeling arise from permeability variations, the effects which they cause or influence have been termed heterogeneity effects. The use of

By T. K. Perkins, O. C. Johnston

Because of the influence of dispersion on miscible-displacement processes, diffusion and dispersion phenomena in parous rocks are of current interest in the oil industry. This paper reviews and summarizes a great deal of pertinent information from the literature. Porous media (both unconsolidated packs and consolidated rocks) can be visualized as a network of flow chambers, having random size and flow conductivity, connected together by openings of smaller size. In such a porous medium, the apparent diffusion coefficient D is less than the molecular diffusion coefficient Do, as measured in the absence of a porous medium. For packs of unconsolidated granular material the ratio D/D, is about 0.6 to 0.7. For all porous rocks, both cemented and unconsolidated, the ratio of diffusion coefficients can also be represented as D/Do = 1/FF where F is the for- mation electrical resistivity' factor and 4 is the porosity. If fluids are flowing through the porous medium, dispersion may be greater than that due to diffusion alone. At moderate flow rates the porous medium will create a slightly asymmetrical mix zone (trailing edge stretched out), with the longitudinal dispersion coefficient approximately proportional to the first power of average fluid velocity (if composition is nearly equalized in pore spaces by diffusion). If the velocity in interstices is large enough, there will be insufficient time for diffusion to equalize concentration within pore spaces. In this region, longitudinal dispersion increases more rapidly than fluid velocity. At low velocities in interstices, transverse dispersion is characterized by a region in which transverse diffusion dominates. If the fluid velocity gets high enough, there will be a transition into a region where there is stream splitting with mass transfer but with insufficient residence time to completely damp-out concentration variations within pore spaces. There are several variables that must be con- trolled to get consistent longitudinal and transverse dispersion results, viz., (1) edge effect in packed tubes, (2) particle size distribution, (3) particle shape, (4) packing or permeability heterogeneities, (5) viscosity ratios, (6) gravity forces, (7) amount of turbulence, and (8) effect of an immobile phase. INTRODUCTION Diffusion and dispersion in porous rocks are of current interest to the oil industry. This interest arises because of the influence of dispersion on miscible-displacement processes. In a recovery process utilizing a zone of miscible fluid, there is the possibility of losing miscibility by dissipating the miscible fluid or by channeling or "fingering" through the miscible zone. Diffusion and dispersion are two of the mechanisms that may lead to mixing and dissipation of the slug. On the other hand, dispersion may tend to damp-out viscous fingers which may be channeling through the miscible slug. 58 Hence, dispersion may be detrimental or beneficial (if it prevents fingering through the miscible zone). Therefore, it is doubly important that we understand these processes. In this paper we review, summarize and interpret a great deal of information from the literature. In particular, we will briefly discuss molecular diffusion in miscible fluids. Then we will discuss what differences to expect for diffusion in a porous rock. If there is movement of the fluid through the rock, then there may be an additional mixing or "dispersion". Furthermore, the dispersion longitudinally (in the direction of gross fluid movement) and transverse to the direction of fluid movement will not be equal. We will discuss both types of dispersion as well as several variables which can affect dispersion (viscosity differences, density differences, turbulence, heterogeneity of media, etc.). This group of variables has sometimes led to difficulty when comparing literature data. DIFFUSION OF MISCIBLE FLUIDS If two miscible fluids are in contact, with an initially sharp interface, they will slowly diffuse into one another. As time passes, the sharp interface between the two fluids will become a diffuse

By G. W. Thomas

A mathematical tnodel of forward combustion in an oil reservoir is treated in this paper. The model describes a radial system having a vertical section of essentially infinite thickness, all of which is permeable to gas flow. Combustion, however, is presumed initiated over a limited thickness of the total vertical section. In the interval supporting cotnbustion, the mechanisms of radial conduction, convection and heat generation are taken into account. Above and below the burning interval, heat transport in the radial direction is by cottduction and convection. Vertical heat losses from the ignited interval are accounted for by conduction alone. A general solution is presented for the temperature distribution caused by radial movement of the combustion front. The results show that no feedback of heat occurs into the ignited interval when convection and conduction are acting in the bounding media. Peak temperatures are also 5 to 10 per cent higher than in the case where heat transport in the bounding media is by conduction alone. We arbitrarily define vertical coverage to be that fraction of the total ignited interval which is at 600F above atnbient, or greater, at any given time. The radial distance at which the vertical coverage becomes zero is the propagation range of the combustion front. It was found that an increase in vertical coverage results when the oxygen concentration, fuel concentration or gas-injection rate is increased. Moreover, the combustion front can be propagated 10 to 15 per cent further than in the case where only conduction is acting above and below the ignited interval. INTRODUCTION In the theoretical treatment of forward combustion in a radial system, one of the problems encountered is the determination of the transient temperature distributions caused by an expanding cylindrical heat source. Bailey and Larkin' and Ramey' simultaneously presented analytical solutions to the problem assuming heat transport by conduction alone. In a subsequent publication, Bailey and Larkin3 included the effects of both conduction and convection while treating linear and radial models. In this latter work, however, vertical heat losses were largely neglected. Selig and Couch' dealt with a radial model in which both conduction and convection were acting. Only a limiting case involving vertical heat losses was considered, however. Namely, temperatures on the boundary of the bed of interest were set equal to zero. Solutions thus obtained were representative of a system having a maximum vertical heat flux. Chu5 recently treated a more general case in which a permeable bed was considered bounded by impermeable media. Conduction and convection took place within the bed, and only conduction outside of the bed. The effects of vertical heat losses were included in his study. Solutions were obtained by numerical techniques. This paper is an extension of the theoretical work of other authors pertaining to forward combustion in a radial system. In particular, a mathematical model of the process is treated in which heat generation occurs over a small vertical interval of a larger permeable section. In the interval supporting heat generation, and above and below this interval, the mechanisms of radial conduction and convection are also presumed acting. Heat losses from the ignited interval are accounted for by vertical conduction. An analytical solution for the temperature distribution caused by radial movement of the burning front is presented. The effects of certain process variables are indicated and comparisons with Chu's results are made. THEORY To render the mechanism of forward combustion tractable to mathematical treatment, we idealize the problem to the extent of assuming continuous reservoir media possessing homogeneous and isotropic properties. The following additional assumptions are implicit in this analysis. 1. The thermal parameters, i.e., heat capacities, thermal conductivities and thermal diffusivities are invariant with temperature and pressure. Moreover, the bounding media possess the same thermal properties as the bed of interest. 2. The temperatures of the porous media and its contained fluids at any point and at any time are equal. 3. The reaction rate between the oxidant gas and the fuel is infinite. This assumption implies that the incoming oxygen concentration instantaneously goes to zero within an infinitesimal distance, i.e., the width of the combustion zone is negligible. 4. The rate of gas injection is constant and corresponds to the average rate throughout the lifetime of the project. 5. The fuel concentration is constant throughout the volume of rock swept out by the burning zone. 6. There is complete burnoff of fuel. This assumption demands that the rate of propagation of the burning front equals the rate of fuel burnoff. In a radial system, with a

By P. M. Blair, D. W. Peaceman

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

By L. T. Stanley

A method is described whereby the behavior of the front is approximated when gas is injected into a thick reservoir having reasonably homogeneous properties. The method is applied to the Abqaiq field of Saudi Arabia. The producing member is the Arab-D, a high-permeability calcarenite which has no continuous betlding planes. Field performance indicates the injected gas to be partially miscible with the reservoir oil because gas volume factors experienced are lower than those predicted by the gas law. Eqlrilibril~m calculations confirnz ihis and also show that the volume factor may vary appreciably, depending on the relative quantities of oil and gas in contact in the reservoir; consequently, volutnetric balances are made for the purpose of determining the behavior of the gas volume factor during injection history. The effect of gravity segregation on position and shape of the gas front is approximated by dividing the producing section into three zones vertically and tracing the movement of the average interface in each zone, with counterflow of oil and gas taking place betweell zones behind the front. This results in a gas front that advances more rapidly along the top of [he formation than along the base, producing an "umbrella" effect in cross section that becomes more pronounced as the front progresses. The calculations cover the gas-injection period of Abqaiq field history through 1958, and the position of the gas front is plotted at one-year intervals. Finally, some comparisotls between the calcrlated fronial position and the gas-oil contact as Measured by neutron logs are made. INTRODUCTION Fig. 1 is a structure contour map of the Abqaiq field. For analytical purpcses, the field has been divided into two areas designated A and B, as shown on the map. Gas injection into Area A was begun in early 1954 into two wells located near the crest of the structure on the north-south axis of the anticline. During the period covered by this paper, the average injection rate has been 122 MMscf/D with a maximum sustained rate of 160 MMscf/D. The reservoir pressure has been maintained at 2,450 psig since start of gas injection. The reservoir originally was undersaturated but, prior to pressure maintenance, the pressure was drawn below the bubble point in most of the Area A section, creating an in-place free-gas saturation of 8 per cent at the crest of the structure. The average critical gas saturation is 13 per cent established by laboratcry core tests. That the entire reservoir remained below critical gas saturation prior to gas injection is borne out by the fact that none of the crestal wells have exhibited high gas-oil ratios during production history. The Area A reservoir is a carbonate section about 200-ft thick having average porosity and permeability of the order of 22 per cent and 500 md, respectively. The productive horizon is the Arab-D member of Jurassic age and is encountered at an average depth of

By R. J. Wygal

A method is presented by which homogeneous, isotropic, and reproducible aggregates of particles such as glass beads, sand, ground glass, etc., may be prepared simply and quickly by using a "particle distributor". Several applications of the packing method to the construction of model reservoirs are described. The flow characteristics of these packings me given. Some factors which limit the usefulness of unconsolidated packs for simulating reservoir rock are overcome through artificial cementation of the aggregates. INTRODUCTION Many laboratory studies of porous media are based on experiments using test packs of sand, glass beads and other granular materials. More confidence can be had in the results of such studies if the test packs are known to have uniform properties throughout, to be mechanically stable and to be accurately reproducible. This paper presents a method by which such packs of a variety of shapes and sizes may be constructed using a particle distributor. Some factors which limit the usefulness of unconsolidated packs as a means for simulating reservoir rock are described. A method of overcoming these limitations, which involves artificial cementation of the packs, is presented. DEVELOPMENT Particles packed by the tamping or vibrating techniques are usually either non-uniform or mechanically unstable. An invention by Currie and Gregoryl overcomes these difficulties to some extent. The invention is based on the use of a stream of falling particles impinging on and passing through a distributor matrix of several layers of closely spaced spheres. The particles are well dispersed by numerous collisions with the spheres. They become homogeneously distributed at some distance below the distributor matrix and at that position are collected in the sample container. This apparatus, however, was found to be ineffective when small particles were used since small particles accumulate on the top surfaces of the spheres. An improved distributor matrix was made, first in the form of a stack of five sieves, and then in the more adaptable form of spaced layers of wire screen. Fig. 1 shows a longitudinal view through a five-layered distributor of eight-mesh wire screen. This device gives excellent results.* Fine particles pass through the distributor with very little hold up. The fine structure of the wire matrix with its many deflecting surfaces ensures good particle distribution. There are no dead areas immediately below the screen as were sometimes observed to occur below a matrix of spheres. Using such a device, 20-30 mesh sands were packed to a porosity of 32.1 per cent at rates varying from 0.17 to 6 lb/rnin. The lowest porosity achieved by tamping small increments of similar sands into a container was 33 per cent. Motion pictures at 3,000 frarnes/sec of the surface of a growing pack of 2-3/4 mm glass beads have helped to show the mechanism by which packs of low porosity are obtained. The beads, although packed at the rate of 5 lb/min, were observed to hit the pack surface singly, for the most part, and at well spaced points of impact. Part of their energy is transferred to surface beads which are thereby

By D. N. Dietz

A method for determining average reservoir pressure is presented, which is simpler to apply than that devised by Matthews. Brons and Hazebroek. For bounded reservoirs, identical results are obtained if stabilized-flow conditions prevail. The present method yields inferior results in the transient state. The method can, with a slight modification, also be used for water-drive reservoirs. INTRODUCTION In the method proposed by Matthews, Brons and Haze-broek' for determining average reservoir pressure in a multi-well reservoir, the cumulative production (the production time of each well) enters twice: once when the buildup is plotted against In (t + t)/At to arrive at p*, and a second time when the correction, p*-p is determined with one of the several formulas for differently shaped drainage areas. Once a steady state has been attained, however, the previous production history should be immaterial. The same pressure distribution could have been arrived at after different cumulative productions of the individual wells. In principle, therefore, it should be possible to determine average pressures without referring to cumulative productions. In the following paragraphs an expression is presented for the difference between the pressure in a producing well and the average pressure of its drainage area. Then the build-up time needed to overcome this difference is indicated. CASE OF A CIRCULAR DRAINAGE AREA AND A CENTRAL WELL PRESSURE DISTRIBUTION BEFORE SHUT-IN The general differential equation of radial flow (see Muskat, Eq. 10.2) may be written In the steady state, the rate of production of the well is equal to the rate of expansion of the fluid contained in the the drainage area; thus

Jan 1, 1966

By D. G. Russell

A simple method has been developed with which flowing bottom-hole pressure data from two-rate flow tests in oil or gas wells can be analyzed to estimate the formation permeability, skin factor and average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized producing rate of the well is changed to another, higher or lower, rate. The new method yields the same information as a conventional pressure buildup analysis, but eliminates the need for closing in the well. The analysis of a two-rate flow test is of the same degree of difficulty and requires about the same engineering time for application as a conventional pressure buildup analysrs. The extended closed-in periods experienced with conventional buildups because of long, low-rate afterproduc-tion periods are eliminated by fIow tests. Other anomalous pressure buildup effects, such as "humping" due to weli-bore phase segregation, can be successfully eliminated with the new method. Generally, flow tests of about 24 hours' duration run with conventional pressure measurement equipment are sufficient for interpretation purposes. Field testing of the two-rate flow test method has established that it is a reliable and economical method which can be used in many instances to complement or even replace conventional pressure buildup methods. INTRODUCTION The principal method for estimating formation permeability and well damage, or skin factor, in a producing oil or gas well is the analysis of shut-in bottom-hole pressure buildup data.' This familiar method has been used quite successfully by reservoir engineers for many years. It is based on the solution of the radial flow equation for constant rate conditions, and requires that the well be closed in for a sufficient period of time to obtain a clearly defined linear portion on the plot of observed t + ?t bottomhole pressure vs log t + ?/?t(where At is shut-in A? time, and t is producing time to the instant of shut-in). From the slope of the plot and other normally obtainable data, the permeability, skin factor, and reservoir pressure at infinite shut-in time (if the reservoir were infinite) can be estimated. Over the years several drawbacks have become apparent in the use of conventional shut-in pressure buildups for determining permeability and skin factor. The conventional pressure buildup interpretation theory assumes that a well is closed in at the sand face and that no production into the well occurs after shut-in. In practice, of course, the well is closed in at the surface, and inflow into the well continues until the well fills sufficiently to transmit the effect of closing in to the formation. This adjustment period is commonly referred to as the "after-production" portion of the pressure buildup. In the tight reservoirs, long, low-rate after-production periods frequently occur, and the well must be shut in for several days or, in some instances, even weeks to obtain an interpretable buildup curve.' Obviously, such long shut-in times can cause loss in current income, both from reduced oil production and from the fact that personnel and pressure measurement equipment are occupied with a single well for too long a time. In other cases, even long shut-in periods do not seem to be of much aid in obtaining an interpretable buildup." If there is considerable phase redistribution (liquid fallout or bubble rise) after a well is shut in, then curves with no interpretable portion are often obtained during the buildup. In addition to instances in which wellbore effects cause trouble, there are also cases in which the major objection to use of the closed-in pressure buildup is simply the fact that the well must be shut in. When there is no proration and when the well has limited producing capacity. closing in the well means loss of income. From the foregoing discussion it is apparent that it ib desirable to have an alternative method of obtaining the same information as that derived from a conventional buildup without the need of closing in the well. One possible solution which has been offered for this problem is the use of a bottom-hole shut-in tool' which isolates the major portion of the flow string from the formation face during the buildup. In this paper an alternative method which is frequently successful in avoiding wellbore effects and which does not require the use of special equipment is presented. A new, simple method has been developed with which the flowing bottom-hole pressure data from flow tests in oil or gas wells can be used to estimate permeability, skin factor and the average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized production rate of the well is changed to another, higher or lower rate. The need for closing in the well is eliminated, and pressure measurement periods of only 18 to 24 hours arc usually sufficient, even in tight reservoirs. Thus, the new

By S. C. Miller, R. D. Carter, H. G. Riley

A merhod based on theoretical considerations is presented for estimating srabilized performance of gas wells from data obtained from short-term flow tests. Calculated results are presented for II gas wells evaluated with four-to five-hour flow tests. For seven wells on which comparison war available, satisfactory agreement between measured and calculated performance was obfained. Although the final results are favorable, diflerences in flow capacity k,h obtained from flow tests conducted at different rates were observed. This inconsistency did not seriously affect the calculated results in the test wells. INTRODUCTION In order to negotiate the sale of gas from newly drilled gas wells and to plan otherwise for the orderly exploitation of gas reserves, it is necessary to predict the long-term withdrawal rates physically possible for each well. Since low permeability gas wells do not stabilize except after a lengthy period of flow, the operator is unable, without wasting large amounts of gas, to measure directly the stabilized performance prior to contracting for sale of the gas and connecting the well to a pipeline. This problem is not new, but is one which has recently received considerable attention because of the rising importance of natural gas in the economy and new emphasis on the development of gas reserves. Thus, a method is needed for estimating stabilized performance from information obtained from short-flow tests in which stabilized flow is not attained. Such short tests would be suitable for wells not connected to a pipeline since the tests require only small amounts of time, labor and production of gas. Other requirements of such a method are short calculation time and use of a model of the gas well and its reservoir drainage area which is as consistent as possible with present theoretical, experimental and field knowledge. The method presented herein was developed in an effort to meet these requirements. It has been named the "two-Row method". In addition to the section on the evelopment and description of the method, a description of field application of the method and a discussion of well test results are in- cluded. Finally, there is a section on an additional consideration prompted by some of the test results. The Appendix contains certain details of the derivation of the two-flow method and a sample calculation for Test Well 1. MATHEMATICAL DEVELOPMENT AND DESCRIPTION OF THE METHOD The method is based upon the following assumptions: (1) The system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial. (2) Flow obeys Darcy's law except very near the well where a quadratic flow law is assumed. (3) The system is also characterized by a laminar flow skin factor. The use of a quadratic or turbulent flow term is more fully justified in Refs. 1 and 2. When such a system is opened at the wellbore to constant rate flow from an initial condition of uniform pressure distribution, the radius or drainage quickly passes through the non-Darcy flow and "skin effect" region imrnediately surrounding the well. After the drainage radius or "transient" has passed beyond this region and before it reaches the external boundary, drawdown can be expressed by'." After the drainage radius reaches the external boundary, stabilized flow exists and drawdown is given by'," Eqs. 1 and 2 form the basis for the method. The validity of these equations will now be examined in the light of some well test data previously available. It has been found by actual measurements on gas wells that the relationship between drawdown and flow rate, that is the performance curve, is closely approximated by' Eq. 3 is a straight line on log-log graph paper. Eq. 2 and, for any lixed value of time, Eq. 1 plot as a curve which is slightly concave upward on log-log graph

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey' have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks' have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.' Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method