Search Documents

By R. A. Fitch, J. D. Griffith

A performance calculation method was used in conjunction with experimental studies to develop means of predicting and interpreting miscible floods and to explore possible methods. of improving their efficiency. The calculation method is based upon the division of reservoirs into multiple independent zones or strata. Comparison was made with some experimental floods and with a field project to obtain some measure of the applicability of the method. Two aspects of miscible flooding were considered: (I) the displacement of a miscible front through a reservoir, and (2) the distribution and utilization of the solvent volume injected to maintain miscibility. Calculations and experihental observations indicate that alternate gas-water injection behind a miscible front significantly improves miscible flood performance, both within a single stratum and in a multistrata reservoir. Calculations were made on the extent to which miscibiliry is maintained by a given solvent volume in a stratified reservoir. Two alternate criteria for determining the optimum volume of solvent for injection are discussed. The preinjection of a small volume of water ahead of the solvent is suggested as a method of obtaining more efficient utilization of solvent. Calculations and experiments were made to investigate the effects of water preinjection. INTRODUCTION Miscible fluid displacement as a possible means of recovering oil has been the subject of extensive research and numerous field tests over the past several years. This work has brought out both favorable and unfavorable features of the miscible recovery methods. Many tests have shown that a miscible fluid can displace all of the oil contacted, leaving no high residual oil saturation, as is characteristic of immiscible displacements. The miscible methods have proven dificult to control, however, as evidenced by early breakthrough of the displacing fluid and poor sweep efficiency in several field projects. It became apparent early in the investigation of miscible displacement that the development of improved techniques and better methods of control would be required to extend the range of applicability to any significant portion of our petroleum reservoirs. The purpose of this study was to consider the problem of designing miscible floods to obtain better performance. A performance calculation was used in conjunction with experimental studies to develop means of predicting and interpreting miscible floods and to explore possible methods of improving their efficiency. Two aspects of miscible floods were considered separately in this study and the results are presented in two parts. The first part concerns the manner in which the miscible front is displaced through the reservoir. In these calculations the condition of miscibility between the reservoir oil and the displacing fluid is assumed. In the second portion the distribution and utilization of the preinjected solvent in maintaining rniscibility is considered. This portion dealing with solvent utilization applies only to the enriched gas or LPG slug processes. CALCULATION METHOD AND EXPERIMENTAL STUDIES COMPUTATIONAL MODEL The calculations for expected reservoir performance are based upon the multi-strata concept of reservoir properties. The method assumes multiple, independent strata (no crossflow between strata) and explicitly includes variations in patterns, injectivity, areal sweep and displacement efficiency. Lateral variation in permeability within the strata (for example, directional permeability) may be considered in the calculations, provided the variations are similar in all strata. The calculations based on this model furnish the rates of injection and rates of production of both oil and the displacing fluid for the composite system of strata. To carry out calculations on the mathematical model described, information is required on the performance of the displacement process under consideration in a single stratum element of the pattern to be used. In particular, for a given pattern and displacement process, data are required relating injectivity and composition of the produced stream to the volume of fluid injected. Data on the variations in areal sweep must also be available if it is desired to track the areal sweeps in the multi-strata model. These data may be obtained from experimental work such as model studies or other suitable methods. This representation of reservoir properties is, of course, an approximation and may be inappropriate in some cases. But the representation is not an arbitrary one since most oil reservoirs are essentially successive layers of sediment with the net producing pay zone often comprising but a fraction of the gross formation thickness. The fact that permeabilities measured in the vertical direction are frequently but a fraction of the horizontal permeability adds some evidence to the validity of this model. This type cal-

Jan 1, 1965

By J. F. Wilson

An experimental investigation has been made of gas-driven slug displacements in a system of high gas saturation to evaluate the process for use in a California reservoir. Fluid compositions, temperature, pressure and core permeability duplicated reservoir conditions as closely as possible. The temperature was above the critical of the slug materials. The slug composition required for 100 per cent oil recovery in the linear flow systems was found to agree with that determined from equilibrium phase studies. A qualitative theory for the slug displacement of two-phase systems is proposed and found to be in general agreement with the observed flow behavior. The theory predicts that the length of the reservoir fluid-slug transition Zone will increase with increased initial gas saturation and decrease with increased oil swelling accompanying the composition change from reservoir oil to the critical composition. The precipitation of small amounts of an asphaltic phase during the formation of the transition zone had no adverse effects on the flow behavior of otherwise miscible systems. INTRODUCTION Despite the large number of recent publications relating to oil recovery by miscible displacement, little information is available concerning the behavior of crude oil-natural gas systems at conditions likely to be encountered in field application. The laboratory data reported here were obtained during an investigation of the possibility of using a gas-driven LPG slug in a partially depleted California reservoir. Seven gas-driven slug displacements employing slugs of different compositions and sizes were made in 26- to 27-ft long sandstone cores. The equilibrium phase behavior of the systems was also determined. The features of primary interest in this study are the following. 1. Data were obtained on the displacement of reservoir fluids at reservoir conditions. 2. A high gas saturation existed in the long cores at the start of the displacements. 3. The temperature of the studies was above the critical of the slug materials. 4. A comparison between the flow behavior and equilibrium phase behavior, as illustrated by a pseudo-ternary phase diagram, is afforded for a system of considerable complexity. THEORY Phase Behavior The use of the pseudo-ternary phase diagram to illustrate the phase behavior in miscible displacements has been discussed by Clark, Schultz and Shearin. Sie-vert, Dew and Conley and Arnold, Stone and Luffel' have presented both phase equilibrium and flow data for true ternary systems. For the true ternary, the application of the theory is straight-forward; for a reservoir system represented as a ternary, the theory undoubtedly gives a qualitatively correct picture. One purpose of this investigation was to see if the pseudo-ternary phase diagram could be used to predict with reasonable accuracy the conditions necessary for miscible displacements with actual reservoir systems. MiscIBLE Displacement of Two-Phase System A large number of papers "I3 have treated the miscible displacement of a single-phase fluid from both an experimental and a theoretical viewpoint. However, the normal field application of the process will be in a partially depleted reservoir. Often the state of the system will be such that the mobility of the gas phase is greater

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

By D. M. Kehn

A method has been devloped for chromatographic analysis Of the vapor and liquid phases Of a a system containing methane to components having 20 or more carbon atoms. The method uses a windowed equilibrium cell in which volumetric phase behavior of the system can be observed accurately and from which small samples of gas or liquid can be withdrawn for analysis. Analyses are made using two chromatographs, one for the lighter and one for the heavier components in a sample. Combination of the two analyses yields a detailed analysis of the gas or liquid sample. The complexity of the condensate heavy ends is evident from the chromatograms of these fractions, and the predominance of the paraffin hydrocarborn serves as a useful marker in interpreting the chromatograms. The K-values obtained in this analytical method are presented for a high-pressure condensate system and predict closely the observed volumetric behavior of the system. INTRODUCTION Quantitative analysis of hydrocarbons from natural gas reservoirs is necessary for several reasons—to calculate the amount of sales gas produced, to calculate the amount of natural gasoline produced, to plan a liquid recovery system, or to calculate the potential economic value of a reservoir produced under one or more of several different conditions. Analysis of natural gas fluids produced to the surface consists of identifying and computing the mol fraction of each component of the mixture. Although methane is the predominant component, varying amounts of ethane, propane, butane, pentanes and heavier components are also present. Materials containing up to 30 carbon atoms occur in amounts which decrease with increasing molecular weight. However, the quantities of components in the 20 to 30 carbon atom range are usually so small that their importance is negligible, and they are undetect-able in natural gas by ordinary analytical methods. All the components up to those having 20 carbon atoms may sigsficantly affect phase behavior, however. Commonly, only the methane-through-pentane fraction is analyzed quantitatively for each component, while components heavier than Pentane are lumped and repored as "hexane- plus". Expensive, tedious techniques are required for analysis of this fraction. Consequently the detailed analyses needed for prediction of reservoir behavior are usually undertaken only when major gas fields are being developed. The need for complete analyses of condensate systems is apparent when it is recalled that most gas fields are produced by pressure depletion. As the pressure declines, some of the heavier hydrocarbons are lost as liquids which are in the reservoir. In many instances the amount of liquid in equilibrium with the gas phase at high pressure constitutes only 1 or 2 mol per cent of the total system. Flash calculations generally must predict the actual amount of liquid with an accuracy of a few per cent in order to be useful. This retrograde condensation has been understood for years, but accurate correlation methods to permit quantitative prediction of phase behavior in the retrograde region are not presently available. The increasing importance of natural gas has made accurate prediction of phase behavior and composition of produced natural gas streams an economic necessity. The work reported here was undertaken to provide a rapid, economical method for obtaining the vapor-liquid equilibrium information needed to predict accurately the composition of the fluids produced from a gas reservoir throughout its life. TO develop this method, a pressure cell equipped with windows was designed and built for observing the volumes of liquid and gas present at reservoir pressures and temperatures. Use was made of established chromatographic methods for rapid and detailed analysis of both phases. This paper describes the equipment and techniques developed for obtaining vapor-liquid equilibrium data, presents the results of analyses of a condensate system, and indicates the usefulness of these data in predicting hydrocarbon phase behavior. DESCRIPTION OF EQUIPMENT USED The equipment used in obtaining the required information on phase behavior and the complete analysis of hydrocarbon mixtures through C, will be described first, followed by a discussion of the operation of the equipment. It will be helpful, however, to consider first a brief outline of the technique used. A sample of separator gas and liquid is charged to the windowed cell (see Figs. 1 and 2) where volumetric equilibrium phase behavior at reservoir pressures and temperatures can be determined. Then samples of the coexisting phases are withdrawn. The methane-through-pentane fraction is analyzed with a chromatograph equipped with a hot-wire detector, and the pentane-plus fraction is analyzed with a second chromatograph equipped

Jan 1, 1965

By R. Wiesenthal

A method is developed for improving the low recovery efficiency which results when viscous oils are flooded by water. Viscous oil has been diluted with a lighter liquid miscible in it in any ratio which produces a mixture possessing a reduced viscosity and which is then displaced by water. Light gasoline was used in these tests to reduce the viscosity of the oil in place. Oils possessing viscosities in the range of 1.28 to 324 cp at the test temperature of 30C were displaced from a cylindrical unconsolidated sand core, with a length of 2 m (6.6 ft) and a diameter of 10 cm (3.9 in.) at a flood of 6.27 cc/sq cm/hr. In two test series up to I PV of light gasoline was injected into the core before water flooding. In the first series the recovery of a low-viscosity oil (1.28 cp) increased as the volume of light gasoline was increased. However, only the oil which could not be recovered by water flooding alone could be replaced by the light gasoline. In the other series, the recovery of a more viscous oil (13.5 cp) was not increased even if large quantities of light gasoline were injected. An oil bank built up ahead of the water front so that the flood water met with equivalent saturation conditions, as would be the case in flooding solely with water. However, the recovery of more viscous oils could be increased substantially when the core was first flooded with water, followed by a slug of light gasoline, and then flooded with water again. Furthermore, an oil bank ,built up ahead of the second water front. The formation of this oil bank is particularly favorable for the displacement process; the greater part of the light gasoline injected was pushed out of the core ahead of the oil and recovered. Using a light oil of low viscosity, it was possible to replace only the residual oil that could not be recovered by water flooding by the light gasoline. As the basic study had shown that highly viscous oils could be recovered by flooding first with water, then with light gasoline, followed by a second injection of water, the study was extended by tests made with reservoir fluids and cores taken from the low-gravity oil reservoir of the Wietze field, Germany. A great increase in recovery over flooding with water alone was obtained in these tests. The results may show a way to recover high-viscosity oils by a comparatively simple flooding procedure. INTRODUCTION In fields where oil recovery has not benefited from natural water drive, increased recovery may be obtained by the injection of water into the reservoir, especially in the case of fields which produce high-gravity, low-viscosity oil. There is a rather close relation between the viscosity of the reservoir fluid and the recovery that can be obtained by water flooding. As the viscosity increases the oil recovery decreases proportionately, so that an oil viscosity about 30 times that of water is generally considered to be the upper limit for an economically successful water flood.' This limitation has been recognized theoretically by Buckley and Leverett,' and has been verified by numerous laboratory flooding tests, such as by Croes and Schwarz,V n the viscosity ratio range of oil to water from 1 to 500. The reason for the inferior recovery of viscous oils by water flooding can best be understood by considering the mechanism of flood tests. In such tests water is introduced into one end of a sand-packed tube which is saturated with oil of the viscosity to be tested, and oil is produced from, the other end. Oil, with a viscosity lower than that of water, is pushed ahead of the advancing water with considerable uniformity, and irregularities in the displacement process tend to be eliminated as the oil saturation decreases and the water saturation increases. However, oil recovery is not complete, for the water forms interfaces with the oil and traps residual oil droplets in certain of the pore spaces. These residual oil droplets are distributed rather uniformly through the sand, and they are not susceptible to recovery regardless of the amount of water passed through the test cylinder. In the case of an oil possessing a viscosity greater than that of water, the injected water tends to move irregularly through the sand in the test cylinder. At the place in the sand where initial water movement takes place, resistance to the movement of water is lowered and the water finger so formed will grow perceptibly; it will extend to the outlet end of the test cylinder so that a breakthrough of water takes place long before much of the oil-saturated volume of the sand has come in contact with water. Once the breakthrough of water takes place, the water-oil ratio is gradually approaching the economic limit, even though much of the oil in the test cylinder remains unaffected by water. Disconnected drops of oil are trapped in the portion of the sand that has been invaded by water, just as when an oil of low viscosity is flooded. The poor recovery of the more viscous oils, therefore, is due to the irregularity of water movement through the test cylinder, and the rate of water movement through restricted sections of the oil-saturated sand increases with the more viscous oils. There are three ways whereby the recovery of viscous oils by water flooding can be improved:' (1) adding materials which will increase the vis-

Jan 1, 1965

By W. E. Howlett, R. L. Slobod

Jan 1, 1965

By B. Habermann

Artificially consolidated sand models, representing one-quarter of a five-spot, have been developed and used to study factors aflecting misciblt. displacrmenr. Sweep efficiency at breakthrough, size of the mixing zone between two miscible liquids and per cent area contacted by drive after breakthrough were determined tor the high mobility ratios encountered in actual resvrvoirs. Quantitative relationships between the degree of viscous fingering and mobility ratio were obtained by measuring the length of the fluid interface. Scaled miscible-slug experiments, supported by field evidence reported in the literature, have shown that when the slug is followed by dry gas the process is less efficient than expected. In addition to low areal sweep efficiencies encountered for high mobility ratio displacements, the effectiveness of a miscible slug is greatly reduced. This is caused by an accelerated growth of the mixing zone between the driving and the displacing fluids. INTRODUCTION Considerable interest in miscible displacement as a secondary recovery process has been shown by the oil industry in recent years. In addition to extensive laboratory investigations, numerous field operations have been initiated so that 39 miscible displacement projects had been reported' in the United States by early 1959. By far the most popular approach in these projects is the use of an LPG slug followed by dry gas. Laboratory tests have indicated that with a small miscible slug, amounting to several per cent of the hydrocarbon pore volume, practically 100 per cent recovery of the oil in place can be achieved and that all of the slug material can be recovered by a subsequellt gas drive. However, these experiments were conducted on small-diameter long cores that can be considered as essentially uni-dimensional reservoirs that do not represent actual field conditions. It soon became apparent that the inherently adverse mobility ratios existing in the reservoir cannot be ignored in any miscible displacement field project. When a less viscous fiuid such as an LPG slug drives a more viscous reservoir oil, fingers develop that reduce the areal sweep efficiency. The slze of the miscible slug to be used is a subject of some controversy. It has been proposed that the size of siug needed be calculated by setting lengths of the mixing zone proportional to the square root of the distance traveled. This proposal is based on experimental evidence' and theoretical considerations"' in which fingering was not a factor. A prerequisite for successful application of the miscible displacement process to field operations is a better understanding of the following factors. L What is the effect of mobility ratio on the rate of deterioration of a miscible slug? 2. How large a slug should be used for given reservoir conditions? 3. What is the sweep efficiency? The purpose of this investigation was to provide some additional laboratory data that should help to answer these questions. More specifically, it was intended to study more quantitatively the effect of viscous fiinger-ing, mobility ratio and slug size on the efficiency of this relatively new secondary recovery process by means of experiments on thin, horizontal, porous medium models.

By S. J. Prison, C. R. Smith

The effect of fluid injection to control water coning in oil and gas wells was investigated. Analytical and model techniques were employed. The factors investigated were the position and length of the completion interval, the point of fluid injection, the viscosity of the injected fluid and the relative thickness of the oil and water sections. The resulting influence of these factors on the net producing water-oil ratio was determined. Several important conclusions can be drawn from the study. In general, it was found that the net producing water-oil ratio can be reduced by fluid injection. The magnitude of this reduction depended on the factors listed above. An important practical consideration is that the injection fluid may be either oil or water. If the injected fluid is less dense than the connate water of the reservoir, the fluid will not be lost. This fact is reassuring when valuable oil is being injected. Efforts to suppress water production were more successful when the injection fluid was more viscous than the reservoir oil, or when a zone of reduced permeability existed in the vicinity of the point of fluid injection. Under test conditions, little benefit was derived through the use of impermeable barriers or cement "pancakes". INTRODUCTION The occurrence of water coning has been known for at least 60 years. In thin oil or gas pay sections, the presence of an oil-water or gas-water contact hinders production and often causes early abandonment of the afflicted well if a completion is even attempted. Even when relatively thick pay sections are found, the encroachment of water when a water drive is present will eventually pose serious water coning problems. This water is often corrosive, expensive to separate from the oil or gas and is costly to dispose of. The theory of water coning has been discussed by a number of authors. 1,2,3,4 Briefly, water coning to the producing interval in a well is due to pressure gradients resulting from the production of fluid from the reservoir. These pressure gradients will cause a water cone to rise toward the bottom of the producing interval if a water-oil or water-gas contact exists. The tendency of the water to cone is offset or partially offset by gravity forces since the water has a higher specific gravity than the oil. A balance then exists between two forces, gravitational forces arising from the difference in specific gravities of the oil and water, and the pressure gradients causing the flow of fluids to the wellbore. If the pressure gradient exceeds the gravitational force, water coning to the wellbore occurs and water production results. Through the years considerable thought has been given to the water coning problem. More than 50 U. S. Patents have been granted to inventors on the subject. A relatively complete literature review of the water coning problem has been made. 5 A number of these patents hold considerable promise for the solution or the partial solution of the water and/or the coning problem. Very little has been written describing field tests of techniques for the suppression of water coning. A notable exception is the paperby West. 6 He reports success in reducing gas coning by a combination of gravel packing and oil injection above the oil-producing interval. He also describes a comparable method to prevent water coning, but provides no field examples. This study experimentally and analytically verifies the benefits of oil injection as a means of partially or completely suppressing the water cone. While the gas coning problem was not treated, it is anticipated that results comparable to those obtained in water suppression could be obtained with reduced oil injection since the viscosity contrast between oil and gas exceeds that between oil and water. For the purpose of this study, the conventional potential flow theory was applied to the water-coning problem. The experimental verification centered on both a radial and a linear model. The model study permitted the investigation of complex flow configurations and the use of fluids of differing densities and viscosities. No analytic expressions are available to permit a solution of the problem as stated (see Fig. 1).

By D. C. Lindley, D. H. Gaucher

The displacement of oil by water in a waterflood project is accomplished by the action of transient viscous, gravitational and capillary forces which drive fluid through interconnecting pore spaces toward production wells. The relative importance of each of these forces to the effectiveness of the displacement and the production history depends on properties of the reservoir itself and the manner in which it is operated. There has been considerable discussion'-3 concerning the importance of these factors and the extent to which they control the performance of individual waterflood projects. No data have been presented, however, to demonstrate the effect on waterflood behavior of variations in rock permeability, water-injection rates and mobility ratios in three-dimensional reservoir systems in which viscous, gravitational and capillary forces are allowed to assume the proper relative influence on fluid flow. Much information relevant to water flooding has been gained from the study of field case histories. But, the complexity of the systems involved, the difficulty of defining geometry and obtaining sufficient data, and the fact that the same reservoir is never flooded twice with the same initial conditions, all obscure the role of process variables such as rate. Further, mathematical calculations of waterflood behavior are restricted at this time by available mathematical techniques to highly simplified situations. Scaled models, however, do offer an approach at the present time to evaluation of the effects of the pertinent variables. Several investigators have predicted prototype reservoir behavior from observations of model waterflood performance under a variety of conditions. Two-dimensional aspects of waterflood pattern efficiency have been investigated by Muskat, Aronofsky, Dyes, et al, raig, et al,' and Rapoport, et al. he effects of gravity segregation, viscous fingering and permeability stratification have also been studied in two-dimensional models. Craig, et al," studied the importance of gravitational forces on frontal displacement in three-dimensional models with wells on a five-spot pattern. This work in- vestigates the tendency for water to segregate in homogeneous reservoirs because of gravity and to under-run the oil for a range of water-inject ion rates. The effects of permeability stratifications were investigated by these authors using miscible fluids to simulate water-channeling through permeable strata in systems of negligible capillary force. In this investigation, a three-dimensional model of a five-spot pattern flood was designed to represent typical prototype conditions, and the sand and liquid properties were selected so that capillary, gravitational and viscous pressure gradients were scaled. It was not, however, possible to achieve exact similitude regarding viscous fingering when the more viscous oil was used. The effects of rock stratification on waterflood behavior were studied by changing both the permeability ratio and the position of the more permeable sand stratum in a two-layered, horizontal model. For each type of stratification, constant-rate floods were performed and the rates were varied, in separate experiments, over a tenfold range. To observe the effect of mobility ratio on flood performance, two series of constant-injection-rate floods were performed on the most highly stratified reservoir at a water-oil viscosity ratio of one.

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst's horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter's edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By H. N. Hall

Various factors must be considered in an engineering evaluation of gravity-drainage reservoirs. Among these are: (1) the effect of producing rate on total oil recovery; (2) the effect upon well productivity and ultimate recovery of the pressure level maintained during the producing life of the reservoir; (3) the economic advantage of full or partial pressure maintenance; and (4) estimate of the rate of gas production and injection and the possible purchase of gas under conditions of full pressure maintenance to ascertain compressor facilities needed. All of these factors can be evaluated only when a reliable method is employed for determining reservoir performance in gravity-drainage reservoirs. The purpose of this paper is to present a general method for calculating the performance of a gravity-drainage reservoir. This method is applicable for conditions of complete pressure maintenance, partial pressure maintenance and normal pressure depletion. Provisions are made to take into account variations throughout the reservoir of reservoir configuration, changes in permeability and fluid composition. Based on the method presented in this paper, an IBM 650 computer program has been developed. The past performance of an actual gravity-drainage reservoir producing under conditions of declining pressure and no gas injection was duplicated using this program. INTRODUCTION In tilted reservoirs the production of oil is influenced by drainage of oil from upstructure to downstructure locations. When this downstructure drainage of oil is sufficient to cause effective segregation of the gas and oil in a reservoir, the reservoir is usually classified as a segregation drive or gravity-drainage reservoir. (Discussion will be restricted to gravity-drainage reservoirs which have no encroachment of edge water.) The important feature in gravity-drainage reservoirs is the density difference between reservoir oil and gas. These phases tend to segregate in the reservoir with the result that in the gas cap the oil saturation is maintained at a higher level by drainage of oil from the gas-cap area. Oil can be produced from the oil zone at a low gas-oil ratio and reservoir energy is thereby conserved. The standard material balance in not adequate for predicting gravity-dramage reservoir performance because it does not take into account the difference in saturation above and below the gas-oil contact. Several authors'.' have presented methods for calculating the performance of gravity-drainage reservoirs in which reservoir pressure is maintained constant by gas injection into the gas cap. Using some simplifying assumptions, these methods can be employed with a desk calculator to give acceptable results. The problem of predicting the performance of gravity-drainage reservoirs under the conditions of declining reservoir pressure is many time more complex than that of constant pressure. fierefore, attempis to develop a method suitable for desk calculation have required excessively simplified assumptions. In the past several years, highspeed digital computers have become more widely available for reservoir engineering problems. These corn puters are well suited to problems such as the prediction of the performance of gravity-drainage reservoirs with pressure decline. Many of the simplifying assumptions necessary for hand computation can be eliminated so that a realistic approach to the gravity-drainage process can be made. CONCEPTUAL PICTURE OF OIL MOVEMENT IN GRAVITY-DRAINAGE RESERVOIRS Before attempting to develop an analytical treatment for conditions occuring in a gravity-drainage reservoir, a concept should be formed concerning the movement of fluids in the reservoir as oil is produced. A review of the literature'.' shows that it is customary to classify gravity-drainage operations into two categories—(1) with complete pressure maintenance, and (2) with declining pressure. The same line of reasoning will be followed in presenting the concept of the movement of fluids in the reservoir because it is easier to visualize the movement of fluids under conditions of complete pressure maintenance. After discussing complete pressure maintenance, an analogy will be made between that and the case of declining pressure. It should be kept in mind throughout that the final aim for the problem of solving gravity-drainage performance with digital computers will be to develop a general program for any kind of gravity-drainage reservoir. COMPLETE PRESSURE MAINTENANCE One feature which is generally common in gravity-drainage reservoirs is a gas cap located at the top of the structure. This is shown in Fig. 1(a). Fig. l(b) shows oil saturations that might occur through the reservoir. In the gas cap, oil saturation is lower than

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By J. W. Lacey, F. H. Brinkman, J. E. Faris

This paper presents an analysis of the high-pressure, gas-driven LPG-slug process, based on fluid flow tests in areal models. Two types of tests were made. One series was made in low-pressure models which permitted observation of fluid movement. Three completely misci-ble analog fluids were used. A second series of tests was made in high-pressure models using methane, propane and a light refined oil saturated with methane at room temperature and 1,550 psig. Under the test conditions of room temperature and a pressure level of 1,550 psig, the phase diagram for the fluids used is similar to those for many of the field systems where the process is considered for use. A method for using these laboratory data to calculate field performance of the process is outlined. As a result of this work, it is concluded that small banks of LPG (5 per cent HV or less) are not effective in increasing oil recovery in horizontal reservoirs. Znstead, where small banks are used, the driving gas quickly penetrates the LPG bank because of fingering and channeling; and from this point on, the process behaves essentially as an immiscible gas-injection project. The validity of this conclusion was substantiated by: (1) laboratory studies of the effect of rate, model size and mobility ratio on miscible displacement in areal models; and (2) calculation of field recovery, which compared closely with actual field recovery. INTRODUCTION Field applications and pilot tests of the gas-driven, LPG-bank, oil-recovery process are on the increase. Most of these tests are employing small banks (2% to 5 per cent hydrocarbon volume) of LPG, which is miscible with both the driving gas and the oil in place, in an effort to attain an effective yet economical miscible displacement of oil by gas. The expectation of miscible displacement with small banks is based on the concepts that: (1) lengths of solvent-oil mixed zones, measured during miscible displacements in long slim cores, are representative of those that will occur in the field; and (2) areal sweep efficiencies5 measured in electrolytic model studies are applicable to miscible displacement in reservoirs. Our experimental evidence indicates that the mixed-zone lengths and sweep efficiencies mentioned are not applicable to miscible displacement in reservoirs. In this paper we present an evaluation of the gas-driven, LPG-bank oil-recovery process based on fluid flow experiments in areal models. These results are used to predict the performance of a field pilot test of this process, and the results are compared with the actual test results. THE EFFECTIVENESS OF LPG BANKS The effectiveness of LPG banks of various sizes in accomplishing miscible displacement of oil by gas was determined by displacement tests in a model representing one-quarter of a confined five-spot pattern. The model, 12 X 12 X 1/4 in. in size, was packed uniformly with glass beads. It was operated at a pressure of 1,550 psig and at room temperature; the fluids used were methane, propane and refined oil saturated with methane at 1,550 psig. The saturated oil had a viscosity of 1.2 cp at room temperature. The results of these tests are shown in Fig. 1, where recovery is plotted as a function of the volume of fluid injected for: (1) an immiscible-gas drive with an oil-to-gas viscosity ratio of 85; (2) three sizes of LPG banks; and (3) two completely miscible displacements with mobility ratios of 85:1 and 10:1. (The miscible displacement data plotted are the averages of several tests.) The results show that a 2½ per cent bank of LPG does not increase oil recovery over that obtained by immiscible-gas drive. The 7 and 17 per cent banks are, respectively, about 30 and 50 per cent effective (with reference to the M = 85 miscible-displacement curve) in increasing oil recovery at the point where 2½ hydrocarbon volumes (HV) of fluid have been injected. The percentage effectiveness of the banks at a given volume of fluid injected is defined as Recovery by Bank — Recovery by Gas Drive Recovery by Miscible Recovery by Flooding Gas Drive We conclude from these results that banks of LPG smaller than about 5 per cent HV in an individual stratum will cause from little to no increase in oil recovery. This finding was substantiated by work in low-pressure models, which permitted visual observation of fluid movements. In these tests, three completely miscible fluids having the same viscosity ratios as those used in

By M. Prats, D. G. Russell

The performance of a well in a bounded, layered reservoir with interlayer crossflow has been investigated mathematically. The system studied comprises a centrally located well in a bounded cylindrical reservoir composed of two layers of contrasting physical properties. The reservoir contains a single fluid of small and constant compressibility. In the case of a well producing at constant pressure, the production rate of the crossflow system declines exponentially after the early transient behavior dies out. For wells which produce with constant rate, the pressure-time relationship is linear following the transient. For very early production time, the production rate of the two-layer crossflow system behaves as that of the two-layer stratified system (the case in which the layers are in communication at the well only). Except for early time, in most practical cases the performance of a two-layer reservoir with cross-flow can be duplicated by that of a single-layer reservoir. The equivalent single-layer reservoir must have the same pore volume and the same drainage and wellbore radii as in the two-layer case. In addition, the equivalent reservoir must have a "kh" (permeability x thickness) product equal to the sum of the "kh" products of the respective layers in the crossflow system. If the porosity and permeability distributions are known, the production and pressure behavior of a well with and without crossflow can be predicted. From a comparison of the observed pressure or production performance of a well and the theoretical behavior, the occurrence of interlayer crossflow can be detected. INTRODUCTION When flow is possible across the bedding plane from a low-permeability layer into an adjacent layer of moderate or high permeability, the earlier depletion of the more permeable layer will cause such flow to occur. If this situation exists, the oil recovery from tight layers can be appreciable; thus, the economic ramifications of interlayer crossflow are apparent. The possible importance and extent of flow between adjacent and connected layers — whether resulting from capillary, gravitational, or viscous forces — have been the subject of considerable speculation for many years. Jacquard1 presented a rigorous analysis of the crossflow due to single-fluid expansion occurring between two adjacent and connected layers when the well produces at constant rate. His presentation was entirely mathematical, however, and no numerical results from evaluation of the mathematical solutions were presented. Our objectives in the preparation and presentation of the present paper are essentially twofold— (1) to analyze the behavior of a well which produces a single compressible fluid from a bounded, layered reservoir with interlayer crossflow, and (2) to present practical numerical results and simplified formulas to facilitate engineering analysis of reservoirs with crossflow. The constant-terminal-pressure case was analyzed rigorously, whereas the constant-terminal-rate case was treated only for semisteady-state conditions. DESCRIPTION OF THE IDEALIZED TWO-LAYER RESERVOIR The reservoir to be studied is assumed to be bounded and cylindrical in shape, with a completely penetrating well at its center. The reservoir is composed of two adjacent horizontal layers between which fluid is free to flow. Each layer is assumed to be homogeneous, isotropic and completely filled with a fluid of small and constant compressibility. The total thickness of the reservoir is small enough

By T. L. Loucks, E. T. Guerrero

Pressure drop characteristics in a system composed of two adjacent concentric regions of different permeability were studied. The differential equations for continuity of mass flow in the two regions were solved using the Laplace transformation and the necessary boundary conditions to give the pressure distribution in the composite reservoir, the resulting equation for pressure drop at the inner boundary was evaluated for a variety of composite reservoirs and compared with the results for a uniform reservoir. From this study it was found that under certain conditions the permeability in both zones, as well as the size of the inner zone, can be determined from the pressure drop curve. INTRODUCTION The theory for the pressure distribution and pressure build-up behavior of a well producing a single, slightly compressible fluid from infinite and finite homogeneous reservoirs was presented by Hornerl and Miller, Dyes and Hutchinson.2 Extensions on this original work to provide improved and extended interpretations and better agreement between theory and observed results have been made by Matthews, et a1,3 van Everdingen,4 Gladfelter, et al5 Stegemeierand Matthews,6 Hurst and Guerrero,7 and Perrine.8 More recently Lefkovits, et al,9 studied pressure build-up behavior in bounded reservoirs composed of stratified layers. Houpeurtl0 has suggested various approaches to the general problem of variable permeability and porosity but presented no analytic solutions for particular permeability variations. Albert, Jaisson and Marionll studied the finite composite reservoir and presented numerical solutions to the unsteady-state case and an analytical solution valid only for large times. They also studied the so-called pseudo-steady state for several examples of radial permeability variations. Very similar examples have been treated in the unsteady state with application to pressure build-up by Loucks in an unpublished manuscript. More recently Hurst12 has presented the complete point-sink solution (valid for all times) for the infinite composite reservoir. He applied these solutions to interference between oil fields along with an even more elegant application of his explicit solution to the material-balance equation including water influx. Mortada 13 approaches the same application by avoiding the point-source limitation but gives the solution for the aquifer region which is valid only for large times. Hopkinson, Natanson and Temple14 have treated both the finite and infinite composite reservoir obtaining the pressure distribution for the inner zone valid for large times. This paper presents a theoretical study of the pressure distribution in an infinite composite reservoir composed of two adjacent concentric regions of different permeability. The object was to determine the manner in which pressure drop at the inner boundary of a composite reservoir depends upon time, the permeability of each zone and the size of the inner zone. Expressions for the pressure distribution in both zones are developed which take into account the radius of the sink and are valid for small times as well as large times. It was felt that an understanding of the pressure drop behavior in various composite reservoirs would be of assistance in the interpretation of some pressure build-up curves which do not behave according to the theory derived for uniform systems. Often the region surrounding the wellbore is either more permeable or less permeable than the reservoir because of the various drilling and completion practices. The effects of reduced permeability due to drilling-fluid invasion and of increased permeability due to fracturing or acidizing need to be more carefully defined. Therefore, an equation for the pressure drop in a composite reservoir was developed, and the effects of both the permeability in each zone and the size of the inner zone were studied. STATEMENT AND SOLUTION OF PROBLEM BOUNDARY CONDITIONS AND ASSUMPTIONS The physical system studied was a composite

By G. J. Heuer, J. N. Dew, G. C. Clark

The effect on well behavior of partial permeability barriers, changes in producing rates and well spacings have been calculated through use of a radial, unsteady-state, two-phase-flow mathematical model. This model allows for variations in permeabilities and porosities with distance from the wellbore. The numerical methods necessary to solve the problem require use of a highspeed digital computer, in this case an IBM 704. In each instance, pressure and saturation gradients, gas-oil ratios and recoveries around a well producing by solution-gas drive have been calculated as a function of time and distance from the wellbore. Comparisons are made to show the effect of changes in producing rate and of varying permeability near the wellbore and out in the formation on the production and pressure history of the well. The effect of different well spacings on producing rate and ultimate recovery is considered. Although the mathematical description of the reservoir is simplified in comparison to an actual reservoir, the results do give some insight into the difficult problem of spacing and proration in a heterogeneous solution-gas-drive reservoir. Results show that for the cases considered spacing has little effect on ultimate recovery and that permeability barriers removed from the well decrease producing rate for a period of time but have only a small effect on ultimate production. INTRODUCTION The effect of spacing and producing rate on the production characteristics and ultimate recovery of a well producing by the solution-gas-drive mechanism have been topics of interest to the oil industry for many years and the subject of many hearings before state regulatory bodies. The problem is very complicated because of the difficulty of simulating in the laboratory a reservoir producing by the solution-gas-drive mechanism under anything like normal field conditions. This paper gives the initial results of a study aimed at determining the effect of these factors on production from reservoirs simulated by a mathematical model. Two-phase, unsteady-state equations are used to calculate pressures, saturations, and rates of oil and gas flow as a function of time and distance from the wellbore. The model is radial, and reservoir properties may be varied with distance from the wellbore. A decrease in permeability at a given distance results in a low-permeability ring concentric with the wellbore through which all the fluids from more distant portions of the reservoir must flow. The current mathematical model allows for variation of permeability, porosity and saturations as a function of distance from the well. Any desired drainage radius for the well may be selected. Drainage radii of 745 and 1,053 ft, which correspond to 40- and 80-acre spacing, were chosen for the calculations which follow. A large amount of the input data has been held constant. No investigation was made into the effect of changes in fluid properties and relative permeability characteristics on oil and gas production behavior. Fluid properties and relative permeabilities are essentially those of the White Mesa portion of the Greater Aneth area in southeastern Utah. A recently concluded hearing before the Utah Oil and Gas Conservation Commission to determine spacing for the field aroused interest in the possible effect of permeability constrictions located some distance from the well on production and recovery. This interest stimulated the present study. In this paper, no attempt has been made to simulate any field production characteristics exactly. Only the solution-gas-drive mechanism has been investigated, and production by liquid expansion has not been considered. This investigation is exploratory in nature but does give some insight into the problems of spacing and proration in a solution-gas-drive field having permeability variations. THE MATHEMATICAL MODEL The mathematical model is similar to that proposed by West, Garvin and Sheldon.' Approximately the same techniques have been used in arriving at a solution. The equations are nonlinear, partial differential equations which have been known for some time. The particular equations used for this study neglect the effects of capillary pressure and gravitational forces. The advent of modern digital computing equipment has made their solution practical. Basically, the method solves for pressure, saturations, and oil and gas flow rates as a function of distance and time. In our model, distance to the

During development mining, ventilation capacity may decrease and the ventilation requirements may depart from initially planned values depending on the mining rate, increases in the methane emission rate, and changes in ventilation airway leakages. Insufficient ventilation air quantities and the occurrences of high methane liberation rates can overwhelm the existing quantity of ventilation air and result in elevated levels of methane gas during development mining. This condition can both limit the development of that section and increase the risk for an ignition that may lead to an explosion. Thus, the prediction of ventilation air requirements prior to mining can enhance worker safety by reducing the likelihood that explosive methane-air mixtures form. This study presents an approach for prediction of methane inflow rates using coalbed methane reservoir modeling. For this purpose, a two-phase coalbed reservoir model of a three-entry type roadway is developed. In the model, grids are dynamically controlled to simulate development of entries. Examples are given for methane inflow and ventilation air requirements as a function of mining rate, development section length, mining height, and the existence of degasification boreholes. This technology can be used to limit the methane concentrations occurring as a result of the influences of various coalbed and operational parameters.

Jan 1, 2008

By E. A. Swedenborg

The unit agreement for the Wertz Dome field, Wyoming, was approved by the Acting Secretary of the United States Department of the Interior on November 4, 1937, effective on December 1, 1937. The stated objectives of the agreement are to conserve and put to beneficial use all oil and gas produced; to make possible a uniform withdrawal of oil and gas in order to maintain equalized reservoir pressures; to provide for an orderly determination of the structural features of the productive horizons; and to permit the injection of gas for pressure maintenance. The purpose of this paper is to show by a discussion of the problems involved and the engineering practices employed by the unit operator, how well the objectives of the unit agreement have been accomplished in the development and production of oil and gas from the field. HISTORY AND DEVELOPMENT The Wertz Dome oil and gas field is located 88 miles southwest of Casper and 38 miles north and westerly from Rawlins, Wyoming. It is served by a state secondary oiled road connecting with U. S. Highway No. 220 at Lamont, a distance of 3 miles from the field. It has been a producing field for 27 years, first as an important producer of gas and more recently of oil. The Producers and Refiners Corp. completed Well No. 1, NE1/4SW1/4 sec. 7, T. 26 N., R. 89 W., the first productive well in the field, in September. 1921, in the Muddy sandstone of Upper Cretaceous age, at a depth of 3435 feet for an initial production of 42,000,000 cu. ft. of gas per day. In 1922 pipe lines were laid to Casper, Wyoming, for supplying fuel to the Standard of Indiana refinery, 'the adjoining town of Mills. and the Producers and Refiners Corp. pump station No. 6 in Mills; and to Parco (now Sinclair), Wyoming, for supplying that refinery with fuel. The gasoline content of the gas was extracted in an absorption plant located near the town of Mills. No compressors were used as the line pressure was sufficiently high to flow the gas through the absorbers. As high as 30,000,000 cu. ft. of gas per day were treated in the plant extracting 6000 gallons of gasoline. In September 1937, the pipe line was discontinued and the plant was dismantled at which time 61,000,000,000 cu. ft. of gas had been produced from the seven Muddy, Cloverly, and Sundance gas wells drilled in the field. Of this volume of gas 800,000,000 cu. ft. were produced from the Frontier formation, 34,200,000,000 cu. ft. from the Muddy sandstone, 24,000,000,000 cu. ft. from the Cloverly, locally called Dakota, and 2,000,000,000 cu. ft. from the Sundance formation. The pipe line was discontinued because of the near depletion of the gas reserves in Wertz and nearby fields, and hecause of the local field requirements demanding the retention of the remaining gas reserves. In December 1936, the Sinclair Wyoming Oil Co. completed Well No. 10, SE1/4NW1/4 sec. 7. T. 26 N., R. 89 W., at a depth of 5886 feet. as the first Tensleep sand producer with an intial of 1700 barrels of oil per day. The well was originally drilled only 14 feet into the sand and in August 1939 was deepened to a depth of 6161 feet, 289 feet in the sand, and the daily produe tion increased to 8350 barrels of oil. In February 1948 Well No. 22, SE1/4 SE1/4 sec. 1. T. 26 N., R. 90 W., was deepened from the Tensleep to the Basal Amsden sand and completed at a depth of 6635 feet for an intial of 600 barrels of oil per day in the Basal Amsden. In April 1948, Well No. 2. SW1/4 NE1/4 sec. 1, T. 26 N., R. 90 W.. was deepened from the Tensleep to the Madison limestone and was completed at a depth of 7193 feet for an intial of 1145 barrels of oil per day in the Madison. Well No. 26, SE1/4SE1/4 sec. 1, T. 26 N.. R. 90 W., was deepened from the Ten-sleep to the Cambrian and in October 1948 it was completed in the Cambrian for an intial of 277 barrels of oil per day. GEOLOGY AND PRODUCTIVE ACREAGE The surface is Steele shale of Upper Cretaceous age, covered in part by gravel beds. The topography of the field is that of a rolling hill with a steep draw or gulch along the northern edge of the field. The Steele shale lies conformably on the lower formations as

Jan 1, 1949