Search Documents

By J. N. Sicking, F. H. Brinkman

A new method for obtaining equilibrium vaporization ratios (K-values) for reservoir fluids has been developed and tested. By application of the method, complex experimental measurements of liquid and vapor phase compositions are eliminated. This simplified technique reduces the cost of experimental equilibrium-ratio data for reservoir studies of condensates and volatile crude-oil systems. The method is designed for systems of constant composition and, therefore, is best suited for depletion studies where compositional changes at high pressures are minor. The basic data required, in addition to the composition of the initial reservoir fluid, are the relative vapor-liquid volumes and densities at reservoir temperature and variom reservoir pressures. Tests demonstrated that the method predicts equilibrium ratios accurately for condensates. A single test on a crude oil was not conclusive; further testing will be necessary before the accuracy of the method can be determined for crude-oil systems. In addition to determining equilibrium ratios, the calculation method provides information on the physical properties of the "plus" component in the vapor and liquid phases. The "plus" component is that mixture of components heavier than the least volatile fraction analyzed. This information is useful in studies of both natural depletion and cycling operations for condensate reservoirs where the heptanes-plus component in the gas phase is produced from the reservoir. INTRODUCTION As more volatile oil and condensate reservoirs are found, the use of phase behavior techniques to predict their performance is increasing in importance. These techniques have long been used for condensate fields and have more recently been applied to crude-oil fields containing oils of medium-to-high volatility. In these phase behavior methods, equilibrium ratios (K-values) are used to predict compositional changes in the reservoir fluids—thereby accounting for the recoverable oil that exists in the gas phase. The reliability of the predictions depends to a large extent on the equilibrium ratios used. These values must be obtained for each component for the entire pressure range being investigated. Unfortunately, because of the complex nature of hydrocarbon mixtures, accurate K-values are hard to obtain. The equilibrium ratios for a particular component will vary not only with the temperature and pressure, but also with over-all composition of the system. The importance of composition is quite critical at elevated pressures, but becomes negligible at pressures below about 300 psia. Therefore, because most phase behavior problems involve the high-pressure region, each fluid system becomes a special case. Experimental programs to determine characteristic K-values are quite difficult and time-consuming. Thus, it is often necessary to resort to approximations of the K-value data. Charts giving K-values for various mixtures and classes of mixtures are available in the literature. However, there are two major difficulties in using them: (1) the K-values of the "plus" component (that mixture of components heavier than the last one analyzed) must be obtained by extrapolation from the K-values of the other components; and (2) the K-values obtained must finally be adjusted by trial and error to agree with observed volumetric data. To eliminate these difficulties, a new method of determining equilibrium ratios was developed. Briefly, after the composition of the system as a whole has been analyzed, the method uses empirical correlations and the gross fluid properties of the system (relative vapor-liquid volumes and densities) to calculate K-values. Because the calculative procedure is long, it is best solved on a digital computer. About one hour of machine time on an IBM 650 computer is required to develop a K-chart for the fluid being examined. DEVELOPMENT OF THE METHOD OF OBTAINING K-VALUES Equilibrium ratios are defined as the ratio of the mole fraction of a component in the vapor phase to its mole fraction in the liquid phase. This statement is expressed in Eq. 1. A typical plot of equilibrium ratios for a particular system is shown in Fig. 1. It should be noted that, at pressures near the saturation pressure, the K-values appear to converge to a common point. This apparent convergence point is called the convergence pressure and is a characteristic of the system involved. Various empirical correlations of K-values have been noted. It has been observed that an isothermal plot of

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 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 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 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 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 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 A. S. Michaels, R. S. Timmins

A study of the effect upon oil recovery from an un-consolidated porous medium of chromatographic transport of selected reverse-wetting additives during water displacement is described. Flow tests were performed on hydrocarbon oil-saturated sand columns (containing connate water) into which small volumes of short-chain aliphatic amine solutions were introduced during water flooding. Concurrently, meusurements of oil-silica contact angles, oil-water interfacial tensions and oil-water distribution coefficients were measured us a function of additive concentration. Certain treatments were found to cause a slight reduction in oil recovery, while others produced a significant increase (up to 50 per cent of the normtzl irreducible rninimum oil saturation). Recovery eficiency was correlated most satisfactorily with the maximum degree of reverse wetting produced during additive transport; conditions which resulted in only mild oil wettability of the sand reduced recovery, while those favoring strong oil wetting increased recovery. A mechanism of reverse-wetting additive action during chromatographic transport is developed which ascribes the results to transient changes in oil-solid contact angle during flow, with consequent release of oil trapped by capillary forces. Implications of these observations with regard to improved waterflood efficiencies are discussed. INTRODUCTION The need for finding improved methods of oil production is amply illustrated by the large volume of oil considered economically unrecoverable by existing production practice. Water flooding has long been an excellent method of secondary recovery. However, even in those areas of a formation which are well swept by a water flood, 20 per cent or more of the oil remains trapped within the formation and is not recovered. This trapping of oil within the pores of the rock structure results primarily from interfacial forces. If these interfacial forces were altered, it is possible that the trapped oil could be released and recovered. The addition of surface-active agents to a water flood to alter the interfacial properties of the oil-water-solid system has received considerable attention as a method of increasing oil recovery. Studies of surfactant treatment have yielded varied and, in some cases, apparently conflicting results. Increased oil recovery using surface-active agents has been reported by a number of investigators. Others have found that use of surfactants in a water flood either has no effect upon oil recovery or decreases it. Since surface-active agents tend to alter the formation wettability and oil-water interfacial tension simultaneously, several investigators have tried to determine the role these two factors play in oil recovery. Data have been presented which indicate that the most desirable condition for the solid surface is strongly water-wet. Other results indicate that, although breakthrough recovery may be less, the greatest ultimate recovery is obtained when neutral wettability prevails.", " The desirability of low interfacial tension seems to depend upon the wettability of the formation. Reports of higher recoveries, when the interfacial tension is high in a water-wet system and low in an oil-wet system, are at odds with reports of increased recovery in water-wet systems when surfactants which decrease the interfacial tension are added to the water flood. A review of previous work studying the effects of using surface-active agents in a water flood shows that these agents produce widely different effects and that opinions as to the role of the various interfacial properties are diverse. The consensus, however, is that the interfacial properties are quite important in the displacement of oil from a porous medium by water. Even in those investigations where appreciable increases in oil recovery resulted from treatment with a surface-active agent, economic considerations generally made practical application of the treatment unattractive because surfactants tend to be adsorbed within the formation. This adsorption would necessitate using large quantities of surfactant if all of the flood water were treated to contain a given surfactant concentration, even if this concentration were quite low. Preston and Calhoun suggested that the principles of chromatography could be applied to the use of surface-active agents as a means of increasing oil recovery. A small volume of surface-active material is injec-