Search Documents

By C. E. Kemp, A. B. Dyes

Results of a field research project on the thermal recovery of oil by movement of a combustion front are presented. This field test was conducted in the South Belridge field, This The war paitern was a Southsing1e 2.5-acre five-spot pattern in the 700-ft deep Tulare sand. Results of this project indicate the technical feasibility of producing high viscosity, low-grallity crude oils by thermal means. It was demonstrated that heavy oils could be moved rapidly over considerable distances by thermal drive, and that high percentage recoveries of oil could be effected. The bulk of the oil recovery came after arrival of the burning front at a production well. However, rapid corrosion and high-temperature failure of conventional pipe and equipment both below and above surface occurred. Stainless steels gave more satisfectory performance under these conditions than did conventional steels. Oxygen was seldom observed in gas from wells ahead of the burning front, and it appears that the front was continuous throughout the life of the experiment. INTRODUCTION During late 1953 and early 1954, industry-wide attention was focused on the thermal recovery of oil by movement of a combustion front through an oil zone. Two separate Oklahoma field tests were described in publications by Kuhn and Koch, and Grant and Szasz. To obtain basic information needed for making engineering and economic appraisals of commercial recovery of heavy oil by this thermal recovery method, a new field research projecta3,4 was started May 31, 1955. Twelve oil companies contributed to this project with General Petroleum Corp. as operator. Since combustion recovery experiments* with heavy oils in other areas had been conducted in small patterns,"' an important purpose of this test was to determine whether the process could be operated in a pattern of near-commercial well spacing. The south Belridge field, located 30 miles north of Taft, Calif., was chosen because its oil sands are similar to many heavy oil sands which typically have low primary recovery and leave a large amount of oil in place. such fields are not particularly suited for common secondary recovery methods. DESCRIPTION OF TEST SITE The test interval was a Tulare sand at a depth of about 700 ft. Nine test wells were drilled, cored and completed from Oct., 1955, through Jan., 1956. The initial pattern shown on the map, Fig. 1, consisted of an injection well (II)', four producing wells (1P-4P), and four observation wells (IT-4T). An old well (81) was also recompleted for use as an observation well. The area enclosed by the original four production wells was 2.5 acres; the distance from the injection well to each production well was 233 ft. During the combustion period, Well 5P was drilled to replace Well 2P, and Well 6P was drilled to replace Well 3P. The new wells were 25 ft beyond the old production wells along a diagonal from the injection well, and resulted in an increase in pattern area to 2.75 acres. Isopachs are shown in Fig. 1. The average thickness of the oil sand in the test pattern was 30 ft. Fig. 2 presents two cross sections through the test pattern. Oil sand is represented by white and claystone by gray. The dashed lines indicate layers of conglomerate grading to claystone. The sand thickness was fairly uniform throughout the test pattern except in the 1P area. The sand body dips gently from the northwest to southeast, about 3 degrees from horizontal. The sand is unconsoli-dated, and grades from very fine material to pebbles as large as 2 in. in diameter. Table 1 presents the average initial properties of the formation. The in-place porosity of 37 per cent was estimated from the porosity of coke-consolidated cores taken after combustion. Original core analyses data were adjusted to 37 per cent, giving pore volume saturations of 60 per cent oil, 37 per cent water, and 3 per

By L. Yarborough, L. R. Smith

This paper presents results of a laboratory study of retrograde condensate recovery by revaporization into dry injection gas. Flow tests were performed in 10.6-ft long sand packs at 100F and 1,500 psi. In three runs methane revaporized the liquid from a n-pentane-methane mixture in the presence of immobile water. TWO of these tests were water-wet, and the third was totally oil-wet. In the three runs n-pentane recovery was complete after 2.5 hydrocarbon PV of injection. There was no significant performance difference between the two wettability extremes. In a fourth experiment, a methane-hydrogen sulfide mixture revaporized a synthetic light, sour condensate. No water saturation was present. Equilibrium compositions and volumetric data were obtained for the four-component condensate. The heavy component, n-beptane, was removed after 6 PV production. Comparison of the effluent fluid compositions with known equilibrium data shows that the flowing fluid was equilibrium vapor and that the mixing zone between equilibrium vapor and dry injection gas was short. Data indicated that complete recovery of retrograde liquid occurred after it was contacted by a sufficient quantity of dry gas. INTRODUCTION When pressure declines below the fluid dew point in a gas condensate reservoir, a liquid phase forms. In this process, referred to as retrograde condensation, the quantity of liquid formed is frequently small enough that the liquid is not a flowing phase. To prevent loss of valuable retrograde liquids, the process of dry gas cycling has been employed for several years as a more or less standard practice. In this procedure the reservoir pressure is maintained above the fluid dew point so that the liquid components may be produced as vapor and then separated at the surface. Although full pressure maintenance by gas cycling seems ideal in terms of preventing liquids loss, several factors can reduce the atuactiveness of such an operation. From a study of a condensate reservoir in Alberta, Canada, Havlena et a1.1 con-cluded that cycling under conditions of declining pressure leads to economic advantages and to a high recovery of hydrocarbon liquids. This study considered effects of volumetric sweep efficiency, retrograde behavior of the original wet gas and revaporization characteristics of the retrograde liquid when contacted by dry gas. The first major work concerning revaporization of liquid in a gas condensate system is that of Standing et a1.2 Calculations based upon the PVT behavior of a recombined gas condensate fluid indicated that all retrograde liquid can be recovered if it is contacted by a sufficient quantity of dry gas. The paper considered the effect of variable permeability upon the recovery of retrograde liquid. Standing et al. concluded that recovery of heavier components in the retrograde liquid is greatest if reservoir pressure is allowed to decline below the dew point prior to dry gas injection. Since the work of Standing et al., several laboratory studies have been reported which show that recovery of hydrocarbon liquids by vaporization into dry injected gas can contribute to increased recovery above that obtained by ordinary production practices. Vaporization from retrograde condensate,3 conventional oil4-7 and volatile oil8 reservoirs has been considered. There is little work that deals with revaporization recovery from condensate reservoirs. This study demonstrates that vapor-liquid equilibrium exists when a dry gas is injected into a porous medium containing retrograde condensate; investigates the effect of wettability on revaporization of retrograde condensate; and provides an indication of the amount of dry gas necessary to

Jan 1, 1969

By W. A. Roper

Since 1950, several papers have been published which have described various methods for studying mobility ratio effects. The methods which have been described for studying mobility ratio effects include X-ray shadowgraph techniques, stepwise potentiometric model procedures, stepwise numerical procedures, gelatin model studies and fluid mapper model studies. The purpose of this paper is to describe a procedure for using a resistance network for studying mobility effects. DESCRIPTION OF APPARATUS The basic unit of the experimental model consisted of a network of 840 resistors fastened to a sheet of black Lucite which was mounted on a wooden frame, Fig. 1. The dimensions of the Lucite sheet were 36 X 36 in. The resistors were fastened to the Lucite sheet at mesh points by means of snaps. The mesh points were arranged in a square pattern, 30 X 30 in., spaced at intervals of l1/2 in. The snap fasteners permitted several resistors to be "stacked" in parallel. By adding additional resistors to the basic unit, the specific resistance on each side of the flood front could be varied. The square network of resistors represented one of four symmetrical elements of a five-spot well pattern. The current input was analogous to an injection well and the current output was analogous to a production well." The power supply was essentially a full wave transformer-rectifier system for converting alternating-current into direct-current. The DC voltage could be varied from 41.43 to 290 v in steps of 41.43 v by selecting the proper outlet jacks. The major components of the power supply were: a Triro power transformer, R 11-A; a dual 5Y3 vacuum tube; two gas filled OD-3 tubes; and a Triplett voltmeter, model 227-1, range 0 to 500 DC v. The voltages at the mesh points were determined with a Keithley electrometer, model 210, No. 185, manufactured by Keithley Instruments, Cleveland, Ohio. A Heathkit vacuum tube volt-

By C. D. Stahl, S. M. Farouq Ali, W. E. Culham

One- and two-dimensional mathematical models have been developed that simulate transient, two-phase flow of hydrocarbon mixtures in porous media in a manner that accounts for interphase mass transfer. Numerical simulations of one-dimensional depletion-drive experiments using a two-component hydrucarbon fluid were used to establish the validity of the mathematical models. In addition, the experimental and numerical data were used to demonstrate that production rate had a relatively insignificant effect on the recovery of individual hydrocarbon components from the experimental system, and that attainment of equilibrium between phases is possible for a wide range of liquid and vapor velocities in reservoirs containing light hydrocarbon fluids. Results of some two-dimensional numerical simulations are also presented. INTRODUCTION This study was undertaken to develop a mathematical model that would simulate transient, two-phase flow of hydrocarbon mixtures in porous media under conditions that result in interphase mass transfer and to test the validity of the assumptions used to set up the model. In addition, the study was designed to determine if production rate influences the recovery of individual hydrocarbon components from reservoirs producing by depletion drive. Two-phase flow in porous media, with interchange of components between the two phases, is important in many petroleum recovery processes. Studiesl-5 conducted within the last 3 years have outlined methods of solving multiphase flow problems incorporating mass transfer. Some of these studies have also indicated the importance of accounting for mass transfer under various producing conditions. An earlier work6 first demonstrated the importance of combining relative permeability data with equilibrium ratios in compositional balance methods. The mathematical model presented in this paper is formulated so that a phase behavior package, as described in previous papers,5,17 is not required as an integral part of the routine employed to solve for the primary dependent variables. The finite difference formulation is designed so that all the primary variables can be solved for simultaneously. This is accomplished by utilizing one basic set of equations. These innovations, which are in contrast to other models2,3,5 but are similar in some respects to the approach used by Taylor,l render the total problem computationally simpler than any of the previously referenced formulations. The mathematical model was developed by combining Darcy's law with a continuity equation for each hydrocarbon component. The principal assumptions invoked in the formulation were (1) that capillary forces and diffusional effects are negligible, and (2) that thermodynamic equilibrium exists in the reservoir at all times. No assumption as to the type of vaporization process was made in formulating this model. Experimental data were required to complete this study. These were generated by conducting several depletion drive experiments. The experimental apparatus consisted of a sandstone core enclosed in a pressurized casing. The apparatus was designed in such a manner that the core could be charged with a liquid hydrocarbon mixture and depleted at different production rates. The experimental tests were designed to determine the effect of production rate on component recovery. In addition, direct comparison of experimental and mathematically predicted

Jan 1, 1970

By J. L. Huitt, IM. L. Slusser, E. E. Glenn, M. L. Slusser

This paper is concerned with: (1) an analysis and interpretation of the filtration characteristics of drilling muds on filter paper, and (2) an interpretation of early stage "filtration" on consolidated synthetic porous media, herein referred to as the "surge loss" period. It is shown that water-loss data can be plotted so as to disclose the existence of a surge period, a non-uni-form cake thickness period, and a constant pressure filtration period. The same periods occur during exposure of consolidated porous media to muds. The data on consolidated porous media indicate that: (1) mud particles bridge pores inside porous media during the surge period, and (2) "filter cake" build-up inside porous media occurs during subsequent filtration. The data also show that the surge loss increases with sample permeability and with pressure differential. INTRODUCTION The phenomenon of mud filtrate invasion into porous sands of a drilling well has been recognized for many years.1,2,3,4 Part of API Code 29 is concerned with the filtration testing of a drilling mud to yield data on this subject. In a number of investigations "plastering properties" of muds were studied,5,6,7,8 "plastering" apparently meaning cake forming properties. In such studies the characteristics (thickness, permeability, and toughness) of the cake and the rate at which the liquid phase of the mud was separated from the solids (filtration rate through the cake, commonly referred to as water loss) were measured or at least considered. From the standpoint of its effect on well productivity, filtrate invasion into "clean" water-wet pay sands (i.e., those containing no clay minerals) may not necessarily be a serious problem. An increase in water saturation in the zone immediately around the wellbore would cause loss of oil permeability. However, in such a system, much of this water may be removed by oil production. If filtration occurs into "dirty" sands (i.e., those containing appreciable clay minerals), clay swelling or other mechanisms may cause losses in oil permeability which cannot be overcome by oil production.'.' The phenomenon of mud particle invasion into the pay sand of a drilling well was considered to be a possible serious factor in the early days of rotary drilling. Opinions varied as to whether drilling fluids "mudded up" and "sealed" oil sands. In 1932 Rubel,9 Gill,10 and Parsons," in separate papers, presented differences of opinion together with some laboratory and field data concerning loss of oil production caused by drilling fluid invasion into pay zones. These papers, and the discussions, furnish a large background on this subject and indicate the very clear thinking of many men at a time when petroleum production technology was in its infancy. Jones and Babson,6 and Silent,' experimentally found that the depth to which mud particles invaded water-saturated sand packs was small (no greater than 1/16 in). Arn-quist12 in 1937 observed loss of water permeability due to invasion of mud particles into water-saturated, consolidated, sandstone samples. Williams and Cannon5 studied the filtra-tion properties of drilling muds and mentioned "sludge" which penetrated a short distance into the water-saturated filtering medium (sand pack). Byck5 concluded that formation permeability had no effect on the plas-

Jan 1, 1958

By J. H. Henderson, L. A. Wilson, R. L. Gergins, R. J. Wygal, D. W. Reed

This paper presents a method of predicting the production history of an underground combustion recovery process. A rigorous solution of the thermodynamics and hydrodynamics involved is beyond the scope of this report. However, a practical scheme to appraise combustion recovery performance has been Worked out. By prudent assumptions, regarding the attainment of steady-state conditions, a trial-and-error solution, which satisfies both material balance and three-phase relative permeability requirements has been evolved and tested. As a combustion zone moves through a formation the interstitial water, the water of combustion and a portion of the oil wil1 be vaporized. These vapors will move downstream to a cooler region of the formation where they will condense. Part of the remaining oil will be displaced by the imposed gas drive and what remains behind will be consumed as fuel. The fluids will distribute themselves downstrean? in a manner which will satisfy the three-phase relative permeability characteristics of the formation. Three distinct zones will exist ahead of the combustion zone: (I) a zone termed a water bank which contains three mobile phases, oil, water and gas; (2) a region containing two mobile phases, oil and gas, called the oil bank; and (3) a zone having the fluid saturution distribution which existed prior to combustion. The mobile water will impose a water flood on the oil in the region containing three flowing phases. The process can be visualized as a simultaneous water .flood and combustion-supported gas drive. It is possible to estimate the saturation distributions in the water bank and oil bank as well as the production history for any combustion temperuture if the porosity, three-phase relative permeability characteristics, oil viscosity and initial saturations are known. This has been done for injection pressures ranging from I to 100 atm, oil viscosities of 10 to 1,000 cp and porosities from 20 to 40 per cent. Several of these calculations have been checked in the laboratory, with the result indicating rather good agreement between the predicted and observed production histories. The analytical and experimental techniques are described and the comparison of predicted and observed performance presented. INTRODUCTION Considerable research by the oil industry during the past 10 years has been directed toward the thermal recovery of crude oil. Of all schemes considered, the one receiving the most attention is in situ combustion. This is a process in which part of the oil is actually burned within the reservoir rock. The uniqueness of the process lies in its potential. Theoretically, this potential generates from the following: (1) it may open the door to vast reserves previously not economically producible, (2) it may induce increased production rates where low rates are now realized and (3) it may give high over-all recovery efficiency. The general idea of the process has received considerable attention in the literature '-' and the following brief description is intended for orientation only. In situ combustion can be visualized by considering a linear system, one end of which represents the injection well and the other the producing well. The formation at the inlet end is raised to ignition temperature by placing a gas or electric heater at the face of the formation or by the ignition of thermite, charcoal or similar material. While heating, the water and part of the oil in the heated region is vaporized and moved out into the reservoir by the injection of gas containing oxygen. In addition, part of the oil is displaced as liquid by the imposed gas drive. The oil not displaced is modified to an immobile, coke-like material (hereafter referred to as coke) by the high temperatures existing immediately ahead of the combustion zone. This coke is used as fuel, and the combustion front advances through the formation at a temperature which is probably in excess of 800°F. As the combustion front advancek it will drive ahead all interstitial oil not rendered immobile and all interstitial water plus the water formed by the combustion of hydrogen. As a result the sand behind the combustion zone is freed of all oil and water. Although the general description of the process has been outlined in the literature, there has been no consideration of the behavior of the fluids moving downstream from the combustion front. It is the purpose of this paper to present a technique

By J. C. Savins

The practical analysis of the hydrodynamics of the wellbore has long been a subject of interest to engineers. This paper presents a simplified solution to the problem of computing the pressure drop for the flow of drilling mud in the annulus of the wellbore. This solution is, however, an exact and rigorous solution under the assumptions which have been imposed. These assumptions are that the drilling fluid is a Bingham plastic fluid* and that the annulus is formed by two concentric, stationary, cylindrical pipes. It is further assumed that the fluid is incompressible and that its motion is isothermal and in a steady state. This problem under the same assumptions has been attacked by previous authors. Beck, Nuss and Dunn' proposed that the equation for the flow of a Bingham plastic fluid in a cylindrical pipe could be applied to an annulus if the pipe radius were replaced in the equation by the hydraulic radius. This equation, known as the Buckingham-Reiner equation' (see Appendix 1), was also used in an approximate form. Van Olphen pointed out that even for a simple or Newtonian- fluid the pipe equation (Poiseuille's law) could not be converted to the Lamb equation escriptive of flow in an annulus (see Appendix 1) by using the hydraulic radius. Van Olphen further attempted to give a solution for the annular flow of a Bingham plastic fluid by introducing approximations similar to those which have been used in the case of the Buckingham-Reiner equation. Other attempts to provide approximate or exact solutions have been made by Grodde' and by Mori and Ototakeq. The present authors some years ago in unpublished work derived the correct expressions relating the pressure drop and flow rate for this problem. It was found that the solution consisted of two simultaneous equations, one of which contained a logarithmic term. Thus, obtaining numerical results for any particular case of interest involves very tedious trial-and-error computations. Very recently Laird presented the correct derivation of the two equations which are given in full detail in Appendix 1. In order to reduce the amount of calculation time which would be involved in providing a complete tabular or graphical solution to the problem, a high-speed electronic digital computer has been utilized. For this purpose the two simultaneous equations were transformed into more compact expressions by introducing reduced variables. These expressions are given in the following theoretical section. A similar procedure in this problem has been developed by Fredrickson and Bird1" Their tabular results, however, are very incomplete in the range of practical interest for problems of wellbore hydrodynamics. We have furthermore been able to express our graphical results in terms of convenient and familiar dimensionless groups. THEORETICAL DEVELOPMENT Use of Reduced Variables In terms of reduced variables the two simultaneous equations just discussed take the following form, The reduced variables q, x, a and z are defined in terms of the various measured quantities, where Q is volumetric flow rate, AP/L is pressure gradient, Dl is OD of inner pipe, D, is ID of outer pipe, is plastic viscosity, and is yield point. Thus, we have a dimensionless volume flux, a dimensionless reciprocal pressure gradient, and the ratio of the pipe diameters Before introducing the fourth reduced variable, z, it is of interest to consider the physical significance of the parameter x. As may 'be seen from the velocity profile of Fig. 1 the Bingham plastic fluid has the interesting property that a portion of the stream flows at a uniform velocity without shearing action. This section of the stream is situated approximately in the center of the conduit and is known as the "plug flow" region. Its existence is due to the fact that the shearing stresses within the region do not exceed the yield point, which is one of the two flow properties characterizing the fluid. The parameter x then turns out to be the ratio of the

By C. S. Matthews, H. C. Lefkovits

The performance of depletion-type reservoirs at the stage in which the gas is at such low pressure that gravity is essentially the sole driving force has been investigated both theoretically and by use of scaled models. A simple theory is developed which completely explains the scaled model experiments. It also agrees with some field data. When the range of production rates to be considered is less than 25, the theory predicts that the production rate decline curve for homogeneous reservoirs in which there is a free surface will be of the hyperbolic type with n = 2: The theory gives a physical significance to the parameter y. INTRODUCTION This investigation was undertaken to determine a method for prediction of performance of reservoirs where gravity is the sole driving agent. Depletion-type reservoirs fall into this category after the pressure has fallen to some low value. Some examples of this type are found at Healdton ('Oklahoma), Oklahoma City (Oklahoma), Midway-Sunset (California), and East Coalinga (California). Generally the production rate per unit thickness is quite low in such reservoirs, with the conse- quence that this stage is sometimes called the stripper stage. However, per well production rates may be of the order of 50 B/D because of the considerable sand thickness open in some of the reservoirs. In the past, predictions of behavior of such reservoirs have been based on empirical extrapolation of the rate of production to the economic limit. Difficulties arise when new wells are drilled in the field or when production is artificially curtailed, since in such cases the decline rate of individual wells or leases changes and no clear trend can be detected. It was hoped that the present study would lead to improved methods for evaluation of such reservoirs. The study was initiated by building a scaled model for a particular field and by running a number of experiments on this model. Later a simple theory was developed which explained the laboratory results completely. The theory led, as shown below, to a hyperbolic-type decline equation which had been used previously' in a completely empirical manner. THEORY The General Case The problem to be treated here is the following: We are given a homogeneous cylindrical porous reservoir with exterior radius r, wellbore radius r., and height h,, containing, from the bottom up to a height hit oil, some gas, and interstitial water (see Fig. 1 A). The portion of the reservoir between K, and h, is filled with gas at low pressure, with residual oil, and with interstitial water. If at a time t = 0, the level of the oil in the wellbore is lowered to a height h,, what is the production rate from the reservoir at a time to The configuration within the reservoir at some instant

Jan 1, 1957

By C. W. Volek, G. Mandl

Steam-injection tests in the field have shoum that heat transport into the oil/water region, ahead of the steam zone, may have a significant effect on the production process. Earlier theoretical work on reservoir heating by hot-fluid injection can be used to describe the growth of the steam zone, but only if heat transport from steam zone into oil/water region is neglected. It therefore became necessary to re-examine and extend the earlier theory. It was found that this theory ceases to be consistent with the physical model of the process at a certain critical time, tc, which depends on reservoir thickness, temperature, and quality of the steam. The critical time marks an important change in the heat flow across the condensation frovt: the heat flow, which is purely conductive during 0 < t < tc&apos; becomes predominantly convective at t = tc. Accordingly, at t = tc, the equation which governs expansion of the steam zone changes. On the basis of the new equation, an approximate description for the steam-zone growth during t > tc is developed by making use of upper and lower bounds for the exact solution of the problem. The boundary conditions for heat and mass flow inside the oil/water region are established, and a simple method of detennining the saturation at the downstream side of the condensation front is presented. Experimental results show that the theory gives an accurate description of steam-zone expansion before the critical time. At times near critical time, and also after critical time, experimental data show theoretical steam-zone volumes to be high. Even so, the accuracy of the calculated steam-zone volumes is sufficient for practical purposes. THEORY OF STEAM-ZONE GROWTH INTRODUCTION The engineering evaluation of steam-drive processes is based mainly on a mathematical description of reservoir heating by hot fluid injection presented several years ago by Marx and Langenheim.l This theory can be used to determine the growth of the steam zone,2 provided the flow of heat from the steam zone into the liquid zone ahead of the condensation front (Fig. 1) is neglected. Since we must expect, however, that heat transport into and inside the liquid region (the region between condensation front and production wells) will affect both the water/oil flow and the growth of the steam zone, it is necessary to re-examine and to generalize the Marx and Langenheim theory as applied to steam drive, which will hereafter be termed the "old" theory. Observations made in several field tests emphasize the need to study the heat flow across the condensation front (CF). On several occasions it was discovered that heat was transported in the liquid zone far ahead of the advancing CF, resulting in an early breakthrough of warm water. This seems to occur most often when high injection pressures are applied, when "wet" kONINKLIJKE/SHELL EXPLORATIE EN PRODUKTIE LABORATORIUM RIJSWIJK, THE NETHERLANDS SHELL DEVELOPMENT CO. HOUSTON, TEX.

Jan 1, 1970

By L. A. Rapoport, W. J. Leas, C. W. Carpenter

A program of scaled flow model experiments has been undertaken to study the performance of five-spot water floods. The modeling procedures are discussed and the construction and operation of the flow models are described. The tests were specifically designed to investigate the effects of injection rate and of oil-to-water viscosity ratio on oil recovery. The areal sweep of five-spot water floods was also studied, and direct comparisons between five-spot and linear flooding behavior have been established to permir better utilizatfon of core analysis results for purposes of field waterflood evaluations. INTRODUCTION For a long time reservoir engineering procedures have been confined to the consideration of idealized, one-dimensional or linear systems.11* However, direct application of the linear waterflood mechanism to the evaluation of actual field performance is justifiable only for the purposes of rough approximation. To obtain a realistic analysis of field behavior, it is necessary to account for reservoir geometry and well distribution, and to take into consideration the configuration of the flow of fluids between injection and producing wells. In the past it has been attempted to account for the complex nature of the fluid movement in the reservoir by associating each well pattern with a certain "sweep-efficiency constant", e.g., a sweep efficiency of 72 per cent for the conventional five-spot pattern.&apos; More recently, methods of waterflood evaluation have been improved by the introduction of variable sweep-efficiency factors, defining the areal coverage as a function of cumulative water injection and viscosity or "mobility" &apos; ratio.""&apos; These variable sweep-efficiency factors, derived from miscible fluid displacement experiments, require the assumption of piston-like displacement of oil by water. This implies that only one fluid, either oil or water, is moving at any one point of the reservoir ht any moment. Such an assumption, however, is incompatible with the bulk of the experimental evidence and also theoretical considerations which indicate that usually in the area invaded by water, both oil and water are flowing in varying proportions, in accordance with the relative permeability characteristics of the porous medium. Another recent approach toward the evaluation of five-spot flooding behavior involves correlations of the areal sweep efficiency (observed in immiscible fluid experiments) with a "mobility ratio," based on the consideration of the "average" relative permeability to the displacing fluid. Combination of these correlations with a modification of the Buckley-Leverett frontal-drive equation leads to a computational procedure for predicting the recovery performance of areal water floods. While of considerable interest, this approach has no complete theoretical justification and, therefore, does not entirely clarify the features of five-spot water floods. Rigorous, theoretical treatment of waterflood behavior in two- or three-dimensional systems as encountered in the reservoir, including proper consideration of relative permeability and capillary-pressure functions, is as yet beyond the scope of analytical or computational methods. The, behavior of such water floods can, however, be studied by means of scaled flow model ekperiments, which is the approach followed in these investigations. The flow model experiments described in this paper simulate a uniform, horizontal formation of constant thickness with fully penetrating wells. The system under consideration is entirely liquid saturated and not influenced by compressibility effects, which is representative of the situation achieved after gas fill-up under practical conditions. Gravity segregation effects are assumed to be negligible, which implies a relatively thin reservoir formation and/or low vertical permeability. THEORETICAL CONSIDERATIONS For a flow system of uniform thickness with no appreciable effects of gravity segregation, the analysis of water-oil displacement can be confined to a two-dimensional treatment based on the consideration of a "slice" of unit thickness. Referring to previously established derivations," the displacement of oil by water in such a system can be described by the following, dimen-sionless equations:

By D. R. Parrish, F. F. Craig

In some respects, forward combustion (alias: fie-flooding, underground combustion, in-situ combustion, dry combustion) is an efficient oil recovery method. For example, its displacement efficiency (as much as 95 percent in reservoirs with high porosity and high oil saturation) is good. In its use of heat, however, forward combustion is grossly inefficient. Most of the heat generated during forward combustion is wasted. About one-fifth of the heat is picked up by the incoming air and re-used (regenerated); the remainder is left behind the combustion zone and ultimately is lost. Compressed air is expensive, and heat generated from compressed air via underground combustion is too valuable to be thrown away. Improvement is needed. Since the early 1950&apos;s, improvements in forward combustion by addition of water injection have been discussed, particularly in the patent literature, with increasing frequency.&apos;-&apos;&apos; These concepts are based on the fact that relatively inexpensive water can be used to pick up much of the heat that is ordinarily wasted during forward combustion, and transport it forward to a more useful place, thereby reducing the amount of air required to move a heat front through the reservoir. Several approaches have been proposed. Merriam and Squires1 taught that a combustible (fuel) gas should be injected along with the air and water. This, however, would give at least a tendency toward reverse combustion instead of forward combustion, and would reduce the velocity of the heat bank, defeating the primary purpose of the water injection. Pelzer2 described primarily a batchwise method of forward combustion and water flooding. He advocated burning part way through a reservoir with forward combustion, then switching from air to water injection alone to push the heat toward the producing wells. Though possible, this method has serious disadvantages, the major one being that, by the time water injection is started, much of the heat has been irretrievably lost to the overlying and underlying formations. We have described what is considered a practical Combination Of Forward Combustion And Water-flooding (COFCAW).3,4 Air and water are injected simultaneously (or alternately) after a small heat bank has been formed by forward combustion. Some of the injected fluid fills up the region behind the heat bank; the rest enters the heat bank. Upon encountering heated rock, the water is converted into steam, which then flows ahead of the combustion zone and displaces much of the oil in place. Only enough air is injected to provide heat for vaporization of the water and to compensate for reservoir heat losses. Since the size of the heat bank is minimized, heat losses also are minimized. Because of more efficient oil displacement ahead of the combustion zone, fuel consumption is reduced. Too, unlike a forward combustion zone that cannot move until all fuel at a point has been consumed, COFCAW leaves some

Jan 1, 1970

By J. R. Kyte, L. A. Rapoport

The material balance equation for partiai or complete 1:trter-drivc reser1,oir.s been re-arranged to include a pressure irrferferencr ternr. This pressure interfernce term was ohtained from the theory established hy Morinda. A .successful progranl deissed to solve the resrclritlg c,qrration in a few hours on a nleelitr .sizeti, digital compuier. No prior knowdelge of the properties of the eqrrifer is required in the solution. 7&apos; ri~rl calculations areatle using various assumed con!binuiions of value of the constant expretsing the qquifer properties correspoinds to the best fit ohtained between observed and calculated field history. Details of the method of calculaion are illus trated by an example. INTKCDUCTICN A method 01&apos; calculating the material balance equation for a partial water-drive reservoir completed in an aquifer has been described in the literature." Pressure interference i&apos;rc,m other pools iocated in the same nquiucr is not taken into account in this solution. When present, the pres- sure interference effect must be included in order to make a complete performance analysis. A correction for this effect has been developed using the theory cstahlished by Mor-tada.&apos; Manual calculation of the material balance for water-drive reservoirs is a formidable task. Such manual ca1culation when interference eflects are includcd would be impossible to completc in any reasonable time period. However, a medium-sized stored program, digltal conlputer can he used satisfactorily. The method is described in the following sections. THECKETICAL CONSIDEKATIOhS The material balance equation for a reservoir producing by partial or total water drive, as proposed by Van Everdingen, This relationship is linear for Z/a from which thc slopc, Z, and the intercept, N, can be obtained by the least-squares method from the data points. It is not necessary lo know the aquifer properties to uktain oplimum values tor 2nd IV. Various assumptions are made tor At,,: where 11 is production time, seconds (oilfield between the equal successive time intervals into which the history is divided. t,, h It, times the number of equal time periods from i.e.. the total dimen-sionless time interval from It is noted that the dinlensionless time interval, At,,, involves all the relevant aquifer properties. Accordingly, any assunlption of is an assumption of the combined aquifer properties which pertain to the solution of the equation. Deviations of the data points from a straight line arc calculated for each assumption of and the value of A, correspond ing to the minimunl deviation is chosen as the optimum solution. This method is described in Ref. 1 and 2. Thc material balance solution is usually made with an initial assumption that the reservoir being studied is located in an infinite aquifer from which there are no other withdrawals. Deviation of the points by amounts greater than can be attributed to inaccuracies in field data may be the result of: (1) boundary effects of a limited aquifer, (2) pressure interference from othcr pools completed in the .sarne aquifer, and (3) a combination of these effects. Differentiation of these effects can be made only on the basis of geologic data. Pressure interference from other fields will result in a pressure drop at the original oil-water boundary having an interference component, The measured instantaneous pressure drops at the oil-water

By Denton R. Wieland, Harvey T. Kennedy

The frequency of bubble formation is measured for oils from the East Texas field and the Slaughter field in cores from these fields. East Texas oil was also tested in a Slaughter core. The data checked previous results in that a definite supersatrtration may exist without the formation of bubb1es in periods of observation totaling 158 hours. Measured frequencies varied from zero at supersaturations below 14 psi to 3.1 bubbles/sec/ cu ft of rock at 40 psi supersaturation INTRODUCTION The formation of bubbles is essential in a solution gas-drive reservoir in the displacement of oil. In formations where the basic pore structure consists of fine pores or fractures connected to larger pores or crevices, bubble formation in the fine pores may be the most effective means by which the oil may be displaced, even if secondary methods are applied. The most efficient oil displacement would occur where at least one bubble would form in every pore. As the pressure decreased, the bubbles would expand, due to the gas diffusing into the bubbles, and displace the oil toward the lower pressure area around the wellbore. Kennedy and Olson&apos; have shown that the total number of bubbles formed in a reservoir is primarily dependent upon the rate of pressure drop, the rate of diffusion of gas through oil, and surface area of the rock. The surface tension of the liquid and other characteristics of the rock may also exert an important effect, which would give a different frequency to different rocks and oils and various combinations of different rocks and oils. The purpose of this paper is to present frequency measurements on two rocks and two oils. A hydrocarbon liquid is at equilibrium with gas when it contains the maximum volume of gas in solution that it can have at the prevailing temperature and pressure. If the temperature is maintained constant and the pressure is reduced below the equilibrium pressure, the liquid is supersaturated to the extent of the difference in the two pressures until the gas comes out of solution as a free phase. As shown by previous work, and confirmed in the present investigation, the number of bubbles formed per second increases with an increase in supersaturation. Kennedy and Olson showed that a tenfold increase in the rate of pressure decline resulted in a tenfold increase in the number of bubbles formed per cubic foot of reservoir rock, e.g., for decline rates of 0.1, 1, and 10 psi/day, the number of bubbles formed per cubic foot would be 40, 400, and 4,000, respectively. For another hypothetical system, Stewart, Hunt, and Berry&apos; estimated that for pressure declines of 0.1, 1, and 10 psi/day, the number of bubbles formed would be approximately 10, 100, and 1,000 cu ft of reservoir rock. Their value of bubbles formed is lower than given by Kennedy and Olson; however, the order of magnitude agreed satisfactorily. The rate of diffusion will also affect the number of bubbles formed. If the diffusion coefficient is high, the gas in solution in the oil surrounding a single bubble will tend to diffuse into the bubble rather than forming other bubbles. On the other hand, the formation of more bubbles will tend to occur if the diffusion coefficient is very low. In his thesis dated Aug., 1953, Wood" reports on the frequency of bubble formation in a system consisting of oil and rock from the Rangely field. In the present work, the apparatus of Wood was redesigned to allow greater precision, especially at low supersaturations, where bubble growth is slower, and detection is more difficult.

Jan 1, 1958

By Irving Fatt, Robert J. S. Brown

INTRODUCTION The wettability of reservoir rocks is recognized as one of the major factors that determines their multiphase flow properties. Multiphase flow properties in turn govern reservoir performance. For this reason, the reservoir engineer is vitally interested in the wettability of reservoir rock. The production research scientist is interested in wettability because the solution of the problem of improving waterflood recovery by use of surface-active additives may hinge on a detailed understanding of the wettability of reservoir rock. Early workers believed that all reservoir rocks were preferentially water wet. Continued study soon uncovered certain rocks which had flow properties that indicated that they were preferentially oil wet. Recent studies show that there is a spectrum of wettability. Some reservoir rocks have multiphase flow properties which indicate that they have an intermediate wettability. The differences in recovery observed among different reservoirs may, in part be the result of small differences in wettability. Differences in wettability amcng reservoir rocks may be the result of differences in the fraction of the total surface area that is preferentially oil wet and preferentially water wet. Nuclear magnetic relaxation measurements may provide a method for determining the fraction of preferentially oil wet and preferentially water wet area in porous media. Measurements to be reported here were made on water-saturated unconsolidated sand packs which were mixtures of preferentially oil wet and preferentially water wet sands. The resulting data show definitely that the fraction of each wettability in these sand packs can be determined by this nondestructive method of measurement. THEORY OF NUCLEAR MAGNETIC RELAXATION The details of the theory of nuclear magnetic relaxation have been presented previously1. In this paper, only a very brief review of the theory will be given to show its application to the measurement of fractional wettability. Nuclear magnetic relaxation is a process of adjustment undergone by certain atomic nuclei when their magnetic environment is changed. They adjust, or relax, to be in equilibrium with a new magnetic field. The reason that some atomic nuclei are affected by a magnetic field is that they themselves behave in some respects like tiny magnets. In particular the proton, the nucleus of the hydrogen atom, behaves as a magnet. The proton spins, somewhat as the earth spins on its axis. The spinning charge produces a magnetic field just as current traveling in a loop of wire produces a magnetic field. No net charge is required, however; for instance, the neutron has spin and magnetic moment. When the spinning proton is subjected to an externally applied magnetic field, its axis tends to line up with the magnetic field, somewhat as a compass needle tends to line up, but the spinning proton cannot do this as easily as the needle can. The proton also has mass, and its spinning mass makes the proton behave somewhat like a gyroscope. As is well known, a force which tends to change the axial direction of a gyroscope merely causes the gyroscope to precess: lhat is, the axis begins to trace out a conical surface. The gyroscope cannot actually assume the position suggested

Jan 1, 1957

By M. King Hubbert, David G. Willis

A theoretical examination of the fracturing of rocks by means of pressure applied in boreholes leads to the conclusion that, regardless of whether the fracturing fluid be of the penetrating or non-penetrating type, the fractures produced should be approximately perpendicular to the axis of least stress. The general state of stress underground is that in which the three principal stresses are unequal. For tectonically relaxed areas characterized by normal faulting, the least stress should be horizontal; the fractures produced should be vertical with the injection pressure less than that of the overburden. In areas of active tectonic compression, the least stress should be vertical and equal to the pressure of the overburden; the fractures should be horizontal with injection pressures equal to or greater than the pressure of the overburden. Horizontal fractures cannot be produced by hydraulic pressures less than the total pressure of the overburden. These conclusions are compatible with field experience in fracturing and with the results of laboratory experimentation. INTRODUCTION The hydraulic-fracturing technique of well stimulation is one of the major developments in petroleum engineering of the last decade. The technique was introduced to the petroleum industry in a paper by J. B. Clark,&apos; of the Stanolind Oil and Gas Co. in 1948, and since then its use has progressively expanded so that by the end of 1955 more than 100,000 individual treatments had been performed. The technique itself is mechanically related to three other phenomena concerning which ah extensive literature had previously developed. These are: (1) pressure parting in water injection wells in secondary-recovery operations, (2) lost circulation during drilling, and (3) the breakdown of formations during squeeze-cementing operations, all of which appear to involve the formation of open fractures by pressure applied in a wellbore. The most popular interpretation of this mechanism has been that the pressure had parted the formation along a bedding plane and lifted the overburden, notwithstanding the fact that in the great majority of cases where pressures were known they were significantly less than those due to the total weight of the overburden as determined from its density. Prior to 1948, this prevalent opinion had already been queried by Dickey and Andresen,2 in a study of pressure parting; and by Walker,"&apos; who, in studies of squeeze cementing, pointed out that the pressures required were mostly less than those of the overburden, and inferred that the fractures should be vertical. J. B. Clark,&apos; in his paper introducing hydraulic fracturing, and later Howard and Fast," and Scott, Bearden, and Howard," also of Stanolind, postulated that the entire weight of the overburden need not be lifted in producing horizontal fractures, but that it was only necessary to- lift an "effective overburden," requiring a correspondingly lower pressure. Hubbert,7 in discussing the paper by Scott and associates, pointed out that the normal state of stress underground is one of unequal principal stresses; and in tectonically relaxed areas, characterized by normal faults, the least stress should be horizontal. Therefore, in most cases, fracturing should be possible with pressure less than that of the overburden and, moreover, such fractures should be vertical. Harrison, Kieschnick, and McGuire,8 also on the expectation that the least principal stress should be horizontal, argued strongly in favor of vertical fracturing. Scott, Bearden, and Howard6 observed that, when

Jan 1, 1958

By J. R. Kyte, L. A. Rapoport, R. J. Stanclift, S. C. Stephan

Experimental studies covering a wide range of core materials and fluid properties have been conducted to determine the mechanism of oil displacement by water in a partially gas-saturated porous medium. In all instances, the presence of a gas phase was found to have a beneficial effect in reducing residual oil saturations. The practical significance of this benefit is discussed, and a simplified procedure is outlined for evaluating the effects of free gas on water flooding by means of short core tests. INTRODUCTION In addition to oil and water, reservoirs subjected to water flooding frequently contain, also, a gas phase. Common engineering procedures account for the presence of this "free gas" only from the viewpoint of volumetric balance, implying that the only role of the gas consists of providing "fill up" space. It is usually visualized that during the initial stages of the water invasion, the oil, moving ahead of the water, displaces part of the gas and that subsequently the remaining portion of the gas phase is totally compressed and dissolved in the advancing oil bank. Thus, consideration of a two-phase, water-oil flooding mechanism supplies an adequate basis for predictions of oil recovery if the pressure build-up caused by the flood is sufficiently great, so as to reduce the free gas saturation to a negligible value in any portion of the reservoir by the time that portion is reached by the advancing flood water. In many in- stances, however, waterflooding operations are carried out when the reservoir pressure is still relatively high SO that the pressure build-up associated with the water flood may not result in complete dissipation of the gas phase. Under such circumstances it is necessary to ascertain whether or not the presence of free gas has an effect on waterflood behavior and, if so, to account for any such effect in the evaluation of field operations. The difficulties encountered in attempting to arrive at a comprehensive picture of the displacement mechanism of oil by water in the presence of free gas stem from the fact that most of the published data on this subject were obtained under conditions that were either not clearly representative of those existing in a reservoir or not sufficiently controlled to permit evaluation of the effects specifically caused by the gas phase. Some of the most precise data found in the literature were obtained in tests using a "static" gas phase created prior to each water flood by circulating oil through partly gas-saturated cores until a "trapped," residual gas saturation was established throughout the system.&apos; Other investigations were conducted under more representative conditions where the gas phase was free to move during the water flood: but in most cases the results of these tests were obscured by auxiliary phenomena which arise from the compressibility and solubility of the gas and preclude quantitative conclusions.2,3,4 Prior to water flooding, a reservoir usually contains oil and gas, distributed in a manner in which both are free to move. Thus, the principal objective of the present studies was to determine the type of displacement mechanism occurring when water invades a porous medium containing oil and a "mobile" gas phase. Further-

Jan 1, 1957

By H. A. Jr. Koch, C. A. Hutchinson

Miscible phase displacement of oil from reservoirs has been emphasized in the past few years. The reason for this emphasis lies in the high oil recovery attainable by this process. Removal of capillary effects in the reservoir leads to recoveries approaching 100 per cent in the area contacted by the miscible phase. The miscible slug process is one means of obtaining a miscible displacement. Here a band or slug of I,PG is injected into the reservoir prior to gas injection. The idea is to maintain the band of LPG "wedged" between the gas and oil phases and thus achieve a miscible phase displacement. A second method lor achieving miscibility is through the injection of a gas which is not miscible with the reservoir fluid but which develops a zone of miscibility in [he reservoir through mass trans-ier with the reservoir oil.&apos; This mass transfer results in either an enrichment of the lean injected gas by intermediates from the oil or an enrichment of the oil by intermediates from a rich injection gas or one that has been enriched on the surface by LPG addition. We are interested here in discussing the process in which miscibility is developed at the displacement front by the evaporation of interrnediatcs from the oil phase into the gas phase. This process "builds up" its own slug of miscible material at the displacement front and therefore does not require the injection of LPG to obtain miscibil-ily. Each process has its own area of applicability. Generally, the high pressure gas process is applicable only with reservoir fluids which con-!ain a high concentration of inter- mediates. If the high pressure gas process is technically feasible at pressures less than 4,500 psi, it is probably more desirable economically than the slug process. The slug process has broad applicability in the shallower reservoirs and with reservoir fluids which contain a relatively low concentration of LPG and natural gasoline constituents. This paper deals with some new concepts of the high pressure gas injection process where it is proposed that flue gas can be substituted for hydrocarbon gas without sacrificing our goal of miscibility. MECHANISM Introduction Considerable effort has been devoted to study of the mechanics ot the high pressure gas injection proc-one generalization result-ing from some of these studies was that the composition of the injected gas is relatively unimportant in establishing the miscibility pressure* for a given reservoir fluid. This generalization is correct for the composition range of gases typically encountered in the field. Two such gases are a gasoline plant tail gas containing 85 per cent methane and 15 per cent ethane, and a field separator gas containing 70 per cent methane and 30 per cent heavier components. The most important factor which sets the miscibility pressure in the operation is the reservoir fluid composition, particularly the concentration of LPG-natural gasoline constituents. The injected gas is the agent by which the LPG-natural gasoline constituents are concentrated at the displacing front to create a miscible displacement. Based on these results, it appeared feasible that some inexpensive gas, such as flue gas, might be substituted for hydrocarbon gas for use in the high pressure gas process. A re-examination of the phase relations of the high pressure gas injection process should clarify the principle behind using flue gas (essentially nitrogen) as an injection gas. Three Component Diagram The phase relations of the high pressure gas injection process have been illustrated by the use of the three component diagrams.&apos;," In Fig. 1 we have arbitrarily represented the multi-component reservoir system by three components; methane, ethane through hexane, and heptanes plus. The solid curve ABC is the phase boundary curve. It represents the locus of compositions which have fixed saturation pressure at a fixed temperature; the lower branch AB shows bubble point compositions, the upper branch BC, the dew point compositions. Point B is the coniposition of the critical mixture at this temperature and pressure. The dashed lines (tie lines) connect vapor and liquid compositions which are in equilibrium. Let us consider Reservoir Fluid D which we wish to displace in 21 miscible manner by gas injection. Let us further restrict the discussion to the case where miscibility between an injection gas and the reservoir fluid at the displacement front is developed by gas enrichment in the reservoir. For this case, any gas whose composition lies between Points C and E on the right side of the three component diagram can be used to give a miscible displacement of Reservoir Fluid D. This is true because the more mobile injected gas moves faster than the displaced oil and is in continuous contact with virgin oil at the displacement front. This leads to a continuing enrichment 01&apos; the gas at the displacement front by evaporation of the C, - C,

By George G. Binder Jr., James W. Lacey, Arthur L. Draper

An experimental investigation of miscible fluid displacement has been made in linear porous media under highly adverse mobility ratio conditions. Various refined oils were displaced at field rates by liquid propane in both horizontal and vertical models 9 to 36 ft in length. In horizontal tests, mixed zones or transition zones form between the two liquids early in the displace-lents. Initially, the zones grow rapidly. III 2-in. diameter models, the mixed regions become nearly stable zones from 6- to 20-ft long. From then on these quasi-stable zones apparently grow only by molecular diffusion. The lenghts of the quasi-stable zones appear to be almost independent of the core material in which displacement;; are made and increase with increasing rates of injection and oil/solvent viscosity ratios. The lengths of the quasi-stable zones. measured in 2-in. diameter models at viscosity ratios and injection rates anticipated in the field indicate that essentially piston-like displacements of oil by solvent would occur. However, the mixed zone lengths increase markedly as the model diameter increases. In vertical displacements controlled by gravity drainage, only short trans-sition zones form. The entire lenghts of these zones can be accounted for by molecular diffusion. Therefore, in reservoirs having characteristics allowing a high degree of gravity drainage, nearly complete oil recovery should be attainable by miscible fluid displacement. INTRODUCTION The oil industry is faced with an ever increasing cost for finding new oil reservoirs. This, coupled with the fact that conventional methods of oil recovery usually leave from 1 to 3 bbl of oil in the ground for every barrel produced, has aroused considerable interest in developing more efficient methods. An integral part of the mechanism of many of these new methods of recovery is a miscible fluid displacement in which oil is displaced by a much less viscous materiai such as liquid propane. A numher of irlvestigators have explored the mechanics of miscible fluid displace-ment,&apos;.2.1.&apos;."." but in the experiments reported many of the important factors have been varied only over a limited range far removed from norma1 oilfield values. This paper presents the results of an experimental study of oil displacement by liquid propane at field rates from linear, natural-sandstone models of various sizes. Both horizontal and vertical flow tests were made. EXPERIMENTAL APPARATUS AND PROCEDURE Refined oils with viscosities varying over a wide range were displaced at several rates by liquid Propane from both horizontal and vertical flow systems. The behavior of the displacement was observed by propane-oil analyses of the effluent and, when possible, by differential pressure measurements along the cores. A schematic flow diagram of the apParatus used to conduct the experirents is shown in Fig. 1. The cores used in the tests were 9-ft lengths cut from natural consolidated sandstone and turned to the desired diameter. These were cast in steel tubing with a low melting point (160°F) alloy which expands on solidification. Flanges welded to the ends of the steel pipes permitted sealing the ends or butting one core to another. Pressure seals were maintained with O-rings. The cores were cleaned by flushing

By W. B. Gogarty, W. C. Tosch

A new recovery process for producing oil under both secondary and tertiary conditions utilizes the unique properties of micellar solutions (also known as microemulsions, swollen micelles, and soluble oils). These solutions, which displace 100 percent of the oil in the reservoir contacted, can be driven through the reservoir with water and are stable in the presence of reservoir water and rock. Basic components of micellar solutions are surfactant, hydrocarbon and water. They may also contain small amounts of electrolytes and co surfactants such as a1cohol.r. The specific reservoir application dictates the type and concentration of each component. A salient feature of [he process is the capability for mobility control. Micellar solution slug mobility, by way of viscosity control, is made equal to or less than the combined oil and water mobility. Mobility control continues with a mobility buffer that prevents drive water from contacting the micellar solution. Laboratory and field flooding have proven that the process is technically feasible and that surfactant losses by adsorption on porous media are small. Introduction projects are under way to recover the maximum amount of oil under the most favorable economic conditions.&apos; : New techniques are being developed to increase oil recovery,3" Polymer solutions are becoming an important means of controlling mobility in a waterflood. Thermal methods such as in-situ combustion and steam injection are being used in reservoirs containing highly viscous crudes. Surfactant flooding is receiving attention as a method of reducing interfacial tension to increase recovery.*&apos;" Exotic recovery processes have been considered primarily for &apos; perations. Economics are unfavorable in most cases for tertiary recovery. studies at the Denver Research Center of the Marathon oil CO. have led to a new oil recovery method.* Micellar solutions (sometimes called microemulsions, swollen micelles, and soluble oils) are used to recover oil by miscible-type waterflooding. Basically, these solutions contain surfactant, hydrocarbon, and water. The method can be used in either secondary or tertiary operations. First, thc concept of thc process is considered in terms of the requirements for an effective miscible waterflood ing operation. Next, micellar solution properties are described including structure, composition, and phase behavior with reservoir fluids. Fluid characteristics are then considered as related to mobility control, and, finally, laboratory and field results are presented to illustrate the efficiency of the process. Concept of the Process Unit displacement efficiency and conformance determine the effectiveness of any oil recovery mechanism. In theory, a miscible waterflood should be capable of a 100-percent unit displacement efficiency with a correspondingly high conformance. Requirements for the slug of a miscible waterflood include (1) 100-percent displacement of oil in the reservoir contacted, (2) controllable mobility, (3) the capability of being driven through the reservoir with water, (4) a low unit cost to enhance economics, and (5) the ability to remain stable in the presence of reservoir water and rock. Micellar solutions satisfy requirements for the slug of a miscible waterflood process. Our discovery that these solutions acted as though they were miscible by displacing all fluids in the reservoir and by being displaced by water solved the miscibility problem. Adequate mobility control is possible by variations in solution viscosity through compositional changes. Economic requirements are met since micellar solution costs below $6/bbl appear possible, Mi cellar solutions stabilize surfactant in the presence of reservoir rock and water, thus reducing the importance of the problem of surfactant loss by adsorption. Fig. 1 illustrates schematically how these solutions are used. Operations start with injection of a micellar solution slug that serves as the oil displacing agent. Next, a mobility buffer of either a water-external emulsion or water solution containing polymer (thickened water) is injected to protect the slug from water invasion. Finally, drive water (water used in a regular waterflood) is injected to propel the slug and mobility buffer through the reservoir. Reservoir oil and water are displaced ahead of the slug, and a stabilized oil and water bank develops as shown in Fig. 1. Stabilized bank saturations are independent of original oil and water saturations. This means that, for a particular type of reservoir, the displacement mechanism is the same under secondary and tertiary recovery conditions. Oil is produced first in a secondary operation. For tertiary conditions, water is produced first. Movement of the slug through the reservoir is stabilized by the mobility buffer. An unfavorable mobility ratio usually exists at the interface between the buffer and drive

Jan 1, 1969

By N. D. Shutler

This paper describes a numerical mathematical model of the steamflood process that depends on fewer restrictive assumptions than models previously reported. The solution, however, is obtained economically. Example calculations are presented that, on comparison with experimental results, tend to validate the model. Results that expose certain process mechanics are discussed. The model describes the simultaneous flow of three phases — oil, water and gas — in one dimension. It includes the effects of three-phase relative permeabilities, capillary pressure, and temperature- and pressure - dependent fluid properties. Interphase mass transfer of water-steam is allowed, but the oil is assumed nonvolatile and the hydrocarbon gas insoluble in the liquid phases. The model allows heat convection in one dimension and two-dimensional heat conduction in a vertical cross-section spanning the oil sand and adjacent strata. The hydrocarbon-steam gas composition is tracked, but the effect of gas composition on water-steam phase behavior is neglected. The model is solved numerically in three separate stages. The three-phase mass balances are solved simultaneously using Newtonian iteration on nonlinearities occurring in the accumulation terms. The energy balance is solved separately by noniterative application of the alternating-direction implicit procedure. Separate solution of the composition balance is accomplished by straightforward solution of the finite difference equations. The method of effecting nonsimultaneous, stable solution of the mass and energy balances is the key to the success of the model. INTRODUCTION Mathematical tools as well as laboratory and field experiments are necessary to help us understand the complex steamflood process. A mathematical model can expose process mechanics and show the relative importance of process variables, but this ability is often limited by restrictive assumptions. Most known models of steam processes,l-5 with the exception of the model of Gottfried,6 are "simplified" in that they involve analytic approximations and require many restrictive assumptions. The primary utility of these methods lies in the routine use as an aid in engineering design. By contrast, the comprehensive model presented here is numerical and requires far fewer restrictive assumptions. It finds its primary utility as a research tool. It serves as an aid in understanding the nature of the process, in interpreting laboratory experiments and in evaluating and developing simpler mathematical models for engineering design. The major reason why previously presented models have been confined to the "simplified" class is evidenced by the one published exception. In 1965, Gottfried6 presented a numerical model for the combustion process of which the steam process is a subset. The result is a comprehensive tool (though it neglects capillarity and two-dimensional heat conduction) that is troubled with convergence problems and that requires 2 to 3 hours of IBM 7094 time to complete a calculation. Though our present model does not simulate combustion, it does consider capillarity and two-dimensional heat conduction and it overcomes the convergence and computer-time problems. MATHEMATICAL DESCRIPTION OF STEAMFLOODING FLUID FLOW The equations employed to describe three-phase fluid flow are of a familiar form. Darcy&apos;s. law provides expressions for the velocities of the three phases (oil, water and gas), which, when combined with oil, water and gas mass balances, give the partial differential equations governing flow of the

Jan 1, 1970