Search Documents

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.' 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" ' ratio.""' 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's, improvements in forward combustion by addition of water injection have been discussed, particularly in the patent literature, with increasing frequency.'-'' 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' 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' calculating the material balance equation for a partial water-drive reservoir completed in an aquifer has been described in the literature." Pressure interference i'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.' 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' 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' 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,' 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,"' 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,' 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.' 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.' 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.'," 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' 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,'.2.1.'."." 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.' : 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.*'" Exotic recovery processes have been considered primarily for ' 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'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' 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.' 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' 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,' 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'. 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 W. H. Somerton

The effects of drilling variable on rotary drilling rates and efficiencies have been studied by a series of laboratory drilling tests. TWO-cone 1.25-in. diameter hits were used to drill vertically upwards into rock samples at corltrolled weights and rates of rotation. Shale, sandstone and specially prepared concrete samples were used in this study. Power input to the ririllirig systeln was measured and drilling chips collected for energy — .size reduction studies. Reasouably good correlations between drilling variables and rates of penetrotion were found. Qrrantities that are difficult to evaluate include rock strength parameters and the effects of bit wear. Effects of bit size and geormetry require further investigation. Analysis of the drilling chips confirmed the premise that, for rocks containing two or more mineral constituents of diflerent strengths, a greater amourrt of rock breakage occurs in the weaker constituent. Drilling conditions which required greater amounts of energy produced finer drilling chips. As bit tooth wear progressed, drilling chips became finer. Efficiency of rotary drilling as a rock breakage mechanism was extremely low. Comparison was made with theoretical energy requirements and with energy requirements for size reduction by comminution methods. INTRODUCTION Many technological developments and innovations have made possible the successful drilling of oil wells by the rotary method to depths exceeding 20,000 ft. Bigger, more powerful rigs, better steels, improved bit design and more careful control of the circulating system have all contributed to this success. Despite these advances comparatively little is known regarding the basic mechanism of rock breakage by the rotary drilling process. Future improvements and, in particular, reduction of drilling costs will undoubtedly require a clearer understanding of the variables controlling effectiveness and efficiency of rotary drilling. Rotary drilling is inherently an inefficient means of producing rock breakage. Loss of energy in transmission from the surface to the cutting mechanism may be large, thus limiting the total amount of energy which may be applied to rock breakage. Increased efficiency may be realized by improved energy transmission methods such as the turbodrill. Conversion of transmitted energy to rock breakage by use of rotating cutter teeth is likewise inefficient. Improvement of cutter efficiency has been obtained essentially by empirical means —changes in bit tooth size, shape and spacing, and tooth deletion to give clean bottom-hole patterns. The present work was undertaken to investigate the factors controlling rates of bit penetration under laboratory drilling conditions. Laboratory drilling apparatus is shown in Fig. 1. The effects of rock strength and bit wear on drilling rates were investigated. Drill cuttings were analyzed to determine the character of breakage and to compare this with other methods of producing rock breakage. VARIABLES INVOLVED IN ROTARY DRILLING Some of the variables controlling bit penetration rates can be studied conveniently by dimensional grouping. R/DN f(F/D2S) or

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

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

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

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

By C. R. Clark, G. Swift, D. S. Roberts

The purpose of this investigation was to measure PVT behavior of various types and combinations of heavy hydrocarbon components from the paraffinic, naphthenic and aromatic classes where the relative proportions Of the various components were selected to approximate those of natural gas or gas condensate systems. Synthetic mixtures were used so that the compositions Of the various components could be measured accurately. n-Heptane, methylcyclohexane and methylbenzene were used to represent the paraffinic, naphthenic and aromatic components of the heavy fraction. In the mixtures studied, the heavy fraction composition was held constant at O. O5-mole lraction, the balance being methane except for one mixture where 0.10-mo1e fraction of an intermediate component, propane, was added. The PVT data for the mixtures were determined in a variable-volume, constant-mass apparatus at psuedo-reducedtemperatures from 1.84 to 2.00, over a pseudo-reduced pressure range from 2.3 to 12.0. The experimental results showed that, regardless of the type of heavy material (paraffinic, naphthenic, aromatic, or combinations thereof) mixed with methane or methane and propane, the compressibility factors at equal values of pseudo-reduced temperature and pressure varied by less than 2.2 percent. INTRODUCTION PVT data is used for natural gas and gas condensate fluids in determining reserves of reservoirs and performance of wells, in metering produced fluids, and in recombining samples for laboratory studies. While there are vast amounts of PVT data reported for natural gas and gas condensate systems from which useful correlations have been developed, ' the compositional analyses for these systems and the resultant correlations typically were made with components analyzed through some arbitrary carbon number, usually C6, with the residue reported as a lumped "heavy" fraction. The heavy fraction was stoichiometrically recombined on the basis of an apparent molecular weight to give the final compositional analysis. It is virtually impossible to make a systematic analysis of the effect that variation of the relative amounts of paraffinic, naphthenic and aromatic constituents of the heavy fraction of natural systems might have on PVT behavior. Because there is considerable latitude in the relative amounts of these constituents, one speaks of crudes or condensate liquids as being "paraffin base", etc. It is of interest to determine if changes in the relative amounts of these three types of hydrocarbons in the heavy fraction cause significant changes in PVT behavior. If so, steps should be taken to describe better the nature of the heavy fraction in PVT correlations. If not, more confidence can be placed in the correlations presently employed. DESIGN OF THE EXPERIMENTAL INVESTIGATION We attempted to determine if changes in the base of the heavy fraction would cause significant changes in the PVT behavior of gas or gas condensate fluids. Synthetic mixtures were used to enable us to systematically vary the base of the heavy fraction and to analyze accurately the fluids. We found that the heavy fraction of naturally occurring condensate systems seldom exceeds 0.05-mole fraction. Thus, the mole fraction of the heavy fraction in the synthetic systems studied was set arbitrarily at 0.05, to maximize whatever deviations in PVT behavior that might occur. Further, we judged that since the C7 hydrocarbons normally are present in greater quantities than higher carbon number hydrocarbons in naturally occurring systems, the use of n-heptane, methyl cyclohexane, and methylbenzene to represent the paraffinic, naphthenic, and aromatic species in the heavy fraction would be most appropriate. Finally, to enable us to compare our

Jan 1, 1970

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

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