Search Documents

By J. C. Martin

Stmplified equations are developed for the flow of fluids in gas drive reservoirs in which the effects of gravity can be neglected. The results show that the pressure distribution is governed by a nonlinear heat flow type of equation, and the saturation distributions are related to the pressure by the same equations QS those developed by Muskat for the average pressure and saturations in gas drive reservoirs. Unrler certain conditions the equation for the pressure can be approximated by the linear form of the heat flow equatlon. This equation tk analogous to the equation for the pressure in single-phase compressible flow The analysis reveals that the equation for the pressure in n~ultiphase flows can be obtained directly from the equation for the pressure in slngle-phase flows by simply replacing the single-phase compressibility by the total compressibility and the single-phase mobility by the sum of the mobilities of all the fluids present. Perrine's method of pressure buildup analysis for multiphme flows is based on these same substitutions. Thus, the development presented in this paper constitutes a theoretical justification for Perrine's method of pressure buildup analysis. One of the most difficult problems associated with the production of oil is the determination of the flow of the fluids within the reservoirs. One purpose of this paper is to present the development of simplified equations for the flow in gas drive reservoirs. The other purpose is to present some of the needed theoretical foundation for the methods currently being used in the analysis of multiphase pressure buildup data. The analysis presented in this paper makes use of a partial linearization of the equations governing multiphase flows. It is felt that the process of linearizing or partially linearizing the nonlinear equations of multiphase flows can lead to useful results since this process has proved extremely valuable in the study of other fields of science and engineering where nonlinear equations are encountered. Pressure buildup data of shut-in wells have been used for many years to estimate the static reservoir pressure. These data have also been used in studying the aver- age permeability both away from and close to the well-bores. The majority of the buildup data are obtained from wells producing two or three fluids. As yet the theoretical foundations for the methods of buildup analysis used in connection with rnultiphase buildup data have not been established. A number of methods of analysis of pressure buildup data have been developed." Most of these methods apply only to single-phase flows, and almost all of the theoretical analysis is confined to single-phase flows. Perrine' developed a method of analysis for rnultiphase flows from the method for single-phase flows of Miller, Dyes, and Hutchinson.' He did this by simply replacing the single-phase compressibility by the multiphase compressibility and replacing the single-phase mobility by the sum of the mobilities of the fluids in the rnultiphase flows. Obviously, such a procedure lacks a certain amount of theoretical justification, and such practices can lead to erroneous results. Fortunately, Perrine's results have a theoretical justification, and this justification is presented in the following discussion. It is assumed that gravitational effects and the cornpressibility of the rock can be neglected. Under these assumptions the equations governing the simultaneous flow of gas, oil and water in porous media areo* The preceding equations can be expanded into the following expressions.

By Wilbur H. Somerton

In order to evaluate the quantitative determination of porosity by various neutron logging systems, four test wells were constructed to simulate some of the possible variab1e.r found in field logging conditions. All of the test holes contain sections of limestone with one section of sandstone included in one well. Each well has a representative section of high, medium, and low porosity. Three borehole sizes are represented by three of the test wells, all containing fresh water. One test well contains saturated brine. Statistical and tune constant checks are made for each system evaluated. Runs and repeat runs are made with the sonde centered and eccentric. Factors evaluated are the effect of borehole size, sonde position, salinity, chemical composition of rock, and statistical variation. A semi-quantitative method of evaluation of these effects is used giving the average maximum per cent error possible for each. Drift and calibration procedures are also discussed. The effect of the various factors appears to be largely a function of the phenomena detected as well as the means of detection. Graphic illustrations of some of these effects are presented. It is concluded, that for accurate log interpretation, the hole size and sonde position should be known. INTRODUCTION The use of neutron logs for the derivation of porosity has become an accepted practice in many areas. In recent years, many different neutron logging systems have been offered to the industry. The neutron sources presently available are of the radium-beryllium or polonium -beryllium type. Detector types used are ionization chambers, Geiger-Mueller-type counters, proportional counters, or scintillation detectors used to detect or measure gamma rays, thermal neutrons, epithermal neutrons, fast neutrons or combinations of these. It is well-known that all systems do not respond alike in all environments. As a result, testing facilities and an evaluation method were developed for use in evaluating the porosity discriminating ability of neutron logging systems using conditions simulating those normally encountered in the field. With such an evaluation, the suitability of a logging system can be determined lor selected field conditions. TEST FACILITIES The testing facilities consist of wells containing rocks of various porosities, constructed to include as many parameters as practicable. Three pits were drilled to a depth of 20 ft and the fourth to a depth of 9 ft, all having 48-in. diameters. Rat holes were dug to a depth of 6 ft in the center of each pit and fitted with sections of pipe having appropriate diameters. The top of this pipe was located approximately 8 in. above the pit bottom. A 6-in, layer of cement was then poured on the bottom and while still wet, 42-in. OD corrugated, galvanized culvert pipe was set on bottom and centered in place. This pipe extends approximately 1 ft above the surface of the ground. The annulus around the culvert pipe was then filled with cement to the top. A 1- or 2-in. layer of sand was placed on the cement bottom to provide an adjustable bed on which to place the series of rocks set in each hole. A cross section of a typical test well containing the porous rocks is shown in Fig. 1.

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

The injection of heat-bearing fluids may offer a wider application to secondary and tertiary recovery from conventional oil reservoirs than inderground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. In all cares, conduction heat losses to overburden and un-derburden will impose an economic limit upon the size of the area which can be swept out from any one injection point for any given set of conditions. The present report describes a method for c.rtin~ating thermal invasion rates, cumulative heated area, and theoretical economic limits for sustained hot-fluid injection at a constant rate into an idealized reservoir. Full allow-ance is made for non-productive reservoir heat losses. Required solutions can be obtained directly from the numerical tables given here, as soon as all operating conditions are specified, withou extensive mathematica1 ~nanipulation on the part of the user. The injection of a heat-bearing fluid into a reservoir is frequently proposed as a means of secondary or tertiary oil production. This form of thermal recovery may offer a wider application to conventional oil reservoirs than underground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. Because of its large gross heat capacity, steam appears to be the most efficient heat injection medium, although mixtures of steam and other gases, hot water, hot oil and hot non-condensable gases also have been employed. Whatever the nature of the injection medium, economic evaluation of such thermal recovery processes will depend upon thermal invasion rates, or rates of productive reservoir heating, at any given time after the start of heat injection. In all cases, conduction heat losses to overburden and underburden will impose an economic limit upon the size of the area which can be swept out from any one injection point, for a given set of reservoir conditions, at any given heat injection rate. The present report describes a method for estimating thermal invasion rates, cumulative heated area and theoretical economic limits for sustained heat injection at a constant rate into an idealized reservoir, making full allowance for non-productive reservoir heat losses. Required solutions can be obtained directly from the nu- merical tables given here, as soon as all operating conditions are specified, without extensive mathematical manipulation. It should be noted that ordinary water flooding may be used to move an annular-heated zone outward, even though the economic areal limit for straight heat injection has been reached. Thus well spacing need not be confined to the areal limits as determined here. However, the present calculations do not cover such supplementary heat scavenging programs. The operating element consists of a radial flow system, concentric about the point of injection, with the idealized step function temperature profile illustrated by the dashed line in Fig. 1. The solid curve shown in Fig. 1 gives a qualitative picture of the true radial temperature distribution. It will be shown later that the expressions given below, which are derived from the simplifying step idealization (dashed line in Fig. I), also can be interpreted in terms of the more realistic temperature distribution shown in Fig. 2. For this case, the thermal invasion radius is simply redefined as the distance from injection well to the mid-point of the temperature distribution.

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By D. R. Parrish, V. J. Berry

Reverse combustion is one thermal method of recovering hydrocarbons from porous undergrortnd formations containing oil or tar. In applying this method, air is introduced via an injection well and the mixture of air and hydrocarbons is ignited in the production well. A combustion zone then recedes toward the injection well, counter-current to the air flow. If the combustion-zorie temperature is sufficiently high, the oil or tar irz place is distilled arid cracked. The hydrocnrhon flows as a vapor to the production well and subsequently is condensed at the surface. Maximum temperature and velocity of morcment arc the two dependent variubles defining the progress of the combustion zone. A theoretical analysis has been made of heat flow in the reverse-comhustion process as.srtrning linear flow in a homogeneous system. The differentirll equatious, which include the oxygen-hydrocar-hon rctrction rate, have been solved numerically. Rescllts indicate that the maximum temperature reached and he cornbustion-zone relocity both increase with an incrense in air-injection rote. Heat loss to surroundings has little effect on the maximum comhustion-zone ten?perature achieved, but it is reflected in a reduced com-bustion-zorzc, velocity. It is also predicted that an increase in the oxygen-hydrocarbon reaction rate results in a reduction in the maximum temperature reached. The calculated results are in agreement with results from reverse-combustion experinments using samples of a tar sand. INTRODUCTION Reverse combustion is one thermal method of recovering hydrocarbons from porous underground forma-tions containing oil or tar.' In applying this method, air is introduced into the underground formation via an in-jection well. In one or more adjacent production wells, the mixture of air and hydrocarbons is ignited at the sand face. The combustion zone thus formed recedes toward the injection well, counter-current to the air flow (Fig. 1). Under certain conditions, the rate of recession of the combustion zone is sufficiently rapid that only a small fraction of the hydrocarbon at any point in the formation is consumed by combustion. If the combustion-zone temperature is sufficiently high, the uncon-sumed oil or tar in place may be partially distilled or cracked. This material then can flow as a vapor with the gaseous combustion products to the production well and subsequently be condensed at the surface. Under other conditions, the combustion-zone temperature may be so low that little, if any, of the hydrocarbon is vaporized. Nevertheless, the viscosity of the hydrocarbon may be reduced enough so that other recovery methods subsequently can be applied to produce oil at attractive rates. The maximum temperature achieved in the combustion zone and the velocity of movement are the two dependent variables defining the progress of the combustion zone. Presumably. knowledge of these factors will permit estimates to be made of oil recovery efficiency, oil and gas producticn rates and average amount of injection air required per barrel of oil recovered. The maximum temperature and combustion-zone vclocity are, in turn, expected to be dependent upon air flow rate, heat losses to surrounding formations, physical properties of the rock and the hydrocarbon. the extent to which the porous formation is saturated with the hydrocarbon, the rate of oxidation of the hydrocarbon and the amount of heat released. It is desirable that the relative importance of each of these factors be understood so that the results of laboratory tests of oil recovery by reverse combustion can be properly interpreted in terms of expected field performance. It is particularly desirable to assess the importance of heat lcsses because these are difficult to control or eliminate

By C. R. Johnson, R. A. Greenkorn, G. W. Widner

This paper describes a short, simple well-testing procedure evolved from unsteady-state concepts to appraise wellbore damage in water-input wells. Input rates are changed at specific intervals over a 30-minute period and the wellhead pressure recorded at the end of each interval. From these data, the pressure drop due to the "skin" can be computed. This pressure drop in terms of psi/B/D of input rate is referred to as the "skin factor". Two of these variable-rate tests are required to determine change of wellbore damage and to evaluate well-treatment techniques. There are several advantages to the variable-rate procedure. The necessary data can be obtained with wellhead gauges and meters. There is little interference with normal operations. The short testing period eliminates many interpretative difficulties that come from interference from adjacent wells. The shut-in period required for standard tests before the plot of pressure vs log of time attains a straight line is often 12 hours or more. Such a time period is incronvenient, especially when a well must be tested repeatedly to evaluate change in wellbore characteristics. This variable-rate procedure has been successfully npplied to water-input wells in three fields. For some of the wells, skin factors were calculated from both the standard and the variable-rate tests. These results show good agreement. THEORY OF THE VARIABLE-RATE PROCEDURE In many waterflood input wells, the pressure in the wellbore differs from the pressure in the reservoir immediately adjacent to the well because of the restriction of flow through the skin of the well. If the true pressure of the reservoir at the well, p,, is that pressure which would exist in the well if no skin were present, and p,' is the actual pressure as measured in the well, this difference, is the skin effect. Fig. 1 shows the change in skin effect with change in flow rate. As the flow is reduced, the skin effect is reduced by an amount equal to the change in wellbore pressure minus the change in reservoir pressure. The equation of pressure change without skin using the point-source solution of the diffusivity equation is (after van Everdingen')

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

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

By S. A. Hovanessian

The method of characteristics which reduces the Buckley-Leverett mnlineur, partial differential equation to two first-order, linear differential equations is applied to the solution of a waterffood recovery problem for multiple lines of producing wells and one line of injection wells. The application of the method is described by a numerical example, and a brief theoretical background is given. From the results of this study, it is concluded that the cumulative water injected and oil produced are the same for single and multiple sets of producing wells. INTRODUCTION In this paper, the method of characteristics is applied to the solution of a waterflood recovery problem with multiple lines of producing wells and one line of injection wells. The application of the method to this type of recovery problem considerably simplifies the procedure of obtaining a numerical solution. It differs from the method described by Kern1 because it eliminates the need for integration in obtaining saturation fronts. In addition, the saturation history at any location within the field is obtained from a single diagram (time vs distance with saturation as a parameter)'. The paper includes a brief theory of the method of characteristics and its application to the solution of oil recovery problems. A numerical example is worked out, and the cumulative oil recovery-vs-time curve for a system with multiple sets of producing wells is obtained. This curve is then compared with the recovery curve of a system with only one set of producing wells. From this study, it is concluded that the cumulative water injected and oil produced are the same for single and multiple lines of producing wells. However, the recovery time is considerably longer for multiple sets of producing wells.

By R. Siy, O. Neiman, B. Hilscher

"A full-scale prototype of “Jameson Downcomer” technology has been installed and tested on a Tailings Oil Recovery (TOR) vessel in Syncrude’s Mildred Lake extraction plant. Jameson Downcomers are a type of mineral flotation equipment used in the coal industry that Syncrude has extended to the recovery of secondary bitumen in the oil sand extraction process. The TOR vessels are currently being operated as single-gravity separation vessels, in series with the Primary Separation Vessels (PSV’s). We have developed a model to differentiate between “floating” and “non-floating” bitumen in the PSV’s, and then applied it to verify that the current TOR performance only captures the “floating” bitumen by simple separation. The existence of “non-floating” bitumen provides an opportunity for recovery enhancement by means of the Jameson Downcomers. To achieve bitumen recovery enhancement, middlings slurry is withdrawn from the TOR vessel, and passed through a number of parallel Jameson downcomers, within which the slurry is accelerated to a high-velocity free jet, and impinged upon the slurry surface inside the downcomer, in the presence of air. The aerated slurry is returned to the middlings zone inside the TOR, resulting in an incremental recovery improvement of bitumen in the TOR. The operation of the downcomers is thought to produce suitably small air bubbles for the flotation of bitumen, either by fluid shearing during the slurry jet impingement, and/or by micro-bubble nucleation. This full-scale prototype design has met or exceeded the recovery expectations of the original Research piloting, recovering at least 40% of the non-floating bitumen. In Syncrude Extraction, this provides high business value by a separation efficiency gain achieved with modest expenditure. As well, the system has demonstrated excellent operability and maintainability.INTRODUCTION Syncrude Overview Syncrude Canada Ltd. is one of the largest producers of crude oil from oil sands. Our crude oil production facility has the capacity to produce over 15% of Canada's total oil requirements. To do this, we surface mine oil sand, extract the raw oil known as bitumen from the sand using water-based processes, and upgrade that bitumen into sweet light crude oil by fluid coking, hydro processing, hydro treating and reblending. Our final product, Syncrude Crude Oil (SCO), is sent by pipeline to three Edmonton area refineries and to pipeline terminals which ship it to other refineries in Canada and the United States. The Syncrude operation is comprised of four major technology areas: Mining, Extraction, Upgrading and Utilities. As well, Syncrude invests more than $40 million annually in science and technology, and is among the top 50 companies in Canada for Research and Development (R&D) investment. Syncrude holds 21 active Canadian and U.S. patents."

Jan 1, 2012