Search Documents

By John C. Calhoun, Raymond G. Loper

Theoretical calculations for the decline of pressure and the variation of instantaneous producing gas-oil ratio with increased cumulative production have been made for reservoir systems under various producing drawdown values and for various spacings. These variables have been included by a modification of previously Published techniques. This consists essentially of utilizing a radial system flowing gas and oil in a succession of steady state calculations. The results of the calculations do not show any variation in the predicted behavior of the internal gas drive reservoir when either the drawdown or the spacing are changed. The prediction of performance on an internal gas drive reservoir has been presented by various authors. ',"' These predictions are in the form of a pressure decline and an instantaneous producing gas-oil ratio as a function of cumulative oil withdrawal. In making these predictions a basic assumption has been that the saturation of the reservoir is uniform throughout at any time, or that the pressure differential imposed on the reservoir system is zero. Because of these assumptions the validity of such predictions has been questioned. Their application to actual systems, where saturations are not uniform and where producing differentials exist, has been a hazardous procedure. The authors cited have used different methods of approach to the derivation of equations for making their predictions and for arriving at a solution of both equations, but basically the same relationships are utilized in each case. In all instances the instantaneous producing gas-oil ratio is calculated from the equation: In utilizing this equation an average reservoir pressure and an average reservoir saturation must be known, the former to define ß, v, s, and uo/ug and the latter to define kg/ko. Predictions can be made, therefore, only for the conditions of a uniform reservoir saturation and a uniform reservoir pressure. This restriction of requiring a uniform oil saturation can be resolved by applying two-phase steady state flow for the radial system to the reservoir's history. Radial steady state flow is defined when the ratio of phases flowing,

Jan 1, 1949

By M. Muskat

Calculations are reported on the differential sensitivity of the updip invasion of oil strata of varving permeability to the driving pressure differential. It is assumed that the water-oil interfaces advance with sharp fronts, or that the microscopic displacement efficiency is independent of rate and capillary pressure effects. It is found that, under such conditions, limitations of withdrawal rates will not greatly inhibit the "fingering" and bypassing tendency through high permeability strata unless the differential fluid head across the original oil column is of the order of or exceeds about 80 per cent of the driving pressure differential. In practice, such restrictions will be feasible only for oil reservoirs and aquifers of high permeability and dip, but may develop automatically under combined gas injection and water drive operations. INTRODUCTION A commonly expressed view regarding the intrusion of edge-waters into oil producing formations is that the geometry of such advance is rate sensitive. Specifically, it has been suggested that high rates of advance, stimulated by rapid withdrawals, will lead to "fingering" and an irregular type of invasion, whereas at low withdrawal rates the water encroachment will be more uniform. Two factors tend to support this view. One is the resultant effect of capillary pressure and gravity gradients in the displacement processes' involved in water intrusion, which are magnified at low rates of flow so as to improve the recovery efficiency. The other is the direct retarding effect of the rising water head as it advances upslope into the producing area. Capillary pressure forces give rise to a rate sensitivity of the oil recovery by water displacement in uniform strata as well as differential effects in non-uniform formations. These are basically independent of the effect of gravity and would also occur in horizontal beds or when the oil and water have the same density. Unfortunately, however, the fundamental equation describing the displacement process under the action of capillary pressures can be solved only by numerical or graphical procedures, and moreover the solutions are sensitive to the character of the capillary pressure and relative permeability curves.* On the other hand, the direct effect of gravity can be given a formal analytical solution if the simplifying assumption is introduced that the water advances with sharp fronts in the individual uniform strata. It depends on the absolute permeabilities of the water invaded zones and hence results in a rate sensitivity of the differential invasion in stratified media. Although a tieatment of this factor alone will not give a complete evaluation of the effect of rate on water intrusion in non-uniform pays, it should provide at least a semiquantitative measure of the sensitivity of "fingering" to withdrawal rate. Such a treatment will be presented in the following sections. GENERAL THEORY The basic element of the theory is the history of water advance updip in a uniform stratum (Fig. 1). In contrast to the assumption of constant flow rate, generally made in the analysis of displacement processes, we shall assume here that

Jan 1, 1951

By M. Muskat

Calculations are reported on the differential sensitivity of the updip invasion of oil strata of varving permeability to the driving pressure differential. It is assumed that the water-oil interfaces advance with sharp fronts, or that the microscopic displacement efficiency is independent of rate and capillary pressure effects. It is found that, under such conditions, limitations of withdrawal rates will not greatly inhibit the "fingering" and bypassing tendency through high permeability strata unless the differential fluid head across the original oil column is of the order of or exceeds about 80 per cent of the driving pressure differential. In practice, such restrictions will be feasible only for oil reservoirs and aquifers of high permeability and dip, but may develop automatically under combined gas injection and water drive operations. INTRODUCTION A commonly expressed view regarding the intrusion of edge-waters into oil producing formations is that the geometry of such advance is rate sensitive. Specifically, it has been suggested that high rates of advance, stimulated by rapid withdrawals, will lead to "fingering" and an irregular type of invasion, whereas at low withdrawal rates the water encroachment will be more uniform. Two factors tend to support this view. One is the resultant effect of capillary pressure and gravity gradients in the displacement processes' involved in water intrusion, which are magnified at low rates of flow so as to improve the recovery efficiency. The other is the direct retarding effect of the rising water head as it advances upslope into the producing area. Capillary pressure forces give rise to a rate sensitivity of the oil recovery by water displacement in uniform strata as well as differential effects in non-uniform formations. These are basically independent of the effect of gravity and would also occur in horizontal beds or when the oil and water have the same density. Unfortunately, however, the fundamental equation describing the displacement process under the action of capillary pressures can be solved only by numerical or graphical procedures, and moreover the solutions are sensitive to the character of the capillary pressure and relative permeability curves.* On the other hand, the direct effect of gravity can be given a formal analytical solution if the simplifying assumption is introduced that the water advances with sharp fronts in the individual uniform strata. It depends on the absolute permeabilities of the water invaded zones and hence results in a rate sensitivity of the differential invasion in stratified media. Although a tieatment of this factor alone will not give a complete evaluation of the effect of rate on water intrusion in non-uniform pays, it should provide at least a semiquantitative measure of the sensitivity of "fingering" to withdrawal rate. Such a treatment will be presented in the following sections. GENERAL THEORY The basic element of the theory is the history of water advance updip in a uniform stratum (Fig. 1). In contrast to the assumption of constant flow rate, generally made in the analysis of displacement processes, we shall assume here that

Jan 1, 1951

By C. C. Miller, C. A. Hutchinson, A. B. Dyes

This paper presents a simple and practical method for the estimation of effective permeability and reservoir pressures from the rise in bottom hole pressure when a well is shut in. The equations on which the method is based consider a stabilized well producing from a uniform formation when flow is shut off at the sand face followed by the unsteady state flow of a compressible liquid from regions far from a well to points closer in during the build-up period. An approximation for extending the method for two-phase flow in the reservoir is also presented. In order to compare the effects of extremes in the conditions existing at the external boundary of the reservoir, two general equations defining the build-up characteristics of the reservoir are shown: first, for a reservoir with constant pressure maintained at the radius of drainage, and second, for a reservoir with no influx over the radius of drainage. These equations define the build-up characteristics for the shut-in well as functions of time, the permeability and porosity of the drainage area, the viscosity and compressibility of the fluid flowing, the production rate prior to shut-in, and the radius of drainage of the well. When a well is shut in at the surface, rather than at the sand face, flow will continue into the well bore after shut-in due to the compressibility of the fluids in the bore. The effect upon the build-up characteristics of this flow into the well after shut-in has been studied by means of the electrical analyzer. From the basic relationships and consideration of the limitations due to continued flow after shut-in, methods have been developed for the estimation of the effective permeability of an area from the plot of the build-up characteristics of the well producing that area. Likewise, equations are presented for the estimation of the reservoir pressure from an early portion of the build-up curve. By means of the electrical analyzer, studies have been made concerning the variations in the build-up characteristics resulting from damaged or improved permeability in the immediate vicinity of the well bore. These studies show that the permeability estimated from the build-up curves represents the average effective permeability of the majority of the drainage area and that the permeability in the immediate vicinity of the well does not affect this estimate. Since the major portion of the drawdown is in the immediate vicinity of the well, variations in the permeability of this area greatly affect the productivity of the well as expressed by the productivity index. From pressure build-up and productivity index tests the ratio of the permeability of some arbitrarily specified area immediately surrounding the well and the effective permeability in the majority of the drainage area can be computed. This ratio is a measure of any damage to the well incurred in drilling or reworking and likewise any benefits as a result of remedial measures. In build-up tests and many other tests that are made on a producing well it is at least desirable, if not absolutely necessary, that the well be stabilized prior to testing. A method is presented here for relating the time required to stabilize a well and the reservoir and fluid properties. An equation is also presented for approximation of the length of shut-in time required to obtain sufficient build-up data for the estimation of permeability and reservoir pressure.

Jan 1, 1950

By Joseph Fry, Ralph H. Espach

In the spring of 1943, when it was evident that the Tensleep bandstone in the Elk Basin Field, Wyoming and Montana, held a large reserve of petroleum, Bureau of Mines engineers obtained samples of oil from the bottom of nine wells and analyzed them for such physical characteristics as the volumes. of gas in solution. saturation pressures or bubble points, shrinkage in volume caused by the release of gas from solution, expansion of the oil with decrease in pressure, and other related properties. The composition of the gas in solution in the oil was studied. The pressures and temperatures existing in the reservoir and the productivity characteristics of the oil wells were determined. The data obtained indicate that the oil in the Tensleep Reservoir of the Elk Basin Field has unusually varying physiral characteristics, such as a saturation pressure of 1,250 psia and 490 cu ft of gas in solultion in a barrel of oil at the crest of the structure and a saturation pressure of 530 psia and 134 cu ft of gas in solution in a barrel of oil low on the flanks. The hydrogen sulfide content of the gas in solution in the oil varies from 18 per cent for oil on the crest to 5 per cent for oil low on the flanks of the structure. Of even greater significance is the fact that these and other variable characteristics of the reservoir oil are related to the position of the oil in the structure. Many geologists and petroleum engineers have considered that all the oil in a petroleum reservoir has rather uniform physical characteristics and that equilibrium conditions prevailed in all underground accumulations of oil and gas; that such is not always so is borne out by the results of the study by the writers. INTRODUCTION The Rocky Mountain region is one in which may be found striking examples of the unusual in oil and gas accumulations, as is evident from the following: The high helium content (7.6 per cent) of the gas in the Ouray-Leadville limestone sequence in the Rattlesnake Field, New Mexico, and gases of similar helium content in other fields; 50" to 55' API gravity distillate in solution in carbon dioxide gas and recoverable through retrograde condensation, in the North McCallum Field, Colorado; the occurrence of gas, oil, or both in closely related structures contrary to the usual concepts of gravimetric segregation: the accumulation of gas and/or oil in structures closely related to other structures that apparently are more favorable but do not contain oil or gas accumulations; the high hydrogen sulfide content (as high as 42 per cent) of the gas associated with oil in some fields in the Big Horn Basin, Wyoming; and the wide range of fluid chararteristics found in the Elk Basin reservoir. Elk Basin, an interesting old oil field that has been producing oil from the Frontier formation since 1915, is situated in a highly eroded basin resulting from the erosion of the crest of an anticline and some of the underlying softer shales. The field came back into national prominence during 1943 when it became known that it was the largest single reserve of new oil discovered in the United States that year. The Tensleep sandstone was found to contain oil in November. 1942, when a well drilled to a depth of 4,538 ft (44 ft into the Tensleep sandstone) flowed oil at the rate of 2,500 B/D. By the end of 1949, 137 oil-producing wells and five dry holes had been drilled, and approximately 32 million bbl of oil had been produced. Approximately 6,000 acres may be considered productive of oil in the Tensleep Reservoir, and estimates of the oil that will be produced average 200 million bbl. The Tensleep Reservoir has further interest because it ha-greater closure than any oil field in the Rocky Mountain region; the closure of the Elk Basin anticline is variously estimated at 5.000 to 10,000 ft. of which the top 2.00 ft of the structure contained oil. SUBSURFACE OIL SAMPLING Fig. 1 is a structural map of the Elk Basin Tensleep Reservoir, on which the nine wells used in this study and the numbers correvponding to the well designations hereafter referred to are shown. Wells 1. 2, 3, 4, and 8 were tested and sampled during October and November. 1943. and Wells 5, 6. 7, and 9 during June and July, 1944. An electromagnetic type sampler developed by the Bureau of Mines and described by Grandone and Cook' was used in obtaining the subsurface oil samples. As the wells were tubed nearly to bottom, the sampler was run as far as possible in the tubing hut never below the top perforations. The following procedure was used in testing and sampling the wells: A well was shut in for at least three days, after

Jan 1, 1951

By Joseph Fry, Ralph H. Espach

In the spring of 1943, when it was evident that the Tensleep bandstone in the Elk Basin Field, Wyoming and Montana, held a large reserve of petroleum, Bureau of Mines engineers obtained samples of oil from the bottom of nine wells and analyzed them for such physical characteristics as the volumes. of gas in solution. saturation pressures or bubble points, shrinkage in volume caused by the release of gas from solution, expansion of the oil with decrease in pressure, and other related properties. The composition of the gas in solution in the oil was studied. The pressures and temperatures existing in the reservoir and the productivity characteristics of the oil wells were determined. The data obtained indicate that the oil in the Tensleep Reservoir of the Elk Basin Field has unusually varying physiral characteristics, such as a saturation pressure of 1,250 psia and 490 cu ft of gas in solultion in a barrel of oil at the crest of the structure and a saturation pressure of 530 psia and 134 cu ft of gas in solution in a barrel of oil low on the flanks. The hydrogen sulfide content of the gas in solution in the oil varies from 18 per cent for oil on the crest to 5 per cent for oil low on the flanks of the structure. Of even greater significance is the fact that these and other variable characteristics of the reservoir oil are related to the position of the oil in the structure. Many geologists and petroleum engineers have considered that all the oil in a petroleum reservoir has rather uniform physical characteristics and that equilibrium conditions prevailed in all underground accumulations of oil and gas; that such is not always so is borne out by the results of the study by the writers. INTRODUCTION The Rocky Mountain region is one in which may be found striking examples of the unusual in oil and gas accumulations, as is evident from the following: The high helium content (7.6 per cent) of the gas in the Ouray-Leadville limestone sequence in the Rattlesnake Field, New Mexico, and gases of similar helium content in other fields; 50" to 55' API gravity distillate in solution in carbon dioxide gas and recoverable through retrograde condensation, in the North McCallum Field, Colorado; the occurrence of gas, oil, or both in closely related structures contrary to the usual concepts of gravimetric segregation: the accumulation of gas and/or oil in structures closely related to other structures that apparently are more favorable but do not contain oil or gas accumulations; the high hydrogen sulfide content (as high as 42 per cent) of the gas associated with oil in some fields in the Big Horn Basin, Wyoming; and the wide range of fluid chararteristics found in the Elk Basin reservoir. Elk Basin, an interesting old oil field that has been producing oil from the Frontier formation since 1915, is situated in a highly eroded basin resulting from the erosion of the crest of an anticline and some of the underlying softer shales. The field came back into national prominence during 1943 when it became known that it was the largest single reserve of new oil discovered in the United States that year. The Tensleep sandstone was found to contain oil in November. 1942, when a well drilled to a depth of 4,538 ft (44 ft into the Tensleep sandstone) flowed oil at the rate of 2,500 B/D. By the end of 1949, 137 oil-producing wells and five dry holes had been drilled, and approximately 32 million bbl of oil had been produced. Approximately 6,000 acres may be considered productive of oil in the Tensleep Reservoir, and estimates of the oil that will be produced average 200 million bbl. The Tensleep Reservoir has further interest because it ha-greater closure than any oil field in the Rocky Mountain region; the closure of the Elk Basin anticline is variously estimated at 5.000 to 10,000 ft. of which the top 2.00 ft of the structure contained oil. SUBSURFACE OIL SAMPLING Fig. 1 is a structural map of the Elk Basin Tensleep Reservoir, on which the nine wells used in this study and the numbers correvponding to the well designations hereafter referred to are shown. Wells 1. 2, 3, 4, and 8 were tested and sampled during October and November. 1943. and Wells 5, 6. 7, and 9 during June and July, 1944. An electromagnetic type sampler developed by the Bureau of Mines and described by Grandone and Cook' was used in obtaining the subsurface oil samples. As the wells were tubed nearly to bottom, the sampler was run as far as possible in the tubing hut never below the top perforations. The following procedure was used in testing and sampling the wells: A well was shut in for at least three days, after

Jan 1, 1951

By C. Özgen Karacan

Degasification and methane control in underground coal mining is an important area augmenting conventional ventilation in order to prevent possible fires and explosions due to excessive methane emissions. Coal bed reservoir engineering has been applied for many years to coalbed methane production. In recent years, modeling the flow and draining methane, not only from coal bed itself but also from the overlying strata, have been recognized as important areas of expertise for supplementing underground ventilation. This paper demonstrates the applications of coal bed reservoir engineering and modeling techniques for optimizing and controlling methane emissions in longwall and continuous miner sections. Practical applications and some considerations are presented for reservoir modeling, borehole performances, and face emission predictions for various coal bed and borehole configurations. A section on field reservoir studies and well bore log analyses for methane control are also presented as a complementary way of controlling methane and optimizing ventilation.

Jan 1, 2009

By J. S. Aronofsky, A. S. Lee

Results and experimental procedures are presented covering a preliminary laboratory investigation of the compaction of reservoir rocks and its effect on porosity and permeability. "Egective" compressibility was measured by subjecting core material to an overburden load and reducing pore pressure. Preliminary results showed relatively high effective compressibilities which varied from 1.8 X 10-' pore volume/pore volume /lb / in.~ for limestones to 7.7 X 10-6 pore volume/pore volume/lb/in.' and even higher for sandstone reservoirs conmining large amounts of fine clay material. It appeared from the initiol results that the effective compressibility of the material could not be correlated directly with porosity and that other factors, such as the amount of clay material present might have to be taken into account. The work to date has indicated that effective compressibility factors will have to be measured for the particular reservoirs heing studied. Oil permeability of core materid was measured under various effective overburden loads from 0 to 8,000 psi. It has been found than reduction of permeability in clean sands was relatively small and the results agreed with published data. For sandstones containing large amounts of clay, preliminary tests indicated very large reductions in permeability with increasing effective overburden pressures. It is planned to continue laboratory work on this phase because of its importance in low-permeability reservoirs. INTRODUCTION During the past few years increasing interest has been shown in the differences obtained between conventional laboratory measurements on cores at atmospheric pressure and those when reservoir conditions are simulated. Since these data are used as the basis for calculations and the prediction of reservoir behavior, their evaluation under reservoir conditions is most important. The effect of the compressibility of the rock upon the reservoir behavior, as fluid pressure declines, is an important subject which has been insufficiently studied. The total pressure on any plane through a reservoir is the result of the grain-to-grain rock pressure due to the overburden and the pressure of the interstitial fluid on the plane. As reservoir fluid pressure declines, the net load on the rock phase must increase, and the result is the compaction of the rock structure which produces a decrease in pore volume. Simultaneously, the rock grains expand into the pore space as the fluid pressure declines. The result of these two effects is defined as the effective compressibility of the rock and is expressed in units of pore volume/original pore volume/psi. Work was started by Imperial in 1954 on the effects of pressure on the porosity of reservoir rock. These studies established effective compressibility factors for some limestones and dolomites. Since that time effective compressibilities have been obtained for other Canadian and some U. S. reservoir rocks. In conjunction with the study of porosity, an investigation has also been initiated to determine the effects of pressure on permeability. One outcrop and several reservoir sandstone cores have been tested to date. ROCK COMPRESSIBILITY Apparatus anD Procedure A diagram of the apparatus used in these studies is shown in Fig. 1. The equipment was constructed so that a uniform pressure, equivalent to the reservoir overburden, could be applied by the oil pump and transmitted on all sides of the core through the Lucite sheath. The sample was completely saturated by circulating brine through the mounted core. Pressure in

By C. H. Stewart, R. A. Fagin

A two-dimensional, three-phasereservoir simulator was programmed for a large memory digital computer. It was designed to provide a practical solution to describing the complex physical relation between the natural forces and the physical properties of a heterogeneous reservoir when subjected to a specific set of conditions. A reservoir study is briefly described to illustrate application of the model. A full volumetric account of three phases (oil, gas and water) is performed simultaneously throughout an integration net representing the reservoir. Absence of one or two of the phases is treated as a special case of the more general situation. Expansion (or contraction) of all phases, including rock expansion, is performed so that the pressure calculation is the general unsteady-state case. To account for the large variations of subsurface elevation encountered in some reseruoirs, and to allow for segregation of the various phases, a gravity head term is included in the basic drive potential. Appropriate fluid and rock properties are used in polynominal surface form (functions of pressure and/or depth) or they can be entered as space variables at each position of the integration net. An unsteady-state water influx calculation, based on the method of van Everdingen and Hurst, l was connected to the boundary of the matrix to simulate aquifers of various sizes. In addition to reservoir calculations, three-phase fluid flow from the producing depth to the wellhead, including provisions for gas lift, was incorporated in the simulator. A workover routine was also built which can automatically switch to a different set of production relations when a gas-oil ratio or water fraction reaches a limit; or it can shut-in the well if prescribed. INTRODUCTION This paper describes a reservoir engineering mathematical simulator used to represent the complex interaction of natural forces and physical properties of a reservoir during natural depletion or with various injection schemes. The simulator, which was programmed for a large memory digital computer, is a two-dimensional calculation which handles three mobile fluid phases simultaneously (oil, gas and water). Basic requisites for the method are individual well production and pressure data, hydrocarbon fluid properties, geological data (producing depth and net sand), capillary pressure data, relative permeability data and permeability and porosity information. Matching the past performance of a combination drive reservoir often has yielded information concerning continuity and the validity of basic data. Detailed predictions of future performance can be made for continuation of current depletion methods (natural depletion) as well as for various types of recovery by gas or water injection. Combination injection cases and pattern studies can also be performed. Workover programs, gas lift and different types of artificial lift programs have been investigated using a technique similar to that described by Kern and Nicholson2 except that conditions of pressure and saturation at the block within which the well is located are used rather than average reservoir conditions. Drilling additional wells to optimize profit was explored, both as to number and location, by placing wells at different spots within the reservoir matrix. Special depletion processes can be examined, such as upstructure drainage and lateral (or strike) waterfloods in thin oil columns. In one case the mathematics of the simulator were modified to calculate the displacement in the vertical plane rather than in the horizontal plane. In this manner specific reservoir problems can be studied, such as coning of gas and/or water around production points, fingering along permeable stringers or, more generally, frontal advances in a heterogeneous section.

Jan 1, 1967

By Mohamed Mortada

A practical analytical method is developed in this paper which provides the practicing reservoir engineer with a handy method for analyzing oilfield interference problems. The procedure employed entails the principle of superposition of individual eflects. A set of charts has been developed and can be used to evaluate the individual effects. The charts represent non-steady state pressure distribution in an extensive aquifer due to a constant rate of production from an oilfield which can be approximated by circular geometry; they are constructed from solutions to the diffusivity equation considering a circular sink as distinguished from a point sink. The solutions are obtained partially from mathenlatical considerations and partially from analog and digital computations. By virtue of the method of images, the charts can be used also to treat the case of production near a fault. In addition, provisions are made for treating oil fields, which are produced by partial water-drive. Treating an actual field case is deferred to a subyequent writing; however, illustrative example is presented to demonstrate the method of application. The pressure distribution in linear and limited radial aquifers is currently under study. INTRODUCTION With the increasing difficulty of discovering new oil reserves, attention is necessarily focused on the efficient development and production of existing reservoirs. Efficient production depends on the availability of reliable means for predicting the pressure-production behavior of oil fields. The prediction is usually performed by analytical methods or computers. Analytical methods are often limited to systems which can be approximated by simple geometry. However, analytical methods have the distinct advantage of being readily performed without special equipment. Analog and digital computers lend themselves to a variety of problems, but are not always locally available and consequently their use is limited. One particular situation in which analytical techniques are needed involves studies of multiple, water-drive oil fields sharing a common aquifer. Such oil fields are in hydrodynamic communication. Therefore, production from any of the oil fields is accompanied by a pressure decline which is transmitted through the aquifer to the other oil fields and manifested as pressure interference. The rate of propagation of the pressure decline is such that the pressure may be significantly reduced many miles away from a producing pool. An example of pressure interference is found in the Woodbine formation in East Texas. There. a dozen or more pools, including several of major size, draw on the Woodbine formation for water drive. The pressure in the water-bearing portion of the formation has been significantly reduced over an area of several thousand square miles. This reduction of pressure has been such that when the Hawkins pool was discovered in 1940, the pressure in it had already been reduced 280 psi below the initial value in the aquifer. This resulted from prior production from other Woodbine sand pools, including, in particular, the East Texas field some 20 miles to the east of Hawkins.' The above example indicates that the interference effect can be significant and should be taken into con-

Jan 1, 1956

By R. N. Hoffman

A new technique for determining capillary pressure curves has been developed and tested. The technique differs from previously reported centrifuge techniques in that the centrifuge is slowly accelerated from zero to the maximum desired speed rather than being held constant at particular, progressively higher speeds. An important advantage of this technique over other methods of determining capillary pressure curves is the short time required to obtain the desired amount of data over the chosen pressure range. For example, to obtain 30 data points between 1.2 and 104 psig with a 1.55-in. long, 3/4-in. diameter core using a brine-air system, 6.6 hours were required with this technique. An equally important development of this paper is an analytic method for the conversion of the data from the centrifuge experiment to capillary pressure curve data. Previously there has been only an approximate conversion available. Although the capillary pressure curves determined by this technique appear to be as accurate as those determined by other techniques, the accuracy could be improved if certain variables, not treated in this experiment, were investigated. Among these are the dynamic distortion of the centrifuge equipment and imperfect initial saturation of the cores. INTRODUCTION Pirson 1 defines capillary pressure in a porous medium as "the differential pressure that exists between two fluid phases at their interfaces when they are distributed under static equilibrium within a porous material". Capillary pressure in rocks is known to be a function of fluid saturation, among other things, and a capillary pressure curve is defined for the purposes of this paper as a plot of the capillary pressure-wetting phase saturation relationship for a particular rock sample. Several methods are used for determining capillary pressure curves for small rock cores. Prominent among these are the semi-permeable barrier, mercury injection, and, to a lesser extent, centrifuge methods. The semi-permeable barrier method is currently the most popular. It features simplicity in both execution and the mathematical conversion of the experimental data into a capillary pressure curve. The main disadvantages of the semi-permeable barrier method are the time required — as long as two months — to obtain several points for the curve, and a fairly low maximum pressure before breakthrough of the non-wetting phase into the barrier occurs, for example, about 32 psig for a brine-air system. It is for these reasons that other methods such as the centrifuge method have been introduced. High accelerations and the absence of a barrier result in quicker attainment of saturation equilibrium at a given pressure. However, the centrifuge method involves much more expensive equipment and more difficult procedures and calculations than does the barrier method. The purpose of this investigation has been to improve the equipment and procedures of the centrifuge method and to develop an analytic method for the conversion of the experimental data into a capillary pressure curve. Hassler and Brunner did the original work in the determination of capillary pressure using a centrifuge. In their work the centrifuge speed was increased in a step-wise manner, each speed being held constant until saturation equilibrium was reached in each core. Saturation equilibrium was indicated when the volume of liquid collected in the graduated pipette of the core holder remained constant. According to Hassler and Brunner, equilibrium was reached in "a few minutes to one-half hour or more". In the centrifuge method, as opposed to the barrier method, the fluid saturation of the core is not a constant throughout the length of the core, but varies with the radius of centrifugation. Also, the capillary pressure cannot be read directly but must be calculated from a knowledge of the centrifuge speed and other parameters. For these calculations, Hassler and Brunner derived the following equations:

By Harold B. Janzen, M. A. Nobles

Campbell and Martin1 described a new exploratory logging tool which was developed as the result of a new concept of the phenomenon that takes place when a porous formation containing mobile hydrocarbons and irreducible interstitial water is invaded by mud filtrate. They postulate that invasion will create a bank of interstitial water between the flushed zone and the undisturbed formation. This postulate was confirmed by Wyllie.' Thus, if R,, is greater than R, and mobile hydrocarbons are present, the resistivity profile with increasing radii from the borehole is high:low:high, as shown in Fig. 1. If mobile hydrocarbons are absent, the bank of interstitial water is indistinguishable from the undisturbed formation and the resistivity profile, when R is greater than R,, is high:low, as shown in Fig. 2. The zone containing the bank of interstitial water and residual oil was named the "low zone". Campbell and Martin conducted extensive electrolytic wedge and resistivity analog model studies to select a logging system that would reflect the presence or absence of the low zone. The system they selected was that of a family of lateral curves. The model studies indicated that the presence of a low zone could be detected in approximately 60 to 70 per cent of the producing sand reservoirs with a family of lateral curves provided that the sonde effectively occupies 90 per cent of the diameter of the hole. These studies were conducted with D, = 5d. The presence or absence of the low zone is best determined by a plot of the apparent resistivity vs spacing as shown in Figs. 3 and 4. A dip or inflection in the profile that repeats at successive depth intervals is considered diagnostic. To reduce the shunting effect of the hole and mud, a sonde of insulating material that is surrounded by a series of flexible rubber bumpers was designed. The bumper efficiency is considered to be the ratio of the bumper diameter to the hole diameter. Field tests conducted by Campbell and Martin on 10 wells with a tool that provided for station readings indicated that dips or inflections could be detected in 19 out of 25 known producing zones. INSTRUMENT DEVELOPMENT Halliburton developed a logging system to simultaneously and continuously record eight lateral curves. The instrumentation is described in the literature.' The original logs made

By A. H. Heim

A method and apparatus for measuring gas porosities of rocks are described. The apparatus can be assembled from commercially available components. In principle, measurements are made by volume substitution at constant pressure. The maximum error is not more than 0.3 porosity per cent. Typical results are given. INTRODUCTION Determining the porosity of rock samples is one of the most important and yet most varied types of measurement in core analysis. Among the many techniques devised are the so-called "gas porosity" methods. An old and well known example is the Washburn-Bunting method.' The U. S. Bureau of Mines2-' described and later improved the apparatus for a now widely used method generally known as the "Boyle's law" method. In the present form of the Washburn-Bunting method,' the volume of air in the pores of a rock sample at atmospheric pressure is extracted and then collected in a graduated burette at atmospheric pressure. The volume of air is read directly as the pore volume of the sample. The absolute error in reading the collected volume of gas is independent of the total volume; thus, the relative error is larger when the volume is small, as it is for rocks of low porosity. In addition, the sample after measurement contains mercury, which limits its use for other analyses. The Bureau of Mines (or Boyle's law) method measures directly the solids volume of a sample from which the pore volume and porosity are derived, using a separate measurement of the bulk volume. Gas at a few atmospheres pressure is introduced into a sample chamber of known volume containing the rock sample. The pressure is accurately measured. Following, the gas is expanded into a burette at 1 atm, and the gas volume is read directly. From the initial pressure p, and the final pressure p2 and volume v,, the initial gas volume v1 is calculated using Boyle's law; that is, p1v1 = p2v2. Volume v, minus the volume of the empty sample chamber is the solids volume of the sample. The accuracy of the method is limited, unless corrections are made, by deviations of the gas from the "ideal" gas-law behavior assumed in the simple form of Boyle's law. The purpose of the present paper is to describe a method for measuring the gas porosity of a rock which avoids many of these difficulties. Gas volumes are measured directly with the same accuracy as the bulk volumes. Pressures of at least an order of magnitude larger than those of previous methods are employed to insure rapid penetration of the gas into the sample. While special equipment may be built to apply the method, the porosimeter may be constructed as well from commercially available components. For simplicity, the apparatus described will be referred to as the "Constant-Pressure gas porosimeter". THE CONSTANT-PRESSURE METHOD Fig. 1 shows schematically the arrangement of components comprising the present Constant-Pressure porosimeter. Briefly, the method is one of volume substitution and may be considered a null measurement. Omitting (for the present) some of the operational details, the method of measurement consists of the following three steps. 1. After evacuation, the volume of the measuring system (a ballast chamber, a manifold, two gauges and their connections) up to the sample chamber is filled with gas to a high pressure (- 1,000 psi). A sample of the gas at this pressure is trapped in one side of a sensitive differential pressure gauge to serve as the reference pressure for subsequent steps. 2. The evacuated sample chamber containing the rock sample is opened to the measuring system. As the gas expands into the chamber, the resulting decrease in pressure unbalances the differential pressure gauge. 3. The pressure is restored by means of a mercury volumetric pump. The volume of mercury injected exactly equals the free or void volume of the sample chamber (volume of empty chamber minus the solids volume of the rock within). From the injected volume and the known empty chamber volume, the solids volume is obtained and the porosity calculated. The pressure and the volume occupied by the gas are the same before and after opening the sample chamber. Expansion and compression of the gas are incidental operations and do not enter into the calculation of porosity. By the pressure balancing or nulling, the free volume of the sample chamber is merely substituted by an equal and measured volume of mercury. Since the measurements are at constant pressure, there are no compressibility corrections necessary for the sample chamber.

By C. S. Matthews, M. J. Fischer

Distortions caused by dip on five-spot sweep pattern were determined by visual experiments on a model composed of parallel plates of glass with holes drilled at the well positions. This model war designed to simulate a reservoir depleted of its primary oil and containing a gas phase. Prior to interference, it was found that the injected fluid took a nearly circular shape about the injection well for the cases studied. However, these circles moved down-dip considerably under the influence of gravity. It was concluded that pattern flooding should be practical for many dipping reservoirs. A quantitative measure is given of how much the well position should be altered to allow for dip. Experimental deternzinarions of injectivity showed this to be about the same in a dipping reservoir as in a horizontal one. INTRODUCTION This investigation was undertaken at the suggestion of the Pacific Coast area, Shell Oil Co., where a pilot water flood in a dipping reservoir was being carried out at the time this paper was in preparation. Prior to the undertaking of this flood there had been considerable speculation as to how gravity would distort a flood pattern in a dipping reservoir. The thought was expressed in some quarters that a five-spot flood would not be applicable to a dipping reservoir because the flood pattern would elongate in the direction of dip and would never form a "seal" around each five-spot. If the flood pattern elongated sufficiently, both the oil and the water might flow largely downdip rather than into the pattern producers. In view of the importance of the California water floods, it was deemed important to determine the effects of gravity on a five-spot pattern in a dipping reservoir. Previous work had also been carried out in this laboratory,1 but only for the single-fluid case. The present investigation extends these results to multifluid cases. Several sets of pictures of flood fronts are presented so that the process can be visualized. From these pictures, curves are obtained which show the correct drilling position for achieving equal times of oil-bank breakthrough. Several injectivity curves are also presented which may aid in prediction of injection rate. ASSUMPTIONS The present results were obtained with a model composed of two parallel plates of glass spaced about 1 mm apart. When the inertial forces are small, as in these experiments, the flow between the plates follows Darcy's law and is analogous to flow through porous media. In using these plates, it is not possible to scale effects of capillarity correctly as regards a transition zone between injected and displaced fluid. Calculations show, however, that effects of such capillarity should not be important in determining the gross flooding pattern in the field. In fact, in this model the capillary forces tend to circularize the pattern to a greater degree than in the field. To preserve a large ratio between the viscous forces and the circularizing capillary forces in the laboratory model, viswus glycerine was used as the flooding liquid and the rate of flooding was very high. The circularizing force is believed negligible in these experiments. In addition to the effect of the surface tension y, which tends to circularize the pattern, another capillary force must be considered. This is the adhesion tension y cos, which acts at right angles to the circularizing term y. If the plates in the model are "wetted" the same all around the periphery, then the adhesion tension,

Jan 1, 1957

By R. E. Collins

The existence of fluid migration across fixed boundaries in oil and gas reservoirs has been known for many years. Several techniques have been developed in the past for estimating The rate of migration across fixed boundaries as an aid in planninng field development and in the valuation of oilfied properties. The principal de/errc7nt to the use of these techniques lies in the rather extensive reservoir and field data required for ca1cula1ion.s of fluid migration. For this reason, a new, simplified procedure has been developed which makes possible the calculation of fluid migration with a minimum of field and resrvoir data. This new technique is based on certain. solutions of the differential equations describing flour in the reservoir which assume, for the portion of the reservoir of interest, that the formation can be approximated by a homogeneous rock of uniforrn thickness, that only a Single mobile fluid phase exists, and {hut fluid production at the well is solely a result of expun.rioti of rile resersoir fluds. The results of the, present work are compiled in a set of curves. These curves can he used to calculate both rule and cummulative fluid migration when the abovc assumptions ore ju.r/ifiril. The only data required for such calcu1ations are the production hisrories of all wells in the field, the permeability and porosity of the reservoir., the, compressibility and viscosity of the fluid. and the shape of the reservoir, INTRODUCTION The migration of fluids across fixed boundaries in oil and gas reservoirs has long been recognized as an engineering and economic problem, but no methods for estimating the extent of migration during the development stage of a reservoir have been published. The investigation reported here has been directed toward providing a simple technique for estimating the extent of such migration and has therefore been restricted to the simple case of a homogeneous reservoir of uniform thickncss containing compressible liquids. It is assumed that only a single mobile fluid phase exists and that fluid production at the wells is solely by expansion of the reservoir fluids. Obviously, few, if any. reservoirs conform to these assumptions during their entire productive life. On the other hand, most reservoirs approximate fluid expansion reservoirs during their initial stages of primary production. Consequently, the present work should be applicable during the field development of most reservoirs, regardless of whether or not they ultimately are gas cap, dissolved-gas, or water-drive reservoirs. The mathematical analysis of this problem has been outlined in the Appendix, while the results of the investigation are presented in the following discussion. DISCUSSION OF CALCULATION PROCEDURE SLNGLE Well, Ini-iimite Reservoir The procedure for calculating fluid migrations is based on solutions of the partial differentia1 equations describing the flow of a homogeneous compressible liquid in a homogeneous reservoir of uniform thickness. This differential equation has the form of the well-known diffusion equation. The equation, as applicable to our problem, is where p is fluid density, is porosity, EL is viscosity, c is compressibility,* and k is permeability; .r, y and t are the space coordinates and time. Derivation of this equation is given in the Appendix. The analytical solution of this equation for a single well of negligible radius producing at a constant rate. q, from a reservoir of infinite areal extent for a time, 1. can he readily obtained. From this solution the rate of flow and also the cumulative flow across a straight line segment of length, 1, at a perpendicular distance. 6. from the well can be computed. One end of the line must be on the perpendicular joining the well to the

By J. M. McDowell

A study has been made to deter~r~irre how the behavior of a water drivee reservoir changes as a function of the permeability of the formation and as a function of the size of the aquifer. The effect of pressure maintenance programs on the rote of nat11rn1 water influx is also studied, (IS well as the influence on the influx of the position of water injection we1ls with reference to the oil-water contact. The results are determined by the use of an electrical reservoir analyzer and are plotted as families of curves of pressure, natural water influx, and water injection vs time or cumulalive production. Two size aquifers are studided. One is limited with a 36-inile radius, and the, other approximates an infinite reservoir. The field remains the same in all respects for the completed. study except that the permeability is changed to agree with that assumed for the aquifer. 1 N TR ODUCTI ON this study was undertaken in -dcr to determine how the pres- sure decline and water influx of a water drive field would change as a function of the rate of fluid withdrawals, as a function of the permeability of the formation, and as a function of the size of the aquifer. The effect of pressure maintenance programs on the rate of natural water influx was also studied as was the position of water injection wells with reference to the oil-water contact. The information has been determined by the use of an electrical reservoir analyzer1, sometimes referred to as a Carter type analyzer. This instrument is an analog type computer which can be set up to represent a water-drive oil field. It can reproduce the past performance of the field in terms of analogous electrical quantities when the proper circuit analog of resistors and condensers has been determined, also, it can determine the future behavior for any method of production which may be selected, or permit the comparison among several different methods. Although most applications of the reservoir analyzer"" have been to thc study of specific field performance, it is intended to use the analyzer here to help understand how various reservoir characteristics effect the final performance of the field and to sec under what conditions water injection will be of benefit and to what extent. Within the oil industry the performance of water-drive reservoirs was first described analytically by LV. Hurst' and M. Muskat2 for certain types of boundary conditions. The solutions of the equations developed involved products of exponential functions and Bessel functions. Later publications by Hurst4 and van Ever-clingen" and the hooks of Muskat6 gave further information on reservoirs with other boundary conditions. They also presented rather extensive curves and tables to help evaluate the analytical expressions involved. A rccent series of publications by Chatas. bascd on the work of Hurst and van Everdingen, gives quite complete charts and tables of data to evaluate the pressure decline or water influx in infinite or limited radial reservoirs. and in linear reservoirs. The curves in all of these publications are given in dimensionless form

Jan 1, 1956

By A. Lubinski, C. R. Stewart, K. A. Blenkarn

Storage porosity has been considered one of the important pore geometry characteristics of heterogeneous-porosity limestones. Storage pores are only containers for fluids, in contrast to flow channel pores which both contain fluids and provide continuity for fluid flow. The concept of another geometry characteristic, "poro-striction", is presented as a second pertinent variable in describing limestone pore space. In simple terms, poro-striction is a measure of the flow resistance between storage and flow channel pores. Alternating-flow core-testing theory provides flow relationships which can be used to divide the pore space of heterogeneous-porosity media into flow channel and storage pores and to measure the "porostriction" of the latter. Experimental application of this theory to naturally occurring heterogeneous limestones shows that "porostriction" and the ratio of storage to flow channel pores can be estimated. Porostriction and porosity ratio are microscopic characteristics which should influence oil-recovery efficiencies during certain types of displacement processes. A knowledge of pore geometry should be valuable in designing or selecting the most effective oil-recovery process for heterogeneous limestones containing storage porosity. INTRODUCTION The primary and secondary recovery of oil from limestone (carbonate) reservoirs has been and will continue to be a major source of the world's supply of petroleum. Oil recovery from this type of formation is strongly influenced by the size, shape and arrangement of the pore spaces which hold the formation fluids. Therefore, it is expected that better knowledge of pore geometry would lead to design of more effective recovery processes. The concept of storage and flow channel pores has been presented as an important characteristic of limestone porosity.1 In this concept, storage pores are only containers for fluids, whereas flow channel pores are those pores which both contain and transmit the fluids. The work reported herein involves the concept of an additional pore characteristic called "porostriction". This characteristic can be visualized as the flow resistance between the storage and flow channel pores. This resistance is an important pore geometry characteristic because it determines how easily gas and oil can replenish fluids removed from the flow channel pores. Fatt and co-workers at the U. of California have studied the effect of storage porosity and a factor similar to porostriction on transient-flow response of models of porous media during a pressure build-up or drawdown test.55 These models consisted of porous sandstone cores with artificially simulated storage pores. It is the purpose of this paper to describe tests which permit the storage porosity and porostriction of actual reservoir core samples to be estimated. The tests used for this purpose involved not transient tests but, rather, an alternating-flow technique in which pressure and flow are varied according to a sine-wave relationship. This technique may offer advantages over transient methods. THE POROSTRICTION CONCEPT Before discussing alternating-flow (A.F.) tests in core samples, it is necessary to define the porostriction parameter. For this purpose attention is focused upon an elemental volume v representative of a porous medium having both flow channel pores and restricted-access storage pores. For such an elemental volume, it is considered that the storage porosity is combined to form a single pore and that the access into the storage porosity may be represented by a single capillary tube connecting the single pore to the flow channels of the elemental volume. Let the radius of the capillary tube be designated as Rc and its length as L,. When fluid flows into or out of the pore, the following equation (Poiseuille's law) relates flow rate q to pressure drop p, and the physical characteristics of the tube and the fluid of viscosity jx.

By Walter J. Karplus

By means of a finite difference expansion a fluid flow field in cylindrical coordinates with axial symmetry, is simulated by a network of electrical resistors. A series of DC analog computing units, connected to the network bounrlary corresponding to the free surface effectively adjust the network parameterr to satisfy simultaneously Laplace's equation and the borrndary conditions. The system automatica1ly "re1axes" to the correct boundary shape in a fraction of a second. A detailed application of the method to the water coning problem in oil wells is presented. INTRODUCTION The water coning problem in petroleum production is a prime example of a class of problems known as free surface problems. A free surface arises in fluid flow regions when the shape of the boundary of a single phase flow region is affected in some manner by the pressure or potential gradients within the region. When there exists an interface of two fluids of different densities, the shape of this interface will be determined by the relative velocities of the two fluids and will therefore comprise a free surface. Other free surface problems of interest to petroleum engineers include the gas coning, water injection, and gravity drainage probIcms. In each case it is desired to predict the shape of the boundary of the oil zone as a function of the flow rate, the well penetration and the physical parameters and dimensions of the system. Such problems do not lend themselves to the classical analytical treatments of boundary value problems because in this case the boundary of the potential field and the potential distribution within the field are interdependent and must be calculated simultaneously. A further difficulty arises from geometric considerations. Since the fluid flow is of the radial type in most cases, such effective methods as conformal transformation and the hodograph are not applicable. It is the purpose of this paper to present a method whereby free surface problems of this type may be solved rapidly to the desired order of accuracy by means of a combination of a network of electrical resistors and conventional analog computer units. The water coning problem is treated in some detail. Then the procedure for extending the technique to other free surface problems is indicated. FORMULATION Most commercial oil reservoirs consist of an oil-saturated sand, bounded above by an impermeable strata and below by a water-saturated sand. Gravimetric separation results in a horizontal oil-water interface. When an oil well is drilled so that it partially penetrates the oil zone, and production is initiated at a constant rate, the oil-water interface will assume a new shape as a result of the pressure gradients in the oil zone. It should be emphasized that this is a steady phenomenon. After an initial transient period the water phase will be completely static and only the oil will be in motion. The new steady-state shape of the interface, as shown in Fig. 1, depends upon the depth of penetration, d, the thickness of the oil zone, I, the production rate, q, the hydraulic permeability zone of the sand, k, the viscosity of the oil, µ, and the densities of the oil and water, p.,. and pw.. The problem is to determine this relationship. Darcy's law + = - (P f g p z).....(1) it.

Jan 1, 1957

By Morris Muskat, J. M. McDowell

In both production and refining operations of the oil industry many processes are controlled by the gas-liquid phase relationships of the hydrocarbon mixtures of interest. The quantitative behavior of fractionating or distillation columns and stills depends on the changes in the equilibrium distributions of the hydrocarbon components between the gas and liquid phases as the pressure or temperature conditions may vary along the length of the column or as the gross operating parameters and compositions are varied. In oil producing operations the nature and amount of gas phase developed within an underground reservoir as it is being depleted and its pressure declines also involve the basic equilibrium gas-liquid phase interactions. The influence of surface conditions of temperature and pressure and various separation procedures on the nature and amount of stock tank oil and natural gas recovered from a given well stream is likewise determined by the same basic phase equilibrium characteristics of the hydrocarbon systems. Such equilibrium separation of the gas and liquid phases in a well stream is of importance both in establishing optimum separator conditions for obtaining maximum stock tank oil yields and in the general problem of crude stabilization. The phenomena involved in these problems and their quantitative aspects can, of course, be established in each individual instance by appropriate laboratory experimentation. It has been found, however, that by associating with the individual hydrocarbon components functional characteristics de- scribing their individual gas - liquid phase equilibrium behavior, the phase properties of the composite mixtures can be predicted by mathematical computation. These characteristics have been termed "equilibrium constants", or "equilibrium ratios"," and are defined by the general equation: where K, is the equilibrium ratio for the jth component, y, is the rnol fraction of that component in the gas phase, and xj is the corresponding mol fraction in the coexisting liquid phase. Considering the Kj as representing available and measurable properties of the individual components of a mixture, and as defined by Eq. (l), it is readily shown that in a coexisting gas-liquid system the rnol fraction concentrations of the individual components in the gas phase, y,, and those in the liquid phase, x,, are given by the expressions: where the nj are the corresponding rnol fraction concentrations in the composite system, and ng is the mol fraction of the whole which is in the gas phase. The basic unknown in Eq. (2) is n,: which gives the gross separation between the gas and liquid phases. This may be determined in principle by imposing the requirement that the sums of the mol fractions in both the gas and liquid phases are unity, i.e. by the equations: The solution of the equivalent Eqs. (3) lying between 0 and 1, when inserted into Eqs. (2), provides a complete description of the compositions of the coexisting gas and liquid phases. While the Eqs. (3) are essentially equivalent to polynomial equations, such transformations are of little value when dealing with multicomponent sys-

Jan 1, 1949

By J. M. McDowell, Morris Muskat

In both production and refining operations of the oil industry many processes are controlled by the gas-liquid phase relationships of the hydrocarbon mixtures of interest. The quantitative behavior of fractionating or distillation columns and stills depends on the changes in the equilibrium distributions of the hydrocarbon components between the gas and liquid phases as the pressure or temperature conditions may vary along the length of the column or as the gross operating parameters and compositions are varied. In oil producing operations the nature and amount of gas phase developed within an underground reservoir as it is being depleted and its pressure declines also involve the basic equilibrium gas-liquid phase interactions. The influence of surface conditions of temperature and pressure and various separation procedures on the nature and amount of stock tank oil and natural gas recovered from a given well stream is likewise determined by the same basic phase equilibrium characteristics of the hydrocarbon systems. Such equilibrium separation of the gas and liquid phases in a well stream is of importance both in establishing optimum separator conditions for obtaining maximum stock tank oil yields and in the general problem of crude stabilization. The phenomena involved in these problems and their quantitative aspects can, of course, be established in each individual instance by appropriate laboratory experimentation. It has been found, however, that by associating with the individual hydrocarbon components functional characteristics de- scribing their individual gas - liquid phase equilibrium behavior, the phase properties of the composite mixtures can be predicted by mathematical computation. These characteristics have been termed "equilibrium constants", or "equilibrium ratios"," and are defined by the general equation: where K, is the equilibrium ratio for the jth component, y, is the rnol fraction of that component in the gas phase, and xj is the corresponding mol fraction in the coexisting liquid phase. Considering the Kj as representing available and measurable properties of the individual components of a mixture, and as defined by Eq. (l), it is readily shown that in a coexisting gas-liquid system the rnol fraction concentrations of the individual components in the gas phase, y,, and those in the liquid phase, x,, are given by the expressions: where the nj are the corresponding rnol fraction concentrations in the composite system, and ng is the mol fraction of the whole which is in the gas phase. The basic unknown in Eq. (2) is n,: which gives the gross separation between the gas and liquid phases. This may be determined in principle by imposing the requirement that the sums of the mol fractions in both the gas and liquid phases are unity, i.e. by the equations: The solution of the equivalent Eqs. (3) lying between 0 and 1, when inserted into Eqs. (2), provides a complete description of the compositions of the coexisting gas and liquid phases. While the Eqs. (3) are essentially equivalent to polynomial equations, such transformations are of little value when dealing with multicomponent sys-

Jan 1, 1949