Search Documents

By A. K. Csazar, L. W. Holm

This paper describes our investigation of a post-water-flood, oil recovery process which consists of injecting a slug of propane followed by water. Also described are the results obtained by applying a modification of the process in which gas was injected ahead of the water. Under the conditions of the latter experiments, misci-bility was not achieved between the propane and gas. Preliminary experiments or) uniform, watered-out sandstone cores showed that an oil bank could be formed and produced by applying this recovery process. However, since reservoirs are not uniform in structure, the process also was applied to porous media containing irregular porosity and to stratified sand systems. As a supplenzerrt to the experinlental work, a mathernatical procedure was developed for calculating the performance of the recovery process in a bounded, layered, porous system with crossflow between layers. As a specific example, the method was applied to predict the perforrnance of the recovery process in a 6-ft long, two-layer, stratified, unconsolidated sand model for comparison with experinlental data. The calculations were programed for the ZBM 704 computer. The equations and calcula-tional procedure presented can be extended to systems containing any number of randomly distributed permeability variations or any number of parallel layers. INTRODUCTION The problem of recovering the oil that remains in a reservoir which has been waterflooded is receiving considerable attention now as an increasing number of water floods reach an economic limit. A large number of the waterflood projects are in shallow reservoirs which are at pressures below 1,000 psi. It has been demonstrated in the laboratory that post-waterflood oil can be recover-ered by miscible displacement, but the LPG-gas, miscible flood and the enriched gas drive cannot be applied effectively at pressures below 1,000 psi. Only a few reports have appeared in the literature2-4 on low pressure, partially miscible recovery methods. However, it is possible to use LPG in a partially miscible displacement process in a reservoir where pressures of 200 to 1,000 psi can be achieved. Under these Pressures and at normal reservoir temperatures, propane is miscible with the oil; but, of course, gas or water used to drive the propane slug would not be miscible with the propane. Because of the lack of complete miscibility, it has generally been concluded that excessive amounts of propane would be required to recover oil and that such a recovery method would not be economical; however, we have found that under conditions present in certain reservoirs, an imrniscible recovery process can be applied effectively. The oil saturation in reservoirs at the economic limit of waterflood projects is usually in the range of 20 to 35 per cent of the pore space." A certain portion of this oil is left trapped by water in various size pores of the rock, but a good part of this so-called "residual" oil can be present in the less permeable lenses or layers of the reservoir rock which were by-passed to some degree by the water. The oil in these permeability traps can be produced only if favorable pressure gradients are formed in the reservoirs between adjacent zones of high and low permeabilities. A low viscosity liquid, miscible with the oil in place, which is driven by water through a stratified or heterogeneous porous system can aid in the development of these favorable pressure gradients. The oil that is released thereby from the permeability traps can be recovered by the subsequent water flood. Studies were made to determine how much oil could be recovered from homogeneous and stratified cores and models, which had been water flooded, by injecting a slug of propane and driving it with water. The effect of injecting a slug of gas ahead of the water was also determined. Most of the work described herein was done with the propane-water combination; unless otherwise specified, no gas was injected. The principal objectives of the investigation were to determine (1) if an oil bank could be formed and (2) what ratio of oil recovered to propane injected would be obtained. A further objective was to develop a method for calculating fluid-flow performance in stratified systems which would account for fluid transfer between zones in hydrodynamic communication but of different permeabilities. THEORETICAL ANALYSIS In a theoretical study of the recovery process, analytical expressions were derived to calculate the pressure distribution, the fluid flux in longitudinal (parallel to layers) and transversal (across the layers) directions, and the fluid distribution at any point in the system. The equations were developed for a two-layer porous system in which it was assumed that the fluids in the system were incompressible and that capillary and gravity effects were

By M. Mortada

Nonsteady-state flow of slightly compressible liquids in porous media of non-uniform properties has been the subject of a number of recent studies. Most of these studies considered one-dimensional flow in systems with cylindrical geometry and with the properties of the porous medium varying with radial distance only. A. Houpeurt' studied the problem in its general form by examining the solutions of the differential equation, In Eq. 1, the permeability k and the porosity + are arbitrary functions of radial distance. Houpeurt concluded that, except for special forms of the functions k and the problem is not well suited for mathematical treatment, but lends itself to high-speed digital computation. Carslaw and Jaeger' discussed the solution of Eq. 1 for the case in which the permeability varies as a power of the radial distance. P. Albert3 studied analytically the pseudo steady-state solutions in systems of finite radial extent and non-uniform properties. He also obtained some numerical solutions for nonsteady-state flow using a high-speed computer. William Hurst4 examined the interference pressure drop due to a point sink located in an infinite system consisting of two permeability regions in series. In this case, the permeability changed at a radial distance r from the point sink. INTERFERENCE BETWEEN OIL FIELDS IN A NON-UNIFORM AQUIFER The interference pressure drop for oil fields located in an extensive aquifer of uniform properties is discussed in a number of places in the literature. This work considers the interference pressure drop for oil fields located in a non-uniform extensive aquifer comprising two regions of different properties as shown in Fig. 1. Region I of the aquifer extends from which is the radius of Oil Field A, to ar,,, where a>1. It has permeability k,, porosity +, and effective compressibility c,. Oil Field B is located in Region II which extends from ar, to infinity and has different properties from those of Region I, but the same formation thickness h. Region II has permeability k2, porosity and effective compressibility c2. The choice of this type of aquifer should not imply that it is of frequent occurrence. It is chosen primarily to illustrate the effects of non-uniform aquifer properties on oilfield interference. It is of interest to point out that a non-uniform radial aquifer in which both the permeability and the porosity are inversely proportional to the radial distance can be treated in a manner similar to a linear aquifer of uniform properties. The situation where only the permeability varies as the inverse of the radial distance is also amenable to analytical treatment. Both aquifers behave quite differently from a radial aquifer with uniform properties so far as the water influx and the interference pressure drop are concerned. METHOD OF SOLUTION The expression for the interference pressure drop in Oil Field B due to a constant rate of water influx in Oil Field A is developed in the Appendix. This is accomplished by the simultaneous solution of two partial differential equations describing the pressure in Regions I and 11. The two differential equations are solved subject to the following conditions: (1) at r,, the pressure gradient is constant; (2) at ar the pressure and

By G. W. Tracy, R. D. Carter

Numerical calculations were made to determine the behavior of reservoirs with high-pressure drawdown and wide well spacing where the initial productivity is low and the wells are completed by hydraulic fracturing. The two-phase flow equations were solved for the flow into a single well. This well was assumed to be producing from a reservoir with hydraulically created horizontal fractures (four different systems with fractures were studied). For comparison purposes, additional two-phase flow calculations were made assuming a reservoir with uniform rock properties. The two-phase flow results were also compared with the conventional calculation methods, which do not include the effect of saturation gradients resulting from a simultaneous flow of oil and gas which are normal to this type reservoir. It was found that the conventional methods predict (1) a high and too optimistic value of ultimate recovery, (2) a high producing rate and a high reservoir pressure at a given oil recovery and (3) a low trend of gas-oil ratio with oil recovery. Included in the two-phase flow calculations were provisions to control the oil production rate by an allowable rate and, also, by a gas-oil ratio penalty rule. For the systems with hydraulic fractures, the producing rate was controlled by the gas-oil ratio penalty rule for most of the life. This is in contrast to the system with uniform rock properties which went "on decline" almost immediately. An unexpected characteristic of the systems which included fractures was the early rise in producing gas-oil ratio from 730 cu ft/bbl to approximately 1,200 cu ft/bbl, followed by a "leveling off" before the normally expected gas-oil ratio rise began. Additional features which are a result of hydraulic fracturing are (I) greater ultimate recovery, (2) higher average producing rates and (3) a lower average reservoir pressure at a given oil recovery. INTRODUCTION Some oil fields discovered during the past few years are producing from certain volumetric ally controlled reservoirs (often referred to as solution or internal gas- drive reservoirs) which are characterized by high-pressure drawdown at the wells. Since the available pressure drawdown at a well is limited by the static reservoir pressure and the producing rate is controlled by the available drawdown, wells completed in this type of reservoir usually produce at a rate less than the allowable from the time of completion. Because of this, this type of reservoir is referred to as a low productivity reservoir. Economic considerations require the use of wide well spacing and well stimulation by hydraulic fracturing to make commercial wells in this type of reservoir. Performance predictions for volumetrically controlled reservoirs have been made using a combination of two standard equations. 1. The "Schilthuis" or "Muskat" type material balance equation is used to relate the average reservoir pressure and the cumulative oil recovery. 2. The results from the material balance equation and the productivity factor as described by Pirson' are used to relate the cumulative recovery with producing rate and time. The material balance equation assumed uniform pressure and liquid saturation conditions throughout a reservoir. The steady-state radial flow formula allows for a pressure gradient toward a well but assumes uniform liquid saturation. These calculation methods are adequate for application to reservoirs wherein the drawdown at the well to realize satisfactory producing rates is small compared to the total pressure. In low productivity, volumetrically controlled reservoirs, the pressure drawdown at the well is large cornpared to the total pressure. Although a precise number cannot be given for the magnitude of a large pressure drawdown, values in excess of 1,000 psi would definitely be included. For practical considerations, this usually occurs when the formation flow capacity is less than about 100 md-ft. However, this limit of formation flow capacity will vary with the well producing rate. The low pressure in the neighborhood of the well which results from a high drawdown causes evolution of large volumes of gas. This causes the gas saturation to be higher near the well than at a greater distance—-hence, a non-uniform gas saturation. Also, the relationship between the relative permeability to oil (K/K) and gas saturation is nonlinear but decreases approximately in an exponential way with increases in gas saturation. Because of this, the following chain reaction is established.

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 T. F. Mueller, P. A. Witherspoon

For the case of an infinite radial system operating at constant terminal rate, the reservoir engineer often uses the "point source" solution of the diffiusivity equation to study pressure interference effects. At early times and at short distances from the inner boundary these solutions are invalid. The amount of error is often not precisely defined because of mathematical difficulties. Work presented here shows that, if dimensionless time is defined appropriately, previous solutions of the pressure equation can be displayed .as a family of curves on one chart. These solutions include the point source solution (referred to in the field of hydrology as the Theis solution) and other solutions obtained with digital computer methods. With these curves, an exact evaluation of the pressure drop within a reservoir or an aquifer can be made by the engineer. Examples of field problem solutions are presented. In most reservoirs the error involved when the Theis solution is employed is often negligible; whereas, in the calculation of interference effects in an aquifer, a substantial error can occur through such an approach. INTRODUCTION Flow equations are used in petroleum engineering to study the behavior of individual wells and reservoirs. In the case of wells, the pressure response at the wellbore face is the major point of interest; whereas, in the case of reservoirs, the pressure response at the interface of the aquifer boundary is sought. To aid in such studies, the flow equations have been solved in terms of the behavior at these two inner boundaries. Only limited work has been published in regard to the pressure conditions away from these points, i.e., within the reservoir or aquifer. Theis and Mortada are among the few who have reported on this problem. The Theis approach employs the exponential integral and is valid for pressure conditions that occur some distance away from the flow disturbance. It is derived from the concept of a point source, as opposed to a flow across some finite area. The Mortada results, on the other hand, are valid at all points within the reservoir or aquifer. They are presented in terms of dimension less ratios of the radius where the pressure is desired to the radius where the flow rate is measured. Their main use, in the past, has been in aquifer studies. The published results are presented in the form of graphs that are limited to a maximum radius ratio of 64. These graphical results are cumbersome to interpolate at non-integral radius ratios, so that one may be forced to utilize a rather involved analytical expression presented by Mortada. BASIC EQUATIONS The solutions of Mortada and Theis are both based on the diffusivity equation as applied to the case of an infinite radial system subject to a constant terminal rate. The equation is obtained by combining the material balance equation with Darcy's flow equation. The assumptions implicit in the use of this equation are as follows: (1) a single fluid is present that occupies the entire pore volume; (2) the reservoir is horizontal, homogeneous, uniform in thickness, and of infinite radial extent; (3) compressibility and viscosity of the fluid remain constant at all pressures; and (4) fluid density obeys the equation Using the diffusivity equation in situations where the above conditions do not hold will result in errors. These errors (not discussed here) but only the errors which arise in the solution of the equation itself. The diffusivity equation for the homogeneous reservoir conditions cited above can be written in cylindrical coordinates, as To obtain a dimensionless equation, so a single solution can be used for applications of different porosity, permeability and fluid properties, the following transformations are usually made: After these transformations are made, Eq. 3 .can be written in dimensionless form, as MORTADA SOLUTION One solution of Eq. 6 has been given by Mortada.

Jan 1, 1966

By A. N. Crownover, H. T. Kennedy, E. P. Miesch, C. H. Bowman

The paper presents correlations of (l) molar volume of gaseous hydrocarbon mixtures with pressure, temperature, composition and properties of the C,-plus fraction; (2) shrinkage of oils during flash and differential liberation of gas, including the calculation of formation volume factor under various conditions; and (3) bubble-point pressure with temperature, composition and characteristics of C, plus. The data on which the correlations are based comprise 1,615 measurements on 900 hydrocarbon systems, including numerous systems containing nitrogen, hydrogen sul-fide and carbon dioxide. In each correlation, the number of data points covered and the accuracy is substantially greater than in previously available work. Thus, the equation yielding molar volumes of gases has an average deviation of 2.04 per cent, applied to mixtures having temperatures up to 313F and pressures up to 9,800 psis, compared to 2.37 per cent for the Benedict-Webb-Rubin equation applied to the same data, and 4.53 per cent for the method based on the law of corresponding states. he equations resented are all explicit in the dependent variable, and require no iteration on the digital cow puter. INTRODUCTION The ease and accuracy of determining the composition of hydrocarbon mixtures, compared to the difficulty of measuring their properties under reservoir conditions, makes it desirable to utilize composition as the key to physical behavior to the greatest possible extent. As a result, there are available correlations between composition, or easily measured characteristics dependent on composition, and practically every important engineering property of reservoir fluids. The task confronting us is one of finding more exact relationships between important variables rather than extending correlations to new properties. This paper describes new correlations of molar volumes of gases, formation volume factors, and bubble-point pressures with composition, temperature and where possible, pressure. Each correlation is obtained by employing a sufficiently large amount of data so the calculated properties are probably as least as accurate as the measurement on which they are based. MOLAR VOLUME OF GASES Although many equations of state have been proposed for pure gases, only a few methods are applicable to hydrocarbons at conditions comparable to those in petroleum reservoirs. Still fewer arc useful in describing the behavior of mixtures, with which the petroleum engineer is largely concerned. The correlation presented here involves procedures similar to those of Alani and Kennedy The van der Waals equation is modified to make a and b functions of temperature instead of being constants for each material.This is a cubic equation, which may have either one or three real roots. The lowest root corresponds to liquid volume, while the highest applies to gas. When the procedure was applied to 703 pressure-temperature points of 164 gases of known composition and volume, the average deviation was 12.08 per cent, and the standard deviation 8.15. The above calculations were made using the constants derived for liquids by Alani and Kennedy and the mixture equation developed by them. A closer approximation to measured molar volumes is obtained by employing different sets of constants for different areas on the pressure-temperature chart, and by changing the relationship between a,,, and b,,, for mixtures and the a, and b, for individual hydrocarbons. The various sets are shown in Table 1, and the areas for which they are recommended are plotted in Figs. 1 through 7. In these figures, the area designated by zero is in the critical region for pure materials, and the values obtained for them may be unreliable. For normal mixtures, however, the area may be added to any of the bounding areas without substantial error. The new relationships employed are

Jan 1, 1966

By A. L. Barnes

Residual oil in watered-out reservoirs is a tremendous reserve which has been unrecoverable by established production methods. A study of the new recovery methods indicated that the forward combustion urocess might- recover oil from such reservoirs; however, no thermal recovery operating experience in a watered-out system was available. The Delaware-Childers pilot thermal test was undertaken to test the feasibility of thermal recovery in this watered-out reservoir. The pilot test consisted of a 2.22-acre, inverted five-spot in the 600-ft deep Bartlesville sand. The reservoir in the pilot area had a porosity of 20.6 per cent, a permeability of 118 md and an average sand thickness of 45 ft. The reservoir contains a 33' API, 6-cp oil. Combustion was started Nov. 22, 1960. The initial air injection capacity was 750 Mscf/D, but it was eventually increased to 2,000 Mscf/D. The test was surrounded by an active water flood; therefore, water production was initially high, but decreased as the heat wave moved toward the producing wells. Oil-bank arrival at an individual well was indicated by a drop in GOR and WOR, and an increase in oil production. Combustion-front arrival was evident at three of the pilot producers, and they were plugged. Cumulative oil production from pilot area wells war over 12,000 bbl. Operational difficulties were negligible and only conventional equipment was necessary. The combustion efficiency of this test averaged over 80 per cent. Results from coring showed that the leading edge of the combustion front tended to be wedge-shaped but a nearly complete sweep of the reservoir was eventually obtained. An isopach map based on evidence from 10 core holes and the existing wells showed that 126 acre-ft had been swept by the heat wave. Using this swept volume, an air requirement of 15.7 MMscf/acre-ft was calculated. It was calculated that 275 ,STB/acre-ft was consumed bv the heat wave. INTRODUCTION There is a large amount of oil remaining in reservoirs that have been water flooded. A study of ways to recover this oil showed that the forward combustion process might be applicable. Results from a number of forward combustion tests have been reported in the literature, but none of these tests were conducted in a watered-out system. The Delaware-Childers pilot thermal test was initiated in 1960 to define the operating characteristics of underground combustion in this watered-out Bartlesville sand reservoir, The purpose of this paper is to present the pilot test results in detail. FORWARD COMBUSTION PROCESS The forward combustion process consists of initiating combustion in the formation surrounding an injection well and driving this heat wave through the formation toward offset producing wells. As the combustion front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase, and recon-densed in the cooler part of the formation. These distilled liquids, water of combustion and gaseous combustion products, form a bank or three-phase region ahead of the burning front. This bank pushes mobile reservoir fluids toward the production wells. The rate of movement of the combustion front is controlled by the rate at which the nondistill able residue which serves as process fuel can be completely burned off the sand. The production performance from a heat wave conducted in a watered-out oil reservoir will differ from one conducted in a dissolved gas-depleted reservoir because of the difference in fluid saturations. The primary depleted reservoir contains connate water saturation, relatively high oil saturation and some gas saturation. The watered-out reservoir contains a highly mobile water saturation, residual oil saturation and little gas saturation. As the combustion front moves forward in either type of depleted reservoir, it drives the distillable oil, the water in place and the water of combustion ahead, and burns the nondistillable portion of the oil. The heat-wave process in a watered-out reservoir differs from the process in a

Jan 1, 1966

By L. A. Wilson, N. H. Harrison, R. R. Clay, R. L. Reed. D. W. Reed

Experiments have been performed in a linear near-adiabatic system for the purpose of extending data on reverse and forward combustion from atmospheric pressure to 1,000 psig. Results obtained from reverse combustion appear to conform qualitatively with the existing description of the process when reaction kinetics are suitably modified to account for the increase in pressure, i.e., increasing tbe pressure decreases the peak temperature and increases the combustion zone velocity. Forward combustion appears to be a fuel dominated process wherein peak temperature and combustion zone velocity are not very sensitive to changes in pressure. The moderate effects of pressure that do exist at low flux virtually disappear at high flux providing all oxygen is consumed. With this provision. increasing the pressure decreases the frontal velocity and increases the peak temperature. Results me shown graphicaIIy which demonstrate the effects of presswe on peak temperature, rate of advance, oil recovery, air/oil ratio, carbon oxides produced and temperature distributions. INTRODUCTION Numerous field tests of oil recovery using the technique of underground combustion are now in progress. Operating pressures used are always substantially in excess of atmospheric pressure. Nevertheless, the literature contains no laboratory data pertaining to the effects of pressure except for those of Martin, Alexander and Dew.1 Unfortunately, these data appear to reflect heat losses which may have obscured some of the effects of increased pressure. The underground combustion processes are exceedingly complex, and general concepts relating to them have necessarily been described in relatively simple terms. Particularly with regard to forward combustion, the interaction between multiphase fluid flow, heat transport and various rate mechanisms is so involved that to date a rigorous treatment is not available. In this paper a linear near-adiabatic system is described which was found capable of operating at pressures up to 1,000 psig and temperatures of at least 1,100 F. Using this equipment, experiments on both forward and reverse combustion were accomplished at elevated pressures and the data presented as a function of air flux or temperature. It will be assumed that the reader is familiar with the papers of Wilson, Wygal, Reed, Gergins and Henderson on forward combustion and that of Reed, Reed and Tracht3 on reverse combustion, so that it is not necessary to repeat the details of the two processes. EQUIPMENT AND PROCEDURE It was decided to study the effects of pressure over the interval 0-1,000 psig using a near-adiabatic system. Since satisfactory operation of such a system requires very thin tubing walls, it proved necessary to insert the entire combustion tube and heater assembly into a high pressure jacket in such a way that the internal tube pressure could be counterbalanced by pressurizing the annulus. This meant that the numerous heater and thermocouple wires had to pass through gas-tight seals that required considerable effort to maintain effective. An unexpected difficulty arose in connection with the nature of the annulus pressuring gas. It turned out that the composition of this gas critically affected the functioning of the adiabatic control system. Operation of the equipment at elevated pressures was not particularly difficult. The problems that did arise were a direct consequence of the large numbers of connections and seals involved. EQUIPMENT The equipment consisted of the combustion tube which fitted within the pressure jacket. This jacket was in turn mounted on a wheeled carriage which facilitated movement and provided for rotation of the

By P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

By D. H. Gray, G. Bergamini, I. Fatt

A method developed by soil scientists for measuring permeability anisotropy of soils has been adapted for use in studies of sandstone permeability aniso-tmpy. Studies on Berea and Bandera sandstone show that permeability of these materials decreases upon application of simulated overburden pressure. Permeability reduction of Berea and Grubb sandstone cores subjected to simulated overburden pressure is shown to be a function of the ratio of radial to axial stress. INTRODUCTION The effects of mechanical stress on reservoir rock parameters are now generally accepted as being important in reservoir engineering calculations. In several oil fields the unique behavior of the reservoir has been attributed to the influence of overburden pressure or internal fluid pressure on rock permeability or fluid content. l Although the gross effects of overburden pressure on permeability, porosity, waterflood efficiency, electrical resistance and sonic velocity have now been established,3-9 there are some very interesting, and perhaps important, details of rock behavior under stress that have not heretofore been studied. In this paper we will show that permeability anistropy of several sandstones is a function of overburden pressure. We will also show that permeability reduction of cylindrical samples of two sandstones, when subjected to mechanical stress, is a function of the ratio of radial to axial stress. Most studies of overburden pressure effects made in the past have attempted to duplicate reservoir conditions by applying hydrostatic external pressure on all surfaces of the sample. Assuming perfect elasticity, that is, no viscoelastic flow during geological time, reservoir rock in situ cannot be under hydrostatic stress as is apparent from the value of 0.25 for Poisson's ratio for most rock. In the absence of tectonic forces the horizontal stress component If there are directional permeability effects, as has been reported from field and laboratory studies, l0,11 then one would expect hydrostatic stress as used in the laboratory to give a different permeability reduction from that occurring under reservoir stress conditions. The influence of stress on directional permeability was measured by a modification of the method first described by Maasland and Kirkham for use on soil samples.12 Their method permits calculation of the permeabilities along each of the three mutually perpendicular axes of a cubic sample. The modification reported in this paper extends the method to observing these permeabilities when the sample is under hydrostatic stress. The Hassler sleeve type of core holder is satisfactory for studying overburden pressure effects on permeability if a known force is applied to the end plugs. Wyble6 used a sleeve core holder but did not load the end plugs. In this paper we combine sleeve pressure and axial force to give permeability changes which we believe are close to those occurring in reservoirs when reservoir pressure changes. The combination of sleeve pressure and axial force assumes, however, that the sleeve pressure simulates both vertical stress and horizontal stress in one direction. METHOD PERMEABILITY ANISOTROPY The Maasland and Kirkhaml2 method of measuring permeability anisotropy uses a cube of porous material which is Open On three faces and has discs with a Center hole for an air inlet cemented to the other three mutually perpendicular faces (Fig. 1). The rate at which air pressure falls in a bottle of known volume, and pressure slightly above atmospheric, when the bottle discharges through each disc, is

By D. Havlena, A. S. Odeh

The material balance equation used by reservoir engineers is arranged algebraically, resulting in an equation of a straight line. The straight line method of analysis imposes an additional necessary condition that a successful solution of the material balance equatiott should meet. In addition, this algebraic arrangement attaches a dynamic ineuning to the otherwise static material balance equation. The straight line method requires the plotting of one variable group vs mother variable group. The sequence of the plotted points as well as the general shape of the resulting plot is of utmost importance. Therefore, one cannot progrm the method entirely on a digital computer ar is usually done in the routine solution of the material balance equation. If this method is applied, then plotting and anaIysis are asential. Only the appropriate equations and the method of analysis and interpremtion with comments and discussion are presented in this paper. Illustrative field examples for the various cases treated are deferred to a subsequent writing. INTRODUCTION One of the fundamental principles utilized in engineering work is the law of conservation of matter. The application of this principle to hydrocarbon reservoirs for the purpose of quantitative deductions and prediction is termed "the material balance method of reservoir analysis". While the construction of the material balance equation (MBE) and the computations that go with its application are not difficult tasks, the criteria that a successful solution of the MBE should fulfill have always been a problem facing the reservoir engineer. True and complete criteria should embody necessary and dcient conditions. The criteria which the reservoir engineer uses possess a few necessary but no sufficient conditions. Because of this, the answers obtained from the MBB are always open to question. However, the degree of their acceptability should increase with the increase in the number of the necessary conditions that they should satisfy. Generally, the necessary conditions commonly used are (1) an unspecified consistency of the results and (2) the agreement between the MBE results and those determined volumetrically. This second criterion is usually overemphasized. Actually, the volumetrically determined results are based on geological and petrophysical data of unknown accuracy. In addition, the oil-in-place obtained by the MBE is that oil which contributes to the pressure-production history,' while the volumetrically calculated oil-in-place refers to the total oil, part of which may not contribute to said history. Because of this difference, the disagreement between the two answers might be of paramount importance, and the concordance between them should not be overemphasized as the measure of correctness of either one. In this paper, a third necessary condition of mathematical as well as physical significance is discussed. It is not subject to any geological or petrophysical interpretation, and as such, it is probably the most important necessary condition. It consists essentially of rearranging the MBE to result in an equation of a straight line. This straight line method of the MBE solution has invalidated a few long time accepted concepts. For instance, it has always been advocated that if a water drive exists, but one neglects to take it into account in the MBE, the calculated oil-in-place should increase with time. The straight line method shows that in some cases, depending on the size of the neglected aquifer, the calculated oil-in-place might decrease with time. The straight line method requires the plotting of a variable group vs another variable group, with the variable group selection depending on the mechanism of production under which the reservoir is producing. The most important aspect of this method of solution is that it attaches a significance to the sequence of the plotted points, the direction in which they plot, and to the shape of the resulting plot. Thus, a dynamic meaning has been introduced into the picture in arriving at the final answer. Since the emphasis of this method is placed on the interpretation of the sequence of the points and the shape of the plot, one cannot completely automate the whole sequence to obtain "the best value" as normally done in the routine application of the MBE. If one uses the straight line method, then plotting and analysis are musts. The straight line method was first recognized by van Everdingen, et al,2 but for some reason it was never fully exploited. The advantages and the elegance of this method can be more appreciated after a few cases are carefully treated and worked out by it. SOLUTION OF THE MATERIAL BALANCE EQUATION SATURATED RESERVOIRS The MBE for saturated reservoirs written in AIME symbols is

By C. Chu

An investigation has been made of the radial heat-wave process using a mathematical model in two-dimensional cylindrical coordinates. This model considers combustion, convection and conduction inside the reservoir, but only conduction in the bounding formations. From a study of the general features of the process, an important phenomenon has been revealed, namely, the feedback of heat into the reservoir on the trailing edge of the heat wave. The effects of various process variables on the performance characteristics of the process have also been investigated. It was found that up to the time when the combustion front reaches a given point, the per cent heat loss, provided it is not higher than 40 per cent, is approximately directly proportional to the square root of thermal conductivity arid fuel content, but inversely proportional to the square root of gas-injection rate and oxygen concentration. The effecr of reservoir thickness is more pronounced, since halving the thickness doubles the per cent heat loss. The most decisive factor in determining the center-plane peak temperature is the fuel content of the reservoir. Within the temperature range investigated, doubling the fuel content doubles the peak temperature in the early stage, but the rate of decline of the peak temperature is high. Reservoir thickness is also a very influential factor. The peak temperature is lowered when the thickness is reduced; however, the effect of thickness becomes less pronounced when the thickness is high. Reduction of oxygen concentration increases the peak temperature in the early stage but lowers it afterwards because of the higher rate of decline of the peak temperature. Increase in gas injection rate or decredse in thermal conductivity geives a higher peak temperature which stays high for a longer period. The propagation range of the heat wave is chiefly governed by the fuel content of the reservoir. An increase of 0.2 1b/cu ft in the fuel content increases the propagation range by 100 per cent. The propagation range is more than doubled by doubling the gas injection rate, or reservoir thickness, or by reducing the thermal conductivity by 50 per cent. Comparatively, oxygen concentration has less effect on the propagation range. INTRODUCTION Several investigators have conducted theoretical studies of a radial heat wave. Vogel and Krueger1 studied a system with a moving cylindrical heat source of constant temper- ature, considering conduction in the radial direction only. Ramey2 included conduction in the vertical direction in his studies. Bailey and Larkin2 attacked a more general problem where initial well heating, vertical heat losses and arbitrary frontal velocities were included. In all these studies, however, conduction was considered to be the only means of heat transfer. Bailey and Larkin in a later paper included the effects of convection in a study involving both linear and radial geometries. Vertical heat losses were neglected in the radial case. Katz5 studied a similar problem in a one-dimensional radial model, using a heat-loss coefficient to account for vertical heat losses. Selig and Couch6 mployed a cylindrical model and investigated two limiting cases. In one case they considered no heat loss from the reservoir whereas in the other they assumed a constant temperature at the interface between the reservoir and its bounding formations. Thomas' studied a more general case but assumed a permeable bounding formation so that the convection effect is not confined to the reservoir. In the present work a more realistic and more generalized model is used. It involves a two-dimensional cylindrical system with combustion, convection and conduction inside the reservoir, but only conduction in the bounding formations. The purpose is to establish the temperature distribution both inside and outside the reservoir, to study the general features of the radial heat wave process, and to investigate the effects of various process variables on the performance characteristics of the process. THEORY We fist consider a circular porous reservoir of thickness H extending vertically from y = — H/2 to y = + H/2. The reservoir extends from a well bore radius r, to an external radius re. A stream of oxygen-containing gas is introduced into the reservoir through the wellbore. The oxygen-containing gas reacts with the fuel contained in the reservoir and forms a combustion front wherever the prevailing temperature can support the combustion. It is here assumed that this combustion front constitutes a cylindrical surface source of heat having an infinitesimal thickness in the radial direction and extending vertically throughout the whole thickness H. This is designated as Region I. We next consider a Region II corresponding to the upper and lower formations bounding the reservoir, extending from y = — - to y = — H/2 and from y = + H/2 to y = + m. Since r, is very small, we may assume that the two bounding formations have the same dimensions, symmetric with respect to the center plane of the reservoir. In this way, we may take the upper half of the system alone into consideration. In contrast with

By V. J. Driscoll

A method is given for checking pore volume reserve estimates and transient water inflzix factors utilizing early well performance. Principal applications are based on the observation that there is a convenient transient time-range in terms of well spacing in which observed pressure drop will give a good indication of fluid-in-place. Primary quantity determined is a porosity-feet value applicable in the vicinity of an observation well. Checks on porosity or net pay thickness and permeability values also are obtained. These are particularly valuable in improving quality of interpretive factors entering into pore-volume reserve determinations for larger fields. In smaller fields, direct determination of total field reserves and approximate location of boundaries may be made in some cases. Well interference and production data plus a build-up curve on the observation well are needed. The build-up curve may be taken at any time. Field applications are given. The method can be applied with the completion of the second well on a new discovery or to older reservoirs by major changes in withdrawal rates from a local area. Appropriate test length permits results to be obtained with equal accuracy on widely spaced wells. Criteria for optimizing test conditions are given. Provision is made for various boundary and well effects including incomplete observation well build-up. INTRODUCTION Due to low quality data, initial pore volume reserve and water influx estimates may often vary by a factor of two or more from results obtained at a later date. Unfortunately, by conventional methods, more reliable answers from performance usually are not available until after most significant economic decisions have been made. Furthermore, in older reservoirs, uncertainties may still exist as to original oil-in-place, recovery factors and residual oil saturations due to incomplete or inaccurate performance data. This paper outlines situations and appropriate arrangements of basic transient theory whereby well interference data can be utilized to improve the quality of estimates under certain conditions. In addition to production data on offset producers, key information required is a measured pressure drop at an observation well, corrected for any incomplete build-up due to previous production, plus a formation capacity value derived from a conventional pressure build-up plot on the observation well. The method, as outlined, is strictly applicable only to a horizontal single-layered system of uniform rock characteristics containing a homogeneous fluid of low and constant compressibility. However, under the set of test conditions devised, requirements can usually be relaxed so that it should be applicable to most low permeability gas, water and undersaturated oil reservoirs. This is an area where more reliable information is certainly needed. BASIC CONSIDERATIONS A, B, C and F as subscripts or in the text refer to producing Well A, observation Well B, a more remote producing Well C and field boundary F. Combinations such as AB refer to distances and dimensionless times in terms of distances between these points, while single letters refer to distances and dimensionless times in terms of the well-bore or field radius. Equations covering definitions of quantities and complex situations are given in the Appendix. Key equations in practical units are given in Table 2. Transient Pressure Behavior Fig. 1 illustrates pressure drop at Point B with time due to initiation of a constant rate of withdrawal from Point A. This is shown as since a quantitative under-

By M. A. Jones

Performance of a pilot flood in eastern Kansas indicates improved recovery and accelerated production resulting from mobility ratio control obtained by adding a high moleculur weight polymer to injected water. Ultimate pilot flood recovery of the 75-cp reservoir oil is anticipated to he 350 bbllacre-ft during the six-year flood life. Cumulative recovery to date totals 242 bbl/acre-ft after 30 months of injection. Reservoir parameters and flood performance are disc~rssed in detail to provide a case history relating to the feasibility of polymer flooding. INTRODUCTION The Vernon Polymer flood of the Brazos Oil & Gas Co., located in Sec. 2, 24 S., 16 E., Woodson County, Kans., was begun in late Oct., 1963. Injection of a 500 ppm solution of Pusher* polymer was initiated in a 15-acre pilot to determine feasibility of improving oil recovery by mobility ratio control. Performance to date and the anticipated recovery of 115,000 bbl from the 329 acre-ft pilot indicate efficient displacement of the viscous reservoir oil. Dependence of flood recovery on the displacing fluid-displaced oil mobility ratio has long been recognized. Displacement efficiency of the reservoir oil decreases with increasing mobility ratio. Therefore, economic reduction of this parameter is an attractive approach to improving oil recovery. particularly from viscous oil reservoirs where mobility ratios are often quite high. Waterflood mobility ratios can be controlled by dilute solutions of certain polymers'; which provide a displacing fluid with considerably reduced mobility. Reservoir parameters and flood performance are emphasized in this paper to provide data relating to polymer flood applicability. THEORY The polymer used in the Vernon pilot exhibits the property of unusually high resistance to flow through porous media.'.' A small concentration in injected water effects a mobility reduction considerably beyond that predicted from increases in bulk viscosity as measured in a capillary or rotating cup viscosimeter. Water mobility in the subject flood is reduced by a factor of 7.55 with the addition of 500 ppm Pusher chemical to injected water while the bulk viscosity is increased from 1.03 to only 1.40 cp at 75F. The exact flow mechanism causing resistance to flow has not been established; however, it is observed only in tortuous flow systems and appears to be influenced by brine properties, rock properties and oil saturation. Consequently, mobility reduction is estimated from flow studies for a given sample of reservoir rock and is expressed in terms of a "resistance factor".' This parameter is defined as the ratio of brine mobility to polymer solution mobility measured at residual oil saturation: RF= ?m/?p= (km/µm)/(kp/µp)M.......(1) where RF = resistance factor (dimensionless) k*. and A, = mobility of water and polymer solution, respectively k,,, and k,, = water and polymer permeability, respectively (darcies) µie and µp = viscosity of water and polymer, respectively (poise). Eq. 1 is similar to the expression commonly used to relate the relative mobility of the water bank at residual oil saturation to oil-bank mobility at irreducible water saturation: where M =mobility ratio of water to oil (dimension-less) A, = mobility of oil k, = permeability to oil (darcies) )i* = oil viscosity (poise). From Eqs. 1 and 2 it is apparent that the polymer-solution-to-oil-mobility ratio can be approximated by dividing the brine-to-oil ratio by the resistance factor RF: Mpn=M10/RF...........(3) where M,-, = mobility ratio of polymer solution to oil (dimensionless). This simplified presentation assumes no dilution of the polymer solution by connate water. Also, the presence of a connate water bank separating the polymer solution from the displaced oil is not considered. Displacement of oil by waterflooding is increasingly dominated by viscous forces as the viscosity of the oil displaced increases. When the oil mobility is less than that of the injected water, the oil-water interface is not a piston-like front.3,' Instead, instabilities occur at the flood front in the form of viscous fingers which tend to penetrate the less mobile oil bank. The magnitude of these protrusions, or the degree of channeling, increases

Jan 1, 1967

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 R. J. Wagner, F. F. Craig, H. G. Riley, J. D. Griffith

Peripheral water injection has been underway in the Sholem Alechem Fault Block "A" Unit, Stephens County, Okla., since 1955. In the engineering planning of the flood, it was recognized that maintenance of pressure in the large, dry gas caps was vitally necessary to prevent the loss of an estimated 27 million bbl of oil through saturation of the gas caps by the oil banks moving updip on the peripheral flood pattern. The development of an abundant, low-cost water supply and the rising cost of makeup-gas purchases indicated that gas-cap control by water injection was desirable from an economic standpoint, provided an eflective water barrier could be placed without premature water breakthrough into downdip oil producers. A laboratory model scaled to the existing reservoir conditions was comtructed to evaluate the feasibility of gas-cap water injection. Model studies indicated that (I) water injected at the edge of the gar-oil contact would form a continuous barrier between the oil and gas zones, and (2) water would advance uniformly downdip without segregation due to gravity effects provided a specified, minimum water-injection rate was maintained. The model test results also indicated that it was not necessary to fill the gas cap completely with water. Water injection at the edge of one of the gas caps was initiated in March, 1958, to field test gas-cap water injection. Water breakthrough occurred at a downstruc-ture producing well in April, 1959, after approximately 2.5 million bbl of water had been injected into the gas cap. This was the volume calculated to be required to form the barrier between the oil zone and gas cap. This performance indicates that prevention of oil encroachment into dry gas caps in a peripheral water flood by water injection at the gas-oil contact appears practical and that the model studies provided a valid prediction of field performance. INTRODUCTION The Sholem Alechem Fault Block "A" Sims Sand Unit was one of the fist major waterflood projects in Oklahoma to utilize a peripheral injection pattern. Design of control measures to minimize oil encroachment into the large crestal gas caps was one of the greatest single problems during engineering planning of this flood. Engineering studies indicated that saturation of the gas caps with oil displaced by the peripheral water flood could result in a loss of about 27 million bbl of oil, or over one-half the estimated waterflood recovery. This estimate was based on a residual oil saturation after water flooding of 35 per cent pore volume, compared to the negligible oil saturation initially in the gas caps. The Sholem Alechem Springer field, located in Stephens County, Okla., is an elongated anticline with structural dips varying from 400 to 1,000 ft/mile. The structure contains numerous faults, some of which are sealing. The portion of the field considered in this paper is Fault Block "A", which is the largest of the five major fault segments in the field. A map of the unit is shown on Fig. 1. The unit is bounded on the east and west by major sealing faults, on the north by a water-oil contact and to the south by a deterioration of sand development. Two domes, each having gas caps, exist in the unit and are separated by a north-south fault as shown in Fig. 1. The sands being flooded are the Sims sand members of the Springer series of early Pennsylvanian age. The Sims sands average over 100 ft in thickness in Fault Block "A" and have a porosity of approximately 19 per cent. Core permeabilities average 160 md, with interstitial water saturation of about 23 per cent pore volume. Three different zones can be correlated over most of the area. However, they are commingled in most of the producing wells. The two upper zones, designated as First and Second Sims, are the most important. Both the First and Second Sims have gas caps in the east dome, while only the First Sims has a gas cap in the west dome. These gas caps, whose limits are also shown on Fig. 1, initially occupied nearly 11 per cent of the total reservoir volume. The major concern of this paper is the control of the two gas caps in the First Sims sand. Fig. 2 is a cross section, to scale, across the east dome of this formation. Notable features are the relatively shallow dip, the deterioration of the First Sims sand on the south side of the gas cap and the water-oil contact on the northern edge. The formation dip at the gas-oil contact in this cross section is approximately 5".

By R. F. Nielsen, R. G. Bentsen

Jan 1, 1966

By R. W. Heins, N. Street

An investigation was made of the Rehbinder-Kuznetsov pendulum as an instrument to aid in the study of rock drillability and of the effect of surface active agents on rock drillability. The "critical weight", i.e., the pendulum weight at greatest sensitivity, was determined for several solids and it was shown that a change in surface environment changed the maximum value of the apparent hardness but not the critical weight, and the critical weight increased with increase of hardness of the solid used. Using silver halides as the solid material, and solutions of either silver nitrate or the appropriate alkali halide, it was shown that there is a hardness maximum at the zero point of charge of the particular silver halide. INTRODUCTION The development of relative hardness measurement by the method of damped oscillations, or what has been called the pendulum sclerometer, is due largely to the Russian workers P. A. Rehbinder and V. D. Kuznetsov and their co-workers.1,2 The apparent simplicity of the method has led to its use in a number of fields, particularly those where relative hardness has to be determined on a brittle solid. As will be shown, the simplicity of the method belies its complexities. THEORY OF THE PENDULUM SCLEROMETER The operation of the pendulum is quite simple. Make a sharp point such as a diamond, hardened steel or tungsten carbide, as the fulcrum of a pendulum, with the center of gravity of the system below the point of support. The attenuation of this pendulum is a function only of air friction, provided no fragmentation, breakage or other work is performed by the point. However, when the solid surface breaks or disintegrates under the point, and the latter penetrates the solid, part of the energy of oscillation is absorbed in the work of creating new surface, and the rate of attenuation increases. In the case of brittle solids, the point of the pendulum produces fragmentary particles, whereas the point produces plastic deformation in a plastic solid; but in both cases damping occurs. This pendulum damping is a function of the hardness or strength of the solid. Whereas a pendulum oscillating on a hard material attenuates primarily by air friction, in the case of a soft material, the points penetrate the surface of the solid, producing new surface, and the pendulum damps much faster. Thus the attenuation rate is a function of the hardness or strength of a material in both plastic and brittle substances. Russian workers initially used a single support for the pendulum, but this created many problems because of oscillation of the pendulum in more than one plane. A revised version used two points, or sometimes a prism as support, successfully confining oscillation to one plane. The version using two points was constructed and used throughout this investigation, making the recording of the attenuation as a function of time relatively easy. That both points should be in contact with the solid during the damping is implied. Rehbinder (c.f. Kuznetsov2) proposed that the hardness H as measured by this method could be expressed by a quantity which is inversely proportional to the relative initial rate of attenuation of amplitude. where H = hardness parameter (arbitrary units) a, = initial amplitude da/dt = rate of attenuation of amplitude with time. In the simplest case of damping with a constant logarithmic decrement, H can be regarded as a constant for the duration of the whole damped oscillation process. That is

Jan 1, 1966

By Caudle B. H., Cooper H. E. Jr., O. K. Kimbler

The results of studies of areal sweepout behavior in nine-spot injection patterns are presented in this paper. The ratio of the production rates of the two types of producing wells in a nine-spot pattern is considered in its effect on breakthrough. A model study was used to develop a method to predict the production histories of both side and corner wells and to develop the production history of the pattern or the field. INTRODUCTION The importance of areal sweepout behavior in injection operations has been recognized for a number of years. We now know that the fraction of a reservoir which will be reached by an injected fluid is a function of the geometric arrangement of the production and injection wells and of the mobility ratio of the fluids. Of these two, the mobility ratio (which is a measure of the relative ease with which the oil and the injected fluid can move through the reservoir) appears to have the greatest effect on recovery. The effect of mobility ratio on production history has been reported for most of the dispersed injection patterns (i. e., five-spot, staggered and straight-line drives).4-6 ere has been very little published on the nine-spot injection pattern. This may be because the two types of producing wells in the nine-spot pattern can—and sometimes should —be produced at different rates. This adds a second variable (the fist being mobility ratio) to be covered in a study of the nine-spot flooding network. Muskat1 has reported the sweep pattern efficiency at breakthrough for a mobility ratio of one at various ratios of producing rates for the two types of wells. For these conditions he found maximum oil production at breakthrough to occur when the corner well (the producing well farthest from the injection well) is produced 10.7 times as fast as the side well. Other investigatorsa2,3 have studied a few of the many possible operating conditions possible with this injection pattern. This paper presents the results of studies which are strictly applicable to a nine-spot injection pattern surrounded by other nine-spots in a homogeneous reservoir and which has no flowable gas present (no oil bank is formed). Although the full ranges of the variables mentioned above have not been covered, the results reported here should allow reasonable estimates of field operations to be made. EXPERIMENTAL The well spacing geometry of the nine-spot pattern is shown in Fig. 1. Here we see that each injection well is surrounded by eight production wells arranged in a square. Four of these producing wells are located at the corners of the square; the other four are in the midpoints of the four sides. These two types of producing wells are called "comer wells" and "side wells" respectively. It is the ratio of the producing rates of these two types of wells that affects the breakthrough recovery, as calculated by Muskat.' This ratio is defined as the total producing rate of a corner well divided by the total producing rate of a side well. The model used in this study represented one-eighth of a nine-spot pattern, as shown in Fig. 1. It contained the equivalent of one-fourth of a side well, one-eighth of a comer well, and one-eighth of an injection well. Because of these proportions, the actual producing rate ratios in the model were one-half of the ratios of producing rates as

Jan 1, 1965

By C. R. Johnson, R. A. Greenkorn, L. K. Shallenbarger

This paper describes a study, based on core data, of the directional permeability of a sandstone reservoir. Directional air permeabilities are explained and correlated with lithology by the tensor theory of permeability, which is extended to the more general case of heterogeneous anisotzopic porous media. In this case, the permeability tensor is made up of two components: (I) anisotropy (variation around a point) which correlates with bedding, and (21 heterogeneity (variation from point-to-point) which correlates with grain size. INTRODUCTION Previous studies of directional permeabilityl.5-8 have concerned themselves only with anisotropy at a point, rather than with both point-to-point differences and local anisotropies. In this paper, we present the existing tensor theory of permeability for anisotropic porous media and then extend it to the more general case of heterogeneous anisotropic porous media. The laboratory measurements of directional permeability are explained in view of this extended theory and the data are discussed in terms of the lithologic factors that correlate with it. The results show what the heterogeneity and anisotropy of the reservoir element are, and that the directional permeability correlates with lithology. Data used in this paper are measurements of air permeability in eight directions, spaced at 45° intervals, on 142 2-in. vertical plugs from 30 cores. The core data are meaningful in terms of directional permeability because the cores were oriented to within 45O during drilling and coring.3 Although we initially thought that only 60 per cent of the core material was reliable, subsequent study showed that almost all of it was reliably oriented. After determining the air permeabilities, the data were reduced to three independent variables for each plug: the major and minor permeability axes, and the direction of the major axis. These were obtained by converting the permeability data to the reciprocal square root of permeability and fitting the transformed data with ellipses according to the tensor theory of permeability.8 The point-to-point areal variation of the minimum permeability axis is related to grain size. The direction of the permeability axes, where the major axis exceeds the minor axis by at least 5 per cent (measured variation is about 4 per cent), correlates with the bedding. The permeability tensor used in this study must be considered as the sum of a scalar and a tensor, with the scalar being the minor axis permeability as a function of position, and the tensor the directional effect, which is additive permeability over the minor axis. In this case, the point-to-point variation or heterogeneity (minor axis) is substantially larger than the variation at a point or anisotropy. It may be important that this separation of "directional permeabilities" be recognized when considering migration of fluids due to permeability variation. Local migration may be due to anisotropy and point-to-point migration may be due to heterogeneity, but the direction and magnitude of these may not be the same. Furthermore, in truly heterogeneous systems, one would expect that anisotropy would be the smaller of the two effects. THEORY TENSOR THEORY OF PERMEABILITY FOR ANISOTROPIC POROUS MEDIA8 Darcy's law for flow in porous media in its usual form is where qis the flow rate vector, k is the permeability, is the viscosity, is the vector differential operator, and P is the Pressure. Consider flow, q, through a unit area in the direction of the unit vector n from elementary flow tubes parallel

Jan 1, 1965