Search Documents

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

A series of laboratory experiments was conducted in which oil was displaced from a porous medium by water-driven slugs of alcohols or similar solvents. The solvents used were soluble to some degree in both oil and water and covered the range of solubilities from complete solubility in oil to complete solubility in water. Displacement experiments were conducted on 2-and 3%-in. in diameter consolidated cores and 1-and 2-in. in diameter unconsolidated sand packs. The cores and sand packs ranged in length from 1 to 30 ft, and they were saturated with brine and crude or refined oils. The solvents used included ethyl alcohol, isopropyl alcohol (IPA), tertiary butyl alcohol (TBAI, secondary butyl alcohol (SBA), n-amyl alcohol (NAA), methyl ethyl ketone, acetone and others. It was found that all of the oil present in a porous medium could be miscibly displaced by injecting a slug of mutual solvent and driving it with water. The oil-recovery efficiency was dependent upon (1) the relative solubilities of the solvent in oil and water, and (2) the distance traversed by the flood. For complete oil recovery from cores, a smaller amount of a preferentially oil-soluble solvent was required, compared to the amount of preferentially water-soluble solvent needed. The size of the solvent slug required varied inversely with the linear flooding-path length raised to the 0.65 power. Water-driven dual solvent combinations (an oil-soluble solvent slug followed by a water-soluble solvent slug) were found to effect complete oil recovery with less total solvent than any single solvent used. In these dual-solvent displacement experiments, the slug size required varied inversely with length raised to the 0.55 power. Based upon the experimental results, a theory was developed to describe the displacement of oil and water by mutual solvents, and equations are presented to predict the production history in a linear system. These equations take into account the properties of the solvents and the porous medium. INTRODUCTION Oil-recovery processes which utilize displacing fluids that are miscible with the reservoir fluids have been studied extensively in recent years. Because of the poor contact efficiency and high pressure requirements of the LPG-gas displacement process there has been considerable interest in the alcohol-water process, and a number of studies have been made on the recovery of oil through the use of solvents which are mutually soluble in oil and water. An investigation by Sievert, Dew and Conleyl indicated that the use of mutual solvents would be limited because the presence of water in a porous medium would cause a phase break in the leading edge of the displacing solvent. Their study also showed that, in consolidated cores containing. only oil, the displacement of oil by a water-driven mutual-solvent slug of tertiary butyl alcohol (TBA) was affected by the viscosity ratios of the fluids involved. Gatlin and slobod3 concluded that an isopropyl alcohol (IPA) slug acts as a miscible piston, completely displacing both oil and water until the alcohol content of the mixing zone falls below the concentration necessary to maintain miscibility. Their study was conducted on uniform unconsolidated sand packs. They concluded further that IPA could be used effectively to recover oil from a watered-out sand. In a paper by Taber, Gamath and Reed4 relating an investigation on sandstone cores, it was stated that the displacement of oil by mutual solvents, particularly IPA, was not a miscible displacement and that no improvement in efficiency could be expected with increase in flooding-path length. However, their analytical analysis of the displacement mechanism using TBA is, in fact, one which indicates that the displacement is controlled by miscible mixing. They suggested that the lack of improvement in efficiency with flooding - path

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 J. B. Agan, R. J. Fernandes

This paper utilizes the layered-system approach, modified to include areal sweep efficiency, to determine the miscible-slug size required and to predict the perfort7zance of a miscible-slug process in a highly stratified reservoir having a permeability variation of 0.76. The slug size in each layer was determined from theoretical calculations and laboratory data which showed that, to maintain miscibility in the reservoir, the slug size had to be 8.2 per cent of the floodable hydrocarbon pore volume. Economic considerations, however, such as the value of the miscible product injected, the cost of injection and the cost of reprocessing the miscible product, limited the amount of miscible product which could be injected into the lower-permeability layers. This resulted in subjecting only 60 per cent of the reservoir to a slug size sufficient to maintain miscibility. Utilizing the relationship of producing gas-oil ratio and areal sweep efficiency as a function of pore volume injected (as obtained from data published by Dyes, Caudle and Erickson2), the theoretical performance was determined. This prediction is compared with actual performance to date. Actual performance reflects an increasing percentage of tniscible product in the produced fluid to a particular level, which remains static for a short period and then increases again to a higher level. This performance suggests that the reservoir has performed as a layered system. It is concluded that the layered-system approach can be used to predict performance in a highly stratified reservoir and that economic consideration is a limiting factor in determining the amount of the reservoir which can be miscibly swept. INTRODUCTION The Paloma zone in the Paloma oil field, Kern County, Calif., is comprised of several sand members separated by siltstone, shale and/or cherty sections. Rapid facies changes within the zone indicate that the sediments were deposited in a "near shore" environment. As a result, the Paloma sands have a high permeability variation, and the Paloma zone is comprised of several sand reservoirs. Since one of the sand reservoirs containing crude oil was not receiving the benefits of gas injection, consideration was given to recompleting certain wells to isolate the sand in order to perform a miscible-phase injection project and, thereby, to increase its ultimate oil recovery. Before proceeding with the project, a review of the literature was made to determine the miscible-slug size that would be required to maintain miscibility throughout the displacement process. Although the review showed that considerable information was available for homogeneous systems, insufficient information was available to determine slug-size requirements and to predict performances for a highly stratified reservoir. This paper has been prepared to describe the method used to determine the miscible-slug-size requirement and to predict the performance of the miscible-phase injection project initiated in the Paloma oil field. Although the Paloma project has only been in operation since Aug., 1958, sufficient production history has been obtained so that a comparison can be made with the predicted performance. DESCRIPTION OF RESERVOIR The reservoir in which the miscible-phase project was initiated in the Paloma field is shown on Figs. 1 and 2. Reservoir properties are as follows: average sand thickness, 43 ft; porosity, 19.1 per cent; water saturation, 47.2 per cent: oil viscosity, 0.328 cp; reservoir temperature, 250°F; reservoir pressure, 3,500 psi: gas saturation, 21.2 per cent at time of miscible product injection; and total net hydrocarbon pore volume, 19,993,000 reservoir bbl. LIMITATIONS The method described herein utilizes a layer-system approach, where injection into and production from each layer is assumed to be proportional to the layer's fractional capacity. The effects of increased fluid mobility on injection and production rates, in layers which have had breakthrough, are not reflected in this method. There has been considerable controversy over the slug size necessary to maintain miscibility throughout the displacement process. In this method, the slug-size

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 S. A. Hovanessian

Analytical solutions are obtained for calculating the pressure distribution in rectangular fields due to injection and/or producing wells located anywhere within the field. The field is assumed to be homogeneous with either constant pressure or no-flow boundaries. These solutions are extended to include the calculation of average reservoir pressure and permeability from pressure build-up curves. Numerical examples are given to illustrate the application of equations to practical problems. INTRODUCTION Pressure distribution analyses are of considerable importance in the field of reservoir mechanics. A number of papers dealing with this subject have appeared in the petroleum literature. perrinel in 1956 presented a summary of the methods available for pressure build-up calculations and discussed the applicability of each in some detail. More recently (1958), Hazebroek, et a1,2 discussed a theoretical method for obtaining pressure distribution in a radial field. Also in 1958, Nisle3 gave a theoretical solution of the pressure distribution in a field extending to infinity with a partially penetrating line source at the origin. Mathews, et a1,4 applied the method of images to existing solutions for pressure distribution in radial fields and obtained expressions for calculating pressures in bounded reservoirs. However, a literature survey revealed that a general analytical solution for the pressure distribution in a rectangular field, with injection and/or production wells located anywhere within the field, was not available. Such a solution could be used to approximate the pressure distribution in a field the shape of which is nearly rectangular. Furthermore, it could be applied to pressure calculations pertaining to a single well taken from a system of wells which were drilled in a rectangular pattern. In this case the boundaries of the well drainage area are considered to be the perpendicular bisectors of the lines joining the well and its nearest neighbors. The basic problem considered in this paper consists of a rectangular field with either constant-pressure 01 no-flow boundaries and with either an injection or production well located anywhere within the field. The physical properties of the rock and fluids present are considered to be constant. Analytical expressions are obtained for the pressure distribution within the field during the production or injection periods. These expressions, which are infinite series, can be used to evaluate pressure at any point within the field as a function of time and position. For fields containing more than one well (production and/or injection), the solutions for each well acting independently are superposed. That is, at any point in the field, the pressures due to each injection 01. production well acting alone can be algebraically added to give the pressure resulting from all wells acting simultaneously. Evaluation of the series solutions given in this paper is readily accomplished by means of an electronic digital computer; and, in most cases, sufficient accuracy is obtained by setting the upper limits of summation at 20. These solutions are extended to include the calculation of average reservoir pressure and permeability from build-up data. A number of numerical examples are given to illustrate application of the solutions to actual field problems. ANALYTICAL SOLUTIONS The partial differential equation representing the pressure distribution in a homogeneous rectangular field (see Fig. I.), having a single source or sink of strength (2, can be written*

By V. J. Berry, J. D. Pendergrass

Well pressure transient tests provide a means for directly obtaining information about formation pressure and reservoir flow capacity. Such tests have also been proposed for determining presence and location of faults or other reservoir closures and for measuring oil in place. For mathematical convenience, most theoretical studies have considered the reservoirs to be homogeneous. Definitive information is not yet available to show whether the actual presence of nonuniformities will make pressure transient behavior different from that of a uniform reservoir. The conclusions reached from actual transient tests are questionable, therefore, insofar as they rely on the original assumption of homogeneity. One type of nonuniformity commonly assumed to exist is that of stratification. In most reservoirs the strata are thought to be in vertical communication. Equations for the transient flow of a single-phase, compressible fluid in a one-well, bounded, circular reservoir have been solved for several situations involving crossflow between multiple strata of various thicknesses and permeabilities. The results show that except for the very early flow period, which usually is too short to be analyzed, the transient performance observed at the well is substantially identical with that of a homogeneous reservoir having the same dimensions and having the same steady-state flow capacity. Thus, stratification does not adversely affect interpretation of well transient tests. This conclusion holds for all commonly encountered combinations of reservoir thickness and external radius. Deviations are observed for unusually thick reservoirs whose outer radii are relatively small. The results of these studies also show that the presence and the amount of stratification cannot be simply diagnosed from reservoir pressure transient data when there is crossflow between strata. INTRODUCTION The last decade has brought wide acceptance of the transient well pressure test for determining reservoir parameters. Following the original work of Hurst and van Everdingen,1 the mathematical theory was thoroughly explored. Numerous authors have suggested how to determine static reservoir pressure, permeability-thickness product, original oil in place and reservoir limits for different reservoir geometry. The same mathematical techniques have been used to predict the transient performance of a reservoir over a long period of time. Most of the theoretical work has been for homogeneous, isotropic systems. Some results have also been presented for a homogeneous, anisotropic reservoir.2 Petroleum reservoirs are not homogeneous. The deposition process seems to favor creation of a stratified formation. This concept is sufficiently well accepted so that the most natural extension of the transient flow theory beyond the homogeneous case is to a stratified formation. Results for a stratified reservoir with no vertical communication between layers can be obtained from the results for a homogeneous reservoir. Lefkovitz, et a1,3 have given a thorough treatment of the two-layer case. Two recent papers4,5 have treated the case of a two-layered reservoir with vertical communication or crossflow, between layers. Russell and Prats4 find that, after a relatively short time, the two-layered reservoir with crossflow exhibits a simple exponential pressure decline. From this time forward, the behavior is not distinguishable from the behavior of a homogeneous reservoir having the same steady-state flow capacity. The results of Katz and Tek 5 are equivalent. Russell and Prats4 also speculate that a multi-layered reservoir with crossflow will behave as a homogeneous system after long enough production time "... providing the contrast in kh between layers is not too great". They also suggest that, at intermediate times, "... the relative positions of the layers with respect to each other will have a great influence on the production behavior and on the time at which the previously mentioned large-time approximation might be valid". Katz and Tek5 remark upon the mathematical difficulty of treating a reservoir having many layers or strata. If relative position of the layers is also

By L. W. Emery

Undergrorlnd combustion operations were initiated in a 60-acre Bartlesville sand "shoe-string" reservoir in Allen Connty, Kans., in 1956. Tests in separate patterns were conducted using various co~nbinations of air and recycle gas to propagate combustion fronts from the injection toward the producing wells. These patterns were made up of 6 injection and 20 prodrrcing wells Gas and liquid prorluctiorz from each pattern was measured on an individual-well basis, and comparisons were made between the three patterns to ascertain the relative effects of injected gas composition on production behavior. Breakthrough of the combustion front at a well was characterized by an increase in water production from the well followed by an increase in bottomhole temperatrrre to approximately 250" F. After burning fronts had broken through at five producing wells, operations were terminated in 1960. From the total project approximately 79,000 bbl of oil were produced during thermal operations at a cumulative produced GOR of 23 Mcf/bbl. No appreciable change in the character of the produced crude was observed. Combustion in the reservoir was maintained with injected gas compositions ranging fronz 6 per cent oxygen in recycle gas to 100 per cent air. lnjectiotz of large quantities of recycle gas resulted in higher producing GOR's from offset wells than were measured from a pattern into ~vhich straight air ~vas injected. The air required to move the combustion front through I acre-ft of reservoir was computed to be 20 MMscf. This valrre was found to be relatively independent of the quantities of recycle gas injected. The recovery efficiency from the swept area was esti~izated to be about 59 per cent. Areas swept were similar in shape to tlzose obtained with a laboratory potentiometric model. Samples of sund taken from behind the burning front by coring indicated almost total oil removal from the sand. Petrographic analysis of the core samples indicated that the sand had been heated to peuk temperature of rlbout 1,200" F. No rignificant difference in peak temperature was forrnd in two areas where compositions of injected gas were quite different. Compression costs for thermal recovery were estimated to be $1,20/bhl of produced oil. INTRODUCTION The use of the "forward combustion" process as an oil recovery method has received a great deal of attention. This method involves ignition of the formation in an injection well, followed by propagation of a combustion front through the reservoir. Combustion is maintained by the injection of an oxygen-containing gas to react with reservoir hydrocarbons. As the flame front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase and recondensed in the cooler part of the formation. In turn, the condensed fluids push oil into the producing wellbores. Completed field tests of the process were first reported by Kuhn and Koch,' and by Grant and Szasz.' Results from other tests have since been reported by Walter,3 by Moss, White and McNeil,' and by Gates and Ramey." ach of these tests essentially utilized a single injection well surrounded by four or more producing wells. Sinclair Research, Inc., elected to do field experimental work using a number of test patterns in a single field in order that comparisons between various operating schemes could be made. The site selected and purchased in 1955 for this experimental work was a 60-acre Bartlesville sand reservoir located in Allen County, Kans. Combustion operations were initiated in mid-1956. Between that time and termination of the project in mid-1960, combustion fronts were propagated from injection wells to producers in three separate well patterns, using different mixtures of air and recycle gas. The test was terminated before sweep of the three patterns was complete so that information about the effect of combustion on the swept areas could be obtained by coring. Results from the test in the form of injection and producing well performance have been carefully recorded, and these form the general basis for this paper. DESCRIPTION OF RESERVOIR The reservoir in which the combustion tests were conducted is a Bartlesville sand "shoe-string", typical of a number of small reservoirs in Southeastern Kansas. Average reservoir characteristics are shown in Table 1. Fig. 1 is an isopachous map of the producing sand showing the reservoir to be approximately 400-ft wide and 2,500-ft long. Maximum net productive sand thickness is 21 ft. Fig. 2 shows a typical core analysis obtained by coring with water-base mud. The reservoir has no appreciable dip and is closed on the sides by degradation of sand into shale. The main body of sand is heavily laminated with shale stringers, which are not continuous between wells. The main reservoir is overlain by 30 to 40 ft of laminated low-permeability sand and shale streaks. No information is available on the original properties of

By N. C. J. Ros

The paper deals with pressure gradients occurring in flowing and gas-lift wells, a knowledge of which can be applied to the determination of optimum flow-string dimensions and to the design of gas-lift installations. The study is based on a pressure-balance equation for the pressure gradient. It appears that a pressure-gradient correlation of general validity must essentially consist of two parts-—one part being a correlation for liquid hold-up and the other part being one for wall friction. Dimensional analysis indicates that both liquid hold-up and wall friction are related to nine dimensionless groups. It is shown that in the field of interest only four groups are really important. On the basis of these four groups a restricted experimental program could be selected that nevertheless covered practically all conditions encountered in oil wells. This experimental program has been carried out in a laboratory installation. Three essentially different flow regimes were found. The pressure gradients in these regious are presented in the form of a set of correlations. Comparison of these correlations with a few available oilfield data showed excellent agreement. INTRODUCTION Prediction of the pressure drop in the flow string of a well is a widely known problem in oilfield practice. Accurate data on the pressure gradient of a simultaneous flow of gas and liquid in a vertical pipe are especially useful for the determination of optimum flow-string dimensions. It is well known that with moderate gas and liquid flows such a vertical string acts as a "negative restriction". The pressure drop decreases (1) when the throughput through a given pipe increases, and (2) when at a given throughput the cross-sectional area is decreased. The reason is that, with increasing velocities, the flow becomes more agitated so that the gas slips relatively more slowly through the liquid. With the resulting increase in gas content in the string, the static head decreases. When the area becomes very small, however, the high velocities entail great wall friction, which causes an increase in pressure drop. For a given flow, therefore, minimal pressure drop is obtained by using a certain cross section. This means that, in principle, each well can be provided with an optimum flow string for minimum pressure drop and, hence, maximum possible production rate. The procedure for the selection of the optimum string has been discussed by Gilbert.' A necessary tool in the procedure, however, is accurate knowledge of the pressure gradient to be expected for various values of the governing variables. Another application of pressure-gradient data lies in the field of gas-lift practice: they provide a means of determining the optimum gas-injection rate, optimum injection pressure and optimum injection depth. Much work has already been done in the study of the pressure gradient of vertical gas-liquid flow. Poett-mann and Carpenter2 presented a pressure-gradient correlation based on measurements in wells. This correlation has been found to provide accurate predictions in high-pressure wells and in high-production wells for flow through both tubing and annuli.2-5 However, when their method is checked on low pressure-low production wells or on wells with viscous crudes, serious discrepancies are found. As we shall see in the next section, this is due to the fact that their correlation factor, representing all irreversible energy losses, is given as a function of only one correlation group. Some important variables, such as gas-liquid ratio and liquid viscosity, are not incorporated in this group so that their specific effects are not accounted for. To study also the mechanism of vertical gas-liquid flow outside the ranges covered by the Poettmann-Carpenter publication and extensions, a laboratory investigation has been carried out. This study is founded on a pressure-gradient equation that is based on a pressure balance. To reduce the number of test runs required, a dimensional analysis has been carried out, followed by a selection of relevant dimensionless groups. These groups guided a subsequent experimental study, and with their aid the experimental program could be minimized while still covering the majority of the situations encountered in oilfield practice. In this paper the choice of a formula for the pressure gradient is discussed first. This is followed by a brief description of the experimental setup. Subsequently, the dimensional analysis is discussed and the relevant dimensionless groups are selected, resulting in the experimental program required. The general relationships of pressure gradient and liquid hold-up are then described; various flow patterns and a certain flow instability (so-called "heading") are discussed and a set of correlations is presented which shows a good agreement with the measurements and a few available field

By R. L. Perrine

This paper shows how stability theory can be used to optimize solvent recovery of oil. Application of the theory leads to definition of the limiting conditions required for stable displacement to occur. One of these conditions is that the size of a solvent "slug" must at least equal a prescribed minimum. The criteria to be satisfied are miscibil-ity and stability. Stability implies that no viscous fingering will occur and that mixing will be caused only by a dispersion process. Miscibility implies that complete recovery will be obtained from the swept regions. Thus, for any specified set of reservoir conditions, an optimum use of solvent is defined. INTRODUCTION The early outlook for solvent flooding as a means to increase oil recovery was very favorable. Some laboratory results indicated that a "slug" containing perhaps 2 to 3 per cent of a hydrocarbon pore volume could be successful. However, other data suggested as much as 30 per cent would be required. The difference is of considerable economic importance. A theoretical explanation for these divergent results has been advanced by the author.1,2 Stability theory defines conditions for two distinct flow regimes. In stable displacement, solvent and oil become mixed by a dispersion process, and the solvent requirement is small. On the other hand, unstable displacement can degenerate into viscous fingering. The practical result is a considerable increase in the extent of mixing and, thus, in the solvent required. The present paper shows how to use stability theory to optimize solvent recovery of oil. Oil is usually considered to be displaced by a small amount, or a slug, of solvent. Gas in turn follows solvent. Any fingering will permit contact of nearly solvent-free gas and oil, leading to immiscibility and reduced recovery. Thus, our optimum is defined by two restrictions — stable flow and miscibility. Although an immiscible solvent process could prove more profitable than any miscible process, such processes lie outside the selected definition of an optimum. The usefulness of computations based on the present work will depend on the validity of the general theory. Because the important idealizations required tend to minimize the predicted solvent requirement, the results should state correctly the limiting conditions required. Therefore, theoretical calculations should have practical value. This paper first reviews several principles that result from stability theory. The following section describes use of these ideas in relatively simple systems which symmetry makes appear one-dimensional. Of particular importance is the step-by-step design of 3 stable minimum-solvent slug. The use of stability theory in cases of multi-dimensional displacement is discussed and, finally, suggestions are given for simplified practical application of the theoretical results. PRINCIPLES OF STABILITY THEORY The use of perturbation methods to derive stability conditions has been shown by the author.l The discussion is presented in terms of the fractional concentration of solvent in the single hydrocarbon phase. Quite generally, any disturbance can be represented as a Fourier series, of which only the least stable term need be considered. Denoting this term at time t = 0 by cn,

By H. D. Yoo, M. R. Tek, D. L. Katz

The behavior of gas-storage reservoirs subject to water drive is investigated through analog simulation on an electronic differential analyzer. The simulation technique developed on an LM-10 computer permits the prediction of reservoir volume or pressure resulting from the movement of water in the surrounding aquifer. The method developed on the analog computer consists of setting up an appropriate transfer-function circuit and feeding the arbitrary time-varying boundary conditions as an input signal. The input may be specified as gas reservoir pressure, pore volume or the water flux. Several cases studied include an isolated gas reservoir on a limited aquifer, interference among three reservoirs adjacent to a common aquifer and the growth of gas-storage volume on an aquifer. It is concluded that the method developed on an electronic differential analyzer provides an excellent technique to simulate and investigate the behavior of gas reservoirs subject to water drive. The agreement between the reservoir performance as predicted from the simulation technique and as measured from actual field data is found to be better than the range usually encountered in predicting water-drive behavior. INTRODUCTION It is generally known that some gas-storage reservoirs are located on top of blanket sands of large extent, saturated with brine called aquifers. Because the volume of the body of water associated with aquifers is usually very large and water is compressible, the cyclic pressure variations encountered in normal storage service inevitably cause unsteady, compressible flow conditions in the adjacent aquifers. The solution to radial diffusivity equation for a limited aquifer for constant terminal conditions has been know111 since the early 1930's. The solution for the constant terminal conditions for an infinite aquifer was puhlished in 1949.2 The proldem of handling the unsteady flow of water through a porous medium with time-varying boundary conditions has been studied by many investigators through the use of graphical,3 digital4 computing techniques or on an R-C type of electrical analog called a reservoir analyzer.5 The digital computer performs a series of arithmetic calculations to superimpose the solutions of constant terminal conditions, while the R-C type of analyzer permits the simulation of the continuous medium by means of "lumped-cell" electrical resistance and capacitance elements based on finite-difference approximations. The graphical method is based on the convolution integral solution of the partial differential equalion. The purpose of the present study was to simulate the behavior of an aquifer on a general-purpose, electronic differential analyzer to obtain the influx of water or the change in pressure at the reservoir continuous1r and quickly. The use of the electronic differential analyzer was based on well known solutions for the constant terminal conditions already available in the literature.2,6-8 THEORETICAL DEVELOPMENT MODEL. The gas reservoir (R) is described in relation to an aquifer in Fig. 1. The permeability and porosity of an aquifer sand, the viscosity of the brine, and the compressibility of the brine and the formation are assumed to be constant. It further is assumed that the blanket sand is not thick enough to give appreciable vertical pressure distribution. The model in consideration is circular with radial geometry. The aquifer radius can be of limited or infinite extent. The gas-water interface is assumed to be fully segregated, and the displacement of the gas by

By P. E. Baker

Approximate solutions are presented for the heat-flow equations in a loss-free linear system with a moving source and with heat transfer by convection and conduction, representing in situ combustion in an oil reservoir. The rate of heat generation and the velocity vf of the source (flame front) are assumed constant. A cold front or convection front caused by the cooling effect of in-jected air develops behind the flame front and moves with velocity vc, which is determined by the thermal properties of the air and rock. Because of thermal conduction the fronts are not sharp; i.e., they do not represent temperature discontinuities. The following solutions apply after a sufficiently long time. where Q is heat generated per unit volume of formation, k is thermal conductivity, CF is volumetric heat capacity, a is thermal diffusivity, and subscripts 1 and 2 refer to values behind and ahead of the flame front, respectively. The parameter CF2 is contained in v . Te is initial reservoir temper- . C2 e ature and air-injection temperature. Other mathematically possible solutions are presented. The different solutions result from the use of different values for various physical parameters. The solution just presented is thought to be most commonly applicable. INTRODUCTION In this paper, approximate analytical expressions are derived for the temperature profile developed in the underground combustion method of oil recovery. A qualitative description of the process may be found in a paper by Tadema.1 The flame front is regarded as a moving source with a constant rate of heat generation. Heat transfer is by conduction and convection. Previous publications, of which a few are cited here, 2-5 have also attacked this problem and have presented solutions pertaining to different methods of approximating or idealizing the process studied. The present paper also treats an idealized case: linear flow is assumed for heat and fluids, and vaporization and condensation effects are neglected. The results obtained are similar to those derived by Bailey and Larkin,3 but are obtained by an entirely different procedure. The method used here should give considerably greater insight into the physics of the in situ combustion process. Four types of cases are considered, distinguished by thermal and fluid-flow properties of the system. A constant rate of heat generation and constant rate of flame-front advance are assumed. Flame-front temperature is determined for each case using a heat balance which involves the entire temperature profile. The applicability of the different solutions to real systems is discussed after they are all presented. Again, this is done for the purpose of understanding better the physics of the process. Consideration of the nonapplicable solutions and of the reasons for their not applying is significant. In a paper by Cooperman,4 the only case solved is one considered here as not applicable. Fluid flow, as such (pressures, saturations, oil recovery, etc.), is not treated. This subject has been covered in a paper by Wilson, et al.6 Also, we are not concerned with requirements for continuance of the process; we only attempt to determine the temperature history for cases where the process has continued for some time. PHYSICAL SYSTEM AND APPROACH TO PROBLEM A linear flow (heat and fluid) system is assumed and the x-axis is taken parallel to the flow direction. The combustion zone (flame front) is vertical,

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing

By G. Thodos, R. B. Grieves

A method has been developed For predicting the critical temperatures and critical pressures of binary mixtures of carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen with the normal paraffin hydrocarbons. For carbon-dioxide and hydrogen-sulfide systems, relations are presented that take into account the peculiar behavior of mixtures with closely boiling components, such as carbon dioxide-ethane and hydrogen sulfide-propane mixtures which exhibit minimum critical temperature points. For hydrogen, nitrogen and carbon-monoxide systems, the extreme critical behavior caused by wide differences in pure component properties is established In addition, those fixed gas-paraffin systems which resemble paraffin-paraffin systems are also accounted for. For a mixture of known composition, the pure component critical temperatures, critical pressures and normal boiling points are all that are required to determine its critical point. Graphical relations are presented relating Tc and PC , of the mixture to the pure component properties. From the treatment of 12 carbon-dioxide and hydrogen-sulfide systems reported in the literature (74 mixtures), the expected error lor the critical temperature is approximately 1.5 per cent, and for the critical pressure, approximately 2 per cent. From the treatment of six hydrogen, nitrogen and carbon-monoxide systems reported in the literature (30 mixtures), the expected error for both the critical temperature and critical pressure is approximately 2.5 per cent. The relationships, which have been developed with only normal paraffins as the hydrocarbon components, may be extended to those isoparaffins and olefins which fall within the allowable volatility ranges. INTRODUCTION Many of the fixed gases — carbon dioxide, hydrogen sulfide, nitrogen, hydrogen, carbon monoxide and oxygen — occur in natural mixtures with hydrocarbons. Carbon dioxide and hydrogen sulfide are frequent components of the fluids produced from underground petroleum reservoirs. Nitrogen, carbon dioxide and hydrogen sulfide are present in varying quantities in most natural gases and gas-condensate well effluents. Hydrogen mixtures are of considerable interest in many phases of refining processes of petroleum. The determination of the critical temperatures and critical pressures of such mixtures is of value in vapor-liquid equilibrium studies, for the prediction of the characteristics of underground reservoirs, and for reduced-state correlations of PVT, transport and thermodynamic properties. The accurate estimation of the critical point for binary mixtures is an important initial step toward a complete analysis for the establishment of the critical temperatures and pressures of multi-component mixtures. Methods for predicting the critical temperatures and critical pressures of binary hydrocarbon systems have already been presented in the literature.8 It is possible to apply these existing methods to fixed gas-paraffin mixtures, but due to their unusual critical behavior, values calculated deviate considerably from experimental values. For systems containing trace quantities of the fixed gases, these methods are acceptable; however, for systems containing more than 5 mol per cent of the fixed gases, these utterly fail to produce reasonable critical values. Consequently, in this study a method has been developed for handling such binary mixtures over the entire composition range. CARBON-DIOXIDE AND HYDROGEN-SULFIDE SYSTEMS The critical behavior and the vapor pressure behavior of mixtures of carbon dioxide and hydrogen sulfide with paraffinic hydrocarbons may be quite similar or quite dissimilar to that of paraffin-paraffin mixtures, depending on the volatilities of the components involved. The critical temperature and normal boiling point of carbon dioxide are very close to the corresponding values for ethane, while its critical pressure is considerably higher than that of ethane. The normal boiling point of

By R. L. Perrine

A stability theory is developed for miscible liquid-liquid displacement within a porous medium. In the usual case considered, a high-density high-viscosity "oil" is displaced downdip by a low-density low-viscosity "solvent'! Perturbation methods are used to find the conditions under which the spreading mechanism changes from the stable dispersion process to unstable viscous fingering. We find that instability is conditional, that there is a dependence on the shape of a disturbance leading to a"diameter" effect and that very difficult experimental scaling problems may result. A useful consequence is the definition of a minimum "slug size" for stable miscible displacement. This should make possible optimum use of the solvent process for oil recovery. The results may apply to many situations in which one fluid displaces another of somewhat different fluid properties within a porous medium. INTRODUCTION The process of miscible liquid-liquid displacement looked very favorable when it first received serious consideration as a means to increase petroleum recovery. Some early laboratory results were interpreted to mean that only a small "slug" of solvent was needed, perhaps 2 to 3 per cent of a hydrocarbon pore volume. This "slug" could recover the oil in an entire reservoir provided the fluids remained miscible.1,2 Only a slow growth of the "mixing zone", or gradual transition from oil to solvent, should occur. This would follow naturally from mixing by a dispersion mechanism such as that described by Scheidegger.3 And, indeed, laboratory cores have shown this kind of behavior provided solvent viscosity was at least as great as that of the oi1.2,4-6 The same kind of behavior has also been observed with solvent viscosity lower than that of the oil after the distance traversed became large. Not all laboratory results have been this favorable, however.7-9 Mixing zone growth may become proportional to the distance traveled rather than the square root of this quantity. Such behavior accompanies high rates in short systems of large diameter, provided the viscosity ratio is adverse. If this more rapid spreading were to persist in a reservoir, the solvent requirement for miscible "slug" displacement could exceed 30 per cent of a pore volume. There is considerable economic importance in the difference between 3 and 30 per cent solvent. A 1ogical conclusion is that two different kinds of flow behavior are possible, with each kind leading to a different result. One kind gives efficient displacement, the other does not. In the efficient case, the solvent bank spreads out only by the mechanism we have termed dispersion. Under these conditions, displacement is as near piston-like as possible. In the second kind of flow, found only with adverse viscosity ratios, viscous fingering occurs. That is, permeability variations cause a small finger of low-viscosity solvent-rich fluid to move ahead of its average position within the mixing zone. This creates a path of low resistance to flow, and an even greater amount of solvent-rich fluid follows. Thus, the process is autocatalytic. Once started, the fingering mechanism rapidly becomes dominant. The question to be answered is this. Under what conditions will each of these two different kinds of flow occur? In particular, what conditions denote the transition from dispersion to viscous fingering as the solvent spreading mechanism? The answer may tell us whether or not the desirable, near piston-like miscible displacement process is practical. An answer to this problem can be obtained from theory by the use of perturbation methods. The procedure is as follows. We first formulate a mathematical representation of the system. Then a small disturbance (or perturbation) in solvent concentration profile is introduced and observed to see what happens. Unstable flow is indicated when a small disturbance will grow larger. This will lead to eventual viscous fingering. If on the other hand the disturbance dies out, the displacement is stable. In this latter case, with dispersion as the spreading mechanism, near piston-like displacement is possible. Thus, this paper presents the development of a stability theory for miscible liquid-liquid displace-

The performance of a two-horizon depletion-type reservoir produced through combination wells is analyzed. By introducing some simplifying approximations, it has been possible to obtain formulas which are easy to handle. Only ordinary differential equations are used, and the development of the analysis can be followed without difficulty by the non-specialist. A rigorous analysis has been made of this problem. It has been found that the approximations introduced in the simplified analysis are fully justified and that errors seldom will be more than 2 per cent. The report shows the effect of the rate of withdrawal upon the relative depletion of the two horizons. Numerical examples are given. INTRODUCTION In its simplest form an oil reservoir consists of a continuous, homogeneous body of porous and permeable rock, enclosed at the top and the bottom by impermeable material. The flow of fluid through such a reservoir has been the subject of numerous papers, and the mathematical analysis of well performance, pressure build-up, etc., has reached a high degree of refinement. One of the basic assump tions introduced into these analyses is the homogeneity of the reservoir. In non-homogeneous reservoirs, these computations still can be carried out satisfactorily, provided small irregularities in the physical properties of the rock are distributed in a random fashion. The flow equations are derived as solutions of the continuity equation, and their solutions are dependent upon the assumption of many constant factors in the reservoir: permeability to oil, compressibility of the reservoir fluid, viscosity saturation, etc. When these axe assumed to be constant, the continuity equation can be solved rigorously. When these assumptions do not hold and corrections are introduced, it becomes very difficult and often impossible to solve analytically. The present analysis is concerned with production from a two-layer depletion-type reservoir in which the layers have different permeabilities. Instead of using the continuity equation, results were derived from some simple physical ideas. In a single-layer reservoir in which the oil is at a pressure above the bubble point, the rate of production often has been assumed to be directly proportional to the pressure drawdown. In the simplified analysis it has also been assumed that the drop in average reservoir pressure is directly proportional with the cumulative withdrawal. These considerations in the past have been the basis of the analysis of reservoir performance and have been applied to a variety of problems, often with marked success. In the problem under investigation these relationships are shown to be compatible with the equation of continuity, but the approach has the advantage of leading to ordinary differential equations in which changes in constants can be introduced without making the analytical solution impossible. It has the additional advantage that it can be understood more readily by the large group of engineers who are not reservoir specialists. For that reason, detailed steps in the mathematical developments have been shown in full. SCOPE OF ANALYSIS The analysis in this paper has been applied to a two-layer depletion-type reservoir in which both horizons are produced simultaneously in the same well. Because of the difference in permeability, it is obvious that the depletion of the two zones will proceed at different rates. For various problems it may be of interest to know the degree of depletion of the two zones at any arbitrary time and to in-

By M. R. Tek, K. H. Coats, D. L. Katz

The nature and the limits of validity of Darcy's law US applied to the flow of natural gas through reservoirs has been considered in order to resolve some controversial aspects of the effect of turbulence on pressure drops. The equivalence between various concepts and viewpoints advanced in the past by several investigators to explain how and why a gas well does not necessarily perform according to Darcy's law is shown. Starting with generalized equations of flow of fluids through porous media, a partial differential equation has been derived which accurately represents the flow at all rates. This equation has been numerically solved using an IBM 704 digital computer. The results permit plots of unsteady radial pressure distribution curves from which specific isochronal backpressure curves may be constructed. These back-pressure curves show the effect of the ,8 factor on the slope of the back-pressure curve. The calculations further indicate that the drainage radius for a gas well in turbulent flow propagates at a rate dependent upon the rate of production at the wellbore. This is quite different from the case with liquid flow or natural-gas flow in laminar regime. Additionally, the effect of reservoir inhomogeneities and crossflow between layers of different permeability on the back-pressure performance of gas wells has been conridered. In light of the current numerical results the significance and linlitation of the rate of flow function Y proposed by Smith" has been discussed. INTRODUCTION The relationship between the pressure drop and flow rate in problems of fluid flow through porous media is known to be affected by the nature of flow through the porous matrix. It has been observed by many that, for a range of flow rates, the pressure drop remains proportional to the rate of flow. When some flow rate is reached, however, it is usually observed that the pressure drop gradually begins to increase more than proportionally to the flow rate. It is well known that this phenomenon was first observed by Osborne Reynolds in 1901 in experimenting with flow through pipes. In his classical experiments, Reynolds made visual observations on the condition of streamlines evidenced by injecting a dye into water flowing through glass tubes. In these experiments, the abrupt transition between steady, "streamline, laminar" flow and unsteady random turbulent flow was found to be a function of the dimensionless group (Dvp/u), now known as the Reynolds number. During these experiments, in addition to observations on the nature of flow regimes, the proportionality between flow rate and pressure drop in laminar flow was contrasted with the nonlinearity between these variables in turbulent flow. Fancher and Lewis' reported data on various consolidated and unconsolidated sands in 1933. Their conclusions were that ". . . the flow of fluids through these porous materials closely resembles that through pipes; that there is a condition of flow in porous systems which resembles viscous flow, another which corresponds to turbulent; that the change from one type to the other takes place at a definite and reproducible condition for each system". In 1947, Brownell and Katz published a method to predict the laminar and turbulent flow behavior from the particle size, bed porosity and the particle sphericity, employing the friction factor-Reynolds number charts for pipes. Several investigators have verified the work of Fancher, Lewis and Barnes and presented their data as friction factor-vs-Reynolds number plots.' The equation which would represent the pressure gradient over the whole range of velocity must have an added term over that represented by Darcy's law. Accordingly, the pressure gradient necessary to sustain flow at the velocity (v) through a porous medium may be represented by the following equation, suggested by Forscheimer." The nature and the range of validity of Darcy's law has been the subject of studies by many investigators over the past years. While everyone seemed to agree on the need for a quadratic correction term to Darcy's law to make it effective over the range of velocities, the concept of inception of turbulence and the use of the term "turbulent flow" remained controversial. Some fluid dynamicists define turbulence as a flow

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

By E. F. Johnson, F. H. Brinkman, H. J. Welge, S. P. Ewing

This paper presents a method for predicting the manrler in which oil will be displaced from a porous body by enriched gas. The calculations apply to a gas rich enough to give a partially, but not a completely, misci-ble displacement. The method — a three-component, two-phase analysis — takes into account condensation of some of the intermediate hydrocarbons from the injected gas into the oil, as well as enhanced volatility of heavier hydrocarbons at elevated pressures and temperatures. The condensation swells the oil and decreases its viscosity, thus aiding in its recovery. The calculations have been extended to apply to actual crude oil-natural gas systems by arranging the components into three groups according to their volatility. As an approximation, each group is then treated as a single component in the analysis. The influence of an angle of dip for an inclined displacement is also taken into account. The recovery predictions are corroborated by experiments which used both consolidated sand cores and un-consolidated glass beads. In some of these tests, actual live crude oil was displaced by a multicomponent gas typical of enriched gases used in oil fields. INTRODUCTION This paper presents a method for predicting the amount of oil that can be displaced from a homogeneous, linear, porous body at various stages during the injection of enriched, or "wet", gas. The porous body can be in either a horizontal or an inclined position. 'This type of displacement is sometimes known as condensing gas drive The method is developed especially for the case in which the injected gas is enriched enough to be partially, but not completely, miscible with the reservoir oil. The need for a calcula-tive procedure for this type of operation is emphasized by the number of field projects where completely miscible drives are not practical, but where near-miscible conditions are feasible. The factors taken into account in the predictive calculations include: (1) the condensation of gas components into the oil, with a resulting increase in oil volume; (2) the lowering of oil viscosity by the addition of lighter ends from the gas; (3) the increase in oil volatility at high temperatures and pressures; and (4) the physical displacement of the oil by the gas. The techniques developed in the paper can be extended to other nonequilibrium displacement processes. Other such processes that we have analyzed include a displacement by lean gas which stripped intermediates from the oil, and a water flood in which the water con. tained in solution a substance somewhat soluble in the oil. ANALYSIS OF ENRICHED-GAS DRIVE GENERAL PRINCIPLES Our method for predicting the amount of oil that can be displaced by an enriched gas uses an analogy between a three-component and a multicomponent system.' The predictive method is based on these assumptions: (1) constant, or nearly constant, pressure; (2) complete equilibrium by diffusion perpendicular to the main direction of flow, but no significant mixing along the direction of flow; (3) constant injection velocity; and (4) flow in a linear porous body. The composition of a liquid or a vapor with respect to three components can be plotted on a three-component, or ternary, diagram like that in Fig. 1. Let Point A represent the composition of the oil originally in place. In this case, Oil A is undersaturated with gas. If Point A lay on the equilibrium Curve BF, the oil would be saturated. In the extreme case where the original oil contained no intermediates or dissolved gas, Point A would lie at the lower left-hand corner of the ternary diagram. In a displacement of Oil A by Gas D, there will be a progressive change in the composition of the oil phase as more and more gas is brought into equilibrium with the oil. The end result of this progressive change is an oil having the composition represented by Point F. This oil is richer in intermediate hydrocarbon and methane than the original oil and, therefore, has a greater forma-

By F. J. Fayers, J. W. Sheldon

A theoretical discussion is presented of the behavior of the interface between two fluids of different physical properties when displacement is occurring along a thin tilted bed. An approximate equation of motion is deduced for the motion of the interface, suitable boundary conditions are derived and the conditions necessary for the approximation to be valid are determined. Steady-state behavior and the circumstances under which a steady-state interface can develop are discussed. Finite-difference solutions of the equation of motion for two characteristic water-drive problems are presenled. The theory of displacement in a thin bed may be regarded as being the porous-medium-equivalent of the classical "shallow-water wave" theory of hydrodynamics. INTRODUCTION The solution of moving interface problems are important in predicting the behavior of oil reservoirs under various conditions of production. Kidder4,5 has obtained analytic solutions for the fingering problem and for another particular case of the motion of the interface between two immiscible fluids of unequal density in a porous medium. In the present analysis, we analyze the problem of the motion of an interface along a slightly dipping bed whose length is long compared with its thickness. The requirement of a thin bed is necessary in order to give a simple differential equation for the motion. The two fluids are of unequal density and mobility. In the derivation of the equation of motion, it is assumed that the lighter fluid is being injected at constant rate at the high end of the porous layer, while heavy fluid is produced at the low end of the layer. The analysis applies equally well, however, to the reverse procedure in which the heavier fluid is injected at the low end of the tilted layer. It will be shown that the behavior of the system is critically dependent on which fluid has the higher mobility and that the interface can in certain circumstances achieve a simple progressive steady-state shape. The fluids may be considered to be immiscible or miscible. It is assumed, however, that the thickness of the interface region is negligible in comparison with other dimensions of the system, and that the strength of the saturation or concentration discontinuity remains constant during the motion. In conjunction with Buckley-Leverett theory and areal sweep-out calculations, the tilted interface problem is one of the fundamental problems of reservoir engineering; it defines the very important underpassing and overpassing phenomena associated with water-drives, gas-drives and secondary recovery operations. The tilted interface problem has been dealt with previously by Dietz.2 In the present paper, we endeavor to derive the pertinent mathematical results more rigorously, and in doing so we obtain modifications in both the differential equation and boundary conditions for those given by Dietz. We also attempt to show the connection between these results and some of the recent work on the stability of interfaces. Finally, we attempt to demonstrate the practical importance of the thin bed theory by showing the difference in behavior of two characteristic waterflood problems. In the favorable mobility case recovery is improved by increasing production rate, while in the unfavorable case undercutting or overcutting is minimized by reducing the displacement rate. PROPERTIES OF THE EQUATION OF MOTION Fig. 1 illustrates a porous layer of constant thickness H and tilt 8 with Fluid 1 of mobility ?l and density pl below Fluid 2 of mobility ?2 and density p2. The fluids are incompressible and flowing with volumetric flux q along the layer. Both the slope a of the interface with respect to the bed and the curvature of the interface are assumed to be small. What is meant by "small" is indicated in the derivation given in Appendix A of an approximate equation of motion for the interface. Other symbols are defined in the "Nomen-

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.