Search Documents

By P. E. Schauer

During the initial nine-year deple-tion of the McKamie-Patton pool gas-condensate reservoir, the pressure behavior and fluid analyses indicated pressure maintenance by gas injection would be attractive relative to continued pressure depletion. As a result, a gas cycling program was inaugurated in Aug., 1949, which has since maintained reservoir pressure. In order to evaluate the benefits of this program, the recoveries for natural depletion and for cycling followed by blowdown were determined by calculating the fractional balance of the molar quantities of vapor and liquid remaining in the reservoir and produced at the surface for various stages of depletion. The results of the investigation show theoretical ultimate recovery due to the cycling program will be 42 per cent greater than would be realized under pressure depletion type operation. INTRODUCTION The benefits of gas injection pressure maintenance programs in condensate reservoirs have been recognized and qualitatively described in several papers1-8 during the past 15 years. Some of these papers also describe methods for calculating the increased recovery realized from cycling; however, as a rule, hypothetical reservoirs were used as examples".'.. This paper presents an evaluation of the pressure maintenance program for the McKamie-Patton condensate reservoir and includes a description of the fractional balance method employed to determine the benefits. Reservoir Development and Characteristics The McKamie-Patton pool, Lafayette County, Ark., was opened in June, 1940, with production obtained from the oolitic Smackover lime at a depth of approximately 9,300 ft. The structure is an anticline 8 miles long by 1 to 2 miles wide, and is one of many productive Smackover structures in Southern Arkansas. Of predominant importance is the condensate-bearing gas cap portion of the reservoir, which covers approximately 4,000 acres, attains a maximum thickness of some 200 ft, and has an average thickness of 90 ft. Development of the cap was accomplished with 21 wells drilled on 160 acre spacing. An isopach map of the cap is shown in Fig. 1. In addition, an oil rim is present along most of the periphery of the pool. Production practice includes field separation of black oil, and a high pressure gathering system which carries total well stream from gas cap wells to a plant where condensate is separated and the gas is processed for recovery of natural gasoline, butane, and propane. Since the gas contains appreciable hydrogen sulfide, it is removed by amine contact, after which free sulfur is produced by the Claus process. Core analyses from 1,767 gas cap formation samples averaged 14.2 per cent porosity and 379 md dry air permeability. Connate water was restored on seven samples of gas cap rock by air-brine capillary techniques, from which the interpreted water saturation averaged 28.17 per cent. During the early life of the field, it was believed the recovery mechanism, based on several of the other regional Smackover reservoirs, might be a natural water drive. However, subsequent production and pressure history revealed that the gas cap would not be produced by water drive, but would behave as a closed reservoir. In this respect, during the nine-year period of natural depletion from initial production to Aug., 1949, the cap production of 84.3 billion

Jan 1, 1958

By H. J. Ramey

A general solution is presented for the transient temperature distribution caused by radial movement of a cylindrical heat source through a homogeneous medium of infinite extent. This problem represents a highly simplified model of the movement of a combustion front during the thermal recovery of oil. Numerical solutions are presented for a heat source moving at a constant velocity, and a velocity inversely proportional to the radial location of the heat source. Numerical solutions are presented for both finite and infinite vertical thickness of the heat source, i.e., with and without vertical heat losses. The numerical solutions are used to estimate the fuel concentration needed to maintain a combustion front during the thermal recovery of oil. The sensible heat carried to the combustion front by the injected gas stream is discussed in regard to the fuel concentration. Numerical solutions for a heat source (combustion front) of finite vertical thickness indicate the field conditions that may be necessary to sustain the temperature of a combustion front above the ignition temperature of the fuel over considerable distances from an injection well. The results also have implications in regard to the quantity of heat required during ignition. The results of this computation apply to a highly idealized model of a thermal recovery process. But the results may be used as a guide in engineering consideration of the thermal recovery process. INTRODUCTIO N The thermal recovery of crude oil has received considerable attention since the publication by Kuhn and Kochl of laboratory and field tests of this method of oil recovery. Other publications2,3,4,5,6, since that time have dealt with various features of the thermal recovery of oil and have illustrated the extremely complex nature of this process. Vogel and Krueger7 described an electric analog computer designed to solve a moving heat source problem which was analogous to movement of a combustion front during the thermal recovery of oil. The heat source was maintained at constant temperature and was a cylindrical source of infinite height (no vertical heat loss) moving radially at either constant velocity, or at a velocity inversely proportional to the radial location of the heat source. Jenkins and Ramey2 presented analytical solutions to a similar heat conduction problem and pointed out the possible importance of vertical heat losses. The following presents a general solution to the transient heat conduction problem introduced by Jenkins and Ramey4 as well as numerical solutions for a wide range of conditions possibly similar to those that might exist during field operation of the analogous thermal recovery process. DESCRIPTION OF PROBLEM The problem considered in this paper is determination of transient temperature distribution caused by a cylindrical heat source of infinite or finite vertical height, moving radially through an isotropic medium of infinite extent. It is further assumed that the heat flux generated at the surface of the moving heat source may be a function of time such that the fuel concentration required may be constant, or permitted to vary, and that heat is generated at the surface of the heat source only. In regard to the thermal recovery of oil the assumption that heat is generated at the surface of the front is equivalent to assuming that the thickness of the combustion zone is infinitely small, or that the reaction rate between fuel and oxidant gas is infinite. The last assumption does not appear stringent. See laboratory information published by Martin, et al3 Under these assumptions the temperature distribution throughout an infinite medium caused by a heat source moving radially with an arbitrary velocity may be described by The heat generation function or heat flux at the surface of the source, +(t), is Note that the fuel concentration may be defined as a function of time. Eq. 1 is derived in the Appendix.

By W. T. Cardwell

Physically absurd, triple-valued saturations appear in the straight-forward solution of the Buckley-Leverell equations for the displacement of oil by water or gas. From an engineering viewpoint, the triple value causen no difficulty. It is well known how it may be compensnted in order to obtain physically meaningful, numerical results. From a scientific viewpoint, the question still arises: What did the triple value mean? This paper explains how and why the triple value arose in non-capillary Buckley-Leverett theory. The discussion should serve as a background for the understanding and use of the modern method of characteristics in di.~pkzc(>ment theory. INTRODUCTION Displacement theory was introduced to petroleum technologists in 1 941 by Buckley and Leverett.' Their first equation, which they wrote down without derivation, was equivalent to the following: ?s/?t= - qr/F ?f/?r (1) where S = saturation of displacing fluid (a fraction) r = time q,. = total volumetric velocity of both the displacing phase and the phase being displaced (usually, but not necessarily, assumed to be a constant) F = porosity (a fraction) f - fraction of the total volumetric velocity that is the volumetric velocity of the displacing phase (assumed to he a function of S only x = distance. Buckley and Leverett called Eq. I "a material balance equation". Actually it is derived from both a continuity equation, or material balance equation, and Darcy's law. Buckley and Leverett jumped from their Eq. 1 immediately to their Eq. 2, which was the equivalent of the following: (?x/?t)x = qr/F df/ds = qr/F f'(s) (2) The jump from Eq. 1 to Eq. 2 involves apparently simple mathematics, but in that apparently simple mathematics lies a subtle point that is part of the key to the meaning of the triple value to be discussed here. Eq. 2 may be readily interpreted as saying that a point of constant saturation (on a saturation-vs-dis-tance curve) moves with a constant velocity that is proportional to the total volumetric rate, inversely proportional to the porosity, and is otherwise a function of the saturation itself. So that if one knows the derivativc with respect to saturation of the fractional flow function, f(S), for each saturation, one knows the velocity of each point of the moving saturation-vs-distance curve. The function f(S) can he determined experimentally and its derivative f'(S) can be calculated. It turns out that the experimentally derived f'(S) is not a monotonic function of S, but instead it has a definite maximum and declines away at both high and low saturations. This in turn means, in accordance with Eq. 2, that on a moving saturation-vs-distance curve, ccrtain intermediate saturations will travel faster than the saturations either higher or lower. The result is indicated in Fig. 1. In Fig. 1, the abscissas represent distances along a column of porous medium, which column is assumed to be uniform in cross-section perpendicular to the ab-scissal direction. The ordinates represent fractions of pore space occupied by the displacing phase, which may be either water or gas. It is assumed that the saturation is uniform over all planes perpendicular to the abscissal direction. The saturation is a function only of time and one space variable. Time variation of the saturation is represented by change in shape of a curve such as curve AGJ, which represents the initial saturation-vs-

By H. A. Koch, R. L. Slobod

This paper discusses a new oil recovery process called the "miscible slug process." This process involves the injection of propane or LPG into the reservoir prior to gas injection. The operating conditions are maintained so that the slug is miscible with both the reservoir oil and the injected gas phase. Thus, a miscible phase displacement is achieved, and high recoveries are obtained from the area of the reservoir contacted. The most striking discovery of the laboratory study which included displacements from cores up to 123 ft in length was that relatively small slugs of LPG are effective over reservoir distances in the miscible displacement of oil. This makes possible the commercial application of the process. The important factors controlling the size of the slug are: (I) reservoir length, (2) reservoir fluid composition, and (3) reservoir pressure at the displacement front. Of lesser importance are the effects of injection rate and porous medium type. Questions which have not been completely answered deal with the effects of gravity and reservoir inhomogeneities on the process. INTRODUCTION One of the major problems facing the oil industry today is recovery of the large fraction of the discovered oil which will remain unrecovered unless some new processes are found. Water flooding has been one successful method used to increase recovery over natural depletion. But even successful water floods may leave 25 to 40 per cent of the pore space filled with residual oil. Recently a process which involves the use of thermal energy has been reported1,'*' which may prove useful in increasing recovery from shallow reservoirs containing highly viscous crudes. This paper deals with a new oil recovery process involving use of miscible phases to displace the oil from the reservoir. Miscible phase displacement oil has been an intriguing idea because the elimination of capillary effects in the reservoir leads to 100 per cent recovery in the areas contacted by the miscible displacing phase. The big deterrent in the past to the practice of miscible phase displacement has been the high cost involved in using large volumes of solvent in the process. The use of high pressure gas to achieve a miscible phase displacement of reservoir oil has been reported4,5 as one economically feasible method. However, this process calls for the maintenance of pressures above 3,000 psi at the dis-placement front. The process described in this paper is not restricted to high reservoir pressures; the only qualification being that the slug of material which is injected into the reservoir must be miscible with both the displaced oil 2nd the displacing phase. Two such cases are the use of an alcohol slug followed by water to achieve a mis-cable water flood, and the use of a low boiling hydrocarbon slug such as propane followed by gas to achieve a miscible gas drive. This paper will deal specifically with the propane dug process. This process will obtain a miscible phase displacement at pressures as low as about 1,200 psi at a reservoir temperature 01 150°F. The process is very simple. The propane displaces oil completely at one front because it is miscible with the oil, and gas displaces the propane completely at a second or trailing front because it is miscible with the propane. A small band of propane should thus remain wedged between the gas and oil phases. This idea has been called the "miscible slug process."6 Unfortunately, the propane slug will become diluted as it flows through the reservoir. This results from a mixing action since the fluids flow through a complex reservoir porous medium rather than simple capillaries. Fig. 1 diagrams a miscible slug displacement. The effects of some of

Jan 1, 1958

By J. S. Aronofsky, H. J. Ramey

Potentiometric model studies of the five-spot water flood to determine the influence of mobility ratio upon injection or production histories are described. The results suggest a method for the correlation of data from field pilot tests for the estimation of an effective mobility ratio and its effect upon field water floods. The effect of the well spacing-to-well radius ratio upon production histories is presented. A method of accounting for this factor is suggested. The influence of injection (or production) rates and injection and production well pressures upon production histories is shown. The five-spot well pattern was investigated for water-to-oil mobility ratios of 0.1 to infinity. INTRODUCTION It has been well established that mobility ratio has an important effect upon flood patterns.1,2,3,4,5 It should be equally apparent that injection and production histories must also be influenced by mobility ratio. However, none of the articles referred to above mention the importance of being able to predict the effect of well pressures (or injection rates) during the entire time history of a waterflood operation. A complete description of production histories (the terms "injection" and "production" are synonymous for the steady-state model under consideration) would provide values for well pressures and interface positions at any time for a specified injection rate. conversely, if the waterflood problem is specified in terms of a constant injection well pressure, then it would be desirable to know the injection rate and interface positions at any time. This paper describes the results of potentiometric model studies which were used to analyze the influence of mobility ratio on injection or production histories. A method is presented for calculating production histories over a range of ratios of well spacing-to-well radius. The basic method described in this paper to evaluate production histories is to: (1) determine interface position as given in Ref. 1, and (2) measure electrical current and voltage drop for each interface position. It should be mentioned that the current state of knowledge on flood patterns as affected by mobility ratio shows that significant discrepancies exist between the results of various investigators2,3,4,5. It is beyond the scope of this paper to discuss the reasons for such discrepancies. However, the method of determining production histories presented in this paper could be applied to any of the swept-area patterns presented in the literature. All results given in this paper are based on a set of swept areas obtained by the potentiometric method. Potentiometric analyzer technique is adequately described in the literature, but the basic assumptions involved in work of this nature bear repetition. It is assumed that: 1. A distinct oil-water interface exists in the reservoir; only water is mobile in the region behind the interface (water zone) and only oil is mobile in the region ahead of the interface (oil zone). 2. There is a fixed mobility value, (k/µ) w associated

Jan 1, 1957

By B. H. Caudle, A. B. Dyes, R. A. Erickson

The study of the influence of fluid mobilities on the sweepout pattern resulting from the injection of gas or water has been extended to cover the production period which follows breakthrough of the injected material. Mobility ratios over the range common in field operations (0.1 to 17) were studied for several pattern floods (five-spot, staggered, and direct line drive). The experimental data required for these studies were obtained by the use of the x-ray shadowgraph technique using miscible oil phases of different viscosities in porous plate mode1s of a reservoir element. From the shadowgraph pictures obtained before and after breakthrough of the injected fluid, flowing ratios at the producing well and cumulative volumes injected were calculated. The method for applying such data in predicting field behavior is illustrated for a water flood of a five-spot. For this case a range of mobility ratios of 0.5 to 5.0 results in: (I) nearly complete (95 to 100 per cent) sweepout pattern efficiencies at abandonment conditions, (2) production after breakthrough being responsible for as much as one-third of the total recovery at the lower mobility ratios, and (3) a twofold variation in the operating life of the reservoir. In determining the merits of a secondary recovery operation, the fraction of the reservoir which is ultimately to be swept by the injected fluid is of major importance. Because of its significance a large amount of research has been applied to the study of areal sweepout efficiency. Muskat, Wyckoff and Botset1,2 have presented the classic analysis of the pattern efficiency at breakthrough for several injection patterns in which the mobility of the fluids ahead of the injected phase front was equal to the mobility of the fluids behind the front. Aronofsky3 as extended this analysis of recovery at breakthrough through numerical computations and by the stepwise use of the potentiometric model to include the effect of different fluid mobilities ahead of and behind the invasion front on the line drive performance. Fay and Prats' have also employed numerical computations to show the influence of one other condition of fluid mobilities on the five spot flood. In a previous paper' a new method of attack on this problem which employed x-ray and porous plates was described. At that time experimental results were also presented which described the effect of fluid mobilities on the reservoir area contacted by the injected fluid at breakthrough for the five spot and direct line drive patterns. Most of the literature pertaining to sweepout pattern efficiency has emphasized the pattern which is obtained at breakthrough of the injected phase. The increase in the sweepout pattern which occurs after breakthrough

Jan 1, 1955

By F. F. Craig

A series of both welter and gas puttern Hoods was made in the laboratory to study the oil recovery performances of .such operations. These tests were contlrrctt~tl on consolidated sandstone models, using oil. water, and gas. The model floods were scaled to reprotlrrce field performance under gas arid water five-spot injection. X-ray shadowgraphs permitted observation of the gross fluid movement within the models. A method was developed for applying the mobility ratio concept to water flooding and dispersed gas drives in a five-spot well pattern. The area1 sweep efficiency. tit breakthrough for dispersed gas drives is much higher than previorsly expected, Iying in the range of 50 to 100 per cent. A method is presented for predictitig the water-oil ratio performance of five-spot pattern water flood., in uniform sands. This method is verified experimentally for the cotidition of 120 free gas initially present and for va1ues of gas saturation normally encountered in fields following depletion operations. Prodrrction perfort?zut~cc for pattern gas injection is also predictable by this method. INTRODUCTION In field pattern hater flooding or dispersed gas injection operations, the displacing fluid (water or gas) is injected and oil is produced through wells which are, area-wise. small openings in a large container or rescrwir of oil. As a result, not all of the area between the injection and producing wells is necessarily contacted by the injected fluid by the time it first reaches the producing wells. The question arises in predicting oil recovery. as to what traction of the pattern area involved is contacted by either the injected water or gas due to the relative position of the injection and producing wells. This is the areal sweep efficiency. It is also desirable to know how well arrangement affects the oil production performance. A solution to this problem has been attempted by many means. Mathematical analysis¹,²,³, and electrolytic models', have been used to establish a breakthrough areal sweep efficiency of 73 per cent for a five-spot pattern in which the ability of the oil to flow ahead of the invading fluid is equal to the ability of the invading fluid to flow in the invaded zone. All of these studies simulate only one of many field conditions. Improved studies of areal sweep efficiency at breakthrough were made using the fluid mapper5, poten-tiometric models" and X-ray shadowgraph techniques. In these it was possible to study the effect of variations in mobility of the injected and displaced fluids, and so to cover a wide range of field conditions. The limitation of these studies is that, with these models, it was not possible to simulate the saturation gradient behind the flood front when oil is displaced by gas or water. In the most recent published technical paperS on the subject, a method was presented for predicting the oil recovery performance after breakthrough in a five-spot pattern water flood. The method involved correlations which were obtained experimentally using miscible fluids. One limitation of this work was that a saturation gradient, usually found in immiscible fluid displacement, was absent in these studies. Also the laboratory studies did not cover the situation of a water flood following oil recovery by depletion drive. as is the situation in most field water floods. The tests reported in this paper were made using oil. water, and gas, thus eliminating the simplifications of former studies. Also reported are laboratory studies on

Jan 1, 1956

By L. D. Wooddy, W D. Moore

Water drives may be natural or artificial. For the case of a natural water-drive reservoir, the volume of water influx corresponding to reservoir pressuve-production history may be calculated by means of the material bal-ance-unsteady state equation. For the artificial water drive or pattern injection water flood, the material balance equation is sufficient to determine the reservoir pressure-production, injection history. However, for both cases to reliably predict reservoir performance, it is necessary to calculate the saturation distribution of the reservoir fluids at different stages of depletion. This paper consolidates the various equations used to determine the performance of water-drive reservoirs into one approach that is applicable to any reservoir that has a natural water drive or is to be produced by water injection. The numerical procedure has been simplified considerably, and the solution is adaptable to either hand calculations or digital computer operation. INTRODUCTION A water-drive reservoir is one in which the predominant source of energy for production of oil results from the encroachment of water into the oil zone. The water drive may be either natural or artificial. Water invasion takes place by expansion of water in the aquifer as a result of the pressure decline transmitted from the oil reservoir. Because water is practically incompressible, a natural water drive requires an extensive aquifer containing a volume of water many times the volume of oil in the reservoir. Provided the permeability of the formation is reasonably high, a moderate decline of reservoir pressure usually causes adequate water influx. If the reservoir pressure decline becomes excessive because of low permeability or a high rate of production, part of the oil production may result from the liberation of dissolved gas which may cause excessive quantities of free gas to be produced along with the oil. In extreme cases, reservoir pressure may decline so rapidly as to cause the oil to be produced predominantly by dissolved gas drive, leading to abandonment of the reservoir before the advancing water has had a chance to flush the entire reservoir. In reservoirs that are not connected to extensive aquifers, an artificial water drive may be created by injecting water into or around the lower edges of a reservoir. Sometimes water is injected within the reservoir through wells interspersed among the producing wells. When injection is carried out in such fashion, the reservoir is in effect subdivided into a number of separate producing segments, in each of which the displacement of oil is simultaneously taking place, and the shapes of which are fixed by the well patterns employed. Because of the natural tendency of oil and water to segregate according to their relative densities, water will tend to migrate downdip regardless of where injected. Accordingly, the injected fluids can be retained within the intended reservoir segments only when conditions are unfavorable to their migration under the force of gravity. Pattern injections are usually confined, therefore, to reservoirs having low structural relief, or low permeability, or when production and injection rates are such that fluid movements in the reservoir are dominated by the pressure gradients rather than by the force of gravity. In the case of natural water-drive reservoirs and artificial water-drive reservoirs in which the water is in-

Jan 1, 1958

By J. C. Martin

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

By C. S. Matthews, W. R. Strickler, M. Prats

To gain an insight into the deviations from normal behavior caused by dip. a mathematical study has been made of single-fluid five-spot floods in dipping reservoirs for flooding techniques likely to be employed in such cases. Streamline patterns and sweep eficiencies are computed for a number of different cases, and the results, where possible, are presented in graphical form. An increase in sweep efficiency can be obtained by moving the central injection well updip, and the opti~num position is found for several cases. The applicnhility of the results is discussed. INTRODUCTION Because of the growing importance of water flooding in California and in other locations where large angles of dip are quite common, it is important to know at least some of the effects of dip on the performance characteristics of water floods in such instances. In the past, only reservoirs having relatively mild dips have been wateflooded. Consequently, little is known today as to the effect of dip on the performance characteristics of a flooding operation from the point of view of field experience; and because of the complexity of the problem, little, if anything, has been done on the problem from the point of view of theory. There are two distinct aspects of the role played by gravity in a flooding operation. The first is that of the preferential segregation of the reservoir fluids at the time when the flood is to commence. The second is the role that it may play in affecting the performance of the reservoir during the course of the flooding operation. In this paper the reservoir is assumed to be filled with only one fluid, so that there is no preferential segregation of reservoir fluids at the commencement of the flood. Moreover, the injected and in-place fluids are assumed to have the same characteristics, so that in effect we consider only a single-fluid problem. Single-fluid results have been used as a guide in the past for predicting and interpreting the behavior of regular five-spot floods in horizontal reservoirs. It is thought that single-fluid results might likewise be used for five-spot patterns in dipping reservoirs as a norm to which the effects of the multifluid problems might he referred.* In confining our attention to the single-fluid case, we therefore do not consider the preferential segregation of injected water in the oil bank, or of the oil bank in the unflooded zone. Instead we consider the deviation from the behavior of a horizontal five-spot which is caused by use of flooding techniques (i.e., the flow potential pattern at injection and production wells) likely to be employed during flooding operations in dipping reservoirs. In other words, we consider effects which are caused by the pattern of injection and production well pressures which we choose to apply. Performance characteristics are given for the five-spot spacing pattern as a function of two parameters; one is the ratio of gravity head to potential difference bctwccn the high updip injection well and the low downdip producer of the same five-spot, denoted by Rg; the second is the orientation of the five-spot pat-

Jan 1, 1956

By Wilbur H. Somerton

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

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

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

By Harvey T. Kennedy, Erasmo T. Guerrero

In this investigation, the effects of the surface tension of brine and the interfacial tension between oil and brine on the recovery of oil by water flooding of Woodbine sand were studied. Variation of the surface forces was obtained by the addition of small amounts of commercially available chemicals, which did not appreciably change the surface tension of the oil. Fourteen surface-active agents, which were found to be effective in lowering the surface and interfacial tensions, were further investigated to determine the extent of their adsorption by the Woodbine Sand cores. It was found that adsorption was sufficient in all cases to prevent any substantial concentration of chemical at the flood front. The flow system was, therefore, brought to equilibrium with definite concentrations in the flood water by circulation before each test. The data thus have no relation to the results obtainable by using the chemicals in full-scale water floods, but should be helpful in predicting the performance of surface-active chemicals which are not lost by adsorption, if such can be found. Flow tests were made on various concentrations 01 six chemicals selected from the 14 investigated for surface activity and adsorption, with and without the presence of gas in the core. It was found that lowering the interfacial tension had a tendency toward reducing the recovery of oil and that the reduction of the surface tension of the brine had little or no effect. A substantial beneficial effect, however, resulted from the presence of gas in the cores. INTRODUCTION The movement and distribution of crude oil in a porous medium are governed by such factors as surfacc and interfacial tensions, viscosity of the oil, pore size configuration, and wetting characteristics of the medium. The knowledge of the part played by each of these factors is limited. If the oil is assumed to be distributed in droplets whose movement is restricted only by the pore openings, the following equation may be applied: ?P= 2r(1/r1-1/r2) .......(1) In this equation ?P is the pressure drop necessary to force an oil droplet through the construction, r is the oil-brine interfacial tension, r1 is the radius of curvature at the forward and smaller end of the deformed droplet, and r2 is the larger radius. This equation neglects the effects of viscosity, as a static or near static condition is assumed. It would apply equally well for a bubble of gas in liquid if r were taken as the surface tension of the liquid. This relation shows that droplets can be more easily displaced if the interfacial tension is decreased. The displacement of oil from the surface of a sand particle may also be considered from the basis of selective wetting assuming that the oil droplets are smaller than the pore constructions and that the sand surfaces may be considered as planes. A relation of the forces involved can be derived from Dupre's equation to give: Wos = row (1 +Cos?)........(2) where Wos = Work of adhesion between the sand and the oil row = Oil-water interfacial tension 0 = Contact angle ros = Oil-solid interfacial tension This equation would indicate that the displacement of oil from a sand surface is more readily obtained when the water-oil interfacial tension is reduced provided the contact angle remains constant. Evidently. however, the values of Cos 0 may vary along with the water-oil interfacial tension due to changes in the other surfaces of the system. The above theoretical considerations suggest that perhaps even the static behavior of a dispersed oil phase in a porous medium is considerably more complex than may be inferred by the simple systems

Jan 1, 1955

By R. L. Slobod, B. H. Caudle, R. A. Erickson

In predicting the sweepout pattern efficiency to be expected in a .secondary recovery operation, the reservoir engineer is often confronted with a situation in which part of the producing formation lies between the last row of wells and the reservoir boundary. To date, there has been little or no reference in the literalure to sweepout pattern efficiencies in such areas. This paper presents the results, of a laboratory study in which the area contacted around the edge of a reservoir injected in a five-spot pattern was determined. The X-ray shadowgraph technique was used in this study. It was found that surprisingly large portions of the reservoir area lying outside the well network are conracted by the injected fluid before, abndonment conditions are reached. In fact, it can he said that at least 00 per cent of the circa lying outside the last row of wells and within a distance of one well spacing of these wells will ultimately be contacted by the injected fluid in the case of most water floods. Frorn the results given in this paper, it is possible to estitilate not only the extent to which encroachment into these areas will occur but also the injection volume necessary to produce the wells to abnndonment conditions. INTRODUCTION While considerable work has hern done in recent years1.?.'.''.' evaluating the fraction of the reservoir contacted as a function of the mobility ratio for several different well patterns, there has been practically no mention of what happens to the oil in that portion of the reservoir which lies between the last row of wells and the reservoir boundary. In fact, the assumption is generally made that little or no recovery can he expected From this area. This assumption is incorrect and in many cases will lead to seriouq errors in calculating reservoir performance. The fraction of the area outside the regular well pattern contacted in a five-spot injection system has been studied using the X-ray technique5 to follow the progress of actual floods in porous models of reservoir elemental areas. In this paper, the method for carrying out such a study and the results showing the fraction of the area in question swept as a function of mobility ratio and the produced fluids are presented. THE REPRESENTATION OF A RESERVOIR BY ITS ELEMENTAL AREAS She portion of the reservoir lying between the edge of a reservoir and a line connecting the outermost wells may be roughly designated as the part of the reservoir outside the regular well pattern. An example of such an area is shown as the shaded portion in Fig. 1. The case of a square reservoir utilizing a five-spot injection program and with one full well spacing between the outermost wells and the boundary (Fig. 1) was selected as the single example which would provide the most insight into

Jan 1, 1956

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

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

By J. Wade Watkins, E. S. Mardock

The accurate evaluation of reservoir-performance characteristics in the secondary recovery of petroleum by water flooding requires use of a water tracer that may be injected into water-input wells and detected at oil-production wells to supplement data obtained from core analyses, wellhead tests, and subsurface measurements. Radioactive iodine has been used successfully as a water tracer in field tests to determine: (I) relative rates and patterns of flow of injected water between water-input and oil-production wells and (2) zones of excessive water entry into oil-production wells. Laboratory evaluations of potential water tracers, previous tracer studies, the value of using a radioactive tracer, general field procedures, and the use of surface and subsurface instruments for the detection of the emitted gamma radiation, are summarized. Data from the field tests are presented graphically and discussed in detail. It is concluded that the radioactive-tracer method, using radioactive iodine, may be used successfully to measure either the relative rates and patterns of flow or zones of excessive water entry into wells under conditions of comparatively rapid transit time between wells. INTRODUCTION Extensive use is made of data obtained from core analyses, wellhead tests, fluid characteristics, and subsurface measurements in evaluating the sweep efficiency of water injected into oil sands for the recovery of oil. Theoretical flow rates and patterns may be calculated from those data using radial-flow formulas, if it is assumed that the physical conditions in the productive formation are homogeneous. Unfortunately, homogeneous conditions rarely, if ever, exist in oil-productive formations. The use of a tracer that may be injected into an oil sand and detected quantitatively, or even qualitatively, at offsetting oil-production wells provides basic data that may be used in determining more accurately the subsurface rates and patterns of flow of injected water between wells than is possible by theoretical calculations based on assumed conditions. Consideration of the data obtained by using a water tracer assists in the application of remedial measures to water-input wells, such as plugging of channels, or selective plugging of highly permeable zones, thereby effecting a more uniform flood and a greater ultimate oil recovery. A water tracer should have the following characteristics: (1) low adsorption on solid reservoir material; (2) high solubility in water; (3) low or negligible solubility in crude petroleum; (4) a wide range of solubility of compounds that may be formed by chemical reactions with ions present in reservoir rocks or waters; (5) high detectability in low concentrations by portable apparatus; and (6) general availability at low cost. Additional desirable characteristics are that the tracer should be nonhazardous in nature under normal working conditions and quantitatively detectable in the well bore, as well as at the surface. Review of Previous Tracer Studies The need for a water tracer has been recognized for many years. As far back as 1906, Dole1 described the use of fluorescein in flow-rate studies made in France in 1901. In 1921, Ambrose" discussed in detail the use of organic and inorganic dyes, chlorides, nitrates, and other anions, lithium salts, and the Slichter electrical

Jan 1, 1955

By Fred H. Poettmann, A. L. Benham

correlation of waters and identification of mixtures of two and three waters is important in secondary recovery evaluation studies and studies involving underground movement of waters. The importance of such studies suggests the need for a method to correlate a great number of water analyses to accurately identify and improve the understanding of mixtures of waters. Geochemical analyses of waters furnish a quantitative measure of each of the ions present in a water. Table L shows a portion of the data on many oilfield water analysis report forms in use today. Under "parts per million" are listed the quantitative amounts of each ion in milligrams per liter. Following the work of Stabler' and Palmer", some laboratories have reported additional data to indicate the reaction capacity, or reacting value of each ion. Table I lists data under two subheads for the "reacting value". Under the first subhead, the reacting value is reported in gram equivalents per million, hereinafter called equivalents per million. This value is determined by multiplying milligrams per liter of the ion by a factor equal to the reciprocal of the equivalent weight of the ion. Under the second subhead (per cent) is listed the percentage of each ion based on the total equivalents per million. Sodium is not usually an analyzed value in routine analyses, but is determined by the difference between total negative equivalent per million values and the sum of all of positive ion values for which analyses were made. This difference is reported as sodium, hereafter labelled "sodium". On this basis, the subtotals of the equivalent per million and per cent reacting values of the positive and negative ions are equal. CORRELATION OF WATER ANALYSES DATA Previously reported methods for correlating oilfield waters have used part of the data in Table 1, and can be classified into two categories. Correlation by Chemical Character Correlation methods using the chemical character of waters, as represented by the per cent reacting values of the ionic constituents, have been used by Palmer2 and Tickell. These methods have the limitation of neglecting the ionic concentration of the waters. This limitation can be serious in areas where waters have the same chemical character, but differ in concentration. Correlation by Concentration of the Ionic Constituents The Reistle' and the Stiff methods use the ionic concentration in parts per million and equivalents per million. These methods are effective in differentiating waters of similar chemical character and varying concentration, as in the highly saline waters of the Mid-Continent area. These methods have the disadvantage that they correlate waters by comparing many graphs of individual water analyses to determine the similarity to, or the variance from "type" waters. PROPOSED METHOD OF CORRELATING WATER ANALYSES DATA Correlation graphs presented in Fig. 1 represent water analyses data as points to facilitate the correlation of waters. The graph on the left of Fig. 1 is a plot of per cent reacting value of magnesium vs calcium. The right

By P. Hazebroek, C. S. Mathews, E. E. Allen, H. C. Lefkovits

A rigorous study was made of the behavior of reservoirs composed of horizontal layers, unconnected except at the well and filled with a compressible fluid. The report is presented in two parts. Part I deals with the practical implications of the results obtained; Part 11 contains the mathematical derivation of the solution. The main portion of Part I is concerned with application of the results to reservoirs composed of two layers. Relative rates of depletion of the layers are studied, and it is shown that differential depletion between layers exists in the transient stage of the reservoir, which generally is of an order of magnitude longer than the transient stage in a single-layer reservoir. In the transient stage, the more permeable layer is depleted faster than the less permeable. However, the reservoir approaches a steady-state situation where both layers contribute equally to production. Theoretical build-up curves from such reservoirs are presented, and the influence of a skin effect and flow into casing and tubing are also evaluated. The build-up curves are predicted to have, first, the familiar logarithmic straight-line section and, subsequently, a rising and flattened section. The results found are compared with an extension of an approximate theory, and it is found that this simplified theory is applicable over most of the life of the reservoir. Extensions to reservoirs containing more than two layers are indicated, and examples of the performance and build-up curves for three-layered reservoirs are presented. The theoretical results on build-up curves were applied to field examples of such reservoirs, and satisfactory results were obtained. From the build-up curve, the pemeability-thickness product, the wellbore damage and the static pressure are obtained. However, it apparently is not possible to determine the properties of the individual layers from the combined build-up curve. A large number of reservoirs have been found where the producing formation is composed of two or more layers of differing physical characteristics, such as permeability and porosity, and where also the thicknesses of the layers differ. Since in general a differential depletion should exist between the layers of such reservoirs, it would be advantageous to know the amount of the differential depletion under a given production history. It also would be valuable to know the pressure behavior, because it could be undesirable to produce a reservoir in such a manner that the pressure in one of the layers would have fallen below the bubble point of the oil at a time when another layer might still be virtually undepleted. The determination of static pressures in shut-in wells with slow build-up is another feature of importance in layered reservoirs. If the behavior of the reservoir differs appreciably from the predicted behavior of a single-layer reservoir, the type of analysis for build-up curves suggested for single-layer reservoirs is not applicable. Several authors have treated the problem of finding the pressure and production characteristics of multilayer reservoirs. Horner treated the problem of n layers in an infinite field, the layers being connected only at the well. Tempelaar-Lietz2 has presented an approximate treatment of the characteristics of a bounded reservoir composed of two layers of different permeabilities and equal thicknesses. Homer gave no information on differential depletion between layers, and the treatment of Tempelaar-Lietz was not applicable for prediction of pressure build-up behavior. In addition, it was thought that a more rigorous treatment than the one by Tempelaar-Lietz would be valuable in determining the range of validity of this solution. The present study, therefore, was made of a bounded reservoir composed of two or more layers,

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

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

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

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