Search Documents

By D. R. Wieland, H. T. Kennedy

The object of the work reported here was to determine the content of elemental sulfur in gaseous methane, carbon dioxide, hydrogen sulfide, and in mixtures of these gases, at pressures and temperatwes encountered in natural gas reservoirs. Sulfur content at equilibrium is reported for pure methane, carbon dioxide, hydrogen sulfide, and on three binary mixtures of each of the three pairs of gases at pressures of 1,000, 2,000, 3,000, 4,000, 5,000 and 6,000 psia and at temperatures of 150°, 200° and 250°F. In addition, the sulfur content of three ternary mixtures at the same temperatures and pressures are reported. The results indicate that the sulfur content is higher in the gases at higher temperatures and pressures. The content is highest in hydrogen sulfide, intermediate in carbon dioxide and lowest in methane. INTRODUCTION At ordinary pressures and temperatures, the concentration of a nonvolatile material, such as sulfur in a gas at equilibrium, is a function of the vapor pressure of the material and is independent of the nature of the gas. As the pressure and temperature increase, the gas assumes some of the properties of liquids, including the power to dissolve other liquids and solids, to an extent dependent on the nature of both the gas and the material dissolved. At equilibrium, the content of sulfur in the several gases considered here may thus be considered as solubilities which are fixed for a given composition of gas, temperature and pressure. The study of the solubility of elemental sulfur in gases is of interest because sulfur is sometimes present in reservoirs producing natural gas and must be present in the vapor phase. Upon reduction of pressure and temperature, the sulfur precipitates from solution in the reservoir and in the tubing and fittings. Due to the fact that the greatest pressure drop in the reservoir is around the wellbore, the volume of sulfur precipitated will be greatest in this locality and can cause a substantial reduction in the permeability of the formation in this area. As the natural gas flows up the tubing string, the pressure and temperature are further decreased, causing further sulfur precipitation. If the volume of free sulfur is large and remedial measures are not taken, complete plugging of the tubing can occur. Knowledge of the content of sulfur as a function of composition, temperature and pressure will enable the operator to determine what changes in pressure and/or temperature are necessary to keep the sulfur in solution. If the sulfur solubility is large and can be controlled, the operator may desire to produce the well at a high wellhead pressure and temperature and reclaim the sulfur at the surface. The high solubility of sulfur in hydrogen sulfide (5.6 per cent by weight at 5,000 psi and 200°F) suggests the feasibility of recovering sulfur by injecting this gas into sulfur-bearing reservoirs and expanding the produced gas to precipitate out the sulfur. The present project was designed to measure the solubility of sulfur in carbon dioxide, methane, hydrogen sulfide and mixtures of these three gases at various temperatures and pressures. These gases were the major constituents of the gas well where plugging by sulfur was encountered. A study of the published literature shows no data on the solubility of sulfur in any gas, although Hannay and Hogarth,' in 1880 reported that it was soluble in carbon disdlide above the critical temperature of carbon disulfide. EQUIPMENT AND PROCEDURE Fig. 1 shows the layout of equipment employed to bring gases to equilibrium with sulfur and to isolate the quantity of sulfur contained in a known volume of gas. Pure gases were measured in the charging bomb at known temperatures and pressures and displaced into the reservoir bomb to make up the desired mixtures. When equilibrium was achieved

By W. Hurst

What is entailed here is the extension of the sinzplified material balance formulas to encompass interference between oil fields. As previously reported, the ex-plicitness as so revealed for the cunzulative pressure drop as a function of all factors contributing to its change in the material balance equation, is now transcribed through the inter-conzmunicating aquifer to effect an increased pressure drop on an adjoining field by interference. Such is performed by mathematical analyses and the application of the Laplace transformations. What is accomplished is that the reiteration problems previously associated with interference studies are nullified, since volumetric changes for the fluids in situ are automatically adjusted by the explicitness so expressed, and what pertains to the superposition principle applies only to the impeded water drive upon a su.bject field, which is likewise incorporated in the over-all pressure drop that results. The mathematics treats with the rigorous solution of the problem, as well as with methods easily amenable to numerical interpretation by the practicing engineer. In all cases, however, what is deduced for areal extents, reveal interference between wells when time becomes large that further substantiate the analyses. Probably of equal significance, all variables and for/nation characteristics are accounted for. This applies to the differences in PVT analyses that occur from field to field and the physical parameters associated with the lithology of the reservoirs. The latter is deduced from a rigorous interpretation of the unsteady-state flow problem for sands of different permeabilities in series. Thus, what purports to be a trial-and-error calculation to include variations in sand conditions within the intercommunicating aquifer to define water drive has little importance compared to the designation of such paranzeters at the reservoir to constitute the essential criteria. INTRODUCTION In a recent publication, the reader has been introduced to the simplified material balance formulas,' It has been shown that the complexities formerly associated for identifying and determining reservoir pressure in the material balance relationship can be resolved by treating with the Laplace transformations, the inversion of which reveals pressure as an explicit function of all factors contributing to its change. For an under saturated oil reservoir, this constitutes an integrated effect from the inception of production; and for a saturated oil reservoir, such represents a survey traverse. However, what is most important, for the first time we are able to express the pressure change in a reservoir directly to the oil, gas and water produced, and to the physical parameters of the formation, in the same manner that reservoir pressure has served in the past to define well interference — which opens up this avenue of endeavor for investigating interference between oil fields. The occurrence of interference is by virtue of the juxtaposition of two or more oil reservoirs producing in the same aquifer. Thus, in any one field, the extended production so incurred creates a pressure lowering in the field, as well as in the aquifer, which can now be transcribed by these simplified material balance formulas to the adjoining reservoir to observe the increased pressure change so induced. These effects are retroactive from field to field, to reveal a pressure drop far in excess of normal depletion to evidence the interference pattern that can result. Such is the purpose for the present undertaking — to develop the mathematical physics treating with this phenomenon of interference between oil fields, employing the simplified material balance formulas. Since Ref. 1 often will be referred to in the text, it is identified as "The Simplification", in the same connotation that another paper has been referred to as the "Laplace Trans-formation. However, the amount of material to be covered in these passages is substantial; therefore, the total emphasis will be placed upon developing these fundamental ideas treating with interference rather than attempting to offer work curves, which always can be prepared. For the practical application of these simplified material balance formulas that form the framework of the present development, the reader will find a comprehensive treatment in the Appendix of Ref. 1, illustrating a factual field example. In essence, what is proposed by this presentation is to make available to the practicing engineer a simplified means to treat with interference between oil fields, independent of the complex machine computations that have attended such undertakings in the past. As the problem unfolds, a common pattern will be discernible, whereby the many reiteration processes that previously required computing devices are eliminated.

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved

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

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

By F. Selig, A. S. Odeh

A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well's shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner's' theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner's widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs' analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner's' method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner's method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner's original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner's approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink

By A. J. Conelius, D. P. Sobocinski

This paper presents a correlation for predicting the behavior of a water cone as it builds from the static water-oil contact to breakthrough conditions. The correlation is partly empirical and involves dimensionless groups of reservoir and fluid properties and of production and well characteristics. The groups were deduced from the scaling criteria for the immiscible displacement of oil by water. The correlation is based on a limited amount of experimenfal data from a laboratory sand-packed model and on results from a computer program for a two-dimensional, incompressible system. Because the correlating groups are dimensionless, they can be used to estimate the performance of water coning cases not specifically considered in the correlation. However, despite its dimensionless nature, the correlation is not completely general and will not provide meaningful estimates of cone behavior in many situations. INTRODUCTION AND BACKGROUND The production of water from oil wells is a common occurrence which increases the cost of producing operations and may reduce the efficiency of the depletion mechanism and the recovery of reserves. We will deal with one cause of this water production, namely, coning. The coning of water into a producing well is caused by pressure gradients established around the wellbore by the production of fluids from the well. These pressure gradients can raise the water-oil contact near the well where gradients are most severe. Gravity forces that arise from fluid density differences counterbalance the flowing pressure gradients and tend to keep the water out of the oil zone. Therefore, at any given time, there is a balance between gravitational and viscous forces at points on and away from the completion interval. When the dynamic forces at the wellbore exceed gravitational forces, a cone of water will ultimately break into the well to produce water along with the oil. We can expand on this basic visualization of coning by introducing the concepts of stable cone, unstable cone and critical production rate. For instance, if a well is produced at a constant rate and the pressure gradients in the drainage system have become constant, a steady-state condition is reached. If, at this condition, the dynamic forces at the well are less than the gravity forces, then the water or gas cone that has formed will not extend to the well. Moreover, the cone will neither advance nor recede, thus establishing what is known as a stable cone. Conversely, if the pressure in the system is in an unsteady-state condition, then an unstable cone will continue to advance until steady-state conditions prevail. If the flowing pressure drop at the well is suffcient to overcome the gravity forces, the unstable cone will grow and ultimately break into the well. It is important to note that in a realistic sense, stable cones may only be "pseudostable" because the drainage system and pressure distribution generally change. For example, with reservoir depletion, the water-oil contact may advance toward the completion interval, thereby increasing chances for coning. As another example, reduced productivity due to well damage requires a corresponding increase in the flowing pressure drop LO maintain a given production rate. This increase in pressure drop may force an otherwise stable cone into a well. The critical production rate, well known in the literature, is the rate above which the flowing pressure gradient at the well causes water (or gas) to cone into the well. It is, therefore, the maximum rate of oil production without concurrent production of the displacing phase by coning. At the critical rate, a built-up cone is stable but is at a position of incipient breakthrough. Numerous papers have been published on critical rates. Some of the better known of these include the work of (1) Muskat and Wyckoff,' who first dealt with the coning problem; (2) Chaney, et al,' who developed expressions similar to those of Muskat but who presented results in a conven ient-to-use graphical form (the "Sun" method); and (3) Meyer and Garder,b hose analysis is based on radial-flow formulas. One assumption in critical production rate analyses is that the cone has built-up to just before its breakthrough into the well. But, these analyses reveal nothing directly about the time it takes for the cone to build up to this incipient breakthrough position. Thus, water-free oil can be produced from a well for prolonged periods at rates above the critical rate before the well reaches the condition to which the critical rate applies. The published literature contains little on the rate of growth of a cone. Experimentally, Meyer and Searcy studied the rate of rise of a cone in a Hele-Shaw model.' Additional related work on water breakthrough and produced water-oil ratios in water driven reservoirs was reported by Muskat," Hutchinson and Kemp," Henley, et al,' and Stevens, et a1.B Theoretically, the basic coning equations for a water-oil system can be developed by applying the conservation of mass to each of the phases, relating flow velocities with pressure by Darcy's law, and relating pressures across water-oil interfaces by capillary pressure. With the usual boundaries at the well and reservoir limits, the solution of the resulting equations for the time behavior of a water-oil interface constitutes a free-surface, boundary-value problem. The shape of the boundary depends on the internal potential distribution, while this

Jan 1, 1966

By R. W. Heins, N. Street

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

Jan 1, 1966

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

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

Jan 1, 1965

By R. W. Parsons

This discussion concerns fitting an ellipse to angle, where the permeability k is measured for a series of plugs cut from varying angles in the horizontal plane. The conclusions and procedures given by Greenkorn, et al. are in no way affected by these comments, since the degree of anisotropy (ratio of maximum to minimum permeabilities) they found is generally quite close to one (see their Table A-1). However, since large anisotropies may arise, I feel that the comments should be made. Both Marcusl and Scheidegger2 have shown that there are two limiting expressions for the permeability of an unbounded anisotropic media. If the permeability is measured in the direction of the pressure gradient, the result is equivalent to the author's Eq. 11. This is the equation used to find kl1 and k22 from the experimental data. If, however, the permeability is measured in the direction of flow, then the resulting equation for the measured permeability is From these equations it follows that plots on angle will be ellipses. The principal axes will correspond with the major and minor permeability axes. The question arises as to what permeability is actually being measured on an anisotropic sample in the laboratory. An overall pressure gradient and a gross flow rate are determined and a permeability calculated by Darcy's law. The details of the macroscopic flow within the sample will be determined by the boundary conditions of the sample and the orientation and magnitude of the major and minor permeability axes. The situation is shown in Fig. 1 for the two-dimensional case. Using the method of conformal mapping, Marcusl has solved the problem of determining the overall apparent (or measured) permeability for this experiment for any condition of W-L ratios and orientation and magnitude of the principle permeability axes. For W-L ratios approaching zero, the streamlines are all forced by the boundary conditions to be in the direction of the measured permeability which is equal to kg. For W-L ratios approaching infinity, the effects of the no-flow bounds are negligible, and the pressure gradient and measured permeatility direction are parall;l. The streamlines are

Jan 1, 1965

By W. B. Gogarty

The reported decreases in water mobility do not seem unusual in view of non-Newtonian fluid properties. Shear stress vs shear rate diagrams have been reported for other solutions of water-soluble polymers. Some of these polymers are similar to the type mentioned by the author. Generally, the shear stress-shear rate is a non-linear function for these solutions. Data for plotting apparent viscosity vs shear rate can be obtained from this function. Apparent viscosity is defined as the ratio of shear stress to shear rate at a given shear rate. When plotted, the apparent viscosity decreases with increasing shear rate. This behavior is typical of a pseudoplastic fluid. For some water-soluble polymer solutions, the apparent viscosity decreases more than 50 times while the shear rate increases 1,000 times. Thus, viscosity of a pseudoplastic fluid only has meaning at a specified shear rate. Results of Fig. 1 could be explained in these terms. Viscosities measured in the Ostwald viscometer represent values at a given shear rate. Some average shear rate is affecting the polymer solutions while flowing through the core. This average value fixes the apparent viscosity as long as the flow rate remains constant. Viscosities measured by the two methods will be equal if shear rates are the same. The results indicate that shear rate in the core is lower (higher apparent viscosity) than in the viscometer. In the paper by Johnson, Bossler and Naumann, the relative permeability is independent of viscosity ratio. Thus, the relative permeability with respect to water flow at residual oil should be independent of the flowing phase viscosity. Polymer solutions will appear as Newtonian fluids The discussion emphasizes the nature of the "resistance factor effect" as discussed in the paper. Repeated anomalies arising in hundreds of experiments led us to the conclusion that non-Newtonian flow is not the only factor. Several of the key anomalies are as follows: 1. Measured viscosities over a range of shear rates from <1 sec-&apos; to 1,000 sec-&apos; do not account for but a minor fraction of the R observed in cores when compared in similar shear-rate ranges. 2. The slope of R vs flow rates in cores is always different from that expected from viscometer shear-rate measurements as shown in Fig. 2. in a core, the level of viscosity being fixed at a given flow rate. With these conditions, the definition of resistance factor R by Eq. 2 is simplified to Since , is constant with rate, R becomes a measure of the apparent viscosity in a core at a given flow rate. Variation in flow rate could easily account for the changes of R shown in Fig. 5. Also, this points to the fallacy of assuming R to be a unique parameter. The constant resistance factors at different flooding velocities appear to be in disagreement with the above discussion. The author furnishes Fig. 2 to support his arguments. As shown, the resistance factors remained substantially constant in the two cores over a considerable range of flooding velocities. However, in the 73-md core, the factor increases at lower rates. This behavior agrees with known characteristics of some pseudoplastic material. These materials act both as Newtonian and as non-Newtonian fluids in different regions of shear rate. Some exhibit first Newtonian, then non-Newtonian, finally, Newtonian character. Others are first non-Newtonian and then Newtonian. This latter type would explain the results with the 73-md core. The Fann-instrument results are not significant since shear rates in the core may be much different than with the viscometer. The higher resistance factor at high rates in the 150-md core is more difficult to explain. The greater resistance at increased flow rates could be attributed to what might be termed temporary bridging. As envisioned, changes in polymer configuration occur at the higher energy associated with the increased flow rate. These changes could cause less effective passage of polymer through the core. Correspondingly, increases in pressure drop will occur. These will be interpreted as higher resistance factors. 3. Most polymer solutions are non-Newtonian and many are more shear-rate sensitive than the polymers in question, yet only a very few polymers demonstrate useful R values. Gogarty&apos;s assumption that viscosities in cores and vis-cometers will be the same if measured at the same shear rate is only valid if non-Newtonian rheology is the only parameter. The experimental evidence does not validate this assumption. The anomalies observed in the equilibrium displacement experiment shown in Fig. 5 are not explained on the basis of varying flow rates since the rates were held constant. M

Jan 1, 1965

By R. L. Dalton, L. A. Rapoport, C. W. Carpenter

Procedures have been developed to study the performance characteristics of unconfined pilot water floods using scaled laboratory models. The effects of operating conditions on pilot performance for variozrs well patterns have been investigated. The limitations of the laboratory model studies are outlined. Relationships between the oil recoveries of pilots and large-scale floods have been established. The applications of model studies to the prediction of large-scale waterflood performance and pilot flood design are discussed. INTRODUCTION Pilot water floods are undertaken to evaluate operational procedures and to gain advance information about the performance of an extensive water flood. In general, however, the performance of a pilot flood is not directly indicative of what can be expected from a large-scale pattern water flood. This difference stems from the fact that the over-all flow configurations are not the same for the two systems. A water flood normally consists of an array of identical well patterns. In such an array, the perimeters of the well patterns are axes of symmetry and act as impermeable boundaries. Thus, an extensive, pattern water flood can be visualized as a repetition of "confined" floods. In contrast, a pilot water flood involving only one or a few well patterns is "unconfined". Here the well patterns are not balanced by other flood units; hence, the perimeter of the pilot area does not act as an effective boundary. Accordingly, only a portion of the fluids injected (and also of those originally contained within the pilot area) is captured by the producers of the pilot. The rest escapes into the surrounding reservoir. This usually causes the amounts of oil and water produced in the pilot to be different from those recovered by a producer in an extensive water flood. Moreover, the recovery observed in the pilot can be expected to be greatly influenced by operating conditions, i.e., by the bottom-hole pressures at the injection and producing wells. The above considerations show that full utilization of pilot test results and reliable prediction of large- scale water floods require a knowledge of the relationship between the recovery performance of confined and unconfined well patterns. To date, only limited basic information is available on the behavior of pilot water floods. Craig, et al,&apos; have investigated the effects of initial gas saturation, oil-to-water mobility ratio, and well damage on the oil recovered from a single, "normal" five-spot pilot pattern (Fig. 1). Rosenbaum and Matthews2 studied the effect of initial gas saturation and mobility ratio on the ratio of production rate to injection rate for various groupings of five-spot patterns. Paulse113 investigated the effect of mobility ratio and distance from the pilot to the edge of the reservoir on the areal sweep in a single "inverted" five-spot (Fig. 1). However, none of these studies considered the effect of operating conditions on pilot flood behavior. Moreover, hardly any of the work to date dealt explicitly with the oil-recovery performance of the pilot floods. The purpose of this paper is to introduce the concept of interrelation between pilot flood behavior and operating conditions, to illustrate this interrelation by means of experimental data and to compare the recovery performances of dif-

By M. A. Torcaso, P. Raimondi

To study mass transport in systems simulating oil recovery processes, different porous media were saturated with a mobile (carrier phase) and a stationary phase. Slugs of carrier phase containing a small amount of solute were displaced with pure carrier phase. By analogy to the chromato-graphic processes, the velocity of the solute can be predicted from a knowledge of the partition coefficient and the saturation provided that equilibrium between the two phases exists. Equilibrium was found to exist for different porous media, solutes and rates. The conditions were varied over the range normally encountered in the laboratory and in the field. The longitudinal dispersion of a solute undergoing interphase mass transfer was also investigated. INTRODUCTION The production of hydrocarbons by gas cycling, enriched gas drive and C02 or alcohol displacement involves, among other factors, relative motion between two phases and compounds, hereafter called solute, which are soluble in both phases. The solute is carried forward by the faster flowing phase at a lower velocity than the average velocity of that phase. Retardation of the solute is caused by chromatographic absorption and desorption in the slower flowing phase and by the degree of departure from equilibrium. At equilibrium the concentration of solute in the two phases can be related by the equation* CSw = K Cso ,............(1) where C, and C are the concentration of solute in the aqueous and oleic phases respectively and K is the equilibrium ratio, or partition coefficient. Displacement theories must contain an explicit assumption with regard to equilibrium, i.e., whether the compositions can be related by Eq. 1. The existance of equilibrium depends, in general, on the relative velocity between the phases. Unfortunately, other factors such as gravity segregation and viscous fingering, also depend on velocity. For this reason, whenever effects of rate on displacement were observed,2-4 it was practically impossible to discern what caused them — lack of equilibrium or the factors mentioned above. Equilibrium between phases has been the subject of extensive studies in fields such as extraction or chroma to graphy. It has received only small attention in flow through the type of porous media encountered in oil production.5 For this reason a method was developed which makes it possible to study the movement of a solute as it is affected by rate, type of porous media, partition coefficient and carrier phase, but in the absence of segregation or fingering. The information obtained enables one to determine when the assumption of equilibrium can be made. Briefly, the method consists of (1) saturating the core with a mobile and an immobile phase, (2) injecting a slug made up of the same fluid as the mobile phase and a small concentration of mutually soluble solute, (3) measuring the lag and the peak height of the slug at arrival and (4) correlating these variables with fluid properties such as partition coefficient and mixing constants of the medium. PROPOSED MECHANISM The principles of chromatography are combined with the equation of longitudinal mixing to predict the velocity of a solute slug compared to the bulk velocity and the peak height of a slug.6-9 The equation so obtained is valid under equilibrium conditions only. Therefore, a comparison between experimental and predicted results will give a measure of departure from equilibrium. This work was done with either the oleic or the aqueous phase being immobile. For simplicity, the following development is based on the case where the oleic phase is immobile. However, the treatment is the same in either case.

Jan 1, 1966

By N. Mungan

Laboratory water floods were performed in oil-wet and waterwet alundum and Torpedo cores, displacing a refined oil with n-hexylamine or Triton X-100 solution. Also, some floods were performed in which a sucrose solution was displaced with n-butyl alcohol. The Purpose ot the tests was to see if the oil recotlery could be increased and to examine what role the inter-facial tension reduction and wettability change play in the recovery mechanism. it was found that the chromatographic transport of amine was influenced by core wettability. In oil-wet cores, the rate of advance of the amine band could be predicted from equilibrium chromato-graphy, while in water-wet cores the amine band moved faster than predicted, indicating a non-equilibrium process. By reducing the interfacial tension to 1.1 dyne/ cm, oil recovery was increased. More oil was recovered from Dri-filmed cores than from water-wet cores. Reversing the wettability of the porous media from oil-wet condition also resulted in some additional oil recovery. Neutral wettability floods did not increase oil recovery. INTRODUCTION Interfacial forces in petroleum reservoirs are responsible for retention of large quantities of residual oil. Increasing exploration costs have created an incentive to attempt recovery of this residual oil by new and novel processes. One novel process involves changing the interfacial forces by introducing a chemical into the reservoir during water flooding. Recently, several investigators have carried out laboratory displacement tests using amines. lt3.4 The amines adsorbed onto initially water-wet core surfaces, changing them to oil wet. Additional oil was recovered when the cores were flooded in a manner reversing the. wettability from oil-wet conditions. In all these studies, the wettability reversal was considered responsible for increased oil recoveries although it was noted that amines also reduced oil-water interfacial tension. In the present study, displacement experiments were performed in water-wet and Dri-filmed cores. The purpose was to separate the effects of lowering the interfacial tension from changing the wettability. LABORATORY STUDY All experimental data given in this paper have been obtained at a temperature of 70F + 1 and prevailing atmospheric pressure. FLUIDS AND MATERIALS The oil used in these experiments* was washed with H2S04 and passed through two 100-200 mesh silica-gel beds to remove unsaturates and surface active impurities. Normal hexylamine** of practical purity was used, and the alundum cores*** were composed of nearly pure A1203. Fused quartz (SiOi) and Lucalox (A12 )t plates were used in the contact angle measurements. DISPLACEMENT TESTS Displacement tests were performed in alundum and Torpedo cores. Before displacement runs, cores were acidized with 0.1 normal HCI, flushed with distilled water and fired at 1,200F to make them water wet. This procedure also reduced water sensitivity of the Torpedo cores. Oil-wet cores were prepared from water-wet cores by flowing several pore volumes of 1 per cent by weight Dri-film SC-87 solution in heptane through the cores, flushing the cores with nitrogen, evacuating, and finally heating the cores at 300F for 4 or 5 hours. Between floods, cores were cleaned by conventional procedures and treated again as above.

Jan 1, 1965

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

Jan 1, 1965

By B. Habermann

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

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By W. D. Kruger

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

By 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&apos;! 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-

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 Frank Van Horn

IT is customary to ascribe two general modes of occurrence to natural gas, namely shale. gas which, as the name indicates, is found in&apos; shale, and reservoir gas, which occurs ill sandstone, conglomerate and dolomitic limestone. Both of these types are found in the Cleveland district. The shale gas has been produced on a small scale for a long period of years, probably dating back as far as 1883. Shale gas is a low-pressure gas, small in volume with variable pressure, and is limited to no definite horizon. It is generally independent of oil, may last for a long time, and is not dependent to any great extent on the structure of the rocks. There are many such wells in the district which range from 400 to 1,840 ft. in depth, although most have been completed at about 800 ft. This indicates that the wells pass through the Cleveland and Chagrin shales into the Huron or Portage shale of the Upper Devonian. During the last 2 years great interest has developed in a deeper-seated reservoir gas. What is commonly called reservoir gas is generally a high-pressure gas occurring in large volume at a definite horizon. Oil may, or may not, be associated with the gas, and in most regions in which reservoir gas is found, the geological structure is of great importance. This type of gas occurs in all large fields such as those of West Virginia which have been supplying Cleveland until the recent discovery of local gas at great depths in commercial amounts.

Jan 1, 1917