Search Documents

By G. W. Tracy, R. D. Carter

A method based on an approach by Hurst' is developed for calculating water influx behavior. Using this method, superposition calculations may be eliminated. The principal diflerence between this method and that of Hurst is that, over finite intervals of time, constant oil-productton rates are assumed by Hurst whereas constant water influx rates are assumed in the presently described method. By adopting the assumption of constant water influx rate for finite time periods, a method is developed for calculating water influx behavior which lends itself to convenient combination with the material balance equation in the Schilthuis form. Above the bubble point, an equation is given for the explicit step-wise calculation of pressure history for a prescribed oil-production history. Use of this method of calculating water influx behavior should result in a considerable saving of time as compared to superposition calculations when the calculations must be made without using a digital computer. The connection is shown between the method described herein and the method developed by van Ever-dingen and Hursr2 (which requires superposition cal-culations). Also, results obtained using the two methods are compared.

By J. E. Smith, P. B. Crawford, B. L. Laudrum

A method is presented by which the displacement efficiency of a steam-injection process may be determined. The method consists of a determination of the position of the steam front at a particular time and the application of the Buckley-Leverett theory to the hot and cold regions of the reservoir with the appropriate hot or cold fluid properties. A uniform, linear flow system is assumed for this study. The method is applied to a number of reservoir conditions. The primary properties which affect the results of a steam-injection process are the oil viscosities at steam and reservoir temperatures, injection rate, pay thickness and oil saturation. The efficiency of steam injection is compared with that of water injection for a 10-acre, linearly flooded reservoir. In reservoirs where the oil viscosity is about 35 times the water viscosity, the steam process yields recoveries about 20 per cent higher than the water flood. Where the oil viscosity is near 500 cp, the steam process may recover from three to four times as much oil as a water flood. INTRODUCTION In a previous paper,' a method was suggested for evaluating a steam-injection process. Emphasis was placed upon the determination of the steam-front position in the reservoir at various times. In the present paper, the displacement process is studied in more detail. The approach is essentially that employed by Jordan, Rayne and Marshall3 in their study of hot-water injection. The steam-injection process is quite complex, and an exact description of the process has not been found; however, a qualitative description can be made. When a quantity of steam is injected into a relatively cold reservoir, the leading edge of the steam front continaally flows into cool sand, and some steam condenses. As more and more steam condenses, a bank of condensed water will form which is augmented by displaced connate water. The reservoir may be divided in three regions—(1) a region containing displaced fluids at reservoir temperatures, (2) a region contain- ing water at reservoir temperature and (3) a region occupied by steam. The reasons for the lack of a region of hot water or oil of any appreciable size have been discussed.' A method has been presented by which the location of the steam front can be determined even when heat losses to the regions surrounding the reservoir are sizable compared with the amount of heat injected. It remains to show the effectiveness of steam injection in increasing the ultimate oil recovery. To illustrate the technique for determining the increase in displacement efficiency and, hence, the increase in ultimate production, two sets of hypothetical reservoir conditions will be studied in detail.

By B. K. Larkin, H. R. Bailey

A model Of heat flow in an underground combustion process is studied. This model includes convection effects and thus is more general than previous studies which considered conduction as the only mechanistn for heat transfer. Both linear (tube run) and radial (field application) geometries are considered. The effects of ignition heaters, vertical losses and finite source width are considered for the linear case. The results are in the form of equations and are presented in graphical form for a number of cases. Convection effects increase frontal temperatures about 25 per cent over those computed for conductive transfer for typical field operating conditions. This increase in temperature is a result of heat being transported toward the front by the injection gas. Even greater temperature increases are realized as the per cent oxygen in the injection gas decreases. It is well known that compression costs are of considerable importance in estimating the economic feasibility of underground combustion. By assuming an ignition temperature for the combustion fuel, predictions of limiting conditions on fuel density and injection rate necessary to sustain the combustion zone are made. For typical field conditions, at least 0.75 lb/cu ft of fuel are needed with air as the injection gas. If the injection gas is 10 per cent oxygen and 90 per cent nitrogen, this figure is 0.69 lb/cu ft. INTRODUCTION The possibility of increasing oil recovery by underground combustion has been considered for many years. Recent field tests indicate that the underground combustion process is technically feasible. The economic feasibility depends to a large extent on the amount of air which must be injected to sustain combustion. Prediction of the success of employing underground combustion in a particular reservoir must be based on existing field tests, laboratory tube runs"' and solutions of mathematical or analog models. Vogel and Krueger5 devised an electrical analog of the heat transfer problem in an underground combus-ion process. They considered the problem of heat conduction from a cylindrical source with increasing radius assuming no vertical losses. Bailey and Larkin and Ramey solved the corresponding problem including vertical losses by considering a mathematical model of the conduction process. Ref. 9, 10 and 11 also present information related to the subject study. The present paper generalizes these results to include both conduction and convection mechanisms for heat flow. A model of the conduction-convection process is described. The partial differential equations governing this model are written and solved for a number of cases of interest. The formulas are evaluated and the results are presented graphically. A CONDUCTION-CONVECTION MODEL Description of the Model A number of models are employed in this paper, each of which approximate a particular physical situation depending on the geometry, i.e., linear or radial flow, the type of vertical heat losses, the nature of the source, etc. The geometry for both the linear and radial case is shown in Fig. 1. A particular linear model is described in this section, and variations from this model are described in later sections of the paper. Consider a homogeneous, semi-infinite porous prism with lateral sides insulated as shown in Fig. 1. The left end, x = 0, is heated for a time t, with a plane heater supplying a constant heat flux and during this time; gas (nitrogen) is passed through the pores of the prism at a constant mass flux. At time t, the heater is turned off and air is injected in place of the nitrogen. This causes combustion to take place, and a uniform source of width 1 is formed which moves at a constant velocity v, to the right. The frontal velocity v, is related to the gas velocity v,, through stoichiometric considerations

By W. L. Penberthy, H. J. Ramey

Experimental work on the combustion oil recovery process has consisted of both laboratory and field studies. Although field experiments are the ultimate test of any oil recovery process, they are costly, time consuming and difficult to analyze quantitatively. Laboratory combustion tube experiments can be opertzted far more rapidly and cheaply, but are subject to scaling and interpretation problems. This paper points out some important design problems, operational criteria and considerations in2portunt to interpretation of results. An analytical hcut model of movement of a burning front axiully along a cylinder with heat loss through an annular insulation was developed. The result was used to identify steady-state temperatrrre distributions both ahead of and behind the burning front, with and without heat loss. Results indicate potential operating limitations on the minimum burning front rielocity (or air flux) which may be used for any given combustion tube. Results also enable estimating the effective thermal diffusivity and over-all heat loss from experimental data and thickness of the brnning zone. Resu lts of operation of a combustion tube constructed recently verify this preliminary theory in the region immediately ahead of and behind the burning front surprisinRly well. INTRODUCTION Many field and laboratory studies of the forward combustion oil recovery process have been conducted since the early publications of Kuhn and Kochl and Grant and Szasz2 in 1953 and 1954. In view of the complex and costly nature of this type of investigation, it is not surprising that no complete theory of the nature of the forward combustion Drocess is yet available. However, gross effects are understood well enough that reasonable design procedures are available for planning field operations. Nelson and McNeil 3,4 have published two comprehensive papers concerning design procedures. One major consideration in planning field operations is the fuel concentration at the burning front. Fuel concentration controls air requirements — an important cost factor in forward combustion. Although fuel concentration can be estimated from field test results by various methods, results are subject to great uncertainty in view of natural limitations on experimental observations. Nelson and hlcNeil4 recommend that fuel concentration be determined from laboratory combustion tube studies. Fuel concentration is only one of many important factors which can be studied by combustion tube experimentation. An obvious goal of imporrance must be development of a comprehensive theory of the forward combustion process. If a theory of this process can be established which matches controlled laboratory experimentation, it should be possible to apply this theory to field operating conditions with some confidence. Laboratory combustion tube studies have already yielded important information concerning the combustion process. However, details concerning the design, construction and operation of combustion tubes are rare. Combustion tubes used by various investigators vary in size, length and mode of operation. Therefore, one purpose of this paper is to present considerations important to design, construction and operation of a combustion tube. In regard to previous combustion tube studies, attention is called to Refs. 1 through 9. These references describe a wide variety of equipment types and present a great deal of pertinent experimental data. In general, combustion tubes usually consist of thin-walled stainless steel tubes containing an oil sand pack mounted within a pressure jacket. Provisions have often been made to heat the tube ex tern a lly by separately controlled heaters to reduce heat losses. This step usually permits operation at low air fluxes (air rate per square foot burning front surface) similar to those encountered in field operations. Burning is usually conducted from the air inlet end of the tube to the outflow end. The tube orientation used has been vertical or horizontal. For vertical tubes, burning has been conducted vertically downwards. As a result of laboratory studies, many of the

Jan 1, 1967

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.'s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,' but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

By C. W. Arnold, H. L. Stone, D. L. Luffel

The purpose of this research is to determine (I) the efficiency of small banks of enriched gar driven by methane in displacing oil from a porous medium and (2) the effects of variation in bank size and composition of that efficiency. Most of the experiments were conducted in a sand-packed tube 20-ft long and 1/2-in. in diameter. The hydrocarbon system generally used was methane, butane and decane at 2,500 psia and 160°F. The results of these experiments indicate that, in the regions contacted by the gas, a small bank of an oil-miscible gas driven by methane can displace all of the oil in a piston-like manner. If the enriched gas is of such composition as to remain immiscible with the oil, displacement of oil is less efficient than for the miscible case, and the gas bank travels through the sand with a velocity less than that of the driving gas. These data along with theories discussed imply that smaller banks and less total gas are required when the enriched gas and oil are miscible. INTRODUCTION Widespread application of enriched-gas drive to the recovery of oil rests upon a key factor — the use of limited quantities, or "banks", of enriched gas. At the present time, the value of liquefied petroleum gas or other enriching agents discourages their use in a continuous injection technique, or even in a large bank, except in a few isolated reservoirs. If small banks of enriched gas driven by methane were as effective in displacing oil as is continuous injection, the enriched-gas drive process might be applied to a larger number of reservoirs. Previous research on the mechanics of the enriched-gas drive process reported by Stone and Crurnpl and by Kehn, Pyndus and Gaskell has utilized continuous injection of enriched gas. This work has shown that two types of displacements occur. With gases containing sufficient intermediates. the oil is displaced misciblv and complete recovery is obtained from the regions swept. When gases are used which contain insufficient intermediate hydrocarbon for miscible displacement, oil is displaced immiscibly. In the latter type, selective solution of the intermediate hydrocarbons causes a swelling and reduction in viscosity of the oil and leads to an increased recovery over that obtained by dry-gas (methane) drive. The size of the enriched-gas bank necessary for efficient displacement of oil is determined by those factors which cause deterioration of the bank. A differentiation may be made between those factors which operate on a microscopic scale and those which act on a macroscopic scale. On the smaller scale, the enriched gas mixes in the direction of flow by diffusion and convection with the fluids immediately preceding and following it. On the larger scale, the gas may by-pass the oil by flowing through permeable streaks, by overriding the oil because of density difference, or by fingering because of unfavorable viscosity ratios. In such cases, the enriching material tends to mix with the oil both laterally and in the direction of flow. The increase in effective area available for diffusion and dispersion of the enriching components leads to a faster degradation of the bank and a need for a larger bank than is necessary for those cases in which no by-passing occurs. The effects of such macroscopic factors in the deterioration of enriched-gas banks have been reported in a separate paper by Blackwell, Terry and Rayne. The present study was confined to the factors which operate on the smaller scale, in particular to the behavior of banks of enriched gas in sands uniformly swept by the gas. Experiments were designed to answer the following questions. 1. Can small banks of enriched gas driven by methane be used to secure oil recoveries comparable to those obtained by continuous injection of enriched gas? 2. What is the optimum bank size (the minimum bank size necessary to obtain a recovery comparable to that obtained by continuous injection of the same enriched gas)? 3. How many total pore volumes of gas must be injected to obtain the maximum recovery when the optimum bank size is used? 4. What is the effect of varying the number of enriching components in a gas bank? This report describes the experimental investigations and discusses the results in terms of their significance to reservoir behavior.

By C. van der Poel

When oil is displaced from a horizontal formation by another fluid of lower density, the latter tends to override the former in the shape of a tongue owing to gravity segregation. This gravity tongue has an adverse influence on the oil recovery. If the fluids are miscible, diffusion (mixing) takes place at the interfacial boundary of the gravity tongue. This mixing should have a favorable effect on oil recovery. The report describes a laboratory study of the magnitude of the mixing zone under various conditions, so as to assess the effect of diffusion on oil recovery both in laboratory experiments and under actual field conditions. The technique used enables visual observation and measurement of the size of the mixing zone in transparent glass-powder packs. The results show that in experiments in small models and cores the width of the mixing zone may well be of the same order of magnitude as the height of the model. In such cases oil recovery is favorably affected by mixing. It can further be concluded that, under conditions prevailing in the field, mixing of the injected fluid with reservoir oil is equal to that caused by molecular diffusion alone, eddy-mixing not taking place to any appreciable extent. A simple calculation, then, shows that molecular diffusion is too small for a beneficial effect to be expected from the injection of miscible fluids in horizontal or nearly horizontal reservoirs unless pay zones are thin. INTRODUCTION This paper gives results of experiments made on the mixing which occurs when a miscible displacement is carried out in a horizontal reservoir. This mixing takes place at the interfacial boundary of the gravity tongue formed when the lighter injected fluid overrides the oil present in the reservoir. The object of the experiments was to simulate the field case where, for example, propane with a viscosity of 0.075 cp under reservoir conditions displaces an oil of viscosity 0.6 cp. In the experiments the lighter fluid had a viscosity of about 1 cp and the heavier one a viscosity of about 8 cp, so as to obtain the same viscosity ratio. In order to enable the results to be compared with those published in the literature, a set of experiments with viscosity ratio equal to one was also performed. EXPERIMENTAL TECHNIQUE The technique developed for the purpose enables the width of the mixing zone to be studied as a function of time and place. A glass-powder pack is saturated with a water-glycerine mixture of suitable viscosity, which represents the reservoir fluid. The pack is then rendered transparent by dissolving sufficient ammonium thiocyanate in the mixture to obtain a solution with a refractive index matching that of the glass powder. Alkaline water to which phenol-phthalein has been added is used as the lighter and less-viscous displacing fluid. In those places where mixing or diffusion of the two liquids occurs, the slightly acidic ammonium-thiocyanate solution neutralizes the alkali in the water, and the glass pack shows up white against the red-colored invading water. If a black cloth is hung over the back of the apparatus, the transparent part of the glass pack appears black. In this way the position of the two phases and the transition zone between them is clearly visible, as shown in Fig. 1 (where the red-colored invading water is the grey-shaded zone lying uppermost). In all experiments the amount of alkali in the water was chosen such that the upper contour corresponded to a concentration of 5 per cent of the dense liquid. The lower contour is determined by the size of the glass grains and the thickness of the pack and, consequently, varies from experiment to experiment. Corresponding concentrations of the dense liquid range from about 97 per cent for packs of small grain size (permeability = 1 darcy) to about 95 per cent for those with large grains

By D. W. Reed, R. L. Reed, Tracht

Laboratory experiments on the reverse combustion of tar sands in a linear adiabatic system have shown that a highly upgraded oil can be produced from an exceedingly viscous, immobile oil. The dependence on the air-injection rate of peak temperature, combustion-zone velocity, oil recovery, air-oil ratio, residual coke and oil, fuel burned and distribution of product gases is shown graphically. Eflects of initial temperature, oxygen concentration, oil saturation and heat loss are discussed. Experiments bearing on the coking properties of heavy oils are mentioned and some results exhibited. Field application of the process hinges on the existence of adequate air permeability and the rate of reaction under reservoir conditions. INTRODUCTION It has been established that oil can be recovered from underground reservoirs by means of at least two fundamentally distinct methods involving in situ combustion of a certain fraction of the oil. Characteristic of both of these known methods is the production of oil from one or more wells by means of hot gases formed when a high-temperature reaction zone is advanced through the reservoir. In both cases, the reaction zone is created by heating certain of the wells to a sufficiently high temperature prior to the introduction of air, and the zone is maintained and advanced through the reservoir by appropriate control of the air-injection rate. In the first of these methods, which is called "forward combustion",' the combustion zone advances in a direction which is generally the same as that of the air flow; whereas in the second method, "reverse combustion",' the combustion zone moves in a direction generally opposite to that of the air flow. Forward combustion, on the one hand, is an ideal combustion process in the sense that a minimum of the most undesirable fraction of the oil is consumed as fuel in the form of coke, a clean sand is left behind and generated heat is used as efficiently as possible. However, the applicability of forward combustion is limited. Since the products of combustion, vaporized oil and connate water must flow into relatively cold regions of the reservoir, there is an upper limit on the viscosity of oil which can be moved by this process in a practical and economical fashion.' On the other hand, it is characteristic of reverse combustion that the vaporized oil and water together with the products of combustion are produced through sand which is already hot and has had its mobile liquid content eliminated. This means there is no upper limit on oil viscosity; indeed, the oil may be an entirely immobile semi-solid. However, fuel for the process is an intermediate fraction of the original oil, and the most undesirable fraction remains on the sand surface as a substantial deposit of coke. Since this coked material is not burned during reverse combustion, it represents energy which is available for the production of oil but is not used for this purpose. It follows that one can expect economics to be somewhat less attractive with reverse combustion than with forward combustion. Nevertheless, it is a process which is designed for reservoirs where forward combustion is impossible and, as such, has become a subject of experimental and theoretical investigation. In this paper, only experiments made with tar sands are discussed. DESCRIPTION OF THE PROCESS We proceed, then, to consider the process of reverse combustion in greater detail. Fig. 1 illustrates a temperature profile defining a combustion zone which moves from right to left when air flows from left to right. In Zone 1, the temperature is the initial reservoir temperature, and the tar sand is as yet unaltered. This statement must be modified to the extent that physical properties of the oil may be changed by low-temperature oxidation at reservoir temperature. As air passes into Zone 2, which has been warmed by conduction, it assists in vaporization of the very light ends (if there are any), and oxidation occurs at a significant rate. In this region, there is almost no production of carbon monoxide or carbon dioxide because predominantly addition-type reactions take place with the formation of oxygenated compounds such as aldehydes and acids together with water. The hydrocarbon-enriched and slightly oxygen-depleted gas stream enters Zone 3 where

I would like to congratulate the authors on an interesting paper and, in particular, on their justification of the use of superposition in treating gas wells. They state, however, that no theory has previously been developed lor treating the variable rate drawdown case. In a paper published a few years ago by my colleagues and myself, a method is described which allows a constant rate drawdown curve to be derived for any well. regardless of real rate variations. This method has the great advantage that it is not necessary to make any assumptions regarding the form or uniformity of the reservoir, and the derivation suggested by the present authors is simply an idealized case of our more general treatment. In the light of several years' experience in the derivation of constant rate drawdown curves: I would like to point out that the method proposed by Odeh and Jones may lead to confusing answers, particularly if applied over an appreciable time interval. The reason is that, if pressures are measured when a well is flowing, non-linearities in plots like Figs. 1 and 2 can arise due to either (1) changes in the skin factor with time and/or flow rate, or (2) an inhomogeneous reservoir. Separation of these two effects is only possible if some of the pressures used in the analysis are measured when the well flow rate is 0 or, in other words. during brief shut-in periods. REFERENCES 1. Arthur, K. B.. Metunomski, Z. G. and Ruche, Ms "Intetet ct methode d'utilisa tion des mesures de pression [)our Petude des gisemonts" Revue de T'I.F.P. (Nov., 1962) XVII, 1299. 2. Thomere, R. and Arthur, K. R.: "Analysis of the Pressure Bet~avior of a Production Well Subject to Interterence and Its Application to Reservoir Models", Paper SPE 1189 presented at SPE :Amund Fall Mecting Ireld in Denwer, Colo. (Oct. 3-6, 1965).

Jan 1, 1966

By J. P. Heller

The purpose of this discussion is to comment briefly on one of the points raised in the paper. The topic is that referred to in the Conclusion 2, in which the authors note a "disproportionate increase in mixing rate" as they proceed from a displacement performed at viscosity ratio 0.998 to one at 1.002. The recovery curve of this displacement, Fig.2A, shows a quite noticeable difference between the two runs in which displaced and displacing liquids were interchanged. This difference is ascribed by the authors to the change in viscosity ratio between the two experiments. This interpretation of the experimental results must certainly be based on rather precise measurements of the fluid viscosities or, at least, of the ratio of their viscosities. The ordinary Ostwald viscometer, employed in the usual way, is only just able to distinguish between liquids whose viscosities differ by two or three parts per thousand. The experimental refinements which then must have been necessary to measure with any precision a viscosity ratio of 1.002 would be worthy of mention. Indeed, one would have expected the authors to devote at least a paragraph to the display of this technique. For the same reason - to emphasize that the details of such unusually critical experimental work should be reported and discussed in the literature — I would have liked to see some comment by the authors on the requirements for temperature control during both the viscosity measurements and the displacement experiments. The temperature coefficient of viscosity of Soltrol C at room temperature has been measured in this laboratory at about - 1.9%/1°C. How much did the coefficients for the fluid mixtures used in the experiment differ from this and from each other Even if the two fluids were identical in this respect, it would appear that a temperature difference of .l°C between them in the flow system would produce a viscosity con- trast equal to that imposed by their composition. Reference to the experimental solutions of these problems would render Conclusion 2 more credible. Published theoretical work on the stability problem contains no indication of such extreme sensitivity to the viscosity ratio in the immediate neighborhood of unity. On the contrary, perturbation calculations such as those made by Chuoke1 and Perrine2 lead to initial finger growth rates proportional to the first power of the viscosity difference. The later instability growth calculated by Scheidegger,3 using the assumption of zero pressure gradient perpendicular to the average flow, is also simply proportional to the viscosity difference. The conflict regarding the sensitivity of the dependence of instability growth or mixing rate on the viscosity contrast needs explicit recognition in the literature. In particular, it would seem that the experimental work reported here should be supported by a more detailed description of the viscosity controls. Further discussion of possible alternative interpretations of the data would also be in order. For instance, is it possible that the unconsolidstted bead pack used changed its characteristics between the runs (due, perhaps, to the handling and inverting involved in the injection of the front) Or, are the data points of Fig. 2A averages from several runs performed in a series in which the two liquids alternated as displaced and displacing fluids? To summarize, it seems to me that the relationship asserted by Conclusion 2 in the paper needs the support of more detailed experimental evidence. REFERENCES 1. Chouke, R. L., van Mews, P. and van der Poel, C.: "The Instability of Slow Immiscible Viscous Liquid-Liquid Displacements in Porous Media", Trans., AIME (1959) Vol. 216, 188.

By K. T. Koone, R. J. Blackwell

One of the more important problems to be solved in designing a miscible flood is related to the size of the solvent bank used. Size of the bank may be critical to economic success. Too large a bank loses money; too small a bank may deteriorate and fail to maintain the miscibility needed for high recovery. An important factor in deterioration of a small bank is permeability channeling. In a highly stratified reservoir, solvent speeds ahead in the more permeable zones and mixes laterally with fluids bypassed in adjacent, low-permeability strata. Numerical solutions have been obtained .for the differential equations that describe the movement of a slug through a two-layer system in which mixing occurs both in the direction of flow and transversely. The solvent slug is assumed to have the same density and viscosity as the resident fluid and the pushing fluid. These solutions have been verified by comparing them with similar concentration profiles obtained in the laboratory in a 36-ft stratified model packed with glass beads. The theoretical study revealed that when the dominant mechanism causing a bank to fail is lateral mixing the bank size needed for a given recovery may increase with length rather than decreasing as the square root of reservoir length, as suggested by one-dimensional mixing theory. From a comprehensive examination of the variables, a generalized correlation is developed that relates strata thicknesses, bank size, fluid velocity, mixing coefficients, system length and simple solvent-resident fluid phase behavior to the area miscibly swept. INTRODUCTION Miscible displacement, or solvent flooding, continues to receive widespread attention as a method for increasing oil recovery over that possible in conventional gas-drive or waterdrive projects. A basic economic requirement in the application of such processes is the use of as little solvent as possible. A basic physical requirement is that enough solvent be used to maintain miscibility. Economics places an upper limit on the size of a solvent slug, and physical considerations establish a lower limit. Consequently, the practicality of any given miscible process requires that the economic limit be greater than the lower limit imposed by physical requirements. Procedures exist for determining the economic limit; however, procedures for determining realistic minimum bank sizes exist for only special reservoir situations. In the past, bank size has usually been selected on the assumption of a piston-like displacement for which only longitudinal mixing is important. This assumption leads to the favorable conclusion that bank size, expressed as per cent of pore volume, is inversely proportional to the square root of length. Collins considered the problem of transverse mixing of solvent with fluids in bypassed zones with the assumptions that no forward mixing occurs and that the concentration is uniform in the permeable stratum of interest. 1 Lauwerier considered a mathematically similar problem in thermal recovery operations.2 Their work suggests the much less favorable conclusion that bank size could be directly proportional to the length or even higher powers of the length. This paper considers the physical behavior of a small solvent bank as it moves through a real reservoir, without imposing many of the restrictive assumptions of past treatments. To facilitate mathematical description, an idealized, two-layer model that permits mixing both laterally and in the direction of flow will be considered. A calculation procedure for solving the descriptive equations will be developed for a displacement involving fluids of equal density and viscosity. In addition, laboratory experiments designed to check the computational results will be described, and the coincidence of calculated and observed results will be discussed. The use of many solutions for a variety of bank sizes, strata thicknesses and characteristic system lengths to develop a usable correlation between

Jan 1, 1966

By W. M. Rutherford

A knowledge of the limits of miscibility between reservoir oil and possible injection fluids is required for selection of the optimum miscible-injection fluid. Limits of miscibility can be estimated from the results of equilibrium phase-behavior experiments. They can also be determined by means of displacement experiments conducted in a high-pressure sandpack. This paper describes the equipment and procedure which have been developed for determining miscibility conditions by stable displacement. A systematic series of displacements of a West Texas reservoir oil was carried out. The results indicate that, at constant pressure, miscibility is a function only of the pseudocritical temperature of the injection gas. This fact, together with improved experimental methods, makes the displacement technique a rapid, reliable means for determining miscibility conditions. In conjunction with the displacement experiments, phase diagrams were constructed for the oil with dry gas and propane and with dry gas and ethane. Phase behavior of the methane-ethane-propane system was determined at 110°F. The experimental work demonstrates the feasibility of using ethane-rich gases to reduce cost and pressure requirements for miscible displacement. INTRODUCTION In recent years, interest in the miscible displacement of oil by light hydrocarbon mixtures has been high. Many pilot and a few field scale projects have been started.1 These projects have made use of various methods for achieving miscibility: (1) the LPG-slug process, (2) the enriched-gas-drive process and (3) the high-pressure gas-drive process. Some field projects have been successful; the results of others are debatable. In general, projects which have performed best have involved the injection of an appreciable fraction of a pore volume of miscible material. Economical application of miscible displacement depends strongly on the cost of the miscible-injection fluid. If an appreciable fraction of a pore volume of material is required for successful application of these methods, a precise knowledge of the minimum requirements for miscibility in terms of composition of injection fluid is essential. Therefore, reliable experimental methods for determining miscibility conditions are important, and a procedure for estimating these conditions from the composition of the reservoir fluid is highly desirable. The subject of this paper is the problem of determining conditions which result in miscible displacement of oil by light hydrocarbon mixtures. Miscibility conditions can be estimated by means of equilibrium experiments conducted in a PVT cell, or they can be determined by means of high-pressure displacement experiments. This paper describes the equipment and procedure which have been developed for the determination of miscibility by high-pressure displacement experiments. These methods have been applied to the displacement of a West Texas reservoir oil with mixtures of dry gas, ethane and propane. In conjunction with the displacement experiments, triangular phase diagrams have been constructed for mixtures of the oil with dry gas and propane and with dry gas and ethane. The effect of injection-gas composition on conditions for miscible displacement in high-pressure sandpacks and cores has been the subject of several published papers.2-5 The experimental methods used in these investigations were such that displacements were unstable, and the effects of fingering and/or gravity layover are clearly evident in the results. Miscibility conditions were probably correct in spite of the instability phenomena, but the experiments evidently were time-consuming, and limited data were reported. Systematic high-pressure flow studies which would support a correlation of miscibility conditions have not been reported; however, Wilson5 has proposed the use of pseudocritical temperature of the injection gas as a parameter and Benham, et al,6 have based a miscibility correlation on observed and calculated equilibrium data.

By A. L. Benham, W. E. Dowden, W. J. Kunzman

The mechanism involved in creating a miscible fluid displacement of typical reservoir fluids through porous media by light hydrocarbon mixtures usually involves a number of mixing and separation stages in the contact area. This mechanism has been illustrated on ternary diagrams for miscible displacement of reservoir fluid by rich hydrocarbon gases. A relatively simple method has been devised for calculating the approximate conditions for a miscible displacement of reservoir fluid by rich, light hydrocarbon gases or by LPG mixtures. Miscibility is favored by increased pressures, decreased temperatures, light reservoir fluids and light hydrocarbon-displacing fluids rich in the heavier components. The conditions for miscible displacements of five typical reservoir fluids by various light hydrocarbon mixtures were determined for a range of temperatures from 70" to 260°F and pressures from 1,000 to 3,000 psia. A correlation of these conditions was obtained using the variables of temperature, pressure, C,+ molecular weight of the reservoir fluid, C,+ molecular weight of the displacing fluid, and mol per cent methane in the displacing fluid. This correlation has been applied to data appearing in the literature with very good success. A partial check of the correlation has been made using a step-wise procedure with a windowed PVT cell. INTRODUCTION The term, "miscible fluid displacement", may be defined as any oil-recovery displacement process where there is an absence of a phase boundary or interface between the displaced and displacing fluids. A water flood would not be a miscible fluid displacement, although cycling in a condensate reservoir would be. Another example is displacement of oil by gasoline. Not so obvious, however, are the conditions under which fluids (such as propane, mixtures of propane and methane, or other similar systems) will give a miscible fluid displacement of oil. There are four standard miscible-fluid-displacement processes being used for the recovery of oil from an "nderground reservoir—(1) high-pressure dry-gas miscible displacement, (2) rich-gas miscible displacement, (3) miscible slug using a slug of LPG followed by a lean dry gas under miscible conditions, and (4) gas driven-slug process using a slug of LPG followed by a lean dry gas under immiscible conditions. The first problem faced by the engineer designing a miscible displacement process is to determine the conditions for miscibility. A number of papers have included discussions of the process of obtaining miscibility in these processes.1-8 However, methods for calculating the conditions for miscibility have not been published. This paper includes a method for calculating the conditions for a miscible displacement of reservoir fluid by a rich gas (LPG diluted with methane or natural gas). In this paper, a rich-gas miscible displacement process is defined as an oil-recovery process in the reservoir. The cost of a rich-gas miscible displacement will be decreased by diluting the LPG with methane or natural gas as much as possible. Therefore, it is important to be able to predict the maximum dilution for miscibility. MECHANISMS The mechanisms for obtaining miscibility between a displacing mixture of light hydrocarbons and displaced reservoir fluid often involve a number of displacing and mixing stages. As stated previously, the cost of a rich-gas miscible displacement is decreased by diluting available LPG with natural gas. The problem is to predict the maximum dilution allowed while maintaining miscibility. To get a picture of the displacing and mixing stages, it is helpful to consider these complex hydrocarbon mixtures as pseu do-ternary systems. Graphical representation of a complex system as a pseu do-ternary system can be thermodynamically rigorous if enough variables are used to describe the character of each pseudo-component in one of the phases. However, even if the characterization of each pseudo-component is not complete, pseudo-ternary representation does serve to illustrate the mechanism for obtaining miscibility and can be used to make some general conclusions concerning the conditions for obtaining miscibility. A phase diagram for the ternary system, methane-

By H. H. Rachford, D. W. Peaceman

A system of partial differential equations describing miscible displacement of fluids in porous media is derived The system takes into account the influence of gravity, the spatial distribution of permeability, diffusion, and fluid viscosities and densities. A numerical procedure for approximating solutions to the differential systems has been tested for a horizontal two-dimensional geometry. In end-to-end displacements of oil with less-viscous solvents, the numerical solutions exhibited fingering qualitatively similar to that observed in laboratory models. Small random spatial variations in permeability about the mean value are sufficient to initiate fingering. Quantitative comparisons of computed results with laboratory data show good agreement. INTRODUCTION Miscible flooding of oil by solvent is receiving increasing consideration for field use.1-4 Unfortunately, such floods potentially present severe problems in loss of recovery through the by-passing of oil by expensive solvent. Consequently, their economic evaluation requires sound techniques for predicting recovery. The purposes of this work are to present a finite-difference method for calculating the multidimensional displacement of oil by solvent and to investigate the validity of the method by comparing results of calculations with data from displacements in laboratory models.5 The formulation of the method to simulate the model experiments treats a case of limited scope in the description of solvent flooding in the reservoir. The model experiments were carried out with fluids that were assumed to form an incompressible, ideal, two-component system with constant diffusivity. Establishing the validity of the method even in its present form provides a major step toward the goal of quantitatively evaluating individual solvent flooding projects. First, it demonstrates the fea- sibility of calculating the course of displacements which are dominated by an inherent macroscopic instability, i.e., viscous fingering. In addition, the method provides almost the only practical means of examining the effect of the size and extent of reservoir inhomogeneities on the development and propagation of the fingers, taking into account the important influences of diffusion and gravitational segregation. It might be inferred that such studies cannot yield results of practical value because in displacements dominated by macroscopic instability the pattern of finger development and the resulting performance should depend critically on small variations in the geometric distribution of reservoir inhomogeneities. The present method offers the capability inherent in computational techniques of predicting performance reproducibly with any arbitrary distribution of reservoir properties and, as such, provides a means of evaluating the sensitivity of behavior to uncertainties in reservoir definition and, thus, assessing the reliability of prediction of performance. THE PHYSICAL PROBLEM DATA USED FOR QUANTITATIVE TESTS OF CALCULATIONS In this work, a calculation for treating the miscible-displacement process is tested by comparing calculated results with experimental data. The experiments chosen for comparison have been described in detail,5 and the quantitative features are summarized in a later section. Briefly, oil was flooded by solvent of equal density from a thin rectangular channel in Lucite packed with uniform Ottawa sand. Under these conditions, two-dimensional geometry is considered adequate to represent the process. Inasmuch as the present work is concerned with testing of the method by comparison with specific two-dimensional data, the system of equations presented in the succeeding section will be oriented to defining the physical system specifically. THE DIFFERENTIAL EQUATIONS Suppose x, y to be a Cartesian co-ordinate system and define h(x, y) to be the height of a point above a horizontal reference plane. It will be assumed

By B. H. Caudle, L. G. Loncarie

Pilot flooding is one method of evaluating a proposed secondary recovery project. However, the amount and rate of oil recovery from an unconfined pilot area is not usually the same as from an equal area in a large-scale flood. This is true because fluids are free to move across the boundary of a pilot area. The result is that some of the oil produced in a pilot flood may come from outside the designated area and some of the displaced oil may leave the pilot area completely. This paper presents the results of fluid-flow model studies on a five-spot pilot flood. Mobility ratios between 0.1 and to have been studied. The effects of changing the ratio of injection to withdrawal rates are shown. INTRODUCTION A small-scale field test, or pilot flood, is an accepted method for evaluating the applicability of a secondary recovery operation and the economic potential of that operation for a reservoir. The pilot flood is carried on in the same manner, with the same injected materials and at the same pressures as would be used were the secondary recovery program to be expanded to a larger area of the field or to the entire field. If the reservoir fluid and reservoir rock properties are the same in the pilot area as in the rest of the reservoir, the production, injection and other pertinent data associated with the pilot test will represent and measure what should be expected from a full-scale development of this secondary recovery program. There is one major exception which, unfortunately, has been overlooked by many and which has been awaiting quantitative definition. This is the sweep-out pattern the reservoir engineer must use in interpreting the pilot flood and in extrapolating it to the full-scale development. The sweep-out area in a pilot flood is not contained within the five-spot pattern. Further, there is oil flow into and/or out of the five-spot during a pilot test. Thus, production data from pilots must be interpreted with this different flow system and sweep pattern in mind. We have studied the sweep-out pattern that will exist for some conditions of operation of pilot flooding a single five-spot (four injection wells and one producer). In our study, reported here, we have assumed that the reservoir is relatively homogeneous and of uniform thickness in the pilot area. The results of this study en-able the engineer to predict the recovery, injection requirements and project life for a full-scale operation from the pilot flooding data. The effects of various injection to withdrawal-rate ratios and mobility ratios in a relatively gas-free (10 per cent pore volume or less) system are shown. Experimental Method The technique used in this study was the same as that for the X-ray shadowgraph studies of sweep-out pattern efficiencies.' The model of the pilot and surrounding area was made from a uniform packing of round-grain sand consolidated with an epoxy resin. Fig. 1 shows the relation between the model and the pilot area. Only one-eighth of the pilot area had to be represented in the model because we assumed equal injec-

By C. C. Templeton, J. C. Rodgers

A key step in feed water treatment for generating wet steam for thermal oil recovery is the removal of calcium and magnesiunt hardness by cation-exchange series softening. Knowing the solubility of any scale forming salts in brines at elevated temperatures is necessary for fixing the level to which the feed water must be softened. Such calcium sulfate solubility data, previously not available above 392F, were determined by the authors in a flow equilibrium apparatus mud will be reported elsewhere. These data were used to develop a method for predicting the solubility of anhydrite in hot water or steam droplets for saturated steam pressures as high as 2,000 psig (637F). (The calcium sulfate solubility product is represented by a combination of two factors, one reflecting the effects of ionic strength and the other accounting for the effects of complex ion formation in either calcium-magnesium-rich or sulfate-rich brines.) The method is applied to a calcium-magnesium-rich brine If moderately high salinity from a pilot hot-water flood, I nd to several sulfate-rich, low-salinity feed waters and l lowdown (cooled steam droplets) samples from steam s ak operations. The predicted calcium hardness levels corresponding to the calcium sulfate solubilities agreed reasonably well with the results of laboratory solubility determinations run on the field samples. Further testing of the method is needed for brines of other composition classes. Existing field cation exchange softeners in the cases tested are performing adequately since all the samples were found to be undersaturated with respect to calcium sulfate at their operating temperatures. Introduction Prevention of scaling caused by precipitation of calcium sulfate (anhydrite) is of considerable concern in connection with thermal recovery processes using wet steam or hot water. To avoid anhydrite precipitation in a heated system, an engineer must keep the product of the calcium and sulfate concentrations in the water or steam droplets below the value of the solubility product of anhydrite for the temperature and brine composition in question. Usually it is most practical to keep the concentration product lower than the solubility product by keeping calcium low in the presence of high sulfate, or by keeping sulfate low in the presence of high calcium. This can be done by a choice of combinations of natural waters and water treatment processes (such as series cation exchange softening to remove calcium). Until recently, few anhydrite solubility data, particularly for solutions containing other salts, were available for temperatures above 392F (211 psig steam); Marshall, Slusher and Jones' studied the CaS0,-NaCI-H,O system up to 392F and surveyed the work of previous investigators. To model natural brines, one needs to study the solubility of anhydrite in aqueous solutions of sodium chloride, sodium sulfate, calcium chloride, magnesium chloride and their mixtures. Since steam pressures as high as 2,000 psig (637F) may be involved in thermal oil recovery projects, a solubility study was conducted between 482 and 617F.' Discussed in this paper is the application of these data to the prediction of anhydrite precipitation in some practical steam soak and hot-water injection projects. Any simple method for predicting the solubility of an inorganic compound over a wide range of temperatures and solution compositions must be based on some assumptions, and therefore must yield approximate results. On the one hand, natural brines contain too many ionic species for all to be included in a simple scheme; on the other hand, there is no adequate theoretical basis for the exact prediction of solubility in even simple solutions of mixed electrolytes. However, it is possible at a given temperature to base a reasonable prediction scheme on two phenomena:'-' the increase in solubility with increasing total concentrations of all ions (as measured by the ionic strength; see the Appendix), and an increase due to formation of cornplexes between calcium ions and sulfate ions and between sulfate and magnesium ions. Stiff and Davis' developed such a scheme for predicting the solubility of gypsum (CaSO, ¦ 2H,O) in brines at temperatures u~ to 212F. However, for the higher temperaturei of present interest, the stable solid phase is anhydrite (CaSO,). This study involves a two-part method for predicting anhydrite solubility products. First, one predicts a value K, for a given ionic strength I and a given temperature T corresponding to mo./msr, = 1 from data of the CaS0,-NaCI-H,O system. Second, one determines a group of factors F .= F(Ca) • F(Mg) • F(SO,), where the individual factors account for increases in the solubility product due to complex formation by high concentrations, respectively, of calcium, magnesium and sulfate. Combining the two parts, one obtains the solubility product in molalities as

Jan 1, 1969

By R. J. Blackwell, J. R. Rayne, J. R. Henderson, W. M. Terry, D. C. Lindley

This paper presents the results of a laboratory investigation of the efficiency of water-solvent mixtures in recovery of oil. These mixtures may have the high displacement efficiencies characteristic of solvent floods and the high sweep efficiencies characteristic of water floods. Thus, the water-solvent process may increase the number of reservoirs in which a miscible-type displacement can be used profitably. The experiments on the use of water-solvent mixtures for recovery of oil were conducted to find the general applicability of the process. These studies demonstrated that, in flowing through sands, water and solvent segregated into a solvent layer on the top and a water layer on the bottom rather than flowing through the sands as a uniform mixture. Calculations based on the simultaneous flow of the water and solvent in layers were used to predict the effective mobility of the mixtures and the optimum operation of the process in steeply dipping, homogeneous reservoirs. As most reservoirs are not suited for the operation of the process under ideal conditions, experimental studies were conducted with sand-packed models scaled to represent more realistic reservoirs. These studies included the effects on recovery of oil of rate of injection, viscosity of oil, variations of permeablity within a formation and variations in water-solvent ratio. For the range of con-ditions studied, higher recoveries of oil were obtained with water-solvent mixtures than with water or practical volumes of solvent alone. INTRODUCTION A group of intriguing—because of their great possibilities—new oil recovery methods at the disposal of the petroleum engineer are the miscible displacement processes. These processes (high-pressure gas drive, en-riched-gas drive and LPG bank driven by methane) displace all of the oil from the portions of the reservoir swept by the injection fluid. The key question confronting the engineer applying one of these techniques to a particular reservoir is, "What fraction of the reservoir can be swept by injection of a practical volume of solvent?"." Intensive laboratory studies have been made during the past several years in seeking answers to the question of the sweep efficiencies which can be expected in solvent floods.' - These studies provided the answer that low-viscosity, low-density solvents channel and by-pass oil in sands with no dip. In horizontal sands, solvent flooding becomes less efficient as the viscosity of the oil increases, the recovery of oil at solvent breakthrough decreases and larger volumes of solvent are required to achieve a given recovery. More efficient displacement of oil by solvent is observed under certain conditions in sands with dip.".' If the permeability and dip of the sand are sufficiently high, gravity segregation of low-density solvent injected updip can prevent channeling. At rates of depletion below a critical rate,"." no channeling occurs. Unfortunately, in many reservoirs, the critical rate is so low that production of oil at rates below this rate is not economically attractive. And at rates over four times the critical rate, channeling is almost as severe in sands with dip as in horizontal sands. These findings pointed out that, for solvent floods to be generally applicable in recovering oil from all types of reservoirs, new methods of improving their sweep efficiencies are needed. Simultaneous injection of water with the miscible fluid was suggested by Caudle and Dyes7 as a method for improving sweep efficiencies. They theorize that water flowing with the solvent would decrease the effective mobility of the solvent and cause it to contact more oil sand. If, indeed, the water and solvent flowed as a uniform mixture, the process should have the advantages of the high displacement efficiencies characteristic of miscible floods and the high sweep efficiencies of water floods. Thus, the method (at least in theory) would greatly increase the number of reservoirs in which miscible-type displacements would be feasible. A research program was conducted to probe the technical feasibility of the water-solvent process. It

By C. S. Matthews, J. L. Mahaffey, W. M. Rutherford

This paper gives results of an experimental study of the sweep efficiency of a miscible displacement in a five-spot. The study was carried out in a parallel-plate glass model in which effects of diffusion were scaled at or near the molecular diffusion level. The experiments show that very early breakthrough (25 to 35 per cent of pore volume (PV) injected) may be expected in miscible floods because of the unfavorable viscosity ratio. However after 1 PV of displacing fluid is injected, the sweep rises to a reasonable value (50 to 60 per cent). Photographs show that small slugs of less than 10 per cent of PV tend to dissipate before breakthrough. A minimum slug size of 15 per cent of PV would appear to be necessary even in a relatively homogeneous formation. Presence of a slug whose viscosity is inter-mediate between that of oil and gas increases the sweep efficiency of the oil-gas system. In a typical system the sweep at breakthrough rose from 26 to 37 per cent of PV for a 25 per cent slug. The increase in sweep brought about by use of a large slug could well pay for the extra deferment cost of the additional slug material. INTRODUCTION Most miscible displacement processes involve the displacement of oil with fluids of much lower viscosity and density. The displacement process at these adverse viscosity and density ratios is dominated by instability phenomena, i.e., viscous fingers and gravity tongues. These phenomena have highly adverse effects on oil recovery. Although a number of laboratory studies have been made to determine the effect of adverse viscosity ratios on five-spot sweep patterns,l,2 the scaling of diffusion effects is uncertain. In the series of scaled model studies reported herein, an attempt was made to scale diffusion. Model studies of miscible displacements in which molecular diffusion predominates are per-mitted by controlling the parallel plate spacing which reduces convective mixing to arbitrarily small levels. To decide how this scaling relates to any particular Geld displacement necessitates an estimation of diffusion effects for the natural rock being considered and conditions under which displacement will be conducted. The approach normally taken is to extrapolate data obtained from stable miscible displacements performed in the laboratory, such as those presented by Brigham et al.3 The validity with which such an extrapolation can be applied to an unstable flow system has yet to be established. If this approach is accepted, a family of oil recovery curves can be generated for a single viscosity ratio based on Brigham's observation that the magnitude of the dispersion coefficient is dependent, among other things, upon specific rock properties. Objective of our test was to define the lower limit of this range by presenting the case where dispersion effects were reduced to the molecular diffusion level in both the transverse and longitudinal directions.

Jan 1, 1967

By A. B. Cook

Gas cycling is generally considered a much less efficient oil recovery mechanism than water flooding. HOWever, recoveries from some fields have been exceptionally high as a result of gas cycling. Recovery from the Pick-ton field, for example, was calculated to be 73.5 perceni of the stock-tank oil originally in place. In evaluating pressure maintenance projects, determining how much of the recovery is due to displacement by gas and determining how much is due to vaporization of the imrnohile oil in the flow path of the cycled gas is very difficrilt. Even though most of the oil is recovered by displacetr~ent, the success of a project may depend on the amount of oil vaporized. A limited number of experiments have heen performed with a rotating model oil reservoir that simulates gas cycling operations and allows a separation of the oil from, tile free gas flowing into the laboratory wellbore at reservoir conditions, thus revealing which is displaced oil and which is vaporized oil. It Iras been determined that the amount of varporizatio'n is .significant if proper conditions exist These experiments show that oil vaporization depends on pressure, temperature, volatility of the oil and amount of gas cycled. Increases in each of these conditions increase the volume of oil vaporized. Data from six experiments affecting vaporization are presented to illustrate reservoir condition that range from favorable to unfavorable. 111 these eaperitnenis recovery by vaporization ranged from 73.6 to 15.3 percent of /he immobile oil (oil not produced by gas displacerrlt). INTRODUCTION Between 1930 and 1950, gas cycling was a popular. oil recovery practice. especially for the deeper reservoirs. Later, with many case history-type studies published for both gas cycling and waterflooding, it was generally believed that waterflooding was far superior to gas cycling, even when gas cycling was conducted as a primary production procedure by complete pressure maintenance. A good example illustrating the advantage of water-flooding over gas cycling is given in a paper by Matthews' on the South Burbank unit where gas injection was followed by waterflooding. The author concluded in part that "Early application of water injection, without the intervening period of gas injection, would have recovered as much total oil as ultimately will be recovered by waterflooding following the gas injection, and total operating life would have been shortened". This appears to be a logical conclusion. However, it should not be applied to all fields. Pressure maintenance with gas in the Pickton field, as reported by McGraw and Lohec;' will result in a much larger percentage of oil recovery than was obtained in the South Burbank unit. The great success in the Pickton field resulted partly from vaporization of the immobile oil in the flow path of the cycled gas. The amount of vaporization is related to the following conditions: volatility of the oil as reflected by the APT gravity of the stock-tank oil; reservoir temperature; reservoir pressure during gas cycling; and the amount of gas cycled. Therefore, the U. S. Bureau of Mines is investigating these effects on vaporization in a research project using a model oil reservoir. Three different stock-tank oils having 22, 35 and 45" API gravities are being used as base stock to synthesize reservoir oils. Experiments are being performcd to determine vaporization at 100, 175 and 250F and at 1,100, 2,600 and 4,100 psia. This is a progress report showing the results from six experiments. Other Bureau of Mines reports"- concerning vaporization are listed. LABORATORY EQUIPMENT AND PROCEDURES The equipment ' consists of an internally chromium-plated steel tube packed with finely sifted Wilcox sand. The tube is approximately 44 in. long and has an ID of 13/4 in. The sand section contains approximately 570 ml of voids, has a porosity of 32 percent, and a permeability to air of 4.3 darcies. A unique feature of the laboratory reservoir (Fig. 1) permits the tube part to rotate at 1 rpm while the outlet and inlet heads are held stationary. The outlet end contains diametrically opposed windows to permit observatlon of the flowing fluids, and two valves, one on the top and the other at the bottom. Oil and free gas. when being produced simultaneously, can be separated by manipulating the two valves to keep a gas-oil interface in view through the windows. Thus, only gas is produced through the top valve and only oil flows through the bottom valve. The laboratory equipment was designed to study vaporization. Therefore, a uniform reservoir was made using dry sifted sand as opposcd to using a consolidated sand core with interstitial water. Furthermore. the reservoir was tilted to minimize fingering of gas. This tilting also in-

By H. S. Price, R. L. Reed, J. E. Warren

The behavior of the reverse-combustion process in a linear adiabatic system is theoretically investigated by means of an idealized physical model. T his model is described by a pair of non-linear equations involving heat and mass transfer which are coupled by a concentration-dependent reaction-rate function of the Ar-rhenius type. The differential equations governing the quasi steady-state temperature and concentration distributions are approximately solved by using either physical or mathematical simplifications. It is shown that the reverse -combustion process can be mechanistically described by simple physical models whose behavior equations can be solved formally. It is further demonstrated that the derived reverse-combustion equations can be solved numerically within the error limits for the experimental data. The theoretically predicted peak temperature-flux-velocity relationship exhibits reasonable agreement with experimental data obtained from Athabasca tar sand. The principal contribution of this work lies in its usefulness as a starting point for developing a more comprehensive theory. A description of the vaporization-coking process must be obtained; then, the theory must be extended to three dimensions to include geometrical effects and heat losses. INTRODUCTION In another paper,'" he physical processes visualized as occurring during reverse combustion have been described. In this paper, the extent to which an elementary theory can account for experimental observations is determined. Theoretical results that can be used to estimate some of the behavior characteristics of importance to this process are also presented. Since it is not possible to include all the features of the combustion mechanism of which we are currently aware, the mathematical formulation refers to the highly idealized system developed herein. This simplified physical model is based on the simultaneous transfer of mass and heat accompanied by a chemical reaction: it is described mathematically by a pair of coupled, nonlinear differential equations. Similar sets of equations have been studied extensively in the field of flame propagation; an excellent review of these investigations has been published by Evans.1 FORMULATION OF THE PROBLEM It will be recalled that the experiments are performed with a thin-walled cylindrical tube that is uniformly packed with sand and oil. The system is maintained substantially adiabatic by means of heat supplied externally in such a way that the radial temperature gradient in any plane normal to the axis of the tube approaches zero. Initially, the tube is at ambient temperature, except at one end where it is rapidly heated to a predetermined temperature. When the prescribed temperature is achieved, air is caused to flow into the cold end of the tube. As the oxygen in the air stream contacts the hot oil, a localized exothermic reaction takes place. The generated heat is conducted and convected away from the reaction zone so that definite temperature and concentration profiles are rapidly established and move uniformly in a direction opposite to that of the air flow (Fig. 1). As a consequence of this mode of oxidation, all of the oxygen in the air is consumed; then, the remaining nitrogen and the gaseous combustion products sweep out all vaporized hydrocarbons and leave an immobile hydrocarbon-coke on the sand.