Search Documents

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.

By G. Thodos, E. M. Phillips

A reduced density correlation for the gaseous and liquid regions of ethane has been developed from all the experimental PVT data available in the literature for this substance. Saturated vapor and liquid densities reported for methane, propane and n-butane were utilized to establish the rectilinear-diameter line for ethane. This line was then used with the saturated-vapor data for ethane to produce the saturated envelope for this substance. Most of the saturated liquid data reported for ethane were found to be inconsistent with the saturated-liquid curve resulting from this study. The correlation covers the regions 0.65 < TR <6.0 and 0 < PR < 15. For regions for which no experimental data have been reported, the isobars of the correlation were extended by the use of the theorem of corresponding states. This correlation is presented both on rectilinear and log-log coordinates. INTRODUCTION At the present time, there is considerable interest in the exact PVT behavior of pure hydrocarbons. Matschke and Thodos12 recently have developed a comprehensive compilation of experimental vapor and liquid density data for methane. Similar studies have been conducted for other pure substances.4, 7 In the present study the density data available in the literature for gaseous and liquid ethane have been critically reviewed and utilized to develop a reduced density correlation for this substance. SATURATED LIQUID AND VAPOR DENSITIES The rectilinear-diameter rule proposed by Cailletet and Mathias5 states that the sum of the saturated liquid and vapor densities of a pure substance is a linear function of temperature which passes through the critical point. This law, although empirical, has been well established by continued verification since its conception in 1886. This rule has been applied to the saturated density data of methane,6 propane 17 and n-butane8 to produce the straight lines presented in Fig. 1 which pass through the critical point TR = 1.00, Pr = 2.00. The rectilinear-diameter relationship for ethane was established from the lines for the other substances in Fig. 1 by plotting the quantity p~, + PR against zc for different reduced temperatures. A critical analysis of the saturated vapor and liquid densities reported for ethane by Kuenen,10 Maass and McIntosh,11 and Sage, Webster and Lacey18 revealed that the values of PR1 + pr resulting from these studies were inconsistent with the rectilinear-diameter line obtained from the reliable density data of the other substances, as shown in Fig. 1. In the analysis of these data, reduced quantities were calculated using the critical con-

By J. H. Moran, E. E. Finklea

The pressure build-up technique is a recognized method of determining permeability from conventional drillstem tests. In this paper an effort is made to extend such techniques to the interpretation of data obtained from the wireline formation tester. Such a study is necessary because of the differences, for this case, in the magnitude of the flow parameters (rate of flow, amount of recovered fluids) and in the flow geometry (flow through a perforation vs flow across the face of the wellbore, etc.) involved in the solution of the equations of flow for compressible fluids. The perforation is replaced by a spherical hole, and the effect of the borehole is neglected, so that the flow can be considered to be radial in a spherical co-ordinate system. Arguments are presented to justify this idealization. Assuming single-phase flow, general relations between pressure and flow rate are developed for a homogeneous medium. The study is then extended to permeable beds of finite thickness. It is shown that the early stages of pressure build-up tend towards spherical flow, while the later stages tend towards cylindrical flow. The thinner the bed, the more quickly flow approaches the cylindrical model. The prevalence of thin beds in practical work makes this analysis quite important. Cases involving permeability anisotropy are treated. INTRODUCTION From wireline formation tester operation, two types of data are obtained: (1) the nature and amount of recovered fluids, and (2) the pressure history recorded during the test. A number of papers have been written dealing with the interpretation of formation production on the basis of the recovered fluids.&apos;.&apos; In general, the methods described have been quite accurate for both high- and low-permeability formations. The present paper will deal with an analysis of the pressures observed. An analysis of the pressure build-up curves obtained in hard-rock country has already been attempted on the basis of the formula proposed by Hor-ner. Although this approach has met with success in many instances, some questions have been raised as to its validity. It is the aim of the present study to place the analysis of pressure build-up in the formation tester on a firmer basis, from which more detailed methods of interpretation can evolve. Because of the great differences between the operation of the wireline formation tester and the conventional drillstem test, modifications are necessary in the interpretation. The major difference relates to the flow geometry. Once the flow geometry has been established other features such as multiphase flow, skin effect, afterflow, etc., well described in the literature, can be introduced. It will be assumed that the mechanical operation of the formation tester is already known to the reader.6 t will suffice here merely to state that the tester provides the means for taking a relatively small sample of the fluid immediately adjacent to the borehole, and for recording the subsequent pressure response. In comparison with conventional drillstem tests, the time required for a satisfactory pressure build-up response is much shorter, because of the relatively small quantity of fluid withdrawn by the wireline tester. This feature is highly desirable in the case of low-permeability formations. For an analysis of the pressure response within the formation, three simple flow geometries are considered— linear, cylindrical and spherical. The spherical and cylindrical flow geometries are most pertinent to the formation tester; therefore, they will receive the major emphasis. Since the configuration of the borehole and the perforation made by the tester complicate the flow geometry, it is necessary to allow for them in the drawdown response. However, because of the volume of formations contributing to the pressure-response, the details of the perforation shape are unimportant in the build-up period. Since relatively small amounts of fluid are withdrawn from the formation, in contrast to a conventional drill-stem test, a study of the "depth of investigation" and the significance of drawdown as well as build-up data will be included. Because the "depth of investigation" will be shown to be rather large, the effect on the build-up curves of the finite thickness of the permeable bed is considered. It is this consideration that leads to the importance of cylindrical flow geometry. Also included is a discussion of permeability anisotropy and its effect on the interpretation of the tester results. The pressure curves recorded by the formation tester will follow two general patterns, depending upon whether the formation is of high or low permeability. Fig. I (a and b) schematically illustrates these two responses. In Fig. 1(a), the high pressure recorded during fill-up of the tool is essentially the pressure differential across the choke in the system. In Fig. l(b), the flow rate is

By R. E. Gilchrist, R. R. Harvey, M. R. Dean

A simple and rapid three-component analysis procedure has been developed for use in fluid-flow studies. The method is based on refractive index measurements, combined with refractive dispersion measured between the Ha and H wave lengths. The procedure requires. calibration curves for the three-fluid system. The only restriction on the fluids employed is that only one be an aromatzc or olefinic compound. Individual sample measurements require only a few minutes. Composition results within the limits of accuracy required for fluid-flow analog studies are obtained. INTRODUCTION The simultaneous analysis of three-component mixtures is a problem often encountered in fluid-flow analog systems. In petroleum production research, the problem is usually encountered in experiments on the displacement of fluid from a porous medium. In the case of immiscible displacement, simple volumetric measurements of the produced fluids are usually sufficient, the analog fluids often being the actual reservoir fluids themselves — oil, water and gas. However, modern secondary recovery of petroleum sometimes involves production under miscible conditions, in which the oil is displaced by liquefied petroleum gas which, in turn, is displaced by dry gas. Since this type of displacement is done under conditions of pressure and temperature that are not easily reproduced in the laboratory, studies are often performed with analog fluids at essentially atmospheric pressure. Since the number of samples taken during one of these experiments is usually quite large, the time required for analysis becomes an important factor. In the past, the difficulty of conventional analysis has been avoided by the addition of various tracers to one or more of the components, or by the use of a tracer fluid as one of the components. With these techniques, a relatively large amount of analytical time is still required. Moreover, in the case of added tracers, the tracer material may migrate in part to one of the other fluids, may be adsorbed on the solid surfaces of the model, or may introduce side effects not characteristic of the prototype fluid. The substitution of a tracer fluid for one of the components can severely limit the choice of analog fluids. Optical refractometry provides a rapid and economical method of analysis of two-fluid systems or of monitoring small changes in composition of an essentially one-component system. Refractometric methods are capable of high degrees of accuracy, require only very small samples and can easily be adapted to automatic methods. They are not generally applicable, however, where three or more fluids are present over wide ranges of concentration. The method of three-component analysis described here has three essential characteristics. 1. It is based on refractometric measurements and thus preserves all of these desirable properties. 2. It does not make any stringent demands on the choice of fluids that can be analyzed. 3. It is capable of accuracy sufficient for fluid-flow analog studies. THEORETICAL The difference in the refractive index of a substance at two different wave lengths of light is known as its dispersivity. When the difference is taken between the specific refractions at the two wave lengths, the quantity is referred to as the specific dispersivity. Specific refraction was defined empirically by Gladstone and Dale1 as r = n-1/d....[1] where n is the index of refraction and d is the density. The definition is given more correctly by the expression r =n2-1/n2+1 1/d..........(2) 1References given at end of paper. MICIFTV OF PETROLEUM ENGINEER JOURNAL.

By H. D. Outmans

In steady vertical flow, the interface of an immiscible liquid-liquid displacement is horizontal for any flow rate below the critical. In nonvertical flow, however, the shape of the interface in the steady state does depend on the flow rate, and the purpose of this paper is to calculate the unsteady interfaces during the transition of one steady state of flow to another. A knowledge of these transient interfaces is of considerable importance in reservoir engineering where the calculation of breakthrough recovery depends on the instant the interface reaches the producing wells and on the shape of the interface at that time. Although the emphasis is put on transient interlaces, which eventually approach stable equilibrium, it is shown that if the displacement exceeds a critical rate no equilibrium is possible. The interface is then unstable and viscous fingers are formed during the displacement. The critical rate and the shape of the transient and equilibrium interfaces are affected by the effective interfacial tension; but since this effective interfacial tension appears in the calculations only in combination with the inverse square of the thickness of the medium, its effect in the reservoir would appear to be negligible compared to its significance in model experiments. INTRODUCTION Stability criteria and the early growth of interfacial disturbances in a plane parallel to the boundaries of a dipping formation in which oil is displaced by an incompressible fluid were described in a previous paper.l This type of instability is significant in thin reservoirs. However, if the reservoir has appreciable thickness, then interfacial stability in vertical planes, normal to the upper and lower boundaries, also becomes important (the displacement is supposed to be parallel to these vertical planes). The difference between the two stability problems is that, in the first case, the intersections of the interface with planes parallel to the boundaries are normal to the direction of the displacement; in the second case, the intersections, this time with vertical planes, are not normal to the displacement. Instead, they are tilted at an angle which depends on the displacement rate. The tilt of steady interfaces was calculated by Dietz2 who also determined the critical rate of displacement for stability in the vertical plane by assuming that this rate would coincide with an interfacial tilt equal to the dip of the formation. The critical rate thus calculated is the same as has been found for thin reservoirs (see Eq. 1.1 of Ref. 1 and of the present paper). Dietz&apos;s calculation of the stable tilt was verified by laboratory experiments and the agreement was found to be fairly good.3 It is doubtful, however, that stable tilts actually exist in the reservoir because a change in production rate is not followed by an instantaneous adjustment of the interface to the new rate but, rather, by a transition period during which the interface changes from one equilibrium tilt to the other. The principal objective of this paper has been to describe these transient interfaces without putting any restrictions on the flow conditions or the shape of the interface, as had been done previously. The second objective was to compute the critical velocity, taking into account capillary effects, and the third was to evaluate, at least qualitatively, the shape of the front at rates above the critical, again without making the simplifying assumptions introduced by previous investigators.2,3 In the following sections two examples are given of the calculation of interfacial motion. The first describes this motion for an initially horizontal interface in a dipping layer, and the second for a vertical interface in a horizontal layer. The mathematical formulation of the problem is nonlinear in the boundary conditions, and this prohibits its solution in closed form. Instead, the solution is obtained in the form of higher-order approximations. 1 Before proceeding to a description of the mathematical model, however, we define two quantities

By T. D. Mueller

Many investigators have used the response of the "dimensionless aquifer" to a unit pressure drop or a unit fluid-withdrawal volume to calculate the performance of an aquifer in supplying water influx to an oil reservoir. In the past, these response functions have been calculated with the aid of the Laplace transform. With the advent of ultra highspeed digital computers, it becomes practical to solve for the response functions with finite-difference techniques. The computer method also permits extension of the dimensionless-aquifer concept to include the nonhomogeneous aquifer wherein the permeability and other properties vary as a function of the space co-ordinates. This paper gives results of calculating the response functions for a series of nonhomogeneous aquifers. Response functions are presented for both linear and radial aquifers whose thickness, permeability-viscosity ratio and porosity-compressibility vary. These functions are new and should prove useful to the petroleum engineer in analyzing the behavior of nonbomogeneous aquifers. Results are presented in the form of charts that can be easily used by the field engineer. INTRODUCTION Aquifers which surround many oil and gas reservoirs have the ability to supply water influx to such reservoirs as oil and gas are withdrawn. This water influx, called natural water drive, provides one of the most effective driving mechanisms for the production of oil and gas. In producing a reservoir, therefore, it behooves one to make the maximum use of natural water drive. To achieve the maximum use, the reservoir engineer must be able to predict the performance of an aquifer under a variety of production schemes that may be proposed for the reservoir. Unfortunately, the physical properties which dictate aquifer behavior often are known only within limits. Seldom do wells penetrate the porous strata of the aquifer. Even when they do, quantitative information regarding porosity, permeability and water compressibility is seldom available. It is known, however, that the water efflux from most aquifer systems is governed by a single, relatively simple, linear, partial differential equation. Also, the general physical location of the aquifer boundaries often are known. A technique originall proposed by Hurst1 and van Everdingen and Hurst2 has been found useful in analyzing reservoirs in this situation. The idea was later expanded by van Everdingen, Timmerman and McMahan3 to include the mathematical technique of least-squares fitting. This latter method will be referred to as the VTM method. The basic assumptions of the VTM method include the following. 1. The location of the physical boundaries of the aquifer are known. 2. The flow conditions at these physical bound-daries are also known. 3. The aquifer is homogeneous; e.g., thickness, permeability, water compressibility and porosity are constant throughout. In the VTM method, a material balance is made on the fluids entering and leaving the reservoir. In the balance equations, the water-influx term is represented as the product of the water influx from an arbitrarily-selected, dimensionless aquifer system times an unknown conversion number. This balance can be formed as many times as there are data points in the history of the reservoir. Each time, the conversion number can be evaluated. If the reservoir engineer has picked the correct dimensionless aquifer to represent the real aquifer, the conversion number remains constant for all balances that have been made over the history period. If such a situation occurs, the reservoir engineer can then predict the performance of the reservoir for any type of production scheme by using the function associated with that particular dimensionless system and the derived conversion number. These functions will be referred to as the response function of the aquifer. If the aquifer is nonhomogeneous (e.g., if the porosity, permeability, thickness, porosity, or

By D. R. Parrish, K. W. Beaver, H. W. Wood, R. W. Rausch

A pilot test of forward combustion in the Shannon pool, Salt Creek field, Wyo., is described. The Shannon sand, 950-ft deep, contains a heavy (25" API), viscous (76 cp) oil. Natural reservoir energy is limited. Primary production, intermittent since 1889, recovered only about 2 per cent of the oil in place. The field is operated by Pan American Petroleum Corp. for the Midwest Oil Corp., the owner. The original pilot was a 1.32-acre five-spot. The expanded pilot has eight producing wells surrounding a roughly triangular area of about five acres with the injection well near the center. A control or comparison well was also recompleted in another part of the field. Operation of the pilot has been little different front an ordinary gas drive. Little special equipment was found to be absolutely necessary. Except for some use of a temperature-resistant cement, all wells were conventionally completed. In spite of poor oxygen consumption, the over-all performance of the pilot has been good. Total oil recovery to June 1, 1961, was 73,971 bbl. The wells of the original pilot alone had produced about 24,000 bbl, equivalent to 50 per cent of the oil in place, when fire breakthrough at the first well occurred. These wells have now produced oil equivalent to more than 74 per cent of the oil in place in the original pilot area and production is continuing. It appears that ultimate recovery will approach theoretical maximums before the wells must be abandoned. Performance of the pilot has been encouraging, and expansion to a fieldwide combustion operation is being investigated. INTRODUCTION The results of both laboratory investigations and field tests of underground combustion have been reported previously.&apos; However, most of the field tests were primarily experimental. More information and experience are needed before forward-combustion operations can be engineered with confidence. The purpose of this paper is to present the results of a successful pilot test of forward combustion. These re- sults should increase confidence in forward combustion as a practical method for commercial oil recovery. HISTORY OF THE SHANNON POOL The Shannon pool is located on the north end of the Salt Creek field in Natrona County, Wyo. It is approximately 50-miles north of the city of Casper. The Shannon pool&apos;s3 discovery well was completed in 1889, making this one of the oldest oil fields in the Rocky Mountain region. Three more wells were drilled in 1890. First production was hauled in wooden barrels by horse and wagon to the railroad in Casper. In 1894 a small refinery, the first in Wyoming, was built in Casper to process the Shannon crude. In the following years the field changed hands several times. It appears that each new owner did some development drilling as several wells were completed in each of the years 1895, 1902, 1905 and 1912, with negligible development in the intervening years. Forty-eight wells were drilled in this period, but many were later abandoned. Discovery of the more prolific Salt Creek field proper ultimately forced suspension of operations at the Shannon pool. After 1915 there was only sporadic production, mostly to supply cheap boiler fuel to drilling rigs in the Salt Creek field. But even this was discontinued in 1931. Since then the pool had been dormant until the recent operations, which are the subject of this paper. The Shannon pool is now owned by the Midwest Oil Corp. Field operations are conducted for Midwest by Pan American Petroleum Corp. THE RESERVOIR Fig. 1 is a map showing subsurface contours of the Shannon pool. The reservoir is on a nose of the Salt Creek anticline dipping to the north at about 500 ft/mile. The trap is provided by a shallow fault on the updip side of the productive area. The downdip limits are bounded by water, but this water has not provided an effective source of reservoir energy. The Shannon sand outcrops at many places in the immediate vicinity, providing good surface indications of the Salt Creek anticline. At the Shannon pool the sand has been lowered by faulting and is overlain by about 900 ft of shale and other sands. The Shannon sand consists of two members. The upper, a water sand, is about 40-ft thick and is separated from the lower member by several feet of sandy shale. The oil pay occurs in the lower mem-

By H. J. Ramey

As fluids move through a wellbore, there is transfer of heat between fluids and the earth due to the diflerence between fluid and geothermal temperatures. This type of heat transmission is involved in drilling and in all producing operations. In certain cases, quantitative knowledge of wellbore heat transmission is very important. This paper presents an approximate solution to the wellbore heat-transmission problem involved in injection of hot or cold fluids. The solution permits estimation of the temperature of fluids, tubing and casing as a function of depth and time. The result is expressed in simple algebraic form suitable for slide-rule calculation. The solution assumes that heat transfer in the wellbore is steady-state, while heat transfer to the earth will be unsteady radial conduction. Allowance is made for heat resistances in the wellbore. The method used may be applied to derivation of other heat problems such as flow through multiple strings in a wellbore. Comparisons of computed and field results are presented to establish the usefulness of the solution. INTRODUCTION During the past few years, considerable interest has been generated in hot-fluid-injection oil-recovery methods. These methods depend upon application of heat to a reservoir by means of a heat-transfer medium heated at the surface. Clearly, heat losses between the surface and the injection interval could be extremely important to this process. Not quite so obvious is the fact that every injection and production operation is accompanied by transmission of heal between wellbore fluids and the earth. Previously, the interpretation of temperature logs&apos;,&apos; has been the main purpose of wellbore heat studies. The only papers dealing specifically with long-time injection operations are those of Moss and White3 and Lesem, et al.&apos; The purpose of the present study is to investigate wellbore heat transmission to provide engineering methods useful in both production and injection operations, and basic techniques useful in all wellbore heat-transmission problems. The approach is similar to that of Moss and White:&apos; DEVELOPMENT The transient heat-transmission problem under consideration is as follows. Let us consider the injection of a fluid down the tubing in a well which is cased to the top of the injection interval. Assuming fluid is injected at known rates and surface temperatures, determine the temperature of the injected fluid as a function of depth anti time. Consideration of the heat transferred from the injected fluid to the formation leads to the following equations. For liquid, Eqs. 1, 1A and 2 are developed in the Appendix. These equations were developed under the assumption that physical and thermal properties of the earth and wellbore fluids do not vary with temperature, that heat will transfer radially in the earth and that heat transmission in the wellbore is rapid compared to heat flow in the formation and. thus, can be represented by steady-state solutions. Special cases of this development have been presented by Nowakl and Moss and White.3 Both references are recommended for excellent background material. Nowak&apos; presents very useful information concerning the effect of a shut-in period on subsequent temperatures. Since one purpose of this paper is to present methods which may be used to derive approximate solutions for heat-transmission problems associated to those specifically considered here, a brief discussion of associated heat problems is also presented in the Appendix. Analysis of the derivation presented in the Appendix will indicate that many terms can be re-defined to modify the solution for application to other problems. Before Eqs. 1, 1A and 2 can be used, it is necessary to consider the significance of the over-all heat-transfer coefficient U and the time function f(t). Thorough discussions of the concept of the over-all heat-transfer coefficient may be found in many references on heat transmission. See McAdams5 or Jakob," for example. Briefly, the over-all coefficient U considers the net- resistance to heat flow offered by fluid inside the tubing, the tubing wall, fluids or solids in the annulus, and the casing wall. The effect of radiant heat transfer from the tubing to the casing and resistance to heat flow caused by scale or wax on the tubing or casing may also be included in the over-all coefficient. According to McAdams, on page 136 of Ref. 5>

By O. R. Wagner, R. O. Leach, H. W. Wood, C. F. Harpke

PAN AMERICAN PETROLEUM CORP.TULSA, OKLA. CASPER, WYO A field test has been made in which additional oil recovery was obtained from a previously-waterflooded "oil-wet" sandstone reservoir. This recovery improvement was accomplished by adjusting the reservoir wettability through chemical treatment of the flood water. The test was made in the Muddy sand of the West Harrisburg Unit, Neb. The chemical used, sodium hydroxide, was injected as a slug of dilute caustic solution through a water-injection well. The natural wettability of the reservoir and the chemical requirements for reversing wettability were determined in the laboratory from contact-angle studies using field water and oil samples. Laboratory flood tests with a synthetic system had shown that reversing the wettability of an oil-wet consolidated core would lead to improved oil recovery. The field performance indicates that the mechanism by which increased oil retovery is obtained in the field is the same as that observed in the lahoratory. Laboratory studies indicate that higher ultimate recoveries and decreascd water-injection requirements result when the adjusting agent is added early in the life of the flood. However, a previously waterflooded area was intentiorzally chosen for the field test so that unambiguous conclusions could be made about the effects of chemical injection on wettability and the extent of oil-recovery improvernent afforded by wettability reversal. Forces existing at the fluid-solid and fluid-fluid interfaces in a porous medium have an important effect on oil recovery during a water flood. Modifying these interfacial forces in the reservoir to improve oil recovery has been the object of much research. Several papers&apos; &apos; have discussed the use of surfactants to lower the water-oil interfacial tension. This paper describes work concerned with improving oil recovery by modifying the forces at the fluid-solid interfaces— that is, by changing the preferential wettability. An earlier paper: presented evidence to show that some reservoir systems could be changed from preferentially oil-wet to preferentially water-wet by the action of simple chemicals added to the water to increase waterflood oil recovery. The process was attractive for further study because it used only inexpensive chemicals and because it called for a gross "wettability reversal" rather than for a precise adjustment. The present paper is an extension of the earlier work. Laboratory flow tests were made using treated water and oil in a consolidated core to determine the amount of additional oil recovery and the producing performance which might be expected from reversing the wettability of an oil-wet reservoir during a flood. A contact-anglc study using water and oil from the Muddy "J" sand of the Harrisburg field, Banner County, Neb., indicated the reservoir to be preferentially oil-wet and susceptible to wettability reversal through chemical injection. Based on these studies, a field trial of wettability-reversal water flooding was initiated. in the Harrisburg field. The primary purpose of the trial was to determine if wettability-reversal flooding would improve oil recovery in an actual field situation. The field trial was also an important test of the laboratory tecnique used to determine the effect on oil recovery of certain wettability manipulation, as well as a test of the contract-angle method of determining reservoir wetra-bility. The paper is presented in two parts. The first part covers the laboratory experiments leading up to the field trial, and the second part covers the field trial. PART I—LABORATORY STUDY The laboratory experiments consisted of displacement tests using an idealized fluid system and of a contact-angle study with crude and water from the Muddy "J" sand of the West Harrisburg Unit, Banner County, Neb. The displacement tests were designed to determine the effects of a particular wettability manipulation, oil-wet to water-wet, on oil displacement in a water flood. The contact-angle measurements were made to determine the natural wettability of the reservoir and to determine if beneficial wetting changes could be brought about by a chemical addition to the flood water. PROCEDURE DISPLACEMENT TESTS The displacement tests were perfornled in a consolidated sandstone core using a refined oil and water. The wetting properties of this system could

By C. W. Carpenter, P. T. Bail, J. E. Bobek

This paper presents the results of model studies made to demonstrate the validity of applying the "basic" waterflooding scaling relationships of Rapoport1 to water-wet reservoirs composed of communicating strata of different permeability. Results demonstrate that these relationships can be validly applied, and also show the combined effect of variables on recovery. This demonstration was needed to confirm that the sizable effects of communicating strata on waterflood behavior, observed in model studies, were actually representative of field behavior. INTRODUCTION Rapoport1 Collins2 and others have presented scaling laws for modeling water-oil displacements. The validity of these laws has been experimentally demonstrated for homogeneous systems only. Recently, however, several investigators have presented the results of multilayer model studies made to determine the effects of viscous, capillary and gravitational forces on oil recovery in systems composed of communicating layers of different permeability. Dyes and Braun,3 using miscible displacements to simulate water flooding, showed that vertical communication adversely affects recovery at unfavorable viscosity ratios in the absence of capillarity. Ogandzhanyants and Egorova, 4 Gaucher and Lindley,5 and Richardson and Perkins6 modeled water-oil displacements in layered systems using immiscible fluids with capillary forces present. These investigators found more efficient behavior than in the miscible system of Dyes and Braun. However, since the scaling laws had not been verified for heterogeneous systems prior to or during these tests, it was necessary to assume that the results were representative of field behavior. This report presents the results of model studies made at Jersey Production Research Co. to find out if the waterflooding scaling relationships are valid when applied to water-oil displacements in heterogeneous reservoirs with vertical communication. Data obtained in the study show the combined effects of viscous, capillary and gravitational forces on oil-recovery behavior. The studies were performed in strongly water-wet systems where strong imbibition forces were present. For systems of different wettability, the effects of capillary imbibition would be somewhat different than those presented. APPROACH The scaling laws studied were those developed by Rapoport1 and designated by him as "basic" scaling laws. These relationships require that the model be scaled to the geometry of the prototype, and that the model and prototype be operated at the same viscosity ratio and have identical relative permeability curves, contact angles and boundary conditions. It is also necessary that the groups** be identical in the two systems, although any of the parameters making up the groups may be different. Finally, the capillary pressure of one system must be directly proportional to that of the second. Studies were conducted in two four-layer models, representing quadrants of five-spots, which were geometrically similar and which had the same glass-bead packing. The fluids used assured the same wettability in both models. Thus, the contact angles and relative permeabilities of comparative sand layers in the two models were identical and the capillary pressures were directly proportional. However, one model had a characteristic length three times that of the other. To make the scaling groups equal and to check the scaling laws, it was necessary to adjust the rate, water-viscosity level, interfacial tension and water-oil density difference to compensate for the difference in model sizes.

By M. H. Gaskell, D. C. Lindley

Small, steeply inclined reservoirs without natural water drives often are found associated with salt domes or other highly faulted structures. Frequently, only one well may be economically justified in these reservoirs, from which oil is produced by an expanding gas-cap drive or by a dissolved-gas drive. Although a large fraction of the oil upstructure from the well may be recovered during primary depletion, only a small fraction of the downdip oil can be produced because of gravity effects. This downdip oil has been defined as cellar oil991 in contradistinction to "attic oil" which is the oil above the structurally highest well in reservoirs where a natural water drive is present. In this study, two methods of producing cellar oil by water flooding a one-well reservoir were investigated by means of a laboratory model. This model was scaled to represent a prototype reservoir having an area of 30 acres, a permeability of 500 md, and a thickness and dip of 30 ft and 30°, respectively. One method considered, the simultaneous injection of water and production of oil from a single well, would require that the lower portion of the productive interval of a well be completed for water injection and the upper portion for oil production. The other method considered, intermittent injection, would involve alternate injection of water into and production of oil from a single well. Injection of water into the formation would repressure the reservoir, causing gas to go back into solution. Since some of the injected water would flow updip from the well, water injection usually would be followed by a waiting period during which the water could drain from above the wellbore. Production of oil by gravity drainage and gas-cap drive could be initiated when water drainage became sufficiently complete. The amount of water produced along with the oil would naturally depend on how soon after water injection ceased that production was started. This process could be repeated in cycles until oil could no longer be produced economically. Of these two methods, the intermittent injection technique would be simpler to put into operation since it does not require a dual completion. Also, it is probably more generally applicable than the simultaneous injection-production method. It is to be expected that in some reservoirs the rock and fluid properties would be such that practical oil-production rates could not be achieved by the simultaneous technique. However, in reservoirs where the simultaneous method will result in practical rates of oil production, it does have the advantage of permitting immediate, continuous oil production; whereas in the intermittent method, oil production is interrupted by periods of water injection and drainage. It was the objective of this study to determine whether these two methods of cellar oil production are practical for the representative prototype investigated and to evaluate their relative merits. DESIGN OF LABORATORY MODEL SCALING REQUIREMENTS The prototype reservoir and fluids were selected to have the dimensions and physical properties shown in Table 1. The laboratory model was then designed in accordance with well known scaling criteria, as discussed by Croes, et al,2 except for one modification. According to Croes, et al, one requirement of scaling a laboratory model to a reservoir prototype is that both the relative permeabilities and dimensionless capillary pressure be the same functions of water saturation for the model and the prototype. Since unconsolidated sand packs are used in the laboratory model to represent the consolidated porous media of the prototype, it is virtually impossible to satisfy this requirement. Much of the difficulty stems from the fact that unconsolidated sands usually exhibit lower residual oil and connate-water saturations than do reservoir rocks. To satisfy the scaling requirement as nearly as possible, Perkins and Collins&apos; have suggested that the relative permeabilities be re-defined as

By F. E. Crane, J. Downie

The general application of Darcy&apos;s law to natural rocks has already been challenged in the literature.1 The evidence shows that the permeability as calculated from the Darcy equation can be a function of the pressure drop and the salt concentration of the water phase. Most explanations for aberrant behavior involve clays and their properties and have been qualitatively satisfactory.Recently, however, Odeh2 revived the theoretical views of Yuster.3 Since Yuster&apos;s concept implies a fundamental error in using Darcy&apos;s relationship for two-phase flow, and not merely that conditions may limit its use, Odeh&apos;s experimental support arouses considerable interest. However, Odeh&apos;s work as presented is thought to to be inadequate. He has omitted important information about his materials and procedures; therefore, acceptance of his conclusions should be withheld. Since a clay effect is possible, Odeh should have been allowed space to present a more detailed account of the rocks which he used and the precautions which he took to avoid confounding the effects of clay and fluid circulation. In addition, it follows from Odeh&apos;s discussion that the oil relative permeability should increase as the viscosity ratio increases. Odeh presents data corroborating this deduction. It also follows that the oil relative permeability should decrease as the viscosity ratio of oil to water decreases. However, our results show that a high relative permeability, when once attained by using a viscous oil, may be maintained when that oil is replaced by an oil of much lower viscosity. EXPERIMENTAL PROCEDURE The results are shown in Fig. 1 and were obtained as follows. Sandstone cores. 3-in. long X l-in, in diameter and which contained clays, were saturated with