Search Documents

By R. G. Jones, W. L. Kinney, R. E. Schilson, R. S. Wilson, G. A. Clark, H. Suralo

This paper describes the field opeturluns at the Fry in situ combustion project in Crawford County, 111. Field operations may be divided into three phases—-preliminary engineering studies of 1960 and 1961, ignition, and, finally, operation of the test. The Fry combustion project was a pilot test of cocur-rent in situ combustion on a 3.3-acre inverted five-spot. In the summer of 1960, air injection tests were conducted using old wells in the area and a new air injection well drilled especially for the test. In the summer of 1961, four new oil wells were drilled on a 3.3-acre pattern for the pilot test. All old wells in the near vicinity of the pilot test were abandoned to preclude any possibility of gas escaping to adjacent formations. A compressor plant was also installed at the test site in 1961. Air injection tests and other engineering studies were made. At the conclusion of these tests, ignition was accomplished with a 40 kw electrical igniter. Operation of the Fry combustion test was little different from the usual oilfield operation with the exception of the high rate of air injection. No difficulty was found with the openhole completion. There were no unusual corrosion or emulsion problems until bottom-hole temperatures exceeded 150F. Four cores were drilled into the burned-out portion of the reservoir with no mechanical difficulties. INTRODUCTION A cocurrent or conventional in situ combustion test was conducted on the Fry unit of the Marathon Oil Co., in Crawford County, Ill. Besides testing combustion in the Fry type of reservoir, the test was intended to evaluate various operating procedures. This paper presents a discussion of the Fry combustion test field operations. The geology of the reservoir and the test performance are discussed in other papers.'.' HISTORY OF THE FRY UNIT The Fry unit consists of 11 leases totaling 337 acres in an area in which oil was originally discovered in 1906. Production from this reservoir, as from all Robinson sand reservoirs in the area, was very high during the first few years. As the dissolved-gas energy was dissipated, early attempts to retard the decline consisted of vacuum, air, and gas repressuring, and recently, most of the area has been subjected to water flooding. Following early air repressuring, the Fry waterflood unit was formed and water flooding initiated in July, 1952. In the area in which the Fry combustion test was conducted, a peripheral water flood of a 160-acre area was developed as shown in Fig. 1. Between 1956 and 1960, it became apparent that the water flood was unsuccessfd. In looking for a suitable site for an in situ combustion test, the Marathon Oil Co. selected the Fry unit as possessing those features most desirable for such a research project. Field operations associated with the Fry combustion test began in the spring of 1960.

Jan 1, 1966

By R. H. McLean

Velocity, kinetic energy and shear in crossflow beneath three-cone jet bits may influence cleaning of the bottom of the borehole and the teeth of the bit. Laboratory investigation shows that each of these parameters is a function of the diameter of the borehole and the product of the volume rate of flow and velocity through the nozzles (QVn). Increasing QV, or decreasing the diameter of the borehole increases each parameter. These functions provide means of predicting the magnitude of each parameter and of scaling the cleaning forces. INTRODUCTION In drilling operations using conventional jet-type rock bits, the impinging jels create an important flow mechanism. Called crossflow, this flow mechanism originates in the impact area of the jets, spreads across the bottom of the hole and supplies the principal source of energy to clean the teeth of the bit and most of the bottom.' Besides providing the means of cleaning cuttings from the bit and the hole bottom, crossflow may also have other, less direct, effects on the rate of penetration. The shear stress generated on the bottom by crossflow will inIluence the thickness and permeability of any filter cake of mud solids or crushed material which forms on the bottom.' These factors may affect the rate of penetration. A previous pulblication introduced some fundamental concepts of crossflow.' Crossllow was shown to occur in a thin layer adjacent to the bottom, and to cover the bottom completely. The maximum velocity in the crossflow above any position on the bottom was found to be directly proporlional to the square root of the product of the volume rate of flow and velocity through the nozzles—the jct QV,,—and inversely proportional to the diameter of the borehole. This information indicated that the effectiveness of crossllow in scavenging the bottom can be improved by maximizing the jet QV. The investigation reported herein amplifies the definition of crossflow. Complete velocity profile data above a representative position on the bottom are analyzed. These data better illustrate control of the capacity of crossflow to scavenge the bottom, and also relate shear stress on the bottom to known, controllable parameters. The conclusion reached in the previous publication that maximization of the jet QV. produces the maximunl cleaning beneath current jet bits is unchanged by these new data; rather, it is strcngthened. Data presented here show that the kinetic energy flux above a representative position on the bottom is maximized by maximizing the jet QV. The shear stress on the bottom will also be shown to be maximized in the same manner. Since the functions relating the jet QV. to velocity, shear stress and kinetic energy also involve the diameter of the borehole, means of equating, or scaling, these quantities in different sizes of boreholes will be illustrated. EXPERIMENTAL EQUIPMENT AND TECHNIQUE JET BIT MODEL Data were recorded from the same laboratory model as used in the aforementioned investigation of the flow around a jet bit.' The model consisted of a 43/4-in. three-cone, jet-type rock bit in a smooth, flat-bottomed borehole constructed of lucite. The bit had a shape and nozzle placement closely resembling larger three-cone jet bits commonly used in field operations. Fig. 1 illustrates the impact area on the bottom of a jet from this bit. Details of the orientation of the jets may be found in the previous publication. TECHNIOU'E OF MEASUREMENT Measurements of the crossflow were made by inserting a very small Pitot tube through the bottom of the simulatetl boreholc. Extreme thinness of the layer of crossflow necessitated accurate measurements of the height of the Pitot tube above the bottom to achieve close definition of the velocity profile. A cathetometer, which could be read to the nearest 0.005 cm, was used to make this measurement.

Jan 1, 1966

By C. D. Stahl, R. Sandrea

This investigation was made to evaluate the comparative effects of the viscosity and the phase behavior of the buffer fluid in the composite solvent displacement of Bradford crude from waterflooded sandstone cores. Buffer slugs of propane, naphtha and other refined hydrocarbons exbihiting increasingly favorable viscosity ratios and decreasing solubility relationships, respectively, with Bradford crude, were used in long Berea sandstone cores. The secondary slug was isopropyl alcohol. The results indicate that higher oil recoveries are obtained for increasingly favorable phase relationships even when these are accompanied by unfavorable viscosity ratios within the range studied. Moreover, when propane is used as a buffer slug with an adverse viscosity ratio of 36, it gives higher oil recoveries than at a similar size slug of an any1 alcohol having a viscosity ratio of 1.09. The investigation was extended to the study of the effect of flooding rate on oil recovery. Residual Bradford crude was displaced from a 6-ft Berea sandstone core at rates varying from 0.3 to 30 ft/ day. The results show that as flooding rates were increased above or decreased below a minimum range of 1 to 2 ft /day, displacement efficiency increased considerably. INTRODUCTION In recent years, the alcohol slug process as a means of tertiary recovery has been the subject of many investigations. Gatlin and Slobod showed the advantages of the combination solvent slugs over the previous single-slug process.1 The combination slugs require less total slug material to obtain similar oil recoveries. The basic combination process involves injecting a primary or buffer solvent which is preferentially oleophilic, followed in turn by a water-driven slug of isopropyl alcohol. This procedure ensures that the buffer slug, which is miscible only with the reservoir oil, displaces and replaces the oil and is in turn displaced, together with the water, by the isopropyl alcohol slug miscible with both fluids. Likewise, for those reservoir oils that show some degree of phase incompatibility with isopropyl alcohol, as does Bradford crude, a proper buffer solvent can artificially make the alcohol process technically feasible. From an economic standpoint, propane or LPG would be the most appropriate primary solvent to be used in the field; however, Taber, et a12 have indicated that the alcohol process becomes highly inefficient in the displacement of light hydrocarbons because of the entrapment of these hydrocarbons during the stabilized bank formation. Moreover, since materials such as propane have such low viscosities, adverse viscosity ratios with most reservoir oils would be penerated, thereby possibly further reducing the efficiency of the buffer solvent. This paper presents the results of an experimental study of the comparative effects of viscosity ratio and phase relationships on the displacement efficiency of Bradford crude by means of dual solvents. The work was extended to include the effect of flooding rate on oil recovery. Rates within the range of field and laboratory values were used. EXPERIMENTAL PROCEDURE APPARATUS The porous media used in this study consisted of two 10-ft and one 6-ft Berea sandstone cores, 3-in. in diameter. These were cast in galvanized steel pipe of 31/2-in. ID and the annulus was filled with Armstrong cement. The 10-ft cores were assembled by connecting two 5-ft units, which were joined together by 8-in. steel flanges fitted to the pipe. Intimate contact between the two sections was obtained by inserting a thin disk of the same sandstone between the core faces. The outlet and inlet ends of the cores were covered with threaded pipe caps. The injected fluid was gravity fed to two variable-rate positive-displacement pumps. The fluids were filtered before entering the core and the effluent was collected in an automatic fraction collector and

Jan 1, 1966

By A. S. Odeh

ABSTRACT A simplified model was employed to develop mathematically equations that describe the unsteady-state behavior of naturally fractured reservoirs. The analysis resulted in an equation of flow of radial symmetry whose solution, for the infinite case, is identical in form and function to that describing the unsteady-state behavior of homogeneous reservoirs. Accepting the assumed model, for all practical purposes one cannot distinguish between fractured and homogeneous reservoirs fram pressure build-up and/or drawdown plots. INTRODUCTION The bulk of reservoir engineering research and techniques has been directed toward homogeneous reservoirs, whose physical characteristics, such as porosity and permeability, are considered, on the average, to be constant. However, many prolific reservoirs, especially in the Middle East, are naturally fractured. These reservoirs consist of two distinct elements, namely fractures and matrix, each of which contains its characteristic porosity and permeability. Because of this, the extension of conventional methods of reservoir engineering analysis to fractured reservoirs without mathematical justification could lead to results of uncertain value. The early reported work on artificially and naturally fractured reservoirs consists mainly of papers by Pollard,l Freeman and Natanson,2 and Samara.3 The most familiar method is that of Pollard. A more recent paper by Warren and Root showed how the Pollard method could lead to erroneous results,4 Warren and Root analyzed a plausible two-dimensional model of fractured reservoirs. They concluded that a Horner-type pressure build-up plot of a well producing from a fractured reservoir may be characterized by two parallel linear segments. These segments form the early and the late portions of the build-up plot and are connected by a transitional curve. In our analysis of pressure build-up and drawdown data obtained on several wells from various fractured reservoirs, two parallel straight lines were not observed. In fact, the build-up and drawdown plots were similar in shape to those obtained on homogeneous reservoirs. Fractured reservoirs, due to their complexity, could be represented by various mathematical models, none of which may be completely descriptive and satisfactory for all systems. This is so because the fractures and matrix blocks can be diverse in pattern, size, and geometry not only between one reservoir and another but also within a single reservoir. Therefore, one mathematical model may lead to a satisfactory solution in one case and fail in another. TO understand the behavior of the pressure buildup and drawdown data that were studied, and to explain the shape of the resulting plots, a fractured reservoir model was employed and analyzed mathematically. The model is based on the following assumptions: 1. The matrix blocks act like sources which feed the fractures with fluid; 2. The net fluid movement toward the wellbore obtains only in the fractures; and 3. The fractures' flow capacity and the degree of fracturing of the reservoir are uniform. By the degree of fracturing is meant the fractures' bulk volume per unit reservoir bulk volume. Assumption 3 does not stipulate that either the fractures or the matrix blocks should possess certain size, uniformity, geometric pattern, spacing, or direction. Moreover, this assumption of uniform flow capacity and degree of fracturing should be taken in the same general sense as one accepts uniform permeability and porosity assumptions in a homogeneous reservoir when deriving the unsteady -state fluid flow equation. Thus, the assumption may not be unreasonable, especially if one considers the evidence obtained from examining samples of fractured outcrops and reservoirs. Such samples show that the matrix usually consists of numerous blocks, all of which are small compared to the reservoir dimensions and well spacings. Therefore, the model could be described to represent a "homogeneously" fractured reservoir. la this paper, the fundamental equation of flow

Jan 1, 1966

By M. M. Mebta, G. W. Dean, W. H. Somerton

With the advent of underground heating operations, interest has developed in the alteration of rock properties by high-temperature treatment. In the present work a number of sandstones were heated to temperatures in the range of 400°C to 800°C under both atmospheric and simulated reservoir pressures. Pertneabilities increased by at least 50 per cent and sonic velocities and breaking strerlgths decreased by an equivalent amount. Differential thermal expansion and other reactions of constituent min-era1 grains are the causes of these alterations. INTRODUCTION In the underground combustion of petroleum reservoirs, temperatures of the order of 600C are reported to have been reached in the combustion zone.' At this temperature rocks are subject to extensive thermal alteration. Temperatures of this magnitude and higher may also occur in subsurface formations when subjected to bottom-hole heating, thermal drilling operations, and underground nuclear explosions. Temperatures of this magnitude might also be generated by conventional rock drilling methods at points of bit-tooth contact. In earlier work, the permanent deformation of rocks resulting from heating was reported. Major structural damage of rocks occurs due to differential thermal expansion of mineral constituents. A number of mineral alterations including crystal inversions, loss of water of crystallization and dissociation, may also contribute to changes in physical structure and properties of rock. In the present work, samples of three typical sandstones were heated to several temperatures up to a maximum of 800C and then allowed to cool to room temperature. Heating was done under both atmospheric pressure and simulated reservoir pressure conditions. Physical properties of the samples were measured before and after heating and comparisons made. Measured properties included permeability, sonic velocity, breaking strength and fracture index. Changes in physical properties were compared to changes in mineralogical characteristics as determined by thin-section. X-ray diffraction and chemical tests. EQUIPMENT AND PROCEDURE Two outcrop sandstones (Bandera and Berea) and one sub-surface sandstone (St. Peter) were selected for the tests. These samples have a wide ranee in composition and physical properties as shown in Table 1. The first series of tests was made on 2-in. diameter by 5-in. long test specimens. Test specimens used in all later work were 3/4-in. diameter by 1 1/8-in. long, this being the specimen size required for heating at simulated reservoir pressures. After careful washing, the cores were oven dried at 100 ± 5C for a minimum of 24 hours before the tests were run. Test specimens were heated in an electric furnace at a constant rate of temperature increase of 3C per minute. When maximum temperature of the run was reached, the sample was allowed to soak for one hour. The furnace was then cooled to room temperature at the same rate of 3C per minute. The entire heating operation was designed for reproducibility without subjecting the test specimens to excessive thermal shock. For samples heated under simulated reservoir pressures, a pressure cell designed by Dean was used (Fig. 1).3 The core sample was inserted into a thin-walled (0.006 to 0.01-in.) copper cup which was then mounted in a high-pressure cell. Provisions were made for the application of internal pore pressure as well as confining pressure. Tests showed that the thin-walled copper cup closed tightly around the core and satisfactorily transmitted confining pressure to the core. The core was heated by placing the entire cell into the electric furnace. The heating program was the same as that used in the atmospheric pressure runs: tempera-ture rise of 3C per minute to maximum temperature of the test, soaking at maximum temperature for one hour, and cooling at a rate of approximately 3C per minute. The cell was designed to withstand 5,000 psi at 1,000C. However, since it was considered likely that repeated heating and cooling would in time weaken the steel, 2,000 psi at 850C was set as a working limit. In the present series of tests, the pore pressure was held constant at 750 psi and the confining pressure at 1.500 psi. The pressure source was a high-pressure nitrogen tank. The two pressures were controlled manually but are accurate well within ± 50 psi.

Jan 1, 1966

By B. A. Eaton, R. H. Jacoby

This paper presents the depletion-performance correlations developed using data from PVT studies and well-completion tests of 27 rich gas-condensate and volatile-oil reservoir fluids. The PVT behavior for each fluid war determined experimentally in the laboratory, and the depletion performance of each fluid was then calculated. Well-completion test data consisting of reservoir pressure and temperature, initial total gas-oil ratio, and stock-tank oil gravity were used with the depletion-performance data to develop general performance correlations. Three correlations are presented. The first equation shows the empirical relation of gas in place per barrel of hydrocarbon pore space with reservoir pressure and temperature, and initial total gas-oil ratio. The second equation is the correlation of stock-tank oil in place per barrel of hydrocarbon pore space. The third equation relates stock-tank oil production (to 500-psia depletion pressure) with pressure, temperature and oil gravity. All three correlations are also presented in convenient graphical form. Gas production to 500-psia depletion pressure is approximately 92.6 per cent of the gas in place for each fluid, so no correlation was developed for gas production. INTRODUCTION The rate of discovery of shallow, low-pressure and low-temperature oil reservoirs in the United States has been declining in recent years. For this reason, the search for oil has gone to greater depths where high pressures and temperatures are encountered. A high-pressure and high-temperature reservoir is the environment in which many rich gas condensates and highly volatile oils are found. Because of this trend, the performance prediction of such reservoirs has become an important task of the reservoir engineer. The need and usefulness of a performance correlation based on field data that are available from the usual well-completion tests have previously been established. These earlier attempts have shown that the phase behavior of rich gas-condensate systems and volatile-oil systems can be correlated with fluid properties, pressure and temperature. However, very little information on gas-condensate depletion performance5 was available at the time of these earlier attempts. The availability of more recent data has made possible the development of an improved correlation. The data, as listed in Tables I and 2, are for 25 gas-condensate and two volatile-oil systems that were studied. The problem was to correlate the data for in-place fluids and the production data with the initial total gas-oil ratio, initial reservoir pressure and temperature, and stock-tank oil gravity. These correlations were obtained by multiple regression analysis, whereby the data were fitted to an equation by the method of least squared error. BASIS OF THE CORRELATIONS 1. The data used in the development of these correlations were obtained from depletion-performance calculations6 and well-test or production measurements. 2. All of the independent variables used in the correlations are available from the usual well-test completion data. 3. The reservoir volume considered is 1 bbl of hydrocarbon pore space (HCPS). 4. The gas-oil ratio used is the total standard cubic feet of gas (at 14.7 psia, 60F) coming off the entire separator train per barrel of stock-tank oil. This minimizes the effect of separator conditions which might vary from system to system. 5. The correlated pressures are the actual reservoir pressures which might or might not be the saturation pressure of the reservoir fluid.

Jan 1, 1966

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By D. E. Hawk, R. F. Burdyn, F. D. Patchen

Earlier work on a mud separator for barite recovery is extended to the design and construct ion of a rugged field unit. Problems associated with scale-up for field use include the me of dilution water and its effect on capacity, power requirements and pressures developed. Equations for these quantities are derived and applied in maximizing the beneficial effects of dilution and minimizing the power required. Tests on the full-size rotor show that its performance can be predicted by theory. The completed unit is compact and lighter than current decanting types. Its capacity is from 10 to 22 gal/min with corresponding barite recovery efficiencies of 92 to 86 per cent. INTRODUC'TION In an earlier publication a new type of centrifugal device was described for solid-liquid separations.l In that paper an operating equation was derived, and its validity demonstrated in laboratory scale experimentation. Some of the features apparent in this new device seem particularly well suited to the conservation of barite in field drilling-mud systems, hence the development was extended to the design and construction of a field unit. In the present paper the design problems associated with scale-up are considered. A full-scale, rugged self-contained field unit is described, and test results discussed which verify that the performance of the large field unit can be described by a simple operating equation. PRINCIPLE OF OPERATION OF SEPARATOR In principle, the new device is simply a perforated cylinder rotating within a body of fluid which is contained within a stationary case, as in Fig. 1. The frictional drag generated by the rotating cylinder causes part of the fluid to rotate. At the surface of the cylinder, assuming no slip occurs, the angular velocity of fluid and of the cylinder are identical. If there are suspended particles present of greater density than the fluid, the centrifugal force acting on these particles causes them to move radially away from the cylinder at some velocity vp . Where the Stokes settling equation holds, this velocity may be expressed as:2 d2 ?p rw2 If the case is provided with openings at both ends, a suspension can be pumped in at one end at a fixed rate. Part of the suspension can then be removed at a fixed rate at the other end of the case as underflow, and the remainder forced through the perforations and hollow shaft of the rotor to appear as effluent, also at a fixed rate. The radial component of the flow velocity of this stream as it reaches the rotor surface, is simply effluent flow rate divided by cylindrical surface area: Because up is a function of particle size, for some particular particle size, the velocities vp and v, will be equal in magnitude. This is defined as the critical particle size and represents the size above which all particles appear in the underflow and below which all particles appear in the effluent. For this equilibrium condition at the rotor surface, Eqs. 1 and 2 combined give

Jan 1, 1966

By H. J. Ramey

The wellbore acts as a storage tank during drawdown and build-up testing and causes the sand-face flow rate to approach the constant surface flow rate as a function of time. This effect is compounded if non-Darcy flow (turbulent flow) exists near a gas wellbore. Non-Darcy flow can be interpreted as a flow-rate dependent skin effect. A method for determining the non-Darcy flow constant using this concept and the usual skin effect equation is described. Field tests of this method have identified several cases where non-Darcy flow was severe enough that gas wells in a fractured region appeared to be moderately damaged. The combination of wellbore storage and non-Darcy flow can result in erroneous estimates of formation flow capacity for short-time gas well tests. Fortunately, the presence of the wellbore storage eflect permits a new analysis which can provide a reasonable estimate of formation flow capacity and the non-Darcy flow constant from a single short-time test. The basis of the Gladfelter, Tracy and Wilsey correction for wellbore storage in pressure build-up was investigated. Results led to extension of the method to drawdown testing. If non-Darcy flow is not important, the method can be used to correct short-time gas well drawdown or build-up data. A method for estimation of the duration of wellbore storage effects was developed. INTRODUCTION In 1953, van Everdingen and Hurst generalized results published in their previous paper3 concerning wellbore storage effects to include a "skin effect", or a region of altered permeability adjacent to the wellbore. Later, Gladfelter. Tracy and Wilsey4 presented a method for correcting observed oilwell pressure build-up data for wellbore storage in the presence of a skin effect. The method depended upon measuring the change in the fluid storage in the wellbore by measuring the rise in liquid level. To the author's knowledge, application of the Gladfelter, Tracy and Wilsey storage correction to gas-well build-up has not been discussed in the literature. It is, however, a rather obvious application. Gas storage in the wellbore is a conlpressibility effect and can be estimated easily from the measured wellbore pressure as a function of time. Several approaches to the wellbore storage problem have been suggested. As summarized by Matthews, it is possible to minimize annulus storage volume by using a packer, and to obtain a near sand-face shut-in by use of down-hole tubing plug devices. Matthews and Perrine have suggested criteiia for determining the time when storage effects become negligible. In 1962, Swift and Kiel' presented a method for determination of the effect of non-Darcy flow (often called turbulent flow) upon gas-well behavior. This paper provided a theoretical basis for peculiar gas-well behavior described previously by Smith. Recently, Carter, Miller and Riley observed disagreement among flow capacity k,,h data determined from gas-well drawdown tests conducted at different flow rates for short periods of time (less than six hours flowing time). In the original preprint of their paper, Carter et al. proposed that the discrepancy in flow capacity was possibly a result of wellbore storage effects. Results of an analytical study of unloading of the wellbore and non-Darcy flow were recorded by carter.14 In the final text of their paper, Carter et al.!' stated that they no longer believed wellbore storage was the reason for discrepancy in their kgh estimates. In view of the preceding, this study was performed to establish the importance of non-Darcy flow and well-bore storage for gas-well testing. In the course of the study. a reinspection of the previous work by van Everdingen' and Hurst' was made, and the basis for the Gladfelter, Tracy and Wilsey' wellbore storage correction was investigated and extended to flow testing. WELLBORE STORAGE THEORY As has been shown by Aronofsky and Jenkins,11-12 Matthews," and others, flow of gas can often be approximated by an equivalent liquid flow system. The following developnlent will use liquid flow nomenclature to simplify the presentation. Application to gas-well cases will be illustrated later. First, we will use the van Everdingen-HursP treatment of wellbore storage in transient flow to establish (1) the duration of wellbore storage effects, and (2) a method to correct flow data for wellbore storage. DURATION OF WELLHORE STORAGE EFFECTS When an oil well is opened to flow. the bottom-hole pressure drops and causes a resulting drop in the liquid level in the annulus. If V. represents the annular volume in cu ft/ft of depth, and p represents the average density of the fluid in the wellbore, the volume of fluid at reservoir conditions produced from the annulus per unit bottom-hole pressure drop is approximately: res bbl-- (V, cu ft/ft) (144 sqin./sq ft) psi -(5.615 cu ft/bbl)(pIb/cuft) ........(I)

Jan 1, 1966

By P. J. Jacquard, C. Jain

Interpreting pressure measurements is one of the fundamental problems in the operation of oil reservoirs. The principal methods and means of study zuhich are currently used are reviewed. Emphasis is placed on the necessity of improving the quality of the information to turn it to the best possible account. One of the important aspects of the problem is the development of a system of automatic interpretation of the pressure measurements. The proposed method aims at this objective: starting from the relation which exists between a local perturbation of permeability and its effect on the velocity potentials in the wells, we determine by successive approximation the modifications to be made to the permeability map by fitting the least squares between the measured pressures and the pressures indicated by the mathematical model. The method has already been successfully employed in the case of heterogeneous reservoirs of radial- circular symmetry. Its application is extended here to the two-dimensional case, and is illustrated by two examples, for which the results are presented. The conditions required for the method to become fully operational are likewise analyzed. I — INTRODUCTION One of the basic problems of production and recovery from oil fields is interpreting pressure measurements made in wells in the single-fluid-phase regime. Simplifying assumptions, for example supposing the reservoir to be homogeneous, make it possible to use more or less simple methods which permit the estimation, from these pressure measurements, of a mean value of the permeability around the wells. Considering the actual heterogeneity of the majority of reservoirs, and the phenomena of interference between wells, these methods are only indicative. Much has been published on the studies of different cases of comparatively simple heterogeneity: some have discussed the radial heterogeneity around a single wel1,l or vertical heterogeneity.2 The results of these works have shown the influence of variations of permeability on the pressure, but unfortunately they are not very useful in practice because they all are concerned with the phenomena around a single well, without interference between wells, and they correspond to relatively simple and theoretical schemes. These studies of reservoirs or parts of reservoirs on a larger scale are made by the electric analyzer or digital model. We define a permeability map and a porosity map, then reproduce the production history of the wells on the wells of the model, and note the corresponding pressures. These pressures are then compared with actuality, and the discrepancies observed are analyzed. From this analysis we deduce the modifications to be made to the maps of permeability and porosity, and repeat until the coincidence between the pressures indicated by the model and the actually observed pressures is considered adequate. We have adjusted the model. Obviously, in this method of investigation, everything depends on the formula which permits the modification of the permeability to approach the adjustment. Considering the complicated nature of the problem, some oversimplified formulas are poorly adapted to reaching a good adjustment. For example, we have shown that the elementary operation of increasing the permeability not far from a producing well to compensate for too low pressures in the well can, in certain cases, express itself by an effect exactly opposite. 3 The project described here is for a system where we begin with the pressure measurements and the several known parameters of the reservoir, and develop the most probable permeability map, taking account of the measurements available and a subdivision of the reservoir fixed a priori. Ideally, this system is conceived as one in which no manual intervention would be necessary. The model simulating the reservoir is only one part of the whole system, and it is neither the longest nor the most complex part. The need to integrate the model into an aggregate of much more extensive processes has led to the adoption of the digital solution, even though a solution combining the digital and the analog is perfectly conceivable. The simulating model is in that case no more

Jan 1, 1966

By K. Leutwyler, H. L. Bigelow

Jan 1, 1966

By D. N. Dietz

A method for determining average reservoir pressure is presented, which is simpler to apply than that devised by Matthews. Brons and Hazebroek. For bounded reservoirs, identical results are obtained if stabilized-flow conditions prevail. The present method yields inferior results in the transient state. The method can, with a slight modification, also be used for water-drive reservoirs. INTRODUCTION In the method proposed by Matthews, Brons and Haze-broek' for determining average reservoir pressure in a multi-well reservoir, the cumulative production (the production time of each well) enters twice: once when the buildup is plotted against In (t + t)/At to arrive at p*, and a second time when the correction, p*-p is determined with one of the several formulas for differently shaped drainage areas. Once a steady state has been attained, however, the previous production history should be immaterial. The same pressure distribution could have been arrived at after different cumulative productions of the individual wells. In principle, therefore, it should be possible to determine average pressures without referring to cumulative productions. In the following paragraphs an expression is presented for the difference between the pressure in a producing well and the average pressure of its drainage area. Then the build-up time needed to overcome this difference is indicated. CASE OF A CIRCULAR DRAINAGE AREA AND A CENTRAL WELL PRESSURE DISTRIBUTION BEFORE SHUT-IN The general differential equation of radial flow (see Muskat, Eq. 10.2) may be written In the steady state, the rate of production of the well is equal to the rate of expansion of the fluid contained in the the drainage area; thus

Jan 1, 1966

By L. J. Health

Synthetic rock with predictable porosity and permeability has been Prepared from mixtures of sand, cement and water. Three series of mixes were investigated primarily for the relation between porosity and permeability for certain grain sizes and proportions. Synthetic rock prepared of 65 per cent large grains, 27 per cent small grains and 8 per cent Portland cement, gave measurable results ranging in porosity from 22.5 to 40 per cent and in permeability from 0.1 darcies to 6 darcies. This variation in porosity and permeability was caused by varying the amount of blending water. Drainage-cycle relative permeability characteristics of the synthetic rock were similar to those of natural reservoir rock. INTRODUCTION The fundamental behavior characteristics of fluids flowing through porous media have been described in the literature. Practical application of these flow characteristics to field conditions is too complicated except where assumptions are overly simplified. The use of dimensionally scaled models to simulate oil reservoirs has been described in the literature. 1-3 These and other papers have presented the theoretical and experimental justification for model design. Others have presented elements of model construction and their operation. In most investigations the porous media have consisted of either u ncon soli da ted sand, glass beads, broken glass or plastic-impregnated granular substances—materials in which the flow behavior is not identical to that in natural reservoir rock. The relative permeability curves for unconsoli-dated sands differ from those for consolidated sandstone. The effect of saturation history on relative permeability measurements is discussed by Geffen, et a1.4 Wyga15 has shown quite conclusively that a process of artificial cementation can be used to render un consolidated packs into synthetic sandstones having properties similar to those of natural rock. Many theoretical and experimental studies have been made in attempts to determine the structure and properties of unconsolidated sand, the most notable being by Naar and Wygal.6 Others have theorized and experimented with the fundamental characteristics of reservoir rocks. This study was conducted to determine if some general relationship could be established between the size of sand grains and the porosity and permeability in consolidated binary packs. This paper presents the results obtained by changing some of the factors which affect the porosity and permeability of synthetically prepared sandstone. In addition, drainage relative permeability curves are presented. EXPERIMENTAL PROCEDURE Mixtures of Portland cement with water and aggregate generally are designed to have certain characteristics, but essentially all are planned to be impervious to water or other liquids. Synthetic sandstone simulating oil reservoir rock, however, must be designed to have a given permeability (sometimes several darcies), a porosity which is primarily the effective porosity but quantitatively similar to natural rock, and other characteristics comparable to reservoir rock, such as wettability, pore geometry, tortuosity, etc. Unconsolidated ternary mixtures of spheres gave both a theoretically computed and an experimentally observed minimum porosity of about 25 per cent.7 By using a particle-distribution system, one-size particle packs had reproducible porosities in the reproducible range of 35 to 37 per cent.6 For model reservoir studies of the prototype system, a synthetic rock having a porosity of 25 per cent or less and a permeability of 2 darcies was required. The rock had to be uniform and competent enough to handle. Synthetic sandstone cores were prepared utilizing the technique developed by Wygal.5 Some slight variations in the procedure were incorporated. The sand was sieved through U. S. Standard sieves. A 18-20 sand denotes that the sand passed the Number

Jan 1, 1966

I would like to congratulate the authors on an interesting paper and, in particular, on their justification of the use of superposition in treating gas wells. They state, however, that no theory has previously been developed lor treating the variable rate drawdown case. In a paper published a few years ago by my colleagues and myself, a method is described which allows a constant rate drawdown curve to be derived for any well. regardless of real rate variations. This method has the great advantage that it is not necessary to make any assumptions regarding the form or uniformity of the reservoir, and the derivation suggested by the present authors is simply an idealized case of our more general treatment. In the light of several years' experience in the derivation of constant rate drawdown curves: I would like to point out that the method proposed by Odeh and Jones may lead to confusing answers, particularly if applied over an appreciable time interval. The reason is that, if pressures are measured when a well is flowing, non-linearities in plots like Figs. 1 and 2 can arise due to either (1) changes in the skin factor with time and/or flow rate, or (2) an inhomogeneous reservoir. Separation of these two effects is only possible if some of the pressures used in the analysis are measured when the well flow rate is 0 or, in other words. during brief shut-in periods. REFERENCES 1. Arthur, K. B.. Metunomski, Z. G. and Ruche, Ms "Intetet ct methode d'utilisa tion des mesures de pression [)our Petude des gisemonts" Revue de T'I.F.P. (Nov., 1962) XVII, 1299. 2. Thomere, R. and Arthur, K. R.: "Analysis of the Pressure Bet~avior of a Production Well Subject to Interterence and Its Application to Reservoir Models", Paper SPE 1189 presented at SPE :Amund Fall Mecting Ireld in Denwer, Colo. (Oct. 3-6, 1965).

Jan 1, 1966

By Paul Soros

Terminal facilities were recently completed at Conneaut, Ohio, on Lake Erie, for the Bessemer and Erie Railroad Co. (a wholly-owned subsidiary of U.S. Steel Corp.) which may become the blue- print for all future coal shipping terminals. Design and engineering of the new facility was by Soros Associates, New York, N. Y. Unlike terminals of conventional design, for example, those at Toledo and Newport News, which load coal directly from the railroad car to the vessel, the new installation passes the coal into temporary storage, loads the vessel from temporary storage and provides alternative reclaiming systems in case of major system failure.

Jan 1, 1966

By A. M. Rowe, I. H. Silberberg

A computer program was written to predict the phase behavior generated by the enriched-gas-drive process. This program is based, in part, on a new concept of convergence pressme, which is then used to select vapor - liquid equilibrium ratios (K-factors) for performing a series of flash calculations. The results of these calculations are the equilibrium vapor and liquid phase compositions which define the phase envelopes. The program was used to predict phase envelopes for 11 different hydrocarbon systems on which published experimental data were available. The predicted and experimental results compare favorably. INTRODUCTION The reservoir engineer is frequently faced with the problem of predicting what will happen if gas is injected into a reservoir. One aspect of this, general problem is predicting the phase changes that will occur when a non-equilibrium gas displaces a reservoir fluid. In particular, a "dry" gas, upon displacing a volatile oil, will pick up intermediate components from the oil. On the other hand, a "wet" gas, containing a high concentration of intermediates, will lose some of these components to a relatively low-gravity, non-eouilibrium1- crude. It is this latter Drocess. occurring in the enriched-gas displacement, which is treated in this paper. In the past, these phase changes have been determined experimentally and the results incorporated into various modifications of the Buckley-Leverett analysis.112 Such experimental work is time consuming, and the results are sensitive to numerous experimental errors. With the wide availability of high-speed digital computing equipment and numerous correlations pertaining to the vapor-liquid equilibria of hydrocarbon systems, it is now practical to calculate such phase behavior. This paper describes a computer program for performing these calculations. THE ENRICHED GAS DISPLACEMENT PROCESS Experimental results have shown that oil recovery can be significantly increased by enriching the displacing gas with intermediate hydrocarbon c0m~onents.3 The essential features of the phase behavior generated by this enriched-gas-drive process are commonly illustrated with ternary diagrams such as Fig. 1.4 In this figure, Gas D, which contains a high concentration of intermediate hydrocarbons with respect to the undersaturated Crude A, is injected into the reservoir. When D contacts A, gas goes into solution until the oil becomes saturated (Point. B). Further contacting of Gas D and saturated Oil B results in a Mixture C which separates into Vapor Y(c) and Liquid X(c). Liquid X(c) is contacted by additional Gas O, resulting in Mixture E which separates into Vapor Y(e) and Liquid X(e). Repeated contacts of the liquid by the injected gas will eventually result in Liquid X(4 of maximum enrichment existing in equilibrium with Gas Y(d). The equilibrium tie-line X(4 Y(4, when extended, passes through the Point D representing the enriched injection gas. For systems of more than three components, the predicted equilibrium states are dependent upon not only reservoir temperature and pressure, but also the compositions of the crude oil and injected gas. If the gas is sufficiently enriched, a miscible displacement is generated. Line If is tangent to the phas,e envelope at the critical point (Point Z) and represents the limiting slope of the tie-lines as the critical state is approached. Point 1 therefore represents the minimum enrichment of injection gas required to generate a miscible displacement. Point G represents the minimum enrichment required for initial miscibility of the injection gas with Crude A. Accra has presented a method to be used for prediction of oil recovery by the enriched gas displacement process.l To develop the phase behavior data needed, he designed the experimental procedure described in the following quotation from his paper: The original liquid was contacted by a volume of displacing gas and allowed to come to equi-

Jan 1, 1966

By T. D. Mueller, H. Dykstra

This paper describes a method for obtaining compositions of gas and oil phases in equilibrium with each other at a given reservoir temperature and pressure as lean or enriched gas is injected into a reservoir oil not in equilibrium with the gas. The method is based on calculation of equilibrium constants, K, from the equation log Kp = A + B F, where p is pressure and F is component characterization factor. A and B are constants that vary with a composition parameter. The paper explains how A and B could be determined with a minimum amount of data. The method consists of calculating equilibrium constants as a function of liquid composition and using the K's to calculate changes in compositiort of injection gas as it successively contacts original resetvoir oil at the gas-oil front and the change in composition of reservoir oil as it is successively contacted by the injection gas in the vicinity of an injection well. An iterative process is required. Two examples of the method are given which show how oil and gas compositions and properties are calculated with the aid of a computer. INTRODUCTION The calculation of recovery from lean-gas or enriched-gas injection into a saturated or under-saturated reservoir oil requires a knowledge of how phase properties change in a reservoir between injection well and producing well. To determine the phase compositions or properties experimentally would be difficult and time consuming, especially if one wanted to look at the effect of various injection gas compositions. For some fields an engineer would want to make preliminary calculations on the economics of gas injection before he makes the decision to spend a corrsiderable sum of money for laboratory work to obtain the necessary fluid properties. A method for calculatine compositions and properties based on a minimum of laboratory data is described. The minimum amount of data required is injection gas and reservoir oil composition and injection gas solubility as a function of pressure for an undersaturated oil, or the compositions of mixtures of reservoir fluid with two or more enriched gases, which mixtures have a bubble-point pressure that is the same as the expected average displacement pressure. On the other hand, if the economics appear to be attractive, more laboratory data can be obtained in the form of equilibrium constants for three or more mixtures of the reservoir fluid and injection gases. The method is based on the calculation of K's from the equation, log Kp = A + B F, where K is equilibrium constant, p is pressure and F is component characterization factor. This method of plotting K data by means of the above equation is described by Hoffman, Crump and Hocott.1 A and B are constants that vary with a composition parameter. It is shown how experimental data are used to determine how A and B vary with composition. BASIS OF METHOD EQUILIBRIUM CONSTANTS Equilibrium constant (K) data are required for this method. The K of each component is obtained from an equation that gives the product of K and pressure as a function of component characterization factor. The characterization factor F = b(l/Tb — l/T), where b is a constant characteristic of a given component, Tb is the boiling point in OR at atmospheric pressure of the given component, and T is the temperature of the system in OR. This method of plotting K data is described by Hoffman, Crump and Hocott.l The straight-line equation of the log Kp - F plot is log hp = log Kop +BFi, where Kop is at F =0, i designates the component, and B is the slope of the line. The factor F can be

Jan 1, 1966

By K. T. Koone, R. J. Blackwell

One of the more important problems to be solved in designing a miscible flood is related to the size of the solvent bank used. Size of the bank may be critical to economic success. Too large a bank loses money; too small a bank may deteriorate and fail to maintain the miscibility needed for high recovery. An important factor in deterioration of a small bank is permeability channeling. In a highly stratified reservoir, solvent speeds ahead in the more permeable zones and mixes laterally with fluids bypassed in adjacent, low-permeability strata. Numerical solutions have been obtained .for the differential equations that describe the movement of a slug through a two-layer system in which mixing occurs both in the direction of flow and transversely. The solvent slug is assumed to have the same density and viscosity as the resident fluid and the pushing fluid. These solutions have been verified by comparing them with similar concentration profiles obtained in the laboratory in a 36-ft stratified model packed with glass beads. The theoretical study revealed that when the dominant mechanism causing a bank to fail is lateral mixing the bank size needed for a given recovery may increase with length rather than decreasing as the square root of reservoir length, as suggested by one-dimensional mixing theory. From a comprehensive examination of the variables, a generalized correlation is developed that relates strata thicknesses, bank size, fluid velocity, mixing coefficients, system length and simple solvent-resident fluid phase behavior to the area miscibly swept. INTRODUCTION Miscible displacement, or solvent flooding, continues to receive widespread attention as a method for increasing oil recovery over that possible in conventional gas-drive or waterdrive projects. A basic economic requirement in the application of such processes is the use of as little solvent as possible. A basic physical requirement is that enough solvent be used to maintain miscibility. Economics places an upper limit on the size of a solvent slug, and physical considerations establish a lower limit. Consequently, the practicality of any given miscible process requires that the economic limit be greater than the lower limit imposed by physical requirements. Procedures exist for determining the economic limit; however, procedures for determining realistic minimum bank sizes exist for only special reservoir situations. In the past, bank size has usually been selected on the assumption of a piston-like displacement for which only longitudinal mixing is important. This assumption leads to the favorable conclusion that bank size, expressed as per cent of pore volume, is inversely proportional to the square root of length. Collins considered the problem of transverse mixing of solvent with fluids in bypassed zones with the assumptions that no forward mixing occurs and that the concentration is uniform in the permeable stratum of interest. 1 Lauwerier considered a mathematically similar problem in thermal recovery operations.2 Their work suggests the much less favorable conclusion that bank size could be directly proportional to the length or even higher powers of the length. This paper considers the physical behavior of a small solvent bank as it moves through a real reservoir, without imposing many of the restrictive assumptions of past treatments. To facilitate mathematical description, an idealized, two-layer model that permits mixing both laterally and in the direction of flow will be considered. A calculation procedure for solving the descriptive equations will be developed for a displacement involving fluids of equal density and viscosity. In addition, laboratory experiments designed to check the computational results will be described, and the coincidence of calculated and observed results will be discussed. The use of many solutions for a variety of bank sizes, strata thicknesses and characteristic system lengths to develop a usable correlation between

Jan 1, 1966

By C. R. Knopp

Gas-oil relative permeability ratio is an important relationship in oil reservoir predictive calculations. A correlation has been developed from 107 gas-flood k/k tests on Venezuelan core samples. The correlating parameter is based on restored-state water saturation tests and' is applicable to both consolidated and poorly consolidated sandstone reservoirs. Data of the correlation show that there are no distinguishah1e differences between the mass-data groupings for the two c1assifications A procedure is recommended for running .sufficient relative. permeability analyses to compute a geometric mean of the sample group. The geometric mean is more representative of the total core, and probably the entire reservoir. For example, while only one in four of the k,,/k,., test curves agreed closely with the resultant correlation of this report, the geometric mean curves of the 16 suites (three samples or more). showed good agreetment ill three cases out of four. INTRODUCTION The gas-oil relative permeability ratio is an important, fundamental relationship in most oil reservoir predictive calculations. Predictive calculations are made to estimate future reservoir production characteristics and ultimate oil recovery. The k1,/k2, relationship is specifically needed to relate the surface gas-oil ratio to the reservoir oil and gas saturation, and to calculate the relative movement of these phases within the reservoir whenever some of the more complex driving mechanisms are present. Laboratory k1/k2, tests are not generally run as a routine analysis. Consequently, k1/k2 data often are not available when needed because the cost of laboratory work could not be justified or the need for such data had not been properly anticipated. When laboratory k1/k2, data are available, they are often very difficult to interpret. For example, wide divergence is sometimes shown in a family of k1,/k1, tests representative of the producing horizon in a single well. With these considerations in mind, a study was made to determine if a relationship might exist between the k1,/k2, curve and some other simple laboratory test criteria. The most probable k1/k2, curve correlation for Venezuela described in this paper is the result of the investigation. The presented correlation defines the most probable gas-flood k,,/k,, curve through the medium of air-water capillary displacement and centrifuge water saturation tests. The laboratory procedures of these tests are. relatively simple, and inexpensive; test data should be. widely available- from routine analysis. DATA AVAILABLE, LABORATORY METHODS The report correlation utilized 107 gas-Hood k1/k2, tests run on sandstone cores of Venezuelan reservoirs. Table 1 is a general tabulation of data pertinent to the tests, while Table 2 summarizes the data. Thetests include 96 from Western Venezuela and 11 from Eastern Venezuela. Eighty-two- of the 107 test samples were sandstones that varied from poorly consolidated to-unconsolidated; 25 were consolidated. The average sample porosity was 26.7 per cent and the average permeability was 1,121 md; these values typify the better sandstone reservoirs of' Venezuela. The Welge gas-flood technique,' based on fundamental Buckley-Leverett frontal displacement theory, was introduced in about 1952 and is widely accepted in the industry. The laboratory procedure is relatively simple, rapid, and can be performed on small core samples. While there have been some minor variations in sample preparation and laboratory procedure in the tests used for the correlation, these tests can be generally summarized as follows. The core sample was first sol vent-extracted and dried. Connate-water saturation was restored by the oil-flushing or evaporation-blow down methods. At the beginning of gas flood the hydrocarbon pore volume was completeiy saturated with the test oil phase. Unsteady-state gas-oil displacement then began with the injection of nitrogen or helium. while the displaced oil and gas phases were incrementally metered at the out-flow face. From the test data, the k,,/k,, curve was calculated by the Welge method.' The individual oil and gas relative permeabilities were also calculated." CORRELATING PROCEDURES In attempting to establish a basis of correlation, we found that broad mid-range sections of 105 of the 107 k,,/k,, test curves could be closely duplicated by a straight line. Only two curves did not show a degree of linearity in this region. Correlation-curve definition parameters were subsequently developed from this observation of consistent mid-range linearity. Possible correlating variables were limited to the physical properties measured on core samples that (1) were widely available as common test data and (2) could be easily and cheaply obtained through future laboratory work. The more obvious possibilities were porosity, permeability and

Jan 1, 1966

Jan 1, 1966