Search Documents

By John N. Dew, William L. Martin, `

This paper describes the results of a laboratory investigation conducted to obtain data for an evaluation of the in situ combustion process as a method of producing crude oil from reservoirs. Air and fuel requirements, rates of advance, com-bustion temperatures, and coke and fluid distributions are presented. The mechanism of oil recovery by in stiu combustion is discussed. Five crude oils ranging in gravity from 10.9 to 34.2"' API were prodrrced from a semiadiabatic, uncon-.solidnted sand pack by in situ combustion. Experimental conditions were varied over a wide range in order to determine the inter-relation-.ships of process variables. The mini-mum air flux requirement for self-sustained combustion was found to he less than 10 scf/hr-ft'l. The rate of advance of a self-sustained combustion zone was found to be nearly proportional to the air flux at the combustion front. The effects of presrure and injected air flux were studied in a series of experiments using a 21.2" API crude. A minimum air require-ment was observed at an air .flux of 20 scf/hr-ft'. The oil saturation con-sumed as fuel averaged 5.5 per cent of pore volume. The effect of air pressure was, found to be .small for experiments having high combustion eficiencies. This study should promote a better understanding of the problems and mechanisms involved in labora-tory investigations and field applications of the in situ combustion pro-cess. The data presented will be use- ful in the interpretation of results of field tests. When tempered with volu-metric sweep efficiencies, the data can be used in making preliminary economic appraisals of the process as applied to reservoirs containing high porosity unconsolidated media. INTRODUCTION The purpose of this work was to obtain laboratory data for an evaluation of the in situ combustion process as a method of producing crude oil from reservoirs. In situ combustion basically consists of (1) injecting air into a reservoir through selected input wells to create an air sweep through the reservoir, (2) igniting the crude at the injection well, and (3) propagating the combustion front through the reservoir by continued air injection. By this means, oil is swept toward producing wells in the area. The fuel for combustion is supplied by heavy residual material (coke) which has been deposited on the sand grains during distillation and cracking of the crude oil ahead of the combustion front. Recovery of petroleum by a combustion or heat wave process is not a new idea. F. A. Howard was granted a patent in 1923 on a process in which air and a combustible gas were pumped into an injection well and ignited.' Russian engineers reported on field experiments with a crude oil gasification process in 1935'. Other known early field tests include those conducted near Bar-tlesville and Ardmore, Okla., in 1942." More recently completed field tests include an inverted seven-spot by Sinclair Oil Co. in the Delaware-Childers field, Nowata County, Okla.'; a three-well and an inverted five-spot test by Magnolia Petroleum Co. in Jefferson County, Okla."."; and a test by California Research Corp. in the Irvine-Furnace field in Kentucky. The Worthington Corp. has completed a test in cooperation with the Forest Oil Corp. in Clark County, 111.l Several field tests are now in progress. Sinclair has acquired a 600-acre lease in the Humboldt-Chanute field, Allen County, Kans., and is under way with a large-scale operation. Three field tests are under way in California. The General Petroleum Corp. is conducting an inverted five-spot pattern test in the South Bel-ridge field in Kern County under a cooperative agreement with 11 other companies, including Continental Oil Co. California Research Corp. is conducting a four-well pattern test in Midway-Sunset field near Maricopa, and Richfield Oil Corp. is under way with a test in the Ojai field in Ventura County. Although the published information is valuable concerning results of field tests, little laboratory data other than that reported by Kuhn and Koch5 are available to aid in appraising the tests from a technical or economic standpoint. Engineering data needed to evaluate the process include the amount of air required per barrel of oil recovered, minimum injection rates which will support combustion, rates of advance, air required per unit volume of reservoir cleaned, amount of oil recovered, and amount of oil consumed as fuel. The present laboratory investigation was undertaken to provide these data for that portion of the reservoir swept by the combustion zone. The first phase of this investigation was exploratory in nature and was conducted at injection pressures less than 100 psig. Data were ob-

By B. L. Landrum, P. B. Crawford

The rise of a new laboratory model for studying tran-sient fluid flow problems, is described. The theory of he model is based on the analogy between the equa-ions which describe the flow of compressible fluids in porous media and the conduction of heal m solids. The "thermal model" is used to study the Iransient Pressure distributions during fluid displacement in five Standard well patterns. The studies indicate quantita-ively how die pressures change with lane at any poit in the reservoir. It was found that, in general, the equipressure contours are more radial about the wells early in the program than after steady-state conditions are established. In addition two studies of fractured five-spot patterns are presented to illustrate the opplicability of the model for a wide range of reservoir conditions. INTRODUCTION In reviewing the oil recovery techniques involving fluid displacement it was found that a very considerable quantity, if not all of the fluid within the pattern, mav be displaced (during a time-dependent or transient period rather than the usually applied steady-state conditions. Since the oil is displaced during a transient period, it was believed desirable to study the nature of the transients for several of the commonly used fluid displacement patterns. A literature review indicated fairly complete studies of the streamlines and isopotentials for steady-state: yet few data were available on the transient displacement process. Several simulated steady-state multiwell displacement programs were carefully mapped by Foster and Lodge on an electrical conducting sheet as early as 1875, hut the more recent work of Muskat provides the best source for steady-state phenomena having direct petroleum production applications. In order to study the transient phenomena a new lab-oratory model was developed which is based on the similarity of the equations that describe the flow of fluids and the conduction of heal. The equations describing heat flow were developed early in the nineteenth century. Regorous solutions tor most transient problems are in gcneral complex and usually limited to fairly restricted houndary conditions.2 the thermal model may be used to obtain approximate solutions for many transient problems without some of the restrictions required for analytical solution. The application of the analogy between heat conduc-. tion in solids. electrical flow and fluid flow, in a porous media has long been recognised. Techniques originally developed for the solution of heat conduction prol?lems have in many instances proved useful in ohtaining solutions to fluid flow prohlcms. Ir is not surpri5ing then that a heat conduction or thermal model might bf usefull in sludying transient fluid flow. As an example ol the problems suitable for the study with the thermal model, transient pressure distributions during fluid displacement programs are studied for the five-spot, line-drive, staggered line-diibe. seven-spot 211 two fractured patterns. A model was constructed of an iron plate of geometry similar to the reservoir patterns under study. The principle of mirror images as suggested by Bryan was employed in shaping the model. A temperature differential was obtained by soldering two sections of copper tubing to the plate at points analogous to the well location in the field. Hot and cold water was passed through the tubing. The temperature distribution, analogous to the pressure distribution in the field. was obtained a a function of time from thermocouplcs distributed over the plate.

By C. C. Mattax, J. E. Bobek, M. O. Denekas

ABSTRACT investigations in recent years have shown that rock wettability can exert a profound influence on the displacement of oil by water from oil producing reservoirs. Core analyses frequently show oil recoveries from preferentially water-wet rock to be significantly greater than those from preferentially oil-wet rock. Thus, an accurate prediction of water-drive or waterflood oil recovery is dependent on the evaluation of reservoir rock wettability. Laboratory and field studies have been undertaken to examine and develop methods which may be used to measure rock wettahililiy, to investigate factors which may alter the wettability of reservoir rock, and to determine the wettability of specific reservoirs. It has been found that measurement of the rate and volume of spontaneous imbibition of the wet-ting phase by a rock is a reliable and reproducible test for semi-quantitarive determination of preferential rock wettability. Laboratory tests have shown that the wettability of reservoir rock may depend on both the crude oil composition and the rock type. Field and laboratory tests have indicated that coring flulids and core handling techniques can cause significant changes in the wettability of rock surfaces. However, a few fluids, brine in particular, appear not to affect core wettability and may be used when coring to determine reservoir wettability. Core handling and packing procedures have been developed which preserve core wettability during storage and laboratory testing. The wetfability characteristics of six reservoirs have been evaluated by wettabilify tests performed in the field on fresh cores cut with brine or special water-base muds. The tests indicated that all six reservoirs are water wet; however, the de-gree differed from one reservoir to another. INTRODUCTION Although surface properties of reservoir rock and interfacial characteristics of rock, crude oil, and water systems have been discussed for many years, the full significance of these properties to the production of oil is still in the process of evaluation. It is known, however, that one surface characteristic—preferential wettability of the reservoir rock—-is of paramount importance when the reservoir is being produced by water-flood or water-drive mechanisms. The importance of rock wettability has been pointed out by several authors'.',%* who have found that pref- erentially water-wet cores flood more efficiently than oil-wet cores; that is, more oil is recovered from water-wet cores in the early flooding stages than from oil-wet cores. Other authors'," indicate that waterflood oil recovery from cores of intermediate wettability may be greater than that from either strongly water-wet or strongly oil-wet cores. These recent investigations emphasize that the true wettability of a reservoir must be known before the performance of waterflood or water-drive reservoirs can be accurately evaluated. When reservoir wettability is to be determined by core analysis, precautions must be taken to maintain the original wettability which can be altered readily during coring5 nd testing.' In the coring operation the core may be partially or completely penetrated by the drilling fluid which, if it contains surface active materials, may drastically change the wettability of the core. Also, improper core handling during stcrage and testing may change its wettability because of evaporation of fluids or exposure to oxygen or surface active contaminants. One purpose of this paper is to discuss methods for obtaining, preserving, and testing cores without altering their original wettability. In addition, the results of wettability tests on fresh and preserved cores from six reservoirs are presented where special precautions were taken to obtain samples without altering their wettability during the coring operation.

By L. V. Pirela, S. M. Farouq Ali

Some interest has been expressed recently in the application of solvents in conjunction with a thermal drive, such as a steamflood. At least one field project of this type has been reported. This paper makes the first attempt to provide some of the basic information necessary for choosing a solvent for such an application. Results of displacements conducted at elevated temperatures are also discussed. This paper shows that, in the case of alcohol-hydrocarbon-water or brine ternary systems, the system miscibility may increase, decrease, or both increase and decrease within the same system when the temperature increases. The distribution coefficient was consistently found to increase in favor of the oleic phase with an increase in temperature, which is advantageous from the standpoint of oil recovery. Results of solvent displacements in a sandstone core showed that even a relatively small temperature effect in the favorable direction (increase in miscibility and distribution coefficient) can lead to a 50 percent increase in oil recovery. On the other hand, if the system miscibility decreases with the temperature, even though the distribution coefficient increases, the oil recovery at the higher temperature may not increase. INTRODUCTION Miscible displacement has been the subject of many investigations. Miscible-phase solvent flooding, employing solvents such as alcohols that are miscible with both oil and water, has received considerable attention. Investigations conducted by Gatlin and Slobod,l Taber, Kamath and Reed,2 Holm and Csaszar,3 and Farouq Ali and Stah14 have helped in understanding the mechanism of solvent flooding. While the use of alcohols and similar solvents as oil recovery agents is economically unattractive, these solvents nevertheless provide a means of understanding the mechanistic aspects of important miscible-phase oil recovery processes such as the Maraflood process (discussed by Gogarty5). As shown by Taber et a1.,2 the phase behavior of a solvent-oil-water system is of .utmost importance in determining the efficiency of the oil displacement. This paper makes the first attempt to obtain the ternary phase behavior data for nine alcohol-hydrocarbon-water/brine systems at elevated temperatures. Results of core tests at elevated temperatures are also presented which corroborate the effect of temperature on oil recovery, as judged on the basis of phase-behavior data. The use of solvents at elevated temperatures has been suggested earlier by Farouq Ali.6 At least one field test of this type has been reportedO7 With the widespread use of thermal recovery techniques, it is possible in some situations (especially in steamflooding) that the use of a suitable solvent in conjunction with the heat carrier may be economically feasible. This paper attempts to provide some of the basic information needed for judging such feasibility. LOCATION OF THE PLAIT POINT The phase behavior of a typical ternary system consisting of a solvent, hydrocarbon and brine (or water) can be represented by a triangular diagram such as the one shown in Fig. 1 (for temperature TI). Point Pon the binodal curve represents the plait point, being the limit of tie lines such as Y1 Y2 The compositions of the coexisting oleic and aqueous phases are given by points Y2 and Y1, respectively, for any mixture composition along Y1Y2. The characteristics of ternary diagrams have been discussed by Findlay8 and others. Taber et aL3 have shown that the position of the plait point in a particular ternary system plays the decisive role in determining whether or not theoleic phase would be displaced in a continuous manner during an alcohol displacement. The area above the binodal curve also is important, since it represents

Jan 1, 1969

By C. W. Volek, J. E. Chappelear

The injection of a hot liquid into initially cool porous media, saturated with the same liquid and surrounded by two impermeable but heat conducting media (cap and base rock), has been studied both experimentally and theoretically. The temperature dependence of the viscosity was included in the theoretical model, but it was assumed that the specilic heats and densities of the various materials were independent of the temperature. Solutions to the theoretical model were approximated by numerical methods. Both theoretical and experimental results indicate that center-line temperatures are significantly higher than boundary temperatures. Comparison of experimental and theoretical results with a cold/hot viscosity ratio 01 19:1 were in reasonable agreement. Theoretical calculations show that the effect of the temperature dependence of viscosity was very significant at ratios of 100:l to 1000:1, which are typical of those that occur when injecting hot water to flood heavy oil reservoirs. INTRODUCTION We consider the problem of prediction of fluid flow and temperature distribution in an initially cold-fluid-filled reservoir on the injection of the same hot liquid by the use of mathematical and physical models. The results reported are for a two-dimensional rectangular section of the reservoir, as shown in Fig. 1. The injection and withdrawal faces are assumed to be equipotentials for fluid flow. The ultimate purpose of such models would be to predict hot-water injection performance. However, we note that in the work presented here, one of the most significant aspects of the problem — the instabilities resulting from two-phase, water and oil flow — is not included. We will not give a historical review, but refer instead to the paper of Spillette and Nielsen,1 which contains a rather complete bibliography and critical discussion. The physical problem is exactly the same as Spillette and Nielsen, except for certain simplifications in our assumptions. The mathematical details are somewhat different, and we will present the details of our method here. The mass flow equation, which is elliptic in character, is handled by successive overrelaxation. The heat flow equation, which is parabolic in character, is handled by a straightforward explicit approximation. Some difficulty arose in the over-all heat balance due to small errors in the solution of the mass flow equation, and we feel that a different formulation of the heat flow equation would be desirable for future work. In addition to the mathematical solutions for the temperature distributions in a porous medium due to the injection of hot liquids, experimental data are presented to check the validity of these solutions. A schematic diagram of the model is shown in Fig. 2. Certain qualitative physical conclusions were obtained from our numerical and experimental results. These are: 1. Assuming high (infinite) conductivity normal to the bedding plane in the reservoir is a poor approximation, and may lead to overestimates of the total heat losses (to cap and base rock) of as much as 50 percent. 2. More heat is retained in the reservoir (per unit of heat injected) for higher viscosity changes. 3. Any particular temperature isotherm moves more rapidly along the center line for higher viscosity changes. Consequently, the approximation of temperature independent viscosity is not suitable for obtaining quantitatively correct results. MATHEMATICAL MODEL Our mathematical model is that of a reservoir of thickness 2h. The problem is idealized from that of a linear hot-water drive, and we imagine that the input face is sufficiently far away from the injection wells that the stream lines enter it normally. Similarly the production wells are far enough from the outflow face that the flow lines leave it normally. The reservoir and surrounding cap and base rock

Jan 1, 1970

By J. J. Taber

The turbulent flow drag coefficients, or friction factors, have been experimentally determined for the cut-tings normally encountered in drilling operations. The gas law and average drag coefficients for characteristically flat particles (limestones and Shales) and for angular to sub-rounded particles (sandstones) were used to extend Newton's equation to a more useful form. The resulting equations are expressions for the slip velocity of a particle as a function of the particle size and shape, the bottom-hole injection Pressure, and the Bottom-hole temperature. The velocities necessary to lift actual rock cuttings as observed in the laboratory were compared with velocities predicted from the derived equations. Results indicate that there is no single correct circulating velocity and that the commonly used 3,000 ft/min linear velocity may he sufficient to lift only very mall particles. The advantage, in terms of horsepower requirements, of cycling high pressure air to obtain lower compression ratios was shown. INTRODUCTION A promising departure from standard rotary hydraulic drilling is air or gas drilling. In many instances the rate of penetration and bit life increases resulting from this method have been very substantial. Further the application of air drilling in areas of lost circulation or water sensitive pay zones has been highly successful. Purpose For the drilling contractor, the practical question remains, "How much air pressure and volume should I have?' lt is necessary to have a reasonable knowledge of the air velocities and pressures required in air drilling operations for the proper design and most advantageous use of surface equipment and for adequate removal of cuttings from the borehole. Eased on experience, most operators agree that at least 3,000 ft/min linear velocity is necessary for satisfactory hole cleaning.' It was the purpose of this investigation to evaluate the drag coefficients, or friction factors, for the cuttings normally encountered in drilling operations, i.e., sandstones, limestones and shales, and to utilize the values obtained to develop an expression for the minimum velocities and pressures (hence the volumes) necessary to lift these cuttings. THEORY Several investigators have done extensive work on the ability of drilling mud to lift bit cuttings. They conclude that the ability of a drilling mud to lift cuttings is dependent on the density difference between the rock being drilled and the drilling mud, the cutting size and shape, the mud flow constants, and the flow state (laminar, transitional or turbulent) of the mud. However, drilling muds are not true fluids but have a variation of viscosity with velocity in laminar flow and a constant apparent viscosity in turbulent flow. A survey of the literature revealed that no experimental work had been done to determine the velocities necessary to lift actual rock cuttings using a compressible fluid such as air or gas. In 1953, Nicolson6 stated that in air drilling fluid mechanics the particles could be assumed spherical in shape and that ; constant drag coefficient of 0.5 could be used due to the highly turbulent flow of the circulating air. By equating the turbulent resistance on the particle to the gravitational force on the particle, he obtained the expression, ^ = 2.67^y ,.......(1) where V, is slip velocity of spherical particle, ft/sec; D is particle diameter, in.; p, is particle density, 1b/ft8; and p, is fluid density, lb/ft3. However, as stated earlier, most operators seem to agree that adequate hole cleaning can be obtained by circulating a sufficient volume of air to give a linear velocity of 3,000 ft/rnin in the annular space. Basic Equation The velocity with which a solid particle freely falls through a fluid will increase until the accelerating force, gravity, is equal to the resisting forces, buoyancy due to the fluid displaced by the particle and friction due to the relative motion of the particle through the fluid. When the accelerating and resisting forces become equal,

By J. W. Marx, R. H. Langenheim

The authors are to he complimented for a timely presentation of useful information concerning application of heat to oil reservoirs to increase the rate and ultimate recovery of oil. The solution for heated area resulting from constant rate injection of a hot fluid, presented by the authors, and the fluid-flow analogy, solved by Carter,' are particularly interesting in that 110 restriction is placed on the direction of development of the heated area. It is not necessary that the heated area grow radially. Thus the solution could provide useful information for heat injection in any type of well pattern with any specified swept area data. However, it was assumed in the development that the heated region would remain at a constant clevatcd temperature. Vertical heat losses were computed on this basis. Physically. the assumption that the heated region will remain at a constant temperature (independent of distance from the injection point) appears to require that the heat injection medium be stcam, or other condensible gases near the boding point at injection pressure. If a hot liquid below the boiling point or a hot gas considerably above saturation temperature were used, the temperature of the heated region would decrease with distance from the injection well. Lauwerier2 has discussed this problem for a linear system. (See also Jenkins and Aronolsky; and McNiel and Nelson.') Thus, the authors' solution applies more nearly to steam injection than to hot water injection or hot, non-condensable gas injection. Assuming steam as the injection medium, it appears that steam would pass through the constant-temperature heated region and condense upon contact with the unheated sand. The tempcrature in the condensing region would probably be governed by the release of heat of vaporization required for phase equilibrium at the existing pressure distribution. Temperatures in the region would also be influenced by vaporization of hydrocarbon. Thus, it is doubtful that an expression for radial temperature distribution based on heat conduction only (such as the authors' Eq. 9) would have practical significance. As pointed out by Marx and Langenheim, the growth of the heated area under constant rate steam drive will eventually be limited economically by vertical heat loss. The possibility of continuing to move heat economically by cold water injection or by increasing the steam injection rate appears to be a very important consideration. (Movement of liquids away from the injection well should increase injectivity.) Fortunately, the solution presented by the authors may be generalized to thc case of a variable heat injection rate. The general solution is:

By J. Weijdema, D. N. Dietz

In the conventional underground combustion process (dry combustion) much heat is left behind in the swept formation and goes to rva.rte. Econonmy can be improved by heat recuperation through water injection. This is most advantageous if done at the earliest opportunity before much heat is dissiputed to cap and base rock. Water injected simultaneously with the air will flash to superheated steam, which passes through the combustion front together with the nitrogen from the air. A condensation front traveling up to three times as fast as the combustion front drives out the oil. In this type of wet combustion, the water evaporates before it reaches the combustion zone. The evaporation front travels more slowly than the combustion zone. If so much water is injected that the evaporation front overrun the combustion front, combustion in that spot will be quenched and some unburned fuel will be left behind. Air reacts with the oil farther down-stream where steam temperatures occur; at steam temperature, the air reacts rapidly with the oil. Velocity of the combustion front is increased thereby and is governed essentially by the water-injection rate. In the extreme case of high water-injection rate, a short heat wave of constant length is driven through the formation by water injection. Once this wave has been established, no more heat need be generated than that required to make up the heat losses from the short heat wave; a relatively low rate of air injection will suffice. The feasibility of partially quenched combustion has been confirmed in tube experiments. A heat wave at steam temperature is observed. Chemical analyses of flue gas indicate preferential burning of hydrogen while a carbonaceous residue is left in the formation. Introduction A disadvantage of so-called dry in situ combustion is that air-compression costs are rather high. An air consumption of about 400 std cu m/cu m (400 scf/cu ft) of formation swept is an accepted figure. This high consumption is mostly wasted since much heat is left behind in the depleted oil sand. Methods were investigated for recuperating as much as possible of the heat left behind. This paper deals only with basic principles and is confined mainly to one-dimen- sional flow without lateral heat losses; experiments were conducted in relatively narrow, well insulated tubes. If some water is injected with the air, it will turn to superheated steam in an evaporation front, which should travel behind the combustion front. The steam having passed the combustion front causes a steam drive by a condensation front that can travel faster than the combustion front. The latter needs to travel only part of the distance covered by the oil-displacing condensation front, and thus consumes less air. The water-air ratio would seem limited to that at which cold water overruns the combustion. This limitation was deliberately exceeded considerably in theory and experiments. It was found that combustion is then indeed quenched, but only locally. Farther downstream, the oxygen finds residual oil at steam temperature, which is suficiently high to ensure rapid oxidation. Thus, the combustion front uses only part of the available fuel because it is chased through the formation faster than its normal velocity. No heat is left behind. This new process is called "partially quenched combustion". At the upper limit of the water-air ratio, a small heat slug is moved through the formation by the flow of water and steam. Only a small flow of air is needed since it has only to generate sufficient heat to make up for the lateral heat losses of the short heat slug. Theory Although many factors complicate underground combustion, the processes will be presented in their simplest form. For this reason, one-dimensional flow without lateral heat losses is assumed. Heat conduction in the direction of flow also is disregarded. Under these conditions, dry combustion causes very high temperatures. The heat-carrying capacity of the gas stream is small. Heat generated by oxidation of a residual oil saturation is retained in the sand. The available fuel determines the air requirement and the temperature obtained. Accepting the often-mentioned air consumption of 400 std cu m/cu m (400 scf/cu ft) formation, we calculate a temperature of the swept sand of 1,200C (2,192F) (Fig, I). If water is injected at a modest rate with the air, it will flash to superheated steam upon contact with the heated sand. One cu m (35.31 cu ft) of hot formation will evaporate about 0.5 cu m (17.66 cu ft) of water, and thereafter will accommodate (at an estimated 0.80 saturation and an assumed 0.40 porosity) another 0.3 cu m (10.59 cu ft) of water in cold condition. As long as less than 0.5 + 0.3 = 0.8 cu m (28.25 cu ft) of water is injected for every 400 std cu m (14,125 scf) of air (water-

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 K. Leutwyler

The key to realistic casing stress analysis in thermal recovery installations is accurate knowledge of the temperatures involved. Much information leading to prediction of heat losses from tubing string to wellbore has been published. This .same information permits calculation of casing temperature under certain conditions. This paper reviews the basis for most of the published work, the line source solution to the diflusivity equation and examines some of its application criteria in the light of unsteady-state conditions. An alternate finite differences method is presented for correlation, and experimental data pertaining to the heat transfer mechanism in the annulus is reported. In the light of the information presented in the paper, it appears that the assumption of a line source for the tubing-wellbore system is valid even for reasonably short periods of heat flow, .such as "steam soak" operations. The empirically modified equation presented earlier matches computer derived results of the finite differences approach lor periods as short as 50 hours. Previously predicted cas ing temperatures and the associated stress levels are thus verified based on these correlative computations and the reported experimental information. Fall out from the computer output resulting in some studies of the influence of insulating cement properties on casing temperature, and experimentally studied effects of nitrogen placed in the annulus at various pressures., are described. INTRODUCTION Generally speaking, predictions of average casing temperatures have been based on an idealized model involving a centralized tubing string at uniform constant temperature radiating energy towards the casing under steady-state conditions.'. ' Heat is transferred away from the casing to the formation by unsteady-state conduction across thermal barriers consisting of cement and the mud sheath. Heat balances written on the casing permit prediction of instantaneous casing temperatures. During early periods of injection, casing temperature varies appreciably as a function of time; however, for long injection times this variation decays as heat flow from the wellbore to the formation stabilizes, and it becomes virtually negligible as conditions surrounding the wellbore become quasi-steady-state (Fig. 1). Carslaw and Jaeger -escribe a solution to the diffusivity equation for a line source (Fig. 2). A constant supply of heat is introduced into the medium from the line placed into the origin of the x-y coordinate system and parallel to the z-axis. Heat spreads outward in the infinite cylindrical solid and the temperature increase on a line going through point (x, y) parallel to the z-axis at a distance r can be predicted by: The following criterion must be considered: the Carslaw and Jaeger equation is obtained on the basis of a large

Jan 1, 1967

By H. R. Froning, R. O. Leach

In wertability alteration flooding, a chemical agent is rnoved through a reservoir by the flood water to increase oil recovery by decreasing the degree of wetting of the rock by the oil. Substantial amounts of the chemical may be lost during movement through the reservoir. The extent of the loss, and therefore the economics of the process, depends in some cases on factors which are difficult to reproduce in the laboratory. Therefore, a short-duration, low-cost field test method is needed to permit evaluation of chemical requirements under actual field conditions. This paper describes a small scale rest conducted at a single well for measuring chemical requirements, thereby giving a more reliable evaluation of this important factor in the applicability and economics of the process. In the rest a small wafer slug containing the chemical agent and a nonadsorbed tracer is displaced into the reservoir by a known volume of wafer. The well is then placed on producrion. Chemical loss per barrel of pore volume contacted is calculated from [he fractional recoveries of the agent tested and the nonadsorbed tracer. The method has been used to determine within the actual reservoirs the chemical requirements for both a sandstone and a dolomite reservoir. Several chemical agents are potentially available for wettability alteration flooding, although none is universally applicable. For some applications of the method, chemical costs per barrel of additional oil recovered can be substantially less than one dollar. INTRODUCTION Wettability alteration flooding provides a means of increasing oil recovery from reservoirs by decreasing the degree of wetting of the rock by the oil and increasing the displacement efficiency of the flood water. Earlier studies demonstrated a relationship between oil recovery during waterflooding and the degree of wetting of a rock surface by an oil. The application of wettability alteration flooding to the Harrisburg field of Nebraska provided a field test' of this recovery process. Subsequently, additional laboratory and field tests have developed additional procedures for evaluating wettability alteration flooding, and have indicated where the process may be applicable. Applicability of this process to specific reservoirs is determined by a progression of tests to determine sus- ceptibility of the reservoir to alteration of its wettability, to indicate the degree of recovery improvement and to estimate the amount of chemical required to process the reservoir. The economics of applying improved oil recovery processes depends not only upon the degree of improvement in oil recovery achievable by the process but also upon the process costs and the timing of the income and the investment. Emphasis in this paper is on the expenditure aspects of the process. The work reported in this paper indicates that the chemical investments required for wettability alteration flooding are substantial. For evaluating the economics of a potential flooding application it is imperative that a sound estimate of the chemical requirements be made for the reservoir. Generally, true reservoir conditions are not adequately simulated in laboratory chemical propagation tests. Because of wide well spacings, many years might be required to obtain chemical propagation data from conventional pilots or inter-well tests. Consequently, n short-duration, low-cost method is needed to determine chemical requirements in the field. The potential applicability of wettability alteration flooding is discussed, as well as the economics of wettability alteration with respect to the inherent and imposed restrictions on the timing of income and investments. DETERMINATION OF CHEMICAL REQUIREMENTS In the process of moving a chemical bank through reservoir rock, some of the chemical agent lags too far behind the flood front to be effective or is otherwise lost to the reservoir system. The extent to which these losses occur overshadows the reductions in chemical concentration due to diffusion and to mixing with the reservoir fluids. Experience indicates that almost without exception, chemicals which induce a wetting change undergo either sorption reactions or chemical reactions with mineral constituents of the pore surfaces. Other reactions may occur between the added chemical and the reservoir oil and water. Even in limiting consideration to reactions of the relatively inexpensive inorganic salts, bases and acids, the reactions may be exceedingly complex. Reservoir pore surfaces consist of more than silica in sandstone reservoirs, and more than calcite or dolomite in carbonate reservoirs. Many mineral species are present, each exhibiting specific tendencies to react with an injected chemical. The reactions which occur can consume enough of the agent to have an important effect on economics. These reactions can cause a change in pH of the chemical bank, or may remove some of the active chemical by precipitation, adsorption or reaction to form a new chemical which may or may not be effective in

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 B. A. Eaton, R. H. Jacoby

This paper presents the depletion-performance correlations developed using data from PVT studies and well-completion tests of 27 rich gas-condensate and volatile-oil reservoir fluids. The PVT behavior for each fluid war determined experimentally in the laboratory, and the depletion performance of each fluid was then calculated. Well-completion test data consisting of reservoir pressure and temperature, initial total gas-oil ratio, and stock-tank oil gravity were used with the depletion-performance data to develop general performance correlations. Three correlations are presented. The first equation shows the empirical relation of gas in place per barrel of hydrocarbon pore space with reservoir pressure and temperature, and initial total gas-oil ratio. The second equation is the correlation of stock-tank oil in place per barrel of hydrocarbon pore space. The third equation relates stock-tank oil production (to 500-psia depletion pressure) with pressure, temperature and oil gravity. All three correlations are also presented in convenient graphical form. Gas production to 500-psia depletion pressure is approximately 92.6 per cent of the gas in place for each fluid, so no correlation was developed for gas production. INTRODUCTION The rate of discovery of shallow, low-pressure and low-temperature oil reservoirs in the United States has been declining in recent years. For this reason, the search for oil has gone to greater depths where high pressures and temperatures are encountered. A high-pressure and high-temperature reservoir is the environment in which many rich gas condensates and highly volatile oils are found. Because of this trend, the performance prediction of such reservoirs has become an important task of the reservoir engineer. The need and usefulness of a performance correlation based on field data that are available from the usual well-completion tests have previously been established. These earlier attempts have shown that the phase behavior of rich gas-condensate systems and volatile-oil systems can be correlated with fluid properties, pressure and temperature. However, very little information on gas-condensate depletion performance5 was available at the time of these earlier attempts. The availability of more recent data has made possible the development of an improved correlation. The data, as listed in Tables I and 2, are for 25 gas-condensate and two volatile-oil systems that were studied. The problem was to correlate the data for in-place fluids and the production data with the initial total gas-oil ratio, initial reservoir pressure and temperature, and stock-tank oil gravity. These correlations were obtained by multiple regression analysis, whereby the data were fitted to an equation by the method of least squared error. BASIS OF THE CORRELATIONS 1. The data used in the development of these correlations were obtained from depletion-performance calculations6 and well-test or production measurements. 2. All of the independent variables used in the correlations are available from the usual well-test completion data. 3. The reservoir volume considered is 1 bbl of hydrocarbon pore space (HCPS). 4. The gas-oil ratio used is the total standard cubic feet of gas (at 14.7 psia, 60F) coming off the entire separator train per barrel of stock-tank oil. This minimizes the effect of separator conditions which might vary from system to system. 5. The correlated pressures are the actual reservoir pressures which might or might not be the saturation pressure of the reservoir fluid.

Jan 1, 1966

By H. O. Jahns

This paper describes the application of regression analysis for obtaining a two-dimensional areal description of heterogeneous reservoirs from short-term pressure-time data such as that obtained in interference tests. The method replaces the time-consuming trial-and-error procedure commonly used to match field data on an electric analyzer or digital computer with a systematic search which is programmed for a computer. The computer program adjusts the properties of a reservoir model automatically until a least-squares fit is obtained between observed and calculated pressure data. The reservoir is simulated by a single-phase, compressible, two-dimensional model. It is divided into a number of homogeneous blocks whose transmissibility (kh/F) and storage (Fch) values are varied to obtain the least-squares fit. The reliability of these values is determined from their standard deviations and correlation coefficients. Although the method is rigorously applicable to single-phase flow only, multiphase flow can be handled provided saturation changes are small during the test. Possibly the method can also be used to obtain a reservoir description from pressure-production history, but this application is outside the scope of this work. The paper includes, in addition to a description of the numerical procedure, a discussion of some of the problems associated with the method. Rules are given to help in selecting the number of homogeneous blocks and deciding upon their arrangement. The uniqueness of a reservoir description is considered. Finally, the use of the method is illustrated by the interpretation of field data from two interference tests. INTRODUCTION Pressure data from short-term transient tests, such as single-well and interference tests, are widely used to obtain reservoir properties. These tests are usually analyzed by assuming a simple reservoir model; very often, a homogeneous one is used. As a result, analysis of the transient data from each well frequently gives different values for reservoir properties. The problem then arises to combine all these differing results into a more detailed picture of the reservoir. One technique is to simulate the reservoir with a digital computer or with an electrical analyzer and to adjust the reservoir parameters by trial and error until the simulated pressure data are in reasonable agreement with the observed pressure data for all wells. Although this method has been used for both transient tests and pressure-history data, it is time-consuming and subjective. A second technique uses regression analysis to replace the trial-and-error procedure with a systematic search that can be programmed for a digital computer. Use of regression analysis in reservoir description was proposed recently by Jacquard and Jain.1 They divided the reservoir into a number of homogeneous blocks whose properties are varied until a least-squares fit is obtained between observed and calculated pressures. However, they did not consider their technique to be operational, mainly because of ".. .(I) the lack of experience in using the method .. . notably for the improvement of convergence; and (2) limitations imposed by the insufficiency of available computers". 1 While the analysis presented in this paper applies the same general principle used by Jacquard and Jain, the specific method is significantly different. Some differences are (1) the regression problem is solved in a different way which requires less computer time in most cases; (2) a stepwise solution, in which the detail in the reservoir description is increased from step to step, is used to improve convergence; and (3) the reliability of the estimated reservoir .properties.. as measured by their standard deviation and correlation coefficient, is estimated'

Jan 1, 1967

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 W. B. Braden

This paper presents a suitable method for predicting gas-free oil viscosities at temperatures up to 500F knowing only the API gravity of the oil at 60F and the viscosity of the oil measured at any relatively low temperature. The API pravity and the one viscosity value are used as parameters to determine the slope of a straight line on the ASTM slanaord viscosity-temperature chart. Then, knowing the slope of the line and one point on the line, the vrscosities at higher temperatures can be determined. The line slope correlations were developed at I00 and 210F since viscosity data are frequently measured at these temperatures. A procedure is given for predicting line slopes from measurements at other tetnperatures. A nomogram is furnished for solving the relationship. The correlation has been evaluated at temperatures up to 5OOF for oils varyzng in gravity from 10 to 33 " API. The correiution is applicable only to Newtonian fluids. Comparison at 500F of true viscosities and those predicted from values at 100F shows an average deviation of 3.0 per cent (maximum deviation of 6.0 per cent). Predictions from the values at 21 0F for the same oils how an average deviation of 1.5 per cent (maximum deviation of 3.4 per cent). INTRODUCTION Correlations have been developed by Beal' and by Chew and Connally' for predicting viscosities of gas-saturated oils at reservoir conditions. Each of these correlations requires a knowledge of the solution gas-oil ratio and the viscosity of the gas-free oil at the reservoir temperature. For temperatures below 350F, measurements of the gas-free oil viscosities can be made easily using commercially available equipment. In thermal recovery processes, however, reservoir temperatures well in excess of 350F are encountered. Viscosity measurements at such conditions are more difficult and time consuming and require modification of existing equipment or the construction of new equipment. Measurements are further complicated by the difficulty of handling highly viscous oils associated with thermal recovery processes. Therefore, it is desirable to have a correlation which allows accurate prediction of viscosities at high temperatures. A commonly used technique for predicting viscosities at high temperatures is to measure the viscosities at two lower temperatures, plot the values on ASTM standard viscosity-temperature charts and extrapolate to the temperatures desired. If either of the values is slightly in error, the extrapolated value can be significantly in error. To justify an extrapolation, three points are actually necessary. This procedure can consume much time, particularly with heavy oils. Considering the cost of viscosity measurements, it would be desirable to eliminate the need for direct measurements by having correlations which would allow viscosity predictions from other physical or chemical properties. Beal1 investigated the possibility of correlating viscosity with oil gravity at temperatures from 100 to 220F. While showing that a general relationship exists, he also found significant deviations. It is possible that correlations will be developed based on oil composition as more information becomes available. While not eliminating the need for viscosity rneasurements, the method presented herein requires that only one viscosity measurement be made. The API gravity must also be known. The theory is based on the fact that the viscosity of paraffins (high gravity) changes less with temperature than does the viscosity of naph-thenes or aromatics (low gravity). The gravity. therefore, is used as a parameter to determine the slope of a straight line on the ASTM standard viscosity-temperature charts. The correlation is applicable only to Newtonian oils, and deviations due to thermal decomposition and nonhomo-geneity cannot be predicted. Oils containing additives have not been evaluated. PROCEDURE Fifteen oils were used in developing the correlation; eight were crudes and seven were processed oils. Oil gravities varied from 9.9" API (naphthene base) to 32.7' API (paraffin base). The temperature range studied was 81 to 516F. Each oil used had a minimum of three viscosity measurements and each plotted essentially as a straight line on the ASTM charts. In all, 91 viscosity measurements were used in the correlation. Saybolt, rolling ball and capillary tube viscometers were used for the measurements. Viscosity data for Samples 1, 2, 4, 7, 10, 11 and 14 were obtained in Texaco, Inc. laboratories. The data for Samples 3, 5, 6, 8, 9, 12 and 15 were from Fortsch and Wilson,3 and data for Sample 13 were from Dean and Lane.' All data points used in the correlation are plotted in Fig. 1. It is seen that some of the viscosity values deviated slightly from the straight-line plots at the higher temperatures. Properties of the oils after exposure to the

Jan 1, 1967

By W. E. Showalter

This paper discusses some of the results of combstion-drive tests which were made in a test cell using a sand bed 10 in. in diameter x 10-ft long. The test method is illustrated and described. The relationship between the API gravity of the in situ oil and the amount of air required for combustion drive is discussed in detail. Other things constant, the air requirement for combustion drive increases as the APl gravity of the in situ oil decreases. If the test results apply to actual reservoirs, the lowest-priced oils may cost the most to recover by this method. Information is shown which indicates that the effect of pressure on the amount of hydrocarbon burned is not large. A method of predicting air requirements from the API gravity of the in situ oil is presented. INTRODUCTION Combustion drive is the term used to identify the process of interstitial or in situ burning as an oil recovery method. Part of the in situ oil is burned to generate the energy needed to produce the remainder of the oil. Combustion drive as an oil recovery mechanism remains an economic uncertainty in spite of all the work that has been done by the industry in both laboratory and field?-6 This paper will show some of the results of tests which were made in a test cell for the purpose of studying the nature of the combustion-drive process. It will present data which indicate that the API gravity of the in situ oil is a significant indicator of the amount of air required to drive a burning front through oil sand. Air requirement varies inversely as the API gravity of the in situ oil. EXPERIMENTAL The tests were performed in a cell which utilized a cylindrical sand section 10 in. in diameter x 10-ft long. The thin-walled metal pipe which held the sand was wound with twenty external electrical resistance heaters which, by means of an automatic controller, maintained adjacent sections of the wall of the pipe at temperatures equal to the temperatures of the contained sand. Each heater covered 6 linear in. of the pipe. By this means lateral heat loss from the sand section was minimized, thereby causing the sand section to simulate more closely a horizontal increment of a combustion-drive reservoir. Fig. 1 shows a schematic diagram of the test assembly. Thermocouples to measure the temperature in the sand were located every 6 in. along the length of the sand section. The pipe containing the sand was enclosed in a cell designed for an operating pressure of 500 psig. The inlet air pressure was controlled at the inlet, and the gas flow rate was controlled and measured at the outlet of the cell. The oil sand used for the tests was prepared by mixing first water and then oil with the noncon-solidated sand using a closed mixer similar to a cement mixer. Table 1 shows a screen analysis of the sand. Ninety percent of the sand was 100 mesh or finer. This sand was a mixture of 80 per cent No. 120 Nevada White Sand and 20 per cent Tennessee Hi-Fusion Moulding Sand No. 3. The Nevada sand was a clean silica sand. The Tennessee sand con-

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 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

By J. N. Sicking, F. H. Brinkman

A new method for obtaining equilibrium vaporization ratios (K-values) for reservoir fluids has been developed and tested. By application of the method, complex experimental measurements of liquid and vapor phase compositions are eliminated. This simplified technique reduces the cost of experimental equilibrium-ratio data for reservoir studies of condensates and volatile crude-oil systems. The method is designed for systems of constant composition and, therefore, is best suited for depletion studies where compositional changes at high pressures are minor. The basic data required, in addition to the composition of the initial reservoir fluid, are the relative vapor-liquid volumes and densities at reservoir temperature and variom reservoir pressures. Tests demonstrated that the method predicts equilibrium ratios accurately for condensates. A single test on a crude oil was not conclusive; further testing will be necessary before the accuracy of the method can be determined for crude-oil systems. In addition to determining equilibrium ratios, the calculation method provides information on the physical properties of the "plus" component in the vapor and liquid phases. The "plus" component is that mixture of components heavier than the least volatile fraction analyzed. This information is useful in studies of both natural depletion and cycling operations for condensate reservoirs where the heptanes-plus component in the gas phase is produced from the reservoir. INTRODUCTION As more volatile oil and condensate reservoirs are found, the use of phase behavior techniques to predict their performance is increasing in importance. These techniques have long been used for condensate fields and have more recently been applied to crude-oil fields containing oils of medium-to-high volatility. In these phase behavior methods, equilibrium ratios (K-values) are used to predict compositional changes in the reservoir fluids—thereby accounting for the recoverable oil that exists in the gas phase. The reliability of the predictions depends to a large extent on the equilibrium ratios used. These values must be obtained for each component for the entire pressure range being investigated. Unfortunately, because of the complex nature of hydrocarbon mixtures, accurate K-values are hard to obtain. The equilibrium ratios for a particular component will vary not only with the temperature and pressure, but also with over-all composition of the system. The importance of composition is quite critical at elevated pressures, but becomes negligible at pressures below about 300 psia. Therefore, because most phase behavior problems involve the high-pressure region, each fluid system becomes a special case. Experimental programs to determine characteristic K-values are quite difficult and time-consuming. Thus, it is often necessary to resort to approximations of the K-value data. Charts giving K-values for various mixtures and classes of mixtures are available in the literature. However, there are two major difficulties in using them: (1) the K-values of the "plus" component (that mixture of components heavier than the last one analyzed) must be obtained by extrapolation from the K-values of the other components; and (2) the K-values obtained must finally be adjusted by trial and error to agree with observed volumetric data. To eliminate these difficulties, a new method of determining equilibrium ratios was developed. Briefly, after the composition of the system as a whole has been analyzed, the method uses empirical correlations and the gross fluid properties of the system (relative vapor-liquid volumes and densities) to calculate K-values. Because the calculative procedure is long, it is best solved on a digital computer. About one hour of machine time on an IBM 650 computer is required to develop a K-chart for the fluid being examined. DEVELOPMENT OF THE METHOD OF OBTAINING K-VALUES Equilibrium ratios are defined as the ratio of the mole fraction of a component in the vapor phase to its mole fraction in the liquid phase. This statement is expressed in Eq. 1. A typical plot of equilibrium ratios for a particular system is shown in Fig. 1. It should be noted that, at pressures near the saturation pressure, the K-values appear to converge to a common point. This apparent convergence point is called the convergence pressure and is a characteristic of the system involved. Various empirical correlations of K-values have been noted. It has been observed that an isothermal plot of