Search Documents

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

By G. Thodos, W. F. Stevens

The point-source function introduced by Horner&apos; us a solution to the general unsteady-state equation for the flow of fluids through porous media has been utilized to calculate pressure profiles for adjacent wells as a function of time. By trial and error, the time at which the pressure waves meet to cause interference can then be calculated. For this calculation, the time of interference is arbitrarily defined as that time when both pressure waves exhibit a pressure decrease of 25 psi ar the same point within the formation. To illustrate the method, an example involving two adjacent wells producing at different rates has been worked out and is presented in detail. INTRODUCTION The initiation of production of a new well is usually accompanied by pressure drawdown measurements in order to establish the extent and characteristics of the formation. The interpretation of such pressure drawdown data is based on mathematical developments involving a single well in an infinite reservoir. If more than one well is present in a reservoir, the possibility of interaction between wells exists. Therefore, care must be taken in interpreting the pressure drawdown data. It becomes important to be able to determine when this interaction, or interference, becomes significant. This paper presents a method for the prediction of the approximate time of interference between two such adjacent wells. The basic differential equation describing the pressure distribution existing during unsteady-state flow of a slightly compressible fluid has been solved by Hor-ner.&apos; for a constant production rate from a circular homogeneous reservoir of infinite extent. The solution utilizes the point-source function to define the pressure at various points within the formation as a function of time. The resulting expression is: q = production rate, STB/D ,;= formation volume factor, dimension-less µ = fluid viscosity, cp k = formation permeability. md H = formation thickness, ft C - fluid compressibility, psi&apos; F = formation porosity, dimensionless R = distance from centerline of well, ft T = production time, hours Ei(—u) = Ei-function of U Eq. 1 may be used to calculate the pressure profiles generated during production from a single well in an infinite reservoir.

By H. J. Ramey, I. S. Bousaid

Experimental results on the oxidation reaction kinetics in the forward combustion oil recovery process are presented. A total of 48 runs were made wherein a stationary thin layer of coked, unconsoli-dated sand was burned isothermally in a combustion cell. Individual runs were made at various temperature levels to permit determination of the effect of temperature upon the reaction. An expression was obtained for the burning rate of carbon as a function of carbon concentration, combustion temperature and oxygen partial pressure. The carbon burning rate for two types of crude oil indicated a first order reaction with respect to both carbon concentration and oxygen partial pressure. The effect of combustion temperature on the reaction rate constant matched the Arrhenius equation. The activation energy was similar for the two crude oils examined. The activation energy decreased for a porous media containing clay. The rate of oxidation of crude oil at reservoir temperature was found to be significant. Other significant findings included information on hydrogen-carbon content of fuel residues, fuel reactivity and the products of combustion. INTRODUCTION The production of crude oil by underground combustion has been studied in the laboratory by many investigators.1-5 Results of laboratory and field experiments have been reported in the literature describing the forward combustion process. But as yet, no qualitative or quantitative study of the kinetics of fuel combustion involved in this process has been reported. The fuel concentration and the rate at which fuel is burned at the front are important factors governing the air requirement in a forward combustion operation. Although the fuel is essentially unrecoverable crude, the air required to burn the fuel is an important economic factor in this process. Because fuel is burned, the heat transport associated with forward combustion is a key and unique feature of this new oil recovery method. Many investiga-tors5-9 have presented information on the heat transmission and fluid mechanics involved in forward combustion. Berry and Parrish10 demonstrated the utility of considering reaction kinetics in reverse burning. From differential thermal analysis, Tadema11 presented a qualitative discussion of the nature of reactions between oil and oxygen in combustion oil recovery. Although little quantitative work has been done on the reaction kinetics involved in forward combustion oil recovery, an extensive literature does exist on combustion of carbons and oils, and carbonaceous residues from cracking catalyst pellets. Dart, et al., 12 studied the combustion rate for oxidation of carbonaceous residues on clay catalyst pellets, and found the reaction to be second-order with respect to carbon concentration, and first-order with respect to oxygen partial pressure for carbon concentrations less than 2 weight percent of the catalyst weight. The reaction appeared to be first-order with respect to carbon concentration for concentrations greater than 2 percent. Metcalfe13 noted that other workers had found that aging of the fuel during the combustion process was responsible for changing coke properties, and it accounted for the apparent second-order carbon concentration effect found by Dart, et al. It appears that burning of residues from cracking pellets is first-order with respect to both carbon concentration and oxygen partial pressure. Dart, et al., also observed that hydrogen in the hydrocarbon residue appeared to react faster than the carbon. Lewis, et al.,14 studied oxidation of charcoal, coke and graphite in a fluidized bed. Gas velocities were high enough to partially lift and circulate the carbon particles. Their results indicated first-order reaction dependency with respect to both carbon concentration and oxygen partial pressure. One

Jan 1, 1969

By J. B. Dorsey, T. W. Brinkley

Performance of a high-pressure, vaporizing gas-drive project in a volatile oil reservoir is discussed. The project has been in operation for over 8 years. Present recovery from the over-all displaced portion of the reservoir is 40 percent, and an ultimate recovery of 60 percent is indicated. Calculations reflect an over-all displacement volume of 50 percent for the entire reservoir; therefore, ultimate oil recovery will be 30 percent of the entire reservoir oil originally in place. Laboratory vaporization tests conducted prior to project initiation also are presented and discussed, along with an analytical method of determining stock-tank oil production and gar injection requirements related to over-all volumetric displacement efficiency from these data. Introduction The Shoats Creek field is located in Beauregard Parish, La., in the extreme western part of the state (Fig. 1). Discovered in 1956, the field is productive in three Cockfield sand zones, ranging in depth from 8,000 to 9.000 ft. The upper two zones, 1st and 2nd Cockfield, are condensate-bearing gas reservoirs. The lowest zone, the 5th Cockfield, is a volatile-oil reservoir. Laboratory work conducted on crude oil from the 5th Cockfield zone shortly after field discovery demonstrated that, when contacted with dry gas, signidcant volumes of the oil could be vaporized at reservoir temperature and at an appropriately high reservoir pressure. Further study led to the conclusion that early initiation of a crude vaporization project in the 5th Cockfield zone would be superior economically to other methods of reservoir depletion. Consequently, the zone was unitized and gas injection operations started late in 1958. This paper discusses laboratory work performed in conjunction with the project, and project performance to the present. History of Development The 5th Cockfield zone in Shoats Creek is composed of a thin, very shaly, low-permeability sandstone. Oil accumulation is controlled by a structural trap bounded by faults along the north and west and by a water contact on the southern boundary (Fig. 2). The zone covers a productive area of 1,550 acres at an average depth of 9,000 ft. The net sand section is thin, having a maximum thickness of 10 ft and an average thickness of 5.5 ft (Fig. 3). Average reservoir rock properties, as determined from core analysis are presented in Table 1. Shaliness of the zone is illustrated by the typical electric log of the 5th Cockfield section in Fig. 4. Because of the thinness and poor quality of the sand, the 5th Cockfield was developed on a wide spacing. Eleven completions were made in the zone between 1956 and 1961 for an average density of 140 productive acres per completion, with the majority of wells being dually completed in the 2nd and 5th Cockfield. Average initial potential for producers in the 5th Cockfield was 130 BOPD. Oil from the 5th Cockfield has a stock-tank gravity of 50.5" API; with a bubble-point pressure of 3,515 psia, the oil was highly undersaturated at the original reservoir pressure of 4,675 psia. Pertinent properties of the oil as determined in the laboratory are presented in Table 2. Because of the volatile nature of the crude, laboratory studies of its susceptibility to vaporization, as well as a vapor-liquid phase study, were conducted in addition to a standard PVT analysis. The vaporization study was conducted in a visual cell; a stepwise illustration of this work and laboratory measurements are presented in Fig. 5. Field

Jan 1, 1969

By I. Fatt

The chromalographic effect refers to the separation of constituents in a moving fluid phase which occurs when the phase is passed over a stationary phase, either solid or liquid, or large areal extent. In performing chemical analyses, the large surface contact areal is provided in the case of gas-solid chromatography by-passing the gas over a porous solid and for gas-liquid chronla tog rap hy by-passing the gas over a porous solid which has been coated with a nonvolatile liquid. The porous medium is usually placed in a cylindrical tube, although it may be a relatively flat sheet as in paper chroma-tography. According to one theory,&apos; the separation results from the difference in time spent by various components in the stationary phahe. A gaseous component which is conlpletely insoluble in the stationary phase, or is not adsorbed on the solid, proceeds unimpeded through the column. The versatility and power of the various types of chromatogruphy for the analysis of complex mixtures are evident from their widespread use in performing chemical analyses and from the large number of symposia, books, and papers devoted to the subject in recent years. Associated developments are the application of chromatography to the determination of heats of solution, heats of adsorption, solubilities orK-values, etc. Gas-liquicl chromatography analyses may be characterized by the manner in which the samples are introduced and displaced from the column&apos;. The three techniques used are known as (1) frontal analysis, (3) displacement analysis, and (3) elution analysis. The frontal analysis is performed by displacing an inert gas from the column by a continuous stream of the sample. As indicated in Fig. l(A), all but the first fraction leaving the column represent mixtures rather than pure components. The displacement analysis is conducted by filling the column with a sample, then displacing it with a more strongly adsorbed or absorbed vapor. The resulting chro-matogram is illustrated by Fig. 1(B). Finally, in the elution analysis a small sample of the mixture to be analyzed is introduced as a "plug" into a Rowing stream of carrier gas to give an elution-type chromato-gram. Fig. l(C). The three techniques may be used to provide models to which Row patterns in other physical situations may be compared. PURPOSE OF THE INVESTIGATION The purpose of this investigation was to determine whether natural reservoir rocks, when coated with crude oil, provided sufficient surface areas to give rise to a chromato-graphic separation of light hydrocarbon constituents. The tests were conducted with the light hydrocar-bons in the gaseous state. SANDS TESTED Two sands of widely different permeabilities have been studied: (1) Ventura core (16 md), and (2) Torpedo sandstone core (578 md).Table 1 gives a description of the cores that were tested. With the exception of one frontal-type flow arrangement which was performed on Core 3, the tests were of theelution-type in which a methane-propane mixture was eluted through the core by means of helium. The Ventura core, when impregnated with a weathered crude from the same field, became so impermeable to gas flow that it was necessary to conduct the tests on the granulated cores (Cores 2 and 3).

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

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

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

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

By R. O. Leach, O. W. Wagner

Because of unfavorable wetting conditions much residual oil is left when a porous material is Pushed by water. Methods suggested to change reservoir wetting to improve oil displncernrnt efficiency are generally expensitlr. The present 1aborator.y study was undertaken to gain an under.standinx of the factors which determine reservoir wettability, arid to find out if oil displacement efficiency might be improved by a wettahility change accomplished at low cost in on oil reservoir. Contact angle measurements were made on mineral surfaces using sevc.r~zl sets of reservoir oil and water samp1es. Results of the contact angle studies suggest that reservoir wetta-hility may he primarily determined by natural surface-active substances present in the reservoir fluids. The effect of changing sa1inity and pH of the water phase was studied. The re.suits suggest that gross changes in preferential wettability might be acc~o~npli.shed by injection of water containing simple chernicnls to alter pH or salinity in the reservoir. Such treatment could he much less expensive than injection of commercial surface-active agents. Waterflood tests have also been made using synthetic cores and oil and water having wening characteristics similar to those of reservoir fluids. Cores initially oil-wet were flooded in such a way that they were made prefermtial1y water-wct by the advancing flood water. This reversal in preferential wettability achieved greater oil displacement efficiency than when either oil-wet or water-wet conditions were maintained throughout the flood. For the systems studied, the higher the oil viscosity the greater the percentage improvement obtained over conventional waterflood recovery. This suggests that a flooding process making use of wettability-reversal may extend the oil viscosity range over which water flooding is attractive. Because a precise adjustment of reservoir wettability does not seem to be required, and because altering the pH or salinity in some reservoirs may be inexpensive, it appears that a waterflooding process employing wet-[ability-reversal could find .succesful field application. I NTRODUCTION The efficiency with which water will displace oil from a porous material is related to the nature of the capillary forces present. These in turn are controlled by the preferential wetting of the solid by the two fluids. Because of unfavorable wetting conditions, 30 per cent or more of the original oil in place may remain unrecovered in that portion of a reservoir flushed by water. This paper is concerned with the possibility of improving waterflood oil displacement efficiency by alterations in the wettability of the porous material. A laboratory study was made to gain a better understanding of the factors which control reservoir wettability, and to determine if the oil displacement efficiency could be improved by some inexpensive means of manipulating wettability of the porous medium. Contact angle measurements were made with several natural and synthetic oil, water and solid systems (1) to obtain a better understanding of how to duplicate reservoir wettability in the laboratory, and (2) to discover possible means for changing preferential wettability of natural reservoir systems. Flooding tests were also made in synthetic systems to determine if oil displacement efficiency could be improved by those wettability manipulations suggested by the contact angle measurements. Based on these studies a possible method for improving waterflood oil displacement efficiency is presented. This method involves causing an originally oil-wet porous material to become preferentially water-wet during the course of a water flood. The purpose of this paper is to present results of the laboratory studies. THEORY Rock surfaces in some oil reservoirs are believed to be covered with a firmly attached bituminous or other organic coating. Such surfaces would be preferentially oil-wet in the presence of both oil and water, regardless of composition of reservoir fluids. Other reservoirs are believed to contain rock surfaces not permanently coated with such materials, and which would be preferentially wet by water in the presence of water and oil free from surface-active substances. However, when the reservoir fluids do contain certain natural surface-active materials in sufficient quantity, rock surfaces acquire a degree of preferential oil wettability caused by adsorption of these natural surface-active materials on the solid. The equilibrium amount of these materials adsorbed per unit surface area is believed to depend upon their concentration in the bulk liquid phases.

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

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

Jan 1, 1969

By C. R. Sandberg, L. S. Gournay, R. F. Sippel

The effect of fluid-flow rate and fluid viscosity on oil-water relative permeability determinations was studied using the "dynamic flow technique." In this work relative per-nleability curves were obtained for homogeneous small core .samples from several sandstorle outcrop form-utions. Radio-tracers were used for the determination of fluid saturation and for the detection of saturation gradients. Cobalt-60 in the form of cobaltous chloride was used as a water-phase tracer in some of the experiments. Iodine-131 in the form of iodobenzene and Mercury-20-3 in the form of mercury diphenyl were used as oil-phase tracers in other experiments. Flow rates for each phase were varied within a range of 2.5 to 140.6 ml/hr. Oil-phase vis-cosities under flowing condition,> were varied from 0.398 to 1.683 cp. The relative permbilities ob-tained were found to he solely a function of saturation and independent of flow rare, provided there was no .Saturation gradient induced in the core sample by "boundary effect." Even though equilibrium with re-sped to flowing conditions was obtained at the lower flow rates, where a saturation gradient exists, this equilibrium is of a "contingent"-type rather than the "steady-state" equi-librium implicit in the relative per-meability concept. The only effect of increasing the oil or non-wetting phase viscosity was to decrease the flow rate required for the elimina-tion of the boundary effect. Fairly good agreement between experimentally determined and calculated values of the boundary effect was obtained when the non-wet-ring oil phase was the only Flowing phase. INTRODUCTION In the characterization of reservoir rock, as well as in the solution of reservoir production problems, it is most desirable to have reliablc rela-tive permeability measurements for the rock and fluids of interest. Many techniques have been developed for the laboratory determination of the relative permeability of both large and small core samples. In varying degree, difficulties attend the use of all of the methods. Each of the methods requires the metered flow of fluids of known viscosity through the core sample under conditions wherein the pressure drop in the individual flowing phases can be measured or closely approximated. Since the relative permeability is a function of the saturation and distribution of the flowing fluids, some means of obtaining such information is also required. The laboratory determination of the relative permeability of a sample should give the same relative permeability saturation relationship that would exist if the sample were in place in the oil reservoir. The attainment of this objective is complicated by the fact that the use of small core samples for laboratory measurements is usually necessary. As a result, errors in measuring the relative permeability arise from boundary effects at either the fluid input or output ends of the sample; no significant boundary effect is anticipated at the sides of the sample normal to direction of fluid flow. At the input end of the sample a boundary effect exists if the flow fluids are not com- mingled as they would be in the reservoir. At the output end of the sample a boundary effect exists as a result of discontinuity of the capillary properties of the system. Here the flowing fluids leave the small flow channels and enter a collecting system having much enlarged flow channels. Thus, there is a tendency for the wetting phase to be retained by the core sample because of the attraction of the core sample for the wetting fluid. No reference is made to these boundary effects in either Darcy&apos;s equation or in the relative permeability concept. As a consequence the saturation effect peculiar to the input end of the core sample may cause an error in the effective permeability measurement, since some portion of the core sample must be utilized for commingling of the flow fluids. At the output portion of the core sample a saturation gradient or a saturation change along the axis of flow results. Since relative permeability is a function of fluid saturation, this saturation gradient may cause a variable permeability measurement, especially since extent and character of the saturation gradient can vary with fluid-flow rate. The most convenient method of minimizing the boundary effect, especially at the output end of the sample, is to flow the fluids of interest through the core sample at such high rates that the capillary forces arc insignificant as compared to the viscous forces involved in the fluid flow. Since recent work of a preliminary nature has indicated that fluid dispersion in the porous medium increases with the flow rate, the high flow minimizes the boundary effect at the input end of the core sample by enhancing the mixing of the fluids.

By A. S. McLatchie, R. A. Hemstock, J. W. Young

Much attention has been given in the past few years to methods of increasing the recovery of oil from proven reserves. Numerous laboratories have made investigations to evaluate the possibilities of increased oil recovery by high-pressure injection of dry gas, injection of a propane slug followed by dry gas, solvent flooding. and by the injection at relatively low pressures of gases enriched with ethane and propane. Several years ago, Whorton and Kie-schnick&apos; published the results of laboratory studies on high-pressure gas injection which showed that, in the case of light oils, recoveries could be considerably increased by sweeping the reservoir with dry gas at pressures in excess of 3,000 psig. The increase was attributed by those authors to the vaporization of oil at the invading gas front and viscosity and solubility effects produced at the front by the invading gas. Later, Stone and Crump&apos; presented the results of a series of displacement experiments in which a light under-saturated crude oil was displaced from sand-packed columns with unusually high recoveries when the displacing gas was a rich condensate or a dry gas enriched with ethane or propane. They also conducted similar experiments at higher pressure on heavy undersa tu rated crude oil, although in this case the increase in recovery was not as great as in the case of the light oils. In both cases, viscosity reduction and swelling of the by-passed oil behind the invading gas front were believed to be responsible for the more favorable recovery of the original oil in place. Within the past year, other investigators", &apos; have presented the results of laboratory studies of miscible slug and solvent flooding recovery processes. This paper describes the laboratory methods developed for evaluating benefits to be obtained by enriched gas drive in specific reservoirs, and presents the results of several displacements of crude oils which possess a wide range of physical properties. The displacements were conducted at reservoir conditions of temperature and pressure. The apparatus used in the following experiments consisted of a sand-packed tube which served as the model reservoir, a cylinder of injection fluid, and a mercury pump. The stainless steel tube, 8 ft long and 0.53 in. I- was packed with a graded quartz sand. The high-pressure pump discharged mercury into the bottom of the injection fluid cylinder at a constant rate as low as 1 cc/hr. The average porosity and permeability of the sand column were 34.6 per cent and 3.25 darcies, respectively. In each case the tube was charged by the displacement of salt water by the reservoir oil under consideration. Recently a project was initiated to evaluate oil recovery by enriched gas drive from three oil reservoirs. Samples from a fourth reservoir (in this case high-viscosity oil) were studied with a view to obtaining recovery information of general applicability to low-grade, high-viscosity crudes. The oils from these four reservoirs exhibited a wide range of physical properties, and the reservoir conditions of pressure and temperature simulated in the laboratory represented several typical field conditions under which enriched gas drive might be employed. One of the reservoirs selected for laboratory displacement experiments produced an intermediate-gravity crude (29.6" API) by solution gas drive and a weak water drive. A low recovery in the range of 10 to 20 per cent was anticipated. Four displacement experiments were made to determine the effect of injection gas composition and initial gas saturation upon recovery of oil from the laboratory model. Two types of reservoir oil were used in these experiments. In the first run, a desaturated crude with a bubble point of 395 psi was charged to the tube. In the three subsequent runs, a simulated original reservoir oil was charged at 2,000 psi and produced by solution gas drive to 300 psi before beginning gas injec-

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

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

Jan 1, 1970

By J. J. Taber

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