Search Documents

By D. L. Archer, W. W. Owens

Conventional waterflooding often is uneconomic in highly fractured reservoirs because of the gross bypassing of the reservoir oil by injected water. Imbibition and pressure pulse flooding have been used in at least one fractured reservoir in an attempt to achieve better oil recovery performance. This paper presents the results of laboratory flow tests conducted on large cores to evaluate the possible applicability of these methods (particularly pressure pulse flooding) to different types of reservoir systems. Test data were obtained on both water-wet and oil-wet systems, and on systems having two widely different levels of compressibility and flow capacity. Results indicate that pressure pulse recoveries from fractured reservoirs will likely not exceed 5 to 10 per cent pore space with maximum response achieved during the first pulse cycle. Improved recovery by this method is possible from both oil-wet and water-wet reservoirs. Comparable saturation distributions during imbibition and pressure pulse production suggest that an initial pressure pulse cycle to speed production response would not interfere with subsequent imbibition flooding in water-wet reservoirs. INTRODUCTION There is a steady increase in the number of fluid injection projects being initiated each year to improve oil and gas recovery over that obtained by primary production mechanisms. Experience gained in different geographical areas and with different recovery methods is teaching that reservoir anatomy is one of the more important factors controlling the success or failure of such projects. Fractured formations appear to be the rule rather than the exception, especially in carbonate reservoirs. Conventional methods of fluid injection, whether it is waterflooding, gas injection or miscible flooding, have limited applicability to highly fractured reservoirs because of the severe bypassing of reservoir fluids. The result is early breakthrough of the injected fluid and rapidly increasing ratios of injected to in-place fluid with an undesirable effect on return on investment. Thus, innovations to conventional methods must be developed if recovery from fractured reservoirs is to be optimized. One new method proposed for application to highly fractured reservoirs is pressure pulse waterflooding. Although pressure pulsing with gas has been suggested' and rested in at least one reservoir; the Sohio Oil Co. is generally credited with the development and testing of pressure pulsing in conjunction with waterflooding. Publications" of Sohio's waterflood operations in the Spra-berry Trend indicate that the idea for pressure pulse flooding evolved from a critical analysis of the disappointing early performance of their flood initiated in a portion of the trend in April, 1961. This flood was planned to take advtantage of imbibiltion and high pressure gradients across the reservoir matrix blocks (unfractured Mocks of the reservoir rock which are visualized to be generally surrounded 'by the fracture system) evaluated in pilot tests by the Atlantic Richfield Co.5 and Humble Oil & Refining Co.' However, the early recovery and pressure performance of this flood indicalted that high pressure gradients induced between the fracture system and the centers of the matrix blocks were forcing water into the periphery of the blocks and temporarily interfering with the counter-current flow of oil due to imbibition. Fill-up in the blocks was occurring with re-solution of the free gas phase as pressures climbed above the original reservoir bubble point. It appeared that cessation of water injection to permit the capillary forces to become dominant and that expansion of the rock and its contained fluids during pressure reduction might aid in expulsion of oil from the rock matrix into the fractures. Subsequent production performance. aflter cessation of injection, proved this hypotheslis correct. Thus, pressure pulse waterflooding was born. Pressure pulse flooding appears to have several advan'tages over imbibition type flooding in highly fractured reservoirs. First, all wells may be used for water injection which should hasten fill-up and achieve a more rapid increase in reservoir pressure. Second, because of increased injection pressures and rates, flooding gradients will force water into the reservoir matrix more rapidly than would be achieved by imbibition alone. Third, during the pressure depletion cycle of the flood, the compressibility of the system, which has been increased by resolution of the free gas phase, provides energy for the displacement of oil at a rate greater than would correspond to countercurrent flow during imbibition. Fourth, during depletion all wells may be put on production, thus contributing to higher oil withdrawal rates on a reservoir-wide basis. One possible disadvantage of pressure pulse flooding as compared to imbibition floodling, however, is that the outer periphery of each matrix block is flooded to a residual or near residual oil saturation during the Water lnjection stage. This zone of reduced permeability to oil m'ay offer greater restriction to the production of

Jan 1, 1967

By W. B. Braden

This paper presents a suitable method for predicting gas-free oil viscosities at temperatures up to 500F knowing only the API gravity of the oil at 60F and the viscosity of the oil measured at any relatively low temperature. The API pravity and the one viscosity value are used as parameters to determine the slope of a straight line on the ASTM slanaord viscosity-temperature chart. Then, knowing the slope of the line and one point on the line, the vrscosities at higher temperatures can be determined. The line slope correlations were developed at I00 and 210F since viscosity data are frequently measured at these temperatures. A procedure is given for predicting line slopes from measurements at other tetnperatures. A nomogram is furnished for solving the relationship. The correlation has been evaluated at temperatures up to 5OOF for oils varyzng in gravity from 10 to 33 " API. The correiution is applicable only to Newtonian fluids. Comparison at 500F of true viscosities and those predicted from values at 100F shows an average deviation of 3.0 per cent (maximum deviation of 6.0 per cent). Predictions from the values at 21 0F for the same oils how an average deviation of 1.5 per cent (maximum deviation of 3.4 per cent). INTRODUCTION Correlations have been developed by Beal' and by Chew and Connally' for predicting viscosities of gas-saturated oils at reservoir conditions. Each of these correlations requires a knowledge of the solution gas-oil ratio and the viscosity of the gas-free oil at the reservoir temperature. For temperatures below 350F, measurements of the gas-free oil viscosities can be made easily using commercially available equipment. In thermal recovery processes, however, reservoir temperatures well in excess of 350F are encountered. Viscosity measurements at such conditions are more difficult and time consuming and require modification of existing equipment or the construction of new equipment. Measurements are further complicated by the difficulty of handling highly viscous oils associated with thermal recovery processes. Therefore, it is desirable to have a correlation which allows accurate prediction of viscosities at high temperatures. A commonly used technique for predicting viscosities at high temperatures is to measure the viscosities at two lower temperatures, plot the values on ASTM standard viscosity-temperature charts and extrapolate to the temperatures desired. If either of the values is slightly in error, the extrapolated value can be significantly in error. To justify an extrapolation, three points are actually necessary. This procedure can consume much time, particularly with heavy oils. Considering the cost of viscosity measurements, it would be desirable to eliminate the need for direct measurements by having correlations which would allow viscosity predictions from other physical or chemical properties. Beal1 investigated the possibility of correlating viscosity with oil gravity at temperatures from 100 to 220F. While showing that a general relationship exists, he also found significant deviations. It is possible that correlations will be developed based on oil composition as more information becomes available. While not eliminating the need for viscosity rneasurements, the method presented herein requires that only one viscosity measurement be made. The API gravity must also be known. The theory is based on the fact that the viscosity of paraffins (high gravity) changes less with temperature than does the viscosity of naph-thenes or aromatics (low gravity). The gravity. therefore, is used as a parameter to determine the slope of a straight line on the ASTM standard viscosity-temperature charts. The correlation is applicable only to Newtonian oils, and deviations due to thermal decomposition and nonhomo-geneity cannot be predicted. Oils containing additives have not been evaluated. PROCEDURE Fifteen oils were used in developing the correlation; eight were crudes and seven were processed oils. Oil gravities varied from 9.9" API (naphthene base) to 32.7' API (paraffin base). The temperature range studied was 81 to 516F. Each oil used had a minimum of three viscosity measurements and each plotted essentially as a straight line on the ASTM charts. In all, 91 viscosity measurements were used in the correlation. Saybolt, rolling ball and capillary tube viscometers were used for the measurements. Viscosity data for Samples 1, 2, 4, 7, 10, 11 and 14 were obtained in Texaco, Inc. laboratories. The data for Samples 3, 5, 6, 8, 9, 12 and 15 were from Fortsch and Wilson,3 and data for Sample 13 were from Dean and Lane.' All data points used in the correlation are plotted in Fig. 1. It is seen that some of the viscosity values deviated slightly from the straight-line plots at the higher temperatures. Properties of the oils after exposure to the

Jan 1, 1967

By L. T. Stanley

A method is described whereby the behavior of the front is approximated when gas is injected into a thick reservoir having reasonably homogeneous properties. The method is applied to the Abqaiq field of Saudi Arabia. The producing member is the Arab-D, a high-permeability calcarenite which has no continuous betlding planes. Field performance indicates the injected gas to be partially miscible with the reservoir oil because gas volume factors experienced are lower than those predicted by the gas law. Eqlrilibril~m calculations confirnz ihis and also show that the volume factor may vary appreciably, depending on the relative quantities of oil and gas in contact in the reservoir; consequently, volutnetric balances are made for the purpose of determining the behavior of the gas volume factor during injection history. The effect of gravity segregation on position and shape of the gas front is approximated by dividing the producing section into three zones vertically and tracing the movement of the average interface in each zone, with counterflow of oil and gas taking place betweell zones behind the front. This results in a gas front that advances more rapidly along the top of [he formation than along the base, producing an "umbrella" effect in cross section that becomes more pronounced as the front progresses. The calculations cover the gas-injection period of Abqaiq field history through 1958, and the position of the gas front is plotted at one-year intervals. Finally, some comparisotls between the calcrlated fronial position and the gas-oil contact as Measured by neutron logs are made. INTRODUCTION Fig. 1 is a structure contour map of the Abqaiq field. For analytical purpcses, the field has been divided into two areas designated A and B, as shown on the map. Gas injection into Area A was begun in early 1954 into two wells located near the crest of the structure on the north-south axis of the anticline. During the period covered by this paper, the average injection rate has been 122 MMscf/D with a maximum sustained rate of 160 MMscf/D. The reservoir pressure has been maintained at 2,450 psig since start of gas injection. The reservoir originally was undersaturated but, prior to pressure maintenance, the pressure was drawn below the bubble point in most of the Area A section, creating an in-place free-gas saturation of 8 per cent at the crest of the structure. The average critical gas saturation is 13 per cent established by laboratcry core tests. That the entire reservoir remained below critical gas saturation prior to gas injection is borne out by the fact that none of the crestal wells have exhibited high gas-oil ratios during production history. The Area A reservoir is a carbonate section about 200-ft thick having average porosity and permeability of the order of 22 per cent and 500 md, respectively. The productive horizon is the Arab-D member of Jurassic age and is encountered at an average depth of

By H. G. Warren, E. W. Hough

Correlation of interfacial tension of the methane-n-pentane and methane-n-decane systems was made by Hough and Stegemeier by the use of the Weinaug and Katz equation. The methane-n-heptane and ethylene-n-heptane systems were investigated by the authors. Additional systems taken from the literature were the me thane-propane system by Weinaug. and Katz, the methane-n-butane system by Pennington and Hough, and the n- butane -carbon dioxide system by Brauer and Hough. These systems were analyzed by the Weinaug-Katz equation and the Sugden equation on the IBM 1620 11 computer using a regression analysis program. Parachor and exponent values for each component and system were determined. The experimental work was close to the critical point of the light component in the n-butane-carbon dioxide and ethylene-n-heptane systems, between the critical points of the light and heavy components in the methane-n-pentane, methane-n-heptane and methane-n-decane systems, and close to the critical point of the heavy component in the methane-pmpane and methane-n-butane systems. An empirical formula was developed to find a value of the constant B in the Sugden equation for binary hydrocarbon systems from the parachor values of the two components. Interfacial tension for binary hydrocarbon systems was found to be a direct function of liquid-vapor phase density difference in the same manner as a single-component system. DISCUSSION Interfacial tension may be defined as the measure of the specific free surface energy between two phases having different compositions.1 The interfacial tension (IFT) for a binary system may be predicted for any temperature and pressure if the

Jan 1, 1967

By A. M. Skov, L. F. Elkins

First response to large-scale water flooding in the fractured very low permeability Spraberry sand has led to a new unique cyclic operation. Capacity water injection is used to restore reservoir pressure. This is followed by marly months production without water irzjection and the cycle repeated. Expansion of the oil, rock and water during pressure decline expels part of the fluids but capillary forces hold much of the injected water in the rock. At least with reservoir pressure restored and with partial water flood development, field performance has proved this cyclic operation is capable of producing oil from the nzatrix rock at least 50 per cent faster and with lower water percentage than is imbibition of water at stable reservoir pressure. INTRODUCTION The Spraberry Field of West Texas presents unusual problems for both primary production and water flooding. Extensive interconnected vertical fractures in the fractional-md sandstone permitted recovery of oil on 160-acre well spacing, but they made capillary end effects dominant. Primary recovery by solution gas drive is less than 10 per cent of oil in place. The concept of displacement of oil from the sand matrix by capillary irnbibition of water has led to field techniques which promise greatly increased oil recovery. Free exchange of laboratory research, reservoir information and results of field pilot tests among the various companies has been very important in development of this technology. Five units covering a total of 170,000 acres have been formed for water flooding, and 10 other areas covering an additional 175,500 acres are in various stages of unitization. Part of the Driver Unit reaching fillup first has demonstrated very unusual waterflood behavior and indicated numerous operating problems that will develop within and among the various units. SPRABERRY ROCK AND PRIMARY PERFORMANCE The Spraberry, discovered in February, 1949, is a 1,000-ft section of sandstones, shales and limestones with two main oil productive members: a 10-15 ft sand near the top and a 10-15 ft sand near the base. In part of the field some thinner intermediate sands are oil productive, and others are water bearing. All sands have permeabilities of 1 md or less and porosities of 8-15 per cent. Ordinary core analysis and electric and radiation logs are ineffective in differentiating between oil productive and nonprcductive sands. Sands capable of containing producible oil are best identified by mercury injection capillary pressure measurement and, in some cases, by core water saturation. About 3,500 wells have been drilled in the 500,000-acre trend. Vertical fractures were observed in practically all Spraberry cores. Continuity and interconnection of fractures were confirmed by pressure interference among wells during early development.' Major fractures trend northeast-southwest as indicated by oriented cores and confirmed by five fluid injection tests, by analysis of the pressure transients observed during development,''' and by three interference tests in the Driver Unit Water Flood reported herein. Fracture spac- ing probably averages inches to a few feet. Spraberry wells typically produced 100-400 BOPD initially after hydrauLic fracture treatments. By 1962 oil production had declined to an average of 12 bbl/well/day, near the economic Limits of operation. Reservoir pressure had declined from 2,300 psi initially in the Upper Spraberry and 2,500 psi in the Lower Spraberry to 500-1,000 psi. Partial closing of the fractures with declining reservoir pressure is believed to be the cause of such low oil production rates at these relatively high reservoir pressures. Cumulative recovery of 208 million bbl of oil is 80 to 90 per cent of that recoverable by primary means. Performance of the entire reservoir is summarized in Fig. 1. IMBIBITION WATER FLOODING By 1952 reservoir performance indicated low primary recoveries. Most engineers, expecting serious channeling of injected fluids through the fractures, held little hope for secondary recovery. With its extensive background of research on the fundamentals of fluid flow within reservoir rocks, Atlantic's Research and Development Division on short notice in 1952 conceived that displacement of oil by capillary imbibition of water into the rock might significantly increase Spraberry recovery. Laboratory data reported by Brownscombe and Dyes scaled to probable reservoir conditions showed potential waterflood recovery equal to or greater than primary recovery with a 10-15 year flood life.= A pilot test using three 40-acre injection wells, one central producing well and 18 surrounding observation wells demonstrated technical feasibility of the process. Injection of 1.5 million bbl of water from November 1952 to August 1955 proved water entered the rock and displaced oil

By D. R. Wieland, H. T. Kennedy

The object of the work reported here was to determine the content of elemental sulfur in gaseous methane, carbon dioxide, hydrogen sulfide, and in mixtures of these gases, at pressures and temperatwes encountered in natural gas reservoirs. Sulfur content at equilibrium is reported for pure methane, carbon dioxide, hydrogen sulfide, and on three binary mixtures of each of the three pairs of gases at pressures of 1,000, 2,000, 3,000, 4,000, 5,000 and 6,000 psia and at temperatures of 150°, 200° and 250°F. In addition, the sulfur content of three ternary mixtures at the same temperatures and pressures are reported. The results indicate that the sulfur content is higher in the gases at higher temperatures and pressures. The content is highest in hydrogen sulfide, intermediate in carbon dioxide and lowest in methane. INTRODUCTION At ordinary pressures and temperatures, the concentration of a nonvolatile material, such as sulfur in a gas at equilibrium, is a function of the vapor pressure of the material and is independent of the nature of the gas. As the pressure and temperature increase, the gas assumes some of the properties of liquids, including the power to dissolve other liquids and solids, to an extent dependent on the nature of both the gas and the material dissolved. At equilibrium, the content of sulfur in the several gases considered here may thus be considered as solubilities which are fixed for a given composition of gas, temperature and pressure. The study of the solubility of elemental sulfur in gases is of interest because sulfur is sometimes present in reservoirs producing natural gas and must be present in the vapor phase. Upon reduction of pressure and temperature, the sulfur precipitates from solution in the reservoir and in the tubing and fittings. Due to the fact that the greatest pressure drop in the reservoir is around the wellbore, the volume of sulfur precipitated will be greatest in this locality and can cause a substantial reduction in the permeability of the formation in this area. As the natural gas flows up the tubing string, the pressure and temperature are further decreased, causing further sulfur precipitation. If the volume of free sulfur is large and remedial measures are not taken, complete plugging of the tubing can occur. Knowledge of the content of sulfur as a function of composition, temperature and pressure will enable the operator to determine what changes in pressure and/or temperature are necessary to keep the sulfur in solution. If the sulfur solubility is large and can be controlled, the operator may desire to produce the well at a high wellhead pressure and temperature and reclaim the sulfur at the surface. The high solubility of sulfur in hydrogen sulfide (5.6 per cent by weight at 5,000 psi and 200°F) suggests the feasibility of recovering sulfur by injecting this gas into sulfur-bearing reservoirs and expanding the produced gas to precipitate out the sulfur. The present project was designed to measure the solubility of sulfur in carbon dioxide, methane, hydrogen sulfide and mixtures of these three gases at various temperatures and pressures. These gases were the major constituents of the gas well where plugging by sulfur was encountered. A study of the published literature shows no data on the solubility of sulfur in any gas, although Hannay and Hogarth,' in 1880 reported that it was soluble in carbon disdlide above the critical temperature of carbon disulfide. EQUIPMENT AND PROCEDURE Fig. 1 shows the layout of equipment employed to bring gases to equilibrium with sulfur and to isolate the quantity of sulfur contained in a known volume of gas. Pure gases were measured in the charging bomb at known temperatures and pressures and displaced into the reservoir bomb to make up the desired mixtures. When equilibrium was achieved

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved

By D. D. Dunlop, J. R. Welker

The growing interest in the use of CO, in crude oil recovery increases the need for data on the effect of CO, on hydrocarbon physical properties. Data are presented on the solubility of CO, in various dead oils, the swelling changes in CO2-oil solutions and the effect of CO, on dead oil viscosity. This latter property shows the most pronounced effect, with viscosity reductions up to 98 per cent of the uncarbonated viscosities. An empirical method of estimating the viscosity of carbonated oils is presented. The apparatus and procedures used are described in sufficient detail to allow others to make similar studies. INTRODUCTION The effect of dissolved carbon dioxide on the swelling and viscosity reduction of specific hydrocarbon oils has been observed and recorded by a number of investigators.'- me object of this paper is to offer a means of predicting these effects for crude oils free from natural gas, using the dead state viscosity and gravity of the crude oils. The CO, solubility and swelling of numerous crude oils were determined in a visual cell at various pressure levels. The viscosity of the oils carbonated to various pressure levels was then determined by measuring the pressure drop across a capillary tube. From these data, the physical properties were correlated empirically. The resulting correlations allow the prediction of CO, solubility, swelling and viscosity reduction if the dead state gravity and viscosity of the oils are known. SOLUBILITY AND SWELLING MEASUREMENT EQUIPMENT AND PROCEDURE A high pressure visual cell was installed in a constant temperature cabinet. A test gauge was attached at the top of the cell for pressure measurement, and a line was run through the cabinet wall to a wet test meter which was used for volumetric measurement of the gas. The first step in making a test run was to put the oil in the cell up to a level about half to two-thirds of the total volume. This required about 50 to 65 ml of oil. carbon dioxide was then bubbled up through the oil for a time during which the pressure of CO2 in the cell was kept above 800 psia. Saturation of the oil with CO2 at this pressure and ambient temperature was confirmed by slowly bleeding CO2 through a valve to the atmosphere. If the oil was completely saturated with CO2, bubbles of gas would form in the oil at the first small decrease in pressure. If the oil was under-saturated, no bubbles formed until the pressure was decreased to the saturation pressure existing in the oil. If this saturation pressure was lower than that desired, more CO2 was bubbled through the oil until the desired level was reached. After saturation at ambient temperature was completed, the cabinet temperature was adjusted to the desired level and the cell was allowed to reach temperature equilibrium. After temperature equilibrium was reached, the pressure was again decreased slightly, and the oil again checked for full CO2 saturation at the cell pressure. The pressure now had changed because of the difference in solubility of the CO, in the oil at higher temperatures and the expansion of CO2 as the temperature increased. The outlet tube from the cell was then connected to the wet test meter and the CO2 was allowed to flow slowly out of the cell and through the wet test meter at ambient temperature and pressure. The water in the wet test meter had previously been saturated with CO2 at ambient temperature and pressure by allowing CO2 to flow continuously through it lor a period of several hours. The gas flow was stopped at several pressures during the run and the cell was allowed to come to equilibrium; this made possible the measurement of solubility and swelling data at the intermediate pressures. The volume of the oil in the cell was recorded at each of the equilibrium pressures in order to obtain swelling data. DATA AND RESULTS The solubility of CO2 in the oil was calculated by the relationship V — V, where R. = solubility of CO, in crude oil, cubic feet of CO, measured at 60F and 1.0 atm/ bbl of dead state oil at the temperature under which solubility was measured, V, = volume of gas released from the cell between the saturation pressure and zero pressure, corrected to 60F and 1.0 atm, cu ft, V, = volume of CO: contained in the gas space above the oil, corrected to 60F and 1.0 atm, cu ft, and V, = volume of the dead oil in the cell in bbl at the temperature of the run. The volumetric data of Sage and Lacey' were used to calculate V., from the volume of CO2 at high pressures. The swelling factor was calculated as where V, is the volume of the C0,-

By A. J. Conelius, D. P. Sobocinski

This paper presents a correlation for predicting the behavior of a water cone as it builds from the static water-oil contact to breakthrough conditions. The correlation is partly empirical and involves dimensionless groups of reservoir and fluid properties and of production and well characteristics. The groups were deduced from the scaling criteria for the immiscible displacement of oil by water. The correlation is based on a limited amount of experimenfal data from a laboratory sand-packed model and on results from a computer program for a two-dimensional, incompressible system. Because the correlating groups are dimensionless, they can be used to estimate the performance of water coning cases not specifically considered in the correlation. However, despite its dimensionless nature, the correlation is not completely general and will not provide meaningful estimates of cone behavior in many situations. INTRODUCTION AND BACKGROUND The production of water from oil wells is a common occurrence which increases the cost of producing operations and may reduce the efficiency of the depletion mechanism and the recovery of reserves. We will deal with one cause of this water production, namely, coning. The coning of water into a producing well is caused by pressure gradients established around the wellbore by the production of fluids from the well. These pressure gradients can raise the water-oil contact near the well where gradients are most severe. Gravity forces that arise from fluid density differences counterbalance the flowing pressure gradients and tend to keep the water out of the oil zone. Therefore, at any given time, there is a balance between gravitational and viscous forces at points on and away from the completion interval. When the dynamic forces at the wellbore exceed gravitational forces, a cone of water will ultimately break into the well to produce water along with the oil. We can expand on this basic visualization of coning by introducing the concepts of stable cone, unstable cone and critical production rate. For instance, if a well is produced at a constant rate and the pressure gradients in the drainage system have become constant, a steady-state condition is reached. If, at this condition, the dynamic forces at the well are less than the gravity forces, then the water or gas cone that has formed will not extend to the well. Moreover, the cone will neither advance nor recede, thus establishing what is known as a stable cone. Conversely, if the pressure in the system is in an unsteady-state condition, then an unstable cone will continue to advance until steady-state conditions prevail. If the flowing pressure drop at the well is suffcient to overcome the gravity forces, the unstable cone will grow and ultimately break into the well. It is important to note that in a realistic sense, stable cones may only be "pseudostable" because the drainage system and pressure distribution generally change. For example, with reservoir depletion, the water-oil contact may advance toward the completion interval, thereby increasing chances for coning. As another example, reduced productivity due to well damage requires a corresponding increase in the flowing pressure drop LO maintain a given production rate. This increase in pressure drop may force an otherwise stable cone into a well. The critical production rate, well known in the literature, is the rate above which the flowing pressure gradient at the well causes water (or gas) to cone into the well. It is, therefore, the maximum rate of oil production without concurrent production of the displacing phase by coning. At the critical rate, a built-up cone is stable but is at a position of incipient breakthrough. Numerous papers have been published on critical rates. Some of the better known of these include the work of (1) Muskat and Wyckoff,' who first dealt with the coning problem; (2) Chaney, et al,' who developed expressions similar to those of Muskat but who presented results in a conven ient-to-use graphical form (the "Sun" method); and (3) Meyer and Garder,b hose analysis is based on radial-flow formulas. One assumption in critical production rate analyses is that the cone has built-up to just before its breakthrough into the well. But, these analyses reveal nothing directly about the time it takes for the cone to build up to this incipient breakthrough position. Thus, water-free oil can be produced from a well for prolonged periods at rates above the critical rate before the well reaches the condition to which the critical rate applies. The published literature contains little on the rate of growth of a cone. Experimentally, Meyer and Searcy studied the rate of rise of a cone in a Hele-Shaw model.' Additional related work on water breakthrough and produced water-oil ratios in water driven reservoirs was reported by Muskat," Hutchinson and Kemp," Henley, et al,' and Stevens, et a1.B Theoretically, the basic coning equations for a water-oil system can be developed by applying the conservation of mass to each of the phases, relating flow velocities with pressure by Darcy's law, and relating pressures across water-oil interfaces by capillary pressure. With the usual boundaries at the well and reservoir limits, the solution of the resulting equations for the time behavior of a water-oil interface constitutes a free-surface, boundary-value problem. The shape of the boundary depends on the internal potential distribution, while this

Jan 1, 1966

By G. L. Chierici

Starting from Muskat's theory of water and gas coning, maximum permissible oil production rates without water and/or free-gas production have been determined, in a broad range of reservoir and well parameters, using the potentiometric model technique. The main assumptions made are as follows: (1) the reservoir rock is homogeneeous (either isotropic or anisotropic); (2) the volume of the aquifer underlying the oil zone is very small, so that it does not contribute to reservoir energy; and (3) the gas cap expands at a very low rate, so that it can be assumed to be in quasi-static conditions. The results obtained are presented in the form of diagrams which can be used for solving two types of problems: (1) given the reservoir and fluids characteristics, as well as the position and length of the perforated interval, determine the maximum oil production rate without water and/or free-gas production; and (2) given the reservoir and fluids characteristics only, determine the position and length of the perforated interval which optimize the maximum permissible oil production rate, without water and/or free-gas production. INTRODUCTION In oil reservoirs where the oil-bearing formation is underlain by an aquifer which does not participate in the production mechanism, water-coning is a limiting factor to the flow rates of producing wells. Production rates are usually kept to a value that will prevent the water from entering the wells. The entry of water into a well lowers its productivity by increasing the weigbt of the fluid column; moreover, the separation of water from the effluent, at the surface, may constitute a very difficult problem in cases of heavy viscous oils. A similar situation is encountered in oil reservoirs with a gas cap overlying the oil-saturated zone; here a downward gas cone is induced by the flow of oil towards the producing wells. Production rates must be low enough to prevent the gas from being produced; producing gas from the gas cap would be a waste of energy. Of course, water-coning and gas-coning phenomena can occur at the same time in the same reservoir if the oil-producing formation is both overlain by a gas zone and underlain by a water zone. Due to its relevant practical importance, the mechanism of coning was studied by many people.2,3,5-8 Defining the conditions for getting the maximum water-free and/or gas-free oil production rate is a difficult problem, often encountered under one of the following aspects: 1. Predict the maximum flow rate that can be assigned to a completed well without the simultaneous production of water and/or free-gas. 2. Define the optimum length and position of the interval to be perforated in a well, in order to obtain the maximum water and gas-free production rate. A systematic study of these problems was made by means of the electrical analog technique. The results of this study are presented here, under the form of a set of curves providing solutions for the above stated problems. These curves are valid only for homogeneous forrnations, either isotropic or anisotropic. Should the formation be non-homogeneous (by horizontal or vertical variation of permeability, shale diaphragms, fractures, etc.), a specific potentiometric study would be required for each specific case. Especially when shale diaphragms of some radial extension are present, the critical rates observed are much larger than would be expected from the diagrams. STATEMENT OF THE PROBLEM In the present study the aquifer is supposed to be of such limited volume that it does not contribute to the energy of the reservoir. Moreover, the gas cap is supposed to expand at such a low rate that the potential gradient in the gas cap is negligible. Under static conditions water-oil and gas-oil interfaces (T1 and T2) are both horizontal. When the reservoir production starts, below each well these interfaces take a cone-like shape (Fig. 1) having as an axis the axis of the well. This shape results from the equilibrium between potential gradients in the oil zone and gravitational forces due to density differences between oil and water and between oil and gas. Assuming the oil-bearing formation to be homogeneous and the oil to be incompressible, the analysis of the problem (see Appendix) shows that the oil-water and gas-oil interfaces are stable only if the oil production rate of the well is not higher than the following values.

Jan 1, 1965

By M. A. Torcaso, P. Raimondi

The distribution of the oil phase in Berea sandstone resulting from increasing and decreasing the water saturation by imbibition was investigated Three types of distribution were recognized: trapped, normal and lagging. The amount of oil in each of these distributions was determined as a function of saturation by carrying out a miscible displacement in the oil phase under steady-state conditions of saturation. These conditions were maintained by flowing water and oil simultaneously in given ratios and by using a displacing solvent having essentially the same density and viscosity as the oil. A correlation shows the amount of trapped oil at any saturation to be directly proportional to the conventional residual oil saturation Sor The factor of proportionality is related to the fractional permeability to the water phase. Part of the oil which was not trapped was displaced in a piston-like manner (normal part) and part was eluted gradually (lagging part). The observed phenomena are more than of mere academic importance. Oil which is trapped may well provide the fuel essential for forward combustion and thus be beneficial. On the contrary, in tertiary recovery operations, it is this trapped oil which seems to make current techniques uneconomic. INTRODUCTION A typical oilfield may initially contain connate water and oil. After a period of primary production water often enters the field either from surrounding aquifers or from surface injection. During primary production evolution and establishment of a free gas saturation usually occurs. The effect and importance of this third phase is fully recognized. However, this investigation is limited to a two-phase system, one wetting phase (water) and one non-wetting phase (oil). The increase in water content of a water-wet system is termed imbibition. In a relative permeability-saturation diagram such as the one shown in Fig. 1, the initial conditions of the field would be represented by a point below a water saturation of about 35 per cent, i.e., where the imbibition and the drainage curves to the non-wetting phase nearly coincide. When water enters the field the relative permeability to oil decreases along the imbibition curve. At watered-out conditions the relative permeability to the oil becomes zero. At this point a considerable amount of oil. called residual oil, (about 35 per cent in Fig. 1) remains unrecovered. Any attempt to produce this oil will require that its saturation be increased. In Fig. 1 this would mean retracing the imbibition curve upwards. In addition, processes like alcohol and fire flooding, which can be employed at any stage of production, involve the complete displacement of connate water and an increase, or imbibition. of water saturation ahead of the displacing front. Thus, in several types of oil production it is the imbibition-relative permeability curve which rules the flow behavior. For this reason a knowledge of the distribution of the non-wetting phase, as obtained through imbibition, whether "coming down" or "going up" on the imbibition curve, is important.

Jan 1, 1965

By G. C. Bernard, L. W. Bernard, L. W. Holm, W. L. Jacobs

The effect of loam on the permeability of porous media to water was studied as a {unction of foaming agent concentration, specific permeability, pressure gradient, length of a porous medium and its oil saturation. At a given fluid saturation in a porous medium, the permeability to water was found to be the same whether loam was present or not. Foam decreases the permeability to water by developing a higher trapped gas saturation than that obtained by water flooding without loam present, increasing the concentration of foaming agent increased the trapped gas saturation and thereby decreased the permeability to water. The presence of oil reduced the capability of most foaming agents to decrease the permeability of a porous medium to water. A few surfactants were found to be effective loaming agents even in the presence of oil. These results are similar to those reported in a previous paper on the effect of loam on the permeability of porous media to gas. The effect of loam was found to persist in long porous media at moderately high reservoir temperatures and during the passage of many pore volumes of surfactant-free water. INTRODUCTION This paper describes part of a study on a novel approach in the use of surfactants for oil recovery; the use of foam rather than water to displace oil. Previously it was found that foam can displace oil which normally is not displaced by water.4 The foam is formed by successively injecting a suitable surfactant solution and gas into a porous medium. Foam appears to have at least two uses in the field: (1) it shows promise as a superior oil recovery agent, and (2) it shows promise as a selective permeability reducing agent. Foam may be very useful in water floods, or in other oil recovery processes, where highly permeable streaks or unfavorable mobility ratios are a problem. A previous paper reported the effect of foam on the permeability of porous media to gas.1 In the present study the effect of foam on the permeability of porous media to water is reported. The specific objectives of the study were to determine: (1) the effect of foam on the permeability to water in porous media of various specific permeabilities, (2) the effect of foam on the permeability to water in the presence of oil, (3) the effect of foam and crude oil on the trapped gas saturation, (4) the effect of foam on permeability to water at trapped gas saturation, (5) the effect of pressure gradient on the permeability to water under foaming conditions, (6) the persistence of foam during the passage of surfactant-free water through the porous medium, and (7) the effect of various foaming agents, length of the porous medium and temperature on the permeability reduction caused by foam. EXPERIMENTAL PROCEDURES EQUIPMENT AND MATERIALS The experimental apparatus consisted of consolidated and unco n soli dated porous media, wet test meters and constant delivery pumps. The porous media consisted of consolidated sandstone cores (6 to 36 in. long), and unconsolidated sand packs (3 to 30 ft long). The consolidated cores had permeabilities of 32 and 1,000 md and porosities of about 20 per cent. The sand packs had permeabilities of 3,500 to 211,000 md and porosities of about 40 per cent. (Throughout this report a term such as 100 md sand" is used. This term means that the porous medium had a dry, nitrogen permeability of 100 md.) Fluids used in the experiments were distilled water, 1 per cent NaCl solution, aqueous solutions of foaming agents, nitrogen gas, air and crude oil.

Jan 1, 1966

By D. M. Waldorf

Steam permeability measurements have been made in the laboratory on several samples of natural reservoir materials. The steam temperatures and pressures were selected to simulate conditions which might exist in a reservoir during the injection of steam. For each sample tested, the experimental permeability to superheated steam was comparable to that measured with air and no evidence of plugging was detected. Some samples were exposed to water at various temperatures and plugging was found to occur in materials which contained significant quantities of monmorillonite clay. Temperature had little effect on the degree of plug-ning between 75 and 325 F. The measured pemeabilities tended to increase slightly with temperature, but the changes were small compared with the initial loss of per~neability on wetting. Sequential pemzeability measurements were made on two samples using air, water, steam, water and air, in that order. Both samples were water-sensitive and plugged extensively after the initial injection of water. Upon exposure to superheated steatm the samples dehydrated and their permenbilities to superheated steam were comparable to those initially measured with air. The remaining measuretnetzts with water and air confirmed that the water plugging was reversible and that the samples were not seriorrsly damaged during the tests. INTRODUCTION The swelling of water-sensitive clays during water floods has long been recognized as a potential source of reservoir damage. The recent extensive application of steam injection and stimulation has compounded this problem since both hot water and steam (as well as fresh water at reservoir temperatures) are, at sume time, in contact with the producing zone adjacent to the bore of a steam injection well. The purpose of this paper is to present data which compare the sensitivity of some natural sedimentary rock samples to water at various temperatures, and to super-heated steam. Some properties of montmorillonite clay are briefly reviewed, and comparisons are drawn between empirical data and the predicted behavior of the montmorillonite known to be present in the samples. PROPERTIES OF MONTMORILLONIT E CLAY Water initially adsorbs on dry N a -montmorillonite clay in discrete layers in the interlaminar space between clal platelets. The platelet spacing, which is 9.6 A (angstroms) for a dehydrated clay, has been observed to expand in discrete steps to 12.4, 15.5, 18.4 and 21.4 A spacings, indicating the formation of four discrete layers of regularly oriented water molecules.' The first two layers are easily formed by hydrating a dry sample to equilibrium in an atmosphere with carefully controlled humidity. The formation of the higher layers is more difficult. The usual X-ray diffraction patterns of the more highly hydrated samples indicate a gradual increase in the average spacing betwcen 15.5 and 19.2 A, followed by a discontinuous expansion to 31 A when the weight ratio of water to dry clay is between 0.5 and 1.2.' Platelet expansion above 31 A proceeds monotonically as the moisture is increased and no regular arrangement of the platelets ib observed. Water-sensitivity in sedimentary rocks is usually associated with Na-montmorillonite clay when it is in the noncrystal-line state. Mering3 found that the average lattice spacing of sodium montmorillonite hydrated at 68 F and 70 per cent relative humidity was 15.5 A, and that the spacing, at 92 per cent humidity was 16.5 A. The water adsorbed at the higher humidity has the same free energy as liquid water at 65.6 F. Kolaian and Low' used a tensiometer to measure the thermodynamic properties of water in diffuse suspensions of montmorillonite clays relative to pure water. They observed that water in suspensions as dilute as 6 per cent clay became partially oriented when left undisturbed. The bonding associated with this orientation was not extensive because the free energy difference between the water in suspension and pure water was only a few millicalories per mole. They also found that the measured free energy difference decreased rapidly with temperature and became negligible above 100 F. This evidence indicates that montmorillonites contained in sedimentary rocks would dehydrate to a crystalline structure when exposed to superheated steam, and that the rock permeability measured with steam would be equivalent to that measured with air. The effect of elevated temperatures on the swelline of montmorillonite clays in aqueous suspensions has not been investigated. The Gouy-Chapman diffuse-ion-layer theory has been used to predict the swelling pressure of clay suspensions in dilute salt solutions at room temperature with reasonable success. theory also correctly predicts the direction of the thermal response of Na-mont-morillonite swelling pressures in dilute salt suspensions, 9 Over the temperature range of 33 to 68 F, an increase in

Jan 1, 1966

By J. R. Kyte, R. L. Dalton, V. O. Naumann, R. L. Ridings, R. W. Greene

This paper presents a one-dimensional numerical method for calculating dissolved-gas-drive behavior for a single well or a linear reservoir. Pressure and saturation distributions can be calculated, including the effects of capillary and gravity forces. The analysis has been programmed for solution on the IBM 7090 digital computer. A comparison of experimental and calculated behavior for two linear models demonstrates that the method is suitable for calculating dissolved-gas-drive behavior under reservoir conditions of flow. Dissolved-gas-drive tests were run, at different production rates and with oils of different viscosities, on short (6.ft) and long (173ft) flow models. For tests approximating reservoir pressure gradients and flow rates, there was excellent agreement between experimental and calculated results. A series of calculations of well behavior for thin, horizontal, homogeneous, dissolved-gas-drive reservoirs shows that ultimate recovery is independent of production rate or well spacing in such cases. These calculations also show that Muskat's method and similar material-balance analyses are not applicable under some conditions particularly for high pressure drawdowns. INTRODUCTION The problem of calculating dissolved-gas-drive reservoir performance has received a great deal of attention throughout the oil-producing indusfry for many years. Really significant advances, however, have been made only in recent years. The rapid evolution of high-speed computing equipment and techniques has opened up new approaches to the problem and led to realistic, quantitative analyses free of many of the assumptions and limitations of earlier methods. The one-dimensional numerical method described in this paper includes the effects of both gravity and capillary forces. A number of test calculations, based on different rates, spacings, and rock and fluid properties, have shown that the method is adequate for systems that can be reasonably represented by a one-dimensional model (linear or radial) and in which fluid segregation perpendicular to the direction of flow is not important. Along with the development of mathematical methods, an experimental program was conducted and the data were used to test the validity of the theoretical basis of the mathematical method. Linear models were used and all experiments were run using 1,000-psig bubble-point oils. The use of long, high-pressure models enabled us to more closely approach reservoir gradients and flow rates. MATHEMATICAL FORMULATION DIFFERENTIAL EQUATIONS The differential equations describing dissolved-gas-drive behavior in one dimension dare summarized here. They include Darcy's law for flow of oil and gas (Eqs. 1 and 2) and the equations for conservation of mass for oil and gas (Eqs. 3 and 4). The oil and gas pressures are related by capillary pressure (Eq. 51.

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By H. Groeneveld, C. A. Connally, P. J. Hoenmans, J. J. Justen, W. L. Mason

A miscible displacement pilot using a slug of LPG driven by separator gas was conducted in the Cardiurn reservoir of the Pembina field. The injection pattern was a 10-acre, inverted, isolated five-spot. Upon completion of the LPG-gar phase, an experiment was conducted using a slug of water followed by gas. Calculated performance of the pilot is compared with actual performance. Equations are developed to calculate the distribution of LPG into zones of varying permeability, to estimate the progress of the flood at different times in the various zones and to estimate gas rates after breakthrough. The analysis indicates that permeability stratification was a dominant factor in controlling oil recovery and that oil was completely displaced from the swept pore volume. The results of the pilot indicated that miscible flooding is a practical means of pressure maintenance in this reservoir. The total recovery from the pilot area was good in spite of the early breakthrough of LPG. The effects of stratification were reduced by injecting a slug of water into the partially swept reservoir. INTRODUCTION The Pembina field,' located in Alberta, is the largest oil field in Canada and one of the largest in the North American continent. The reservoir is a stratigraphic trap producing from the Cardium sand. Neither bottom water nor free gas has been found. The recovery of oil by the natural depletion mechanism has been estimated at 12.5 per cent. Pressure maintenance studies of various areas have indicated that the recovery can be increased 21/2 times by water flooding, and a large area of the field is presently under water flood. However, reservoir studies of the North Pembina area indicated that miscible flooding might be competitive with water flooding. A pilot test was conducted to evaluate the performance of a miscible flood. A 10-acre, inverted, isolated, five-spot pattern was selected for the pilot. The pattern area was large enough to minimize wellbore fracturing effects and contained sufficient oil to provide significant working numbers. The performance of each of the four producers could be evaluated individually and compared. In the event of breakthrough in one direction, the effect would be isolated from the other producers. The use of a single injector minimized the volume of LPG required, and, because of the high mobility of gas, one well was sufficient to inject the necessary daily volume to replace the high rate of production. With four producers, the test could be completed in time for results to be evaluated, additional engineering studies to be made and a unit to be formed before the reservoir pressure in the North Pembina area declined below the bubble point. The pilot was located in an area developed on staggered, 80-acre spacing. The injection well was drilled at a regular location, while the four producers were drilled 467-ft north, east, south and west of the injector. Each quadrant and its associated producer were identified according to their direction from the injector— that is, north, east, south or west. The eight surrounding producers on 80-acre spacing were shut in to isolate the pilot area and provide for reservoir pressure observation. The pilot wells were completed using permanent-type completion techniques. After coring, casing was run through the pay section and cemented. Inside 51/2-in. casing, 2 1/2-in. tubing was hung. The wells were perforated opposite the Upper Cardium sand and lightly fractured. Fracturing volumes, rates and pressures were low to minimize the extent of the fractures. The fracturing treatments average 1,000 lb of 20-40 mesh sand in 700 gal of a low fluid-loss sand-carrying agent. Feed rates and wellhead fracturing pressures averaged 5.5 bbl/min at 2,535 psig, respectively. After fracturing, the productivity index was measured in each of the five pilot wells. The average PI of the four producers was 0.41 BOPD/psig drawdown. The measured PI'S were approximately the same as PI'S calculated from core analysis data, indicating that the fracturing treatments were just sufficient to overcome

By D. H. Henley, F. F. Craig, W. W. Owens

The oil recovery performance of systems producing entirely by bottom-water encroachment has been experimentally determined in a series of scaled laboratory-model tests. The effects of well spacing, fluid mobilities, rate of production, capillary and gravity forces, well penetration and well completion techniques on the oil recovery performance have been investigated. The laboratory tests were performed using two uniform, un-consolidated sand-pack models. The models have ratios of the interwell distance to the formation thickness of 12 and 2, respectively. Tests at constant total fluid production rate were performed simulating a range of uniform reservoir characteristics and operating conditions encountered in field operations. The performance was determined by material balance and by observation of the encroachment of dyed fluids into the models. The results of the model tests agreed with those obtained mathematically when the conditions previously considered in theoretical studies were simulated, that is, when the oil and water are of equal density and no capillary forces exist. The model study of bottom-water drive indicated that certain variables can aoect the oil recovery performance to a greater degree than can be predicted by present analytical methods. In one comparison, the oil recovery at a water-oil ratio of 20 (obtained at a wide well spacing) varied as much as threefold, depending upon the system's properties and the production rate. Lesser effect of mobility ratio and no eflect of capillary forces over the range studied were observed. The test 'results also showed that the deeper the well penetration into the oil column, the greater the total water production to a producing WOR of 20. However, the ultimate sweep efficiency, and so the oil recovery to this level of WOR, did not vary significantly with well penetration. Horizontal fractures at the top of the formation did not significantly change the sweep characteristics of the reservoir models when values of radius and fracture capacity encountered in actual reservoirs were used. Impermeable pancakes at the bottom of the well moderately increased the oil recovery efficiency both at water breakthrough and at high water-oil ratios. A method is outlined by which the oil recovery performance of other uniform bottom-water drive systems can be estimated from the information obtained in these model tests. INTRODUCTION When oil is produced from a well which partially penetrates an oil zone completely underlain by water, the water rises directly beneath the well in a symmetrical cone when the system is uniform. Two different flow mechanisms can cause the water cone to form—coning and bottom-water drive. In coning, the aquifer is relatively inactive and the cone is formed beneath the well by the pressure gradients associated with the oil flow to the well. The oil can be produced by a solution-gas drive, an edge-water drive or other driving forces in the interwell area. In a bottom-water drive, the driving force for oil production comes from an upward encroachment of the underlying active aquifer. Two papers have analyzed the theoretical performance characteristics of bottom-water drive reservoirs. In the initial mathematical investigation, Muskat' established the equations which determine the pressure distribution in this type of reservoir and solved these equations for certain conditions. Specifically, it was assumed that the water and oil had equal mobilities and equal densities, there were no capillary forces, the pressure throughout the oil zone remained above the bubble-point pressure, a constant pressure existed at the initial water-oil contact and the oil was completely displaced by the encroaching water. These assumptions were used in obtaining analytical solutions. In general, Muskat found that the sweep efficiency to initial water breakthrough to the well was larger for the thicker oil zones, the closer well spacings, the lower ratios of vertical to horizontal permeabilities, the smaller the penetration of the well into the oil zone and the smaller the bore size of the well. The production history after water breakthrough was expressed as a volumetric sweep efficiency at a given producing water-oil ratio. The results indicated that cumulative oil production at producing water-oil ratios of 10 is less affected by the well spacing than is the water-free production history. Muskat studied well spacing which today would be regarded as close. The maximum value of his dimensionless well spacing (ratio of interwell distance to formation thickness) was 4.3. This would require the development of a 50-ft-thick oil sand on less than 10-acre spacing if the vertical and horizontal permeabilities were equal, with

By J. H. Henderson, J. Naar

The displacement of a wetting fluid from a porous medium by a non-wetting fluid (drainage) is now reasonably well understood. A complete explanation has yet to be found for the analogous case of a wetting fluid being spontaneously imbibed and the non-wetting phase displaced (imbibition). During the displacement of oil or gas by water in a water-wet sand, the porous medium ordinarily imbibes water. The amount of oil recovered, the cost of recovery and the production history seem then to be controlled mainly by pore geometry. The influence of pore geometry is reflected in drainage and imbibition capillary-pressure curves and relative permeability curves. Relative permeability curves for a particular consolidated sand show that at any given saturation the permeability to oil during imbibition is smaller than during drainage. ' Low imbibition permeabilities suggest that the non-wetting phase, oil or gas, is gradually trapped by the advancing water. This paper describes a mathematical image (model) of consolidated porous rock based on the concept of the trapping of the non-wetting phase during the imbibition process. The following items have been derived from the model. 1. A direct relation between the relative permeability characteristics during imbibition and those observed during drainage. 2. A theoretical limit for the fractional amount of oil or gas recoverable by imbibition. 3. An expression for the resistivity index which can be used in connection with the formula for wetting-phase relative permeability to check the consistency of the model. 4. The Iimits of flow performance for a given rock dictated by complete wetting by either oil or water. 5. The factors controlling oil recovery by imbibition in the presence of free gas. The complexity of a porous medium is such that drastic simplifications must be introduced to obtain a model amenable to mathematical treatment. Many parameters have been introduced by others in Iprogressing" from the parallel-capillary model to the randomly interconnected capillary models independently proposed by Wyllie and Gardner3 and Marshall4. To these a further complication must be added since an imbibition model must trap part of the non-wetting phase during imbibition of the wetting phase. Like so many of the previously introduced complications, this fluid-block was introduced to make the model performance fit the observed imbibition flow behavior. THE IMBIBITION MODEL ASSUMPTIONS The following assumptions are common to many two-phase models previously described. 1. The wetting characteristics of the medium are clearly defined. 2. When two fluids co-exist in the medium, the wetting phase will fill the smaller pores and the non-wetting fluid the bigger ones. 3. The pores have a circular cross section. 4. Apparent pore radii are related to the capillary pressure (drainage) by

By E. W. Hough, G. L. Stegemeier

Empirical equations for surface tension of propane and normal butane as functions of reduced temperature are obtained from experimental data. Another correlation relating surface tension to enthalpy of vaporization is given for these two compounds. In addition, new paracbor numbers are calculated for the normal paraffin hydrocarbons. These numbers are utilized for the calculation of interfacial tension of two-component systems as functions of pressure and temperature, using a modified form of Weinaug-Katz equation. The experimental data for two binary systems are approximated by the correlation. From these results it is found that the interfacial tension in the high-pressure region remains extremely low at large pressure decrements below the critical pressure. Thus, it appears that condensate systems may have flow characteristics almost like single-phase conditions even though the pressure is within the two-phase region. Experimental data have shown that interfacial tension divided by density difference approaches zero as the critical pressure is approached. A calculation of wetting-phase saturations indicates that the saturation gradient at the two-phase contact becomes increasingly abrupt as the critical pressure is approached. DISCUSSION Prediction of the surface and interfacial tension of the light hydrocarbons and of two-component hydrocarbon mixtures at various temperatures and pressures may be made if other physical properties are known. Extensive experimental work on single-component and binary systemsl is the basis for the correlations outlined in this paper. Interfacial tension is defined as the specific surface-free energy between two phases of unlike fractional composition, while surface tension is defined as the specific surface-free energy between two phase:; of the same fractional composition. The usual definitions relating interfacial tension to a liquid-liquid interface and surface tension to a gas-liquid interface are not clearly defined when the critical region is included, and there is no sharp distinction between a gas and a liquid phase. Interfacial tension is probably the most important single force that makes one-half to one-third of the total oil actually in place in a reservoir rock unrecoverable by conventional gas-drive or waterflood methods. A rough estimate of this figure for the United States is 100 billion bbl. Interfacial tension presently is used by petroleum engineers in the estimation of saturation gradients at the gas-oil contact and at the oil-water contact. The data in this paper should prove useful for estimates of reserves involving gas-oil contacts. Relative permeability undoubtedly is influenced by interfacial tension, for sufficiently small values. These data should be useful in determining how small the values are In addition, these data should eventually add to our fundamental knowledge of surfaces. At the critical point, all surface excesses approach zero and the thickness becomes very large. SINGLE-COMPONENT SYSTEMS It has been observed1 that the following relationships are good approximations to the physical properties of propane and n-butane.

By N. Marusov, D. K. Lowe, J. C. Karp

This paper considers, from an engineering viewpoint, several factors involved in creating, designing and locating horizontal barriers for controlling water coning. This is an effort to consolidate new concepts with previous information so that a reasonable selection of barrier materials, dimensions and vertical position can be made. Coning theories previously developed are briefly reviewed and an effort i,s made to reduce the results of coning-theory calculations to a point where routine calculations can be made with a desk calculator. It is expected that these simplified calculations will give usable predictions of the amounts of improvement attainable with barriers of various dimensions. Apparatus and procedures used for testing the suitability of commercial cements are described and test results are presented. INTRODUCTION Just about as soon as the phenomenon of water coning was recognized as an oilfield problem, there were suggestions that the production of water coning could be controlled or completely suppressed by means of horizontal barriers. It was also suggested that natural barriers such as shale streaks were helpful in restricting bottom-water production. The implication was that wells not having the benefit of continuous shale streaks to suppress bottom waters should, in some way, be supplied with an artificial barrier. With regard to the selection of barrier materials, there are many references which give the merits of liquid barriers.',' These liquids include surface-active agents," recipitates' and emu1sions. There is no real need to review the merits of these materials again. This paper deals with the use of solid barriers which could be made from various cements. The two major restrictions placed on these barrier cements are that (1) they must be commercially available in large quantities, and (2) they must be applicable by means of standard tools and equipment with techniques that have been fully developed and are in general use at the present time. A schematic drawing of a horizontal barrier in a reservoir is shown in Fig. 1. Two maximum stable water cones are illustrated-—-one for the wellbore diameter and the other for a barrier. A maximum stable cone is one that has attained its maximum volume at the critical production rate and is on the threshold of producing water. The volume displaced by each water cone and the rise in water table are indications of the amount of water-free oil produced. Corrections for porosity, connate water and residual oil, of course, must be applied. The water-oil interface which defines a maximum stable cone is called a free surface. It is a surface over which the pressure is constant. This free surface has the shape of a limiting streamline at and below which there is no movement or transport. The shape of the water cone describes the pressure drop in the reservoir. The effect of a completely impermeable barrier on the cone shape is essentially the same as extending the wellbore out to the barrier radius. The placement of a horizontal barrier requires that a horizontal fracture be created using a single-point entry technique. Two methods may be employed for filling this fracture. One method consists of propping the fracture before filling with cement. The other consists of fracturing with fracture fluid, then immediately following with cement