Search Documents

By N. L. Hancox

Previous researchers have used the cutting and rejoining model to represent a porous body when attempting to place Archie's law on a theoretical basis and calculate rock permeabilities (F = Ø-m, where F = formation factor, = - porosity and m = a constant). In using the model, a porous solid is assumed to consist of sections with parallel sides containing straight. continuous pores of total porosity 4 (pore to total volume or, for a homogeneous body, pore to total area). The joining of pores at an interface between two sections is assumed to be random which results in a common pore area Ø2. This reduced area Ø2 is then taken as that area which limits electrical conduction or fluid flow in the body. Agreement between experimental results and results calculated using a model is not always good Because of this, various modifications of the model are presented — for instance, the random nature of the joining is questioned. Formation factor measurements on specially constructed porous specimens are described in this paper. it is concluded that the common pore area limits flow completely only if sections containing straight pores are very thin, and that, in most natural materials, the common pore area at an interface must be greater than that given by random joining. However, the model is still useful in permeability calculations, provided that the fractional area limiting flow is determined from formation factor measurements. INTRODUCTION Two important properties of porous media are permeability k, which describes the ease with which a fluid can flow through a porous medium, and formation factor F, defined as the ratio of the resistivity of a saturated specimen to the resistivity of the saturating liquid which describes the electrical properties of the saturated material. k is defined by Darcy's law,l k= µqL/P Ae . . . . . . . . . . . . . . . [1] where µ = liquid viscosity, q = flow rate, L = length of specimen, p = pressure across the specimen and = - the effective area over which flow occurs, F is given experimentally by Archie's law,2 F = Ø-m ,..............(2) where 4 is porosity defined as the pore to total volume or, for a homogeneous body, the pore to total area; and m lies between 1.3 for loose sand and 2 for cemented sand. Consider a porous body of total cross-sectional area A and length L, saturated with a liquid of resistivity R and with an effective area for conduction of Ae. Then, from the definition of F given previously, F=A/Ae , ................(3) or, from Eqs. 2 and 3, Ae = A Øm ..............(4) Thus, in any attempt to calculate k or F, some assumptions must be made concerning the pore structure of the body and a value derived for A,. To determine Ae, Childs and George3 and Wyllie and Gardner4 suggested an elementary statistical representation often referred to as the cutting and rejoining model of a porous body. The body is assumed to be composed of two or more parallel-sided sections of porosity 4, containing straight, continuous pores with a common pore area of Ø2 at an interface between two sections. The value is given by a random joining of the pores in either section. Two such sections,

Jan 1, 1969

By E. F. Johnson, V. O. Naumann, D. P. Bossler

A method is presented for calculating individual gas and oil or water and oil relative permeabilities from data obtained during a gas drive or a waterflood experiment performed on a linear porous body. The method has been tested and found both rapid and reliable for normal-sized core samples. INTRODUCTION Individual oil and gas or oil and water relative permeabilities are required for a number of reservoir engineering applications. Chief among these is the evaluation of oil displacement under conditions where gravitational effects are significant, such as a water drive or crestal gas injection in a steeply dipping oil reservoir. Numerous proposed methods of obtaining relative permeability data on reservoir core samples have been too tedious and time consuming for practical use, or have yielded questionable and sometimes inconsistent results. A method has been developed by which the individual relative permeability curves can be calculated from data collected during a displacement test. The method is based on sound theoretical considerations. Using this method, with a properly designed experimental procedure, relative permeability curves can be obtained using core samples of normal size (i.e., 2 to 3 in. in length and 1 to 2 in. in diameter) within a few days after receipt of the core. In a recent publication D. A. Efros1 describes an approach to the calculation of individual relative permeabilities that is based on the same theoretical considerations. We believe the approach described in the present paper is more adaptable to practical application than the method implied by Efros. In addition, comparisons with independently determined relative permeabilities are furnished to substantiate the reliability of the new method. DERIVATION Previously the theory of Buckley and Leverett as extended by Welge2 as been used to calculate the ratio of relative permeabilities. In the derivation which follows, this theory is further extended to permit calculation of the individual relative permeabilities. The theory assumes two conditions which must be achieved before the method is applicable. They are that the flow velocity be high enough to achieve what has been termed stabilized displacement,3 and that the flow velocity is constant at all cross sections of the linear porous body. In stabilized displacement the flowing pressure gradient is high compared with the capillary pressure difference between the flowing phases. The high pressure drop insures that the portion of the core in which capillary effects predominate will be compressed to a negligibly small fraction of the total pore space. The assumption of constant flow velocity at all

By R. E. Gilchrist, R. F. Nielsen

When a gas is displaced by mother in a porous medium, and there is a relative immobile liquid present. there is a transition zone in which the gas composition varies from essentially that of the original gas to that of the injected gas. The change in extent of this zone with time may be governed by a mixing process in the direction of flow and by the rate of establishment of gas-liquid equilibrium, if the gases are appreciably soluble in the liquid. The relative importance of these two factors in deterri~ining the growth of the transition zone is discussed, with appropriate equations and experimental data. INTRODUCTION The various solution methods of oil recovery involve a transfer of the components of a mixture, whether within a single phase or between a liquid and gas. The similarity of a natural porous medium to systems used in separation and chromatographic processes makes it seem reasonable that certain concepts used by chemical engineers in "mass transfer" might be applicable to petroleum reservoirs. One approach, the HETP (height equivalent to a theoretical plate), involves a distance or length of time over which equilibrium might be con. sidered as established. Another approach involves direct use of a rate of mixing or rate of establishment of equilibrium. The latter approach has been preferred in relation to production problems. DISCUSSION AND PROCEDURE When a gas or liquid in a porous medium is displaced by one of different composition but miscible with it there will be set up a transition zone in which the composition varies from cssentially that of the original fluid to essentially that of the injected fluid. If the fluid being displaced is a gas, and there is a relatively immobile liquid present, the transition zone may be

By R. I. Lahring, C. G. Inks

A surfactant waterflood, using a block copolymer of ethylene and propylene oxide (PLURONIC L64*), wa., started in 1956 in Kansas to determine the commercial feasibility of such a process. The production formation was the Bartlesville sand at a depth of 1,400 ft. Leares in the test area had attained both primary and secondary production peaks, and were on a normal decline when the tests were initiated. The injection wells were taking water at a reasonable pressure and none of the operating conditions were particularly unusual. The leases were selected primarily because of normal operating conditions, and not because of any serious problems demanding remedie>,. Oil production and water injection data were recorded during the 10-year test period, and their compilation indicates that the surfactant resulted in an increase of about 9 percent in secondary oil production. Other benefits derived from the surfactant injection were lower injection pressure that resulted in reduced power consumption and reduced erosion of the pumps; and lower calcium scale buildup. These fringe benefits more than offset the cost of the surfactant used. Introduction During the past several decades, considerable research has been carried out on secondary recovery of crude oil that remains in the producing formation despite the most efficient, current primary production practices. Methods are constantly sought of increasing recovery of oil, while at the same time improving the economic efficiency of operations. One method that has received intensive study over the years is the use of a surfactant as a waterflood additive. Laboratory tests show conclusively that dilute solutions of surfactants remove considerably more oil from uncon-solidated sand cores than does water alone.'- The loss of surfactant by adsorption at solid-liquid interfaces has caused serious doubt that surfactants can be used economically as waterflooding additives. This is especially true with anionic and cationic surfactants, and to a lesser extent with nonionic surfactants. This work was undertaken to determine the effectiveness in increasing oil production by injecting water containing small quantities of a nonionic surfactant. The surfactant selected was PLURONIC L64*, a block copolymer based on ethylene and propylene oxide, having the following structure The selection of the test site, a typical shoe-string sandstone type of field located at Batesville pool, Toronto District of Woodson County, Kans., was made because of its normal primary and secondary recovery operating conditions, and not because of any serious operating problems. When the tests were begun, the field was experiencing the customary post-peak production decline of a waterflood. Using for the tests a field that was on a normal production decline, while maintaining constant all variables other than the addition of the surfactant. was of varamount importance if success or failure was to be attributed to a particular cause. It was reasoned that if a surfactant was added to the injection waters and nothing else was done to influence the oil production or injection of water, the effects of the surfactant on the flood could be determined precisely. Production History The first well in the tested field was drilled in 1935. Gas repressurization began immediately, but it was not until Aug., 1938, that the major portion of the field was developed. During the period of repressurization, 1,700 MMcf of gas was injected. As can be seen from the curves of Fig. 1, the production peak from gas injection was attained in 1937 and a normal production decline followed. Early in 1942, water injection wells were drilled and injection of water Was initiated. The amount of water injected increased as more injection wells were drilled until a peak injection rate of 2,700 BID was reached in

Jan 1, 1969

By A. R. Gregory, A. L. Podio, K. E. Gray

Dynamic elastic properties of dry and water-saturated Green River shale samples were computed from compressional- and shear-wave velocity measurements. P- and 5-wave velocity measurements were made in three mutually perpendicular directions with respect to the bedding planes. Measurements were also made in several different directions by varying the angle between the bedding planes and the direction of propagation of the wave for angles of 0, 30, 45, 60 and 909 The oriented samples were subjected to both confining pressure and axial loads, in excess of the confining stress, in the direction of propagation. In general, P- and S-wave velocities increased with increasing stress levels, with a corresponding increase in Young's modulus. Water saturation caused the P-wave velocity to increase and the S-wave velocity to decrease. Elastic moduli decreased upon saturation, except for Pois son's ratio, which increased, indicating some degree of weakening of the material. The samples showed a moderate degree of anisotropy; this was to be expected from the laminated nature and shallow occurrence of Green River shale. INTRODUCTION This paper presents some results of an experimental determination of the elastic coefficients of anisotropic materials (in particular, finely layered rocks and minerals such as Green River shale) from measurements of dilatational and shear-ultrasonic-wave velocities. Ultrasonic techniques have been used extensively in nondestructive testing. Several methods have been proposed by Mcskimmin, 1 and some of these have been used successfully to measure ultrasonic velocities in rocks. Hughes and Cross, 2 Wyllie et al., 3 and Birch 4,5 developed pulse first-arrival techniques for the measurement of dilatational and shear velocities. Williams and Lamb 6 used the method of cancellation of a traveling wave, which was later modified by Myers et al.' and perfected by McSkimmin. 8 Although this method is highly accurate, it has not been used as widely as the pulse-transmission methods recently reported by Jamieson and Hoskins 9 King,10 and Mattaboni and Schreiber. 11 It has been common practice to use some form of crystal transducers, either quartz or ceramic, that has been cut or polarized in different directions in order to generate either compressional or shear waves. However, accurate determination of shear wave velocities has been difficult due to problems that arise in obtaining a pure shear wave from cross-polarized crystals, which usually also generate a small amount of compressional energy. As reported by Gregory, 12 this energy can be seen as a long precursor preceding the sharp break of the shear first arrival. The need for generating pure shear waves led to interest in mode-conversion techniques, which are based upon conversion of the mode of vibration through wave reflection or refraction at a discontinuity. Arenberg13 showed that for certain materials and for certain angles of incidence it is possible to generate pure shear modes by reflection at a boundary. Jamieson and Hoskins 9 used a pyrex glass-air interface for generating pure shear waves, and King 10 used this method successfully for measuring shear-wave velocities in rocks. Gregory l4 arrived at a similar result by refraction of a wave at an aluminum-oil interface. A plane compressional wave, traveling in the oil phase, is incident on the aluminum at an angle larger than the critical angle for compressional waves, and thereby generates a purely transverse, plane-polarized wave in the aluminum. During the last few years methods have been developed that allow the simultaneous determination of shear and compressional velocities in solids.

Jan 1, 1969

By D. O. Wyble

Five water-driver fields, for which pressure predictions had been made within a few years after discovery, were re-studied when four to five years of additional pressure and production data became available. It was fortnd that forecasts of the pressure in water-dri~fe reservoirs may be subject to considerable unccrtainty. Predictions car be erroneous when there is insufficient presswe history to define the cllaracteristics of the aquifer, and when there are errors in the pressure, production or geologic data. INTRODUCTION The efficient operation of an oil or gas reservoir requires that the important characteristics of the reservoir be determined soon after its discovery. It is especially desirable to evaluate the pressure-production relationship quickly, for this relationship provides a basis for deducing the type of production mechanism that is operative, and for making pressure predictions. The pressure predictions, in turn, provide a basis for calculating efficient rates of production, determining whether water or gas should be injected, etc. The importance of pressure prediction is such that forecasts should be made as accurate as possible, and the limits of error well defined. Historically, pressure forecasts for some fields have proved very accurate, while for others the forecasts were so erroneous that they resulted in the adoption of wrong, or at best, unnecessary production practices. Pressure predictions made for water-drive fields are especially subject to error. Accurate forecasting for these fields requires that much be known not only about the reservoir proper but about the adjacent water-filled rock as well. The water-containing rock, or aquifer, may extend for many miles, and often there is only meager data regarding its thickness, permeability and porosity. Much of the information must be inferred from the previous pressure behavior of the reservoir. A brief pressure history may provide little information about the characteristics of the aquifer. Hence, pressure predictions made early in the life of water-drive fields can be incorrect. The work described in this paper was done to determine the degree of reliability which can be expected of pressure predictions made within a few years after water-drive fields are discovered. Five fields, for which the original pressure predictions had been based on less than seven years' history, were re-examined with the aid of four to five years of additional pressure data. Where the predictions were inaccurate, the reasons for failure were ascertained. The studies were accomplished by use of the electric reservoir analyzer.' The analyzer was first adjusted to repeat the original pressure-production match using the same values for the various parameters as were used in the original study. Then, the actual production rates experienced during the time elapsed since the initial study were set into the analyzer, and the pressure behavior indicated by the machine for this period was compared with that which actually OCcurred. Where the pressure profile traced out by the analyzer deviated from the observed field behavior, the analyzer was readjusted to obtain a better over-all match. The old and the new solutions were compared to determine the reasons for any differences. COMPARISON OF PREDICTED RESERVOIR PRESSURE BEHAVIOR WITH ACTUAL PERFORMANCE Field A The original pressure-production match is entirely adequate as evidenced by the pressure match of the field's history subsequent to Jan., 1952 (see Fig. 1)'.

By L. Yarborough, L. R. Smith

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

Jan 1, 1969

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

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

By 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 D. R. Parrish, F. F. Craig

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

Jan 1, 1970

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

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

By 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 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 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"&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 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.