Search Documents

By F. E. Crane, J. Downie

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

By C. C. Templeton, J. C. Rodgers

A key step in feed water treatment for generating wet steam for thermal oil recovery is the removal of calcium and magnesiunt hardness by cation-exchange series softening. Knowing the solubility of any scale forming salts in brines at elevated temperatures is necessary for fixing the level to which the feed water must be softened. Such calcium sulfate solubility data, previously not available above 392F, were determined by the authors in a flow equilibrium apparatus mud will be reported elsewhere. These data were used to develop a method for predicting the solubility of anhydrite in hot water or steam droplets for saturated steam pressures as high as 2,000 psig (637F). (The calcium sulfate solubility product is represented by a combination of two factors, one reflecting the effects of ionic strength and the other accounting for the effects of complex ion formation in either calcium-magnesium-rich or sulfate-rich brines.) The method is applied to a calcium-magnesium-rich brine If moderately high salinity from a pilot hot-water flood, I nd to several sulfate-rich, low-salinity feed waters and l lowdown (cooled steam droplets) samples from steam s ak operations. The predicted calcium hardness levels corresponding to the calcium sulfate solubilities agreed reasonably well with the results of laboratory solubility determinations run on the field samples. Further testing of the method is needed for brines of other composition classes. Existing field cation exchange softeners in the cases tested are performing adequately since all the samples were found to be undersaturated with respect to calcium sulfate at their operating temperatures. Introduction Prevention of scaling caused by precipitation of calcium sulfate (anhydrite) is of considerable concern in connection with thermal recovery processes using wet steam or hot water. To avoid anhydrite precipitation in a heated system, an engineer must keep the product of the calcium and sulfate concentrations in the water or steam droplets below the value of the solubility product of anhydrite for the temperature and brine composition in question. Usually it is most practical to keep the concentration product lower than the solubility product by keeping calcium low in the presence of high sulfate, or by keeping sulfate low in the presence of high calcium. This can be done by a choice of combinations of natural waters and water treatment processes (such as series cation exchange softening to remove calcium). Until recently, few anhydrite solubility data, particularly for solutions containing other salts, were available for temperatures above 392F (211 psig steam); Marshall, Slusher and Jones' studied the CaS0,-NaCI-H,O system up to 392F and surveyed the work of previous investigators. To model natural brines, one needs to study the solubility of anhydrite in aqueous solutions of sodium chloride, sodium sulfate, calcium chloride, magnesium chloride and their mixtures. Since steam pressures as high as 2,000 psig (637F) may be involved in thermal oil recovery projects, a solubility study was conducted between 482 and 617F.' Discussed in this paper is the application of these data to the prediction of anhydrite precipitation in some practical steam soak and hot-water injection projects. Any simple method for predicting the solubility of an inorganic compound over a wide range of temperatures and solution compositions must be based on some assumptions, and therefore must yield approximate results. On the one hand, natural brines contain too many ionic species for all to be included in a simple scheme; on the other hand, there is no adequate theoretical basis for the exact prediction of solubility in even simple solutions of mixed electrolytes. However, it is possible at a given temperature to base a reasonable prediction scheme on two phenomena:'-' the increase in solubility with increasing total concentrations of all ions (as measured by the ionic strength; see the Appendix), and an increase due to formation of cornplexes between calcium ions and sulfate ions and between sulfate and magnesium ions. Stiff and Davis' developed such a scheme for predicting the solubility of gypsum (CaSO, ¦ 2H,O) in brines at temperatures u~ to 212F. However, for the higher temperaturei of present interest, the stable solid phase is anhydrite (CaSO,). This study involves a two-part method for predicting anhydrite solubility products. First, one predicts a value K, for a given ionic strength I and a given temperature T corresponding to mo./msr, = 1 from data of the CaS0,-NaCI-H,O system. Second, one determines a group of factors F .= F(Ca) • F(Mg) • F(SO,), where the individual factors account for increases in the solubility product due to complex formation by high concentrations, respectively, of calcium, magnesium and sulfate. Combining the two parts, one obtains the solubility product in molalities as

Jan 1, 1969

By C. Perez-Rosales

A statistical method for determining simultaneously some of the basic geometrical characteristics of porous media such as porsity, specific surface* mean pore width, mean grain thickness and absolute permeability is presented. The proposed method is characterized by its simplicity and the minimum amount of laboratory equipment that it requires. The experimental procedure used to evaluate the various geometrical characteristics is based 'pan the 'surface analysis Of a a In view of this fact the app1icability Of the method is limited to homogeneous and isotropic materials. INTRODUCTION The behavior of a fluid flowing through a porous medium depends, among other things, upon the internal geometry of the medium; hence the importance of developing efficient methods to determine the geometrical properties of porous materials. This paper describes a statistical method to determine simultaneously some of the basic geometrical properties of porous materials such as porosity, specific surface, mean pore width, mean grain thickness and absolute permeability. The mathematical formulation has been developed that the data required to calculate the various geometrical characteristics of a sample can be easily measured by analyzing a section of the sample with an evenly spaced grid. THEORY A simple way of obtaining information about the internal geometry of a porous material is to throw a point at random over a cross-section of a sample. l-Through this procedure the porosity can be detemined.l,2 By considering that the material is homogeneous and isotropic, and that the point is dropped many times, the porosity 4 is given by where N is the total number of times the point is thrown and n is he number of times the point falls within pore areas. Another simple manner of analyzing the structural characteristics of a porous material is to superimpose an arbitrarily long line on a cross-section so that its length is evenly distributed over the surface area. Through this procedure the porosity can also be determined. If L is the total length of the line, and 1 is the sum of the lengths of the line segments within void spaces, the porosity is given by 4 = —.1.............. (2) The advantage of this type of analysis is that geometrical characteristics other than porosity can be obtained. Thus, if c represents the number of intersections between the line and the perimeter of pores, it can be shown3 that the specific surface, defined as he surface area of pores per unit bulk volume, is given by To give a proper description of the internal geometry of a porous medium, it is necessary to define quantities that characterize the notions of pore size and grain size. Unfortunately, because of the complexity of porous materials, it is difficult to give exact geometrical definitions of what is meant by the ConCePtS of "poresl' and "grains", especially when dealing with consolidated materials. Never theless, one often talks about the "size of pores" and the "size of grains" without defining accurately the terms. A possible solution to this problem follows. Assume that an arbitrarily long line is placed on a section of a sample so that the line is evenly distributed over the surface area (see Fig. 1). 1n a system like this, each line segment within a void space will be, by definition, the pore width in a given location and direction. Consequently, if the

Jan 1, 1970

By S. J. Prison, C. R. Smith

The effect of fluid injection to control water coning in oil and gas wells was investigated. Analytical and model techniques were employed. The factors investigated were the position and length of the completion interval, the point of fluid injection, the viscosity of the injected fluid and the relative thickness of the oil and water sections. The resulting influence of these factors on the net producing water-oil ratio was determined. Several important conclusions can be drawn from the study. In general, it was found that the net producing water-oil ratio can be reduced by fluid injection. The magnitude of this reduction depended on the factors listed above. An important practical consideration is that the injection fluid may be either oil or water. If the injected fluid is less dense than the connate water of the reservoir, the fluid will not be lost. This fact is reassuring when valuable oil is being injected. Efforts to suppress water production were more successful when the injection fluid was more viscous than the reservoir oil, or when a zone of reduced permeability existed in the vicinity of the point of fluid injection. Under test conditions, little benefit was derived through the use of impermeable barriers or cement "pancakes". INTRODUCTION The occurrence of water coning has been known for at least 60 years. In thin oil or gas pay sections, the presence of an oil-water or gas-water contact hinders production and often causes early abandonment of the afflicted well if a completion is even attempted. Even when relatively thick pay sections are found, the encroachment of water when a water drive is present will eventually pose serious water coning problems. This water is often corrosive, expensive to separate from the oil or gas and is costly to dispose of. The theory of water coning has been discussed by a number of authors. 1,2,3,4 Briefly, water coning to the producing interval in a well is due to pressure gradients resulting from the production of fluid from the reservoir. These pressure gradients will cause a water cone to rise toward the bottom of the producing interval if a water-oil or water-gas contact exists. The tendency of the water to cone is offset or partially offset by gravity forces since the water has a higher specific gravity than the oil. A balance then exists between two forces, gravitational forces arising from the difference in specific gravities of the oil and water, and the pressure gradients causing the flow of fluids to the wellbore. If the pressure gradient exceeds the gravitational force, water coning to the wellbore occurs and water production results. Through the years considerable thought has been given to the water coning problem. More than 50 U. S. Patents have been granted to inventors on the subject. A relatively complete literature review of the water coning problem has been made. 5 A number of these patents hold considerable promise for the solution or the partial solution of the water and/or the coning problem. Very little has been written describing field tests of techniques for the suppression of water coning. A notable exception is the paperby West. 6 He reports success in reducing gas coning by a combination of gravel packing and oil injection above the oil-producing interval. He also describes a comparable method to prevent water coning, but provides no field examples. This study experimentally and analytically verifies the benefits of oil injection as a means of partially or completely suppressing the water cone. While the gas coning problem was not treated, it is anticipated that results comparable to those obtained in water suppression could be obtained with reduced oil injection since the viscosity contrast between oil and gas exceeds that between oil and water. For the purpose of this study, the conventional potential flow theory was applied to the water-coning problem. The experimental verification centered on both a radial and a linear model. The model study permitted the investigation of complex flow configurations and the use of fluids of differing densities and viscosities. No analytic expressions are available to permit a solution of the problem as stated (see Fig. 1).

By B. H. Caudle, W. J. Bernard

Factors which influence the success or failure of a waterflood can seldom be determined in the laboratory. For this reason pilot waterfloods are initiated in a repreventative portion of the oil reservoir in question. For a pilot flood to predict quantitatively the recovery to be expected in a field-wide waterflood operation, the pilot area must behave as though it were confined (surrounded by similar areas). In this study, laboratory fluid-flow models were used to determine the simplest pilot pattern, for particular conditions of mobility ratio and initial gas saturation, that would behave as though it were confined. Pilot patterns studied ranged in complexity from a single inverted five-spot to a grouping of nine regular five-spots. Only the innermost producing well in each pattern was studied. Model results showed that the optimum number of wells in the pattern depends upon the oil-water mobility ratio and the expected oil-bank size. Unfavorable mobility ratios will, in general, require more wells in the pilot pattern than will favorable mobility ratios. Pilot patterns in reservoirs which contain a dispersed, flowable, free gas saturation will require fewer wells than for the under-saturated case. The single inverted five-spot pattern was found to be unsatisfactory for predicting behavior of fully developed waterfloods. In particular, it is possible that, in reservoirs which contain a flowable, dispersed gas phase, the oil bank will never be observed at the producers due to the large amounts of free gas which continue to be produced with the oil. INTRODUCTION One method which has been used to predict the performance of a waterflood is the pilot flood. The pilot waterflood is a flood which involves only a small cluster of the reservoir wells and is located in a small, representative portion of the reservoir. The object is that oil produced from the pilot can, in some way, be related to the oil recovery to be expected from a field-wide expansion of the waterflood. However, these field pilot waterfloods have often been unreliable in the prediction of oil recovery in a fully developed waterflood. This unreliability has also been demonstrated in several laboratory studies of pilot floods. Some of the investigators have shown that there are situations in which the pilot flood oil production is far too optimistic with respect to the oil recovery in the fully de- veloped flood. Others4-G have shown that the pilot results can also be pessimistic, especially if the pilot waterflood is initiated in an oil reservoir which has been depleted by primary recovery and is at very low pressure. The major reason for this unreliability of pilot water-floods is the migration of fluids into or out of the pilot area. By the well-known method of images, if straight lines can be drawn to represent vertical planes of symmetry in a porous medium which contains pressure sources and sinks (injectors and producers), then these lines are invariant streamlines, or lines across which there is no potential gradient, and therefore no flow. In an actual reservoir, these lines of symmetry can never be established exactly because of reservoir inhomogeneities and irregular reservoir boundaries. However, if the reservoir is relatively large and contains wells in repetitive patterns, these lines of symmetry are commonly assumed to exist for the pattern units sufficiently far removed from the reservoir boundary. Lines of symmetry for the five-spot injection pattern are shown in Fig. 1. Each five-spot unit in this figure can be considered confined with respect to flow across its boundary. In pilot floods this is not the case. The lines of symmetry for the pilot patterns investigated in this study are shown in Fig. 2. It is obvious that the fluid within these pilots is not confined and is therefore able to migrate into or out of the pilot area. Intuitively, one can see that, if more wells are added to the pilot, the innermost unit tends to behave more and more like the confined pattern. However, there is a practical limit to the number of wells which should be placed in the pilot. This limit is usually determined by economic factors. It was the purpose of this study to use laboratory fluid-flow models to determine which of the previously mentioned pilot patterns will force the innermost producing well to behave as it would in a fully developed waterflood. Since fluid migration is influenced by initial saturation conditions and the mobility ratio, these factors were included in the study. The ultimate objective of this study was to develop data which would allow the operator to choose a pilot pattern and operating conditions that will yield a production history which can be applied directly as an estimate of the performance of each production well of the fully developed waterflood. BASIS FOR THE STUDY The basic problem of field pilot floods is the migration of the reservoir fluids into or out of the pilot area. This problem has been the subject of previously reported model studies on pilot floods. These studies have been concerned mainly with the development of arbitrary "correction factors" to be applied to the simple, unconfined pilot systems such as the single five-spot. The correction factors were intended to adjust the production history of the un-

By T. S. Hutchinson, C. E. Kemp

The bottom water drive analysis presented by Muskat has been extended to include fields with wider well spacings. Curves are presented from which volumetric sweepout, water-oil ratio, and rate relations may be calculated. The effects of non-unity ratios on sweepout eficiency and productivity decline have been studied by an approximate method. Also, the effect of gravity on volumetric sweepout has been approximated for a particurlar geometry. INTRODUCTION The purpose of this paper is to report work done in applying the analysis of bottom water drive reservoirs given by Muskat1 to a larger class of reservoirs, to extend this analysis to include the effects of gravity and non-unity mobility ratios in an approximate manner, and to investigate the rate of decline of productive capacity. The bottom water drive type reservoir is one where the oil is underlain by a water zone of sufficient thickness that water movement is essentially vertical. This system is defined in more detail in the appendix and a one well element is shown diagrammatically in Fig. 1. Muskat obtained solutions to the equations describing the volumetric displacement efficiency of idealized bottom water drive systems by assuming that gravity effects due to differences in density betwecn oil and water were negligible and that the mobilities of oil and water were the same. If either of these assumptions is not made, the complexity of the system makes an analytical solution well nigh impossible at the present stage of mathematical development. Coles and Stade' developed an approximation to consider these factors. The authors have refined approximations, starting with the same basic assumption, which should give a somewhat quantitative picture of bottom water drive performance when mobilities are not equal and the effects of gravity are not negligible. The results of these approximate calculations can then be used to obtain the decline in well productivity as the field is depleted.

Jan 1, 1957

By Adrian C. Dorenfeld, Theodore Balberyszski, Strathmore R. B. Cooke

This paper reports results of studies of sulfidiza-tion of base-metal oxides and silicates with gaseous sulfur, hydrogen sulfide gas and pyrite and of their subsequent flotation with xanthate collectors. It has been established that a 100% conversion of the oxides of lead, zinc, copper and manganese into their respective sulfides is possible and is a function of reaction time and temperature. Sulfidization of chrysocolla (with malachite and azurite inclusions) resulted in a product containing covellite, chalcocite and digenite in proportions depending on conditions of sulfidization. Reacting cassiterite with gaseous sulfur, H2S or pyrite at various reaction temperatures resulted in its partial conversion to a tin sulfide. Flotation of the products of sulfidization indicated that their behavior is similar to that of natural sulfide minerals, the recovery being a function of the con-version of the oxides to their respective sulfides. The gradual depletion of high-grade sulfide mineral deposits has turned the attention of the mineral industry to the recovery of metals from the oxides and silicates. Anionic (fatty acid) and cationic collectors have in some cases made flotation of nonsulfides possible, although, in general, flotation of oxide ores is not as yet an economic process. Consequently, many attempts have been made to sulfidize the surface of the relevant oxide mineral and thus make it amenable to to In recent years, Bautista and sollenberger4 have investigated the conversion of a number of oxides to sulfides by reacting them with molten sulfur. This paper presents the results of .investigations carried out at the School of Mineral and Metallurgical Engineering at the University of Minnesota to study the flotation characteristics of some sulfidized minerals. It is divided into two parts: The first describes the sulfidization of a number of oxides and a silicate, and the second presents the flotation behavior of the products of sulfidization. SULFIDIZATION OF METAL OXIDES AND SILICATES Experimental Procedure: SULFIDIZATION WITH HYDROGEN SULFIDE GAS — Sulfidization with H2S gas was carried out in a horizontal tube furnace (Fig. 1). The temperature of the furnace was regulated by means of a Variac and determined with a chromel-alumel thermocouple. The constant temperature zone in the furnace was about 6 in. in length, and it was within this zone that the experiments were carried out. The gases employed were supplied from standard compressed gas cylinders, and passed through a glass-wool filtering tower before entering the Pyrex reaction tube. Gaseous products were bubbled through a water tower (used as a flow indicator) and excess H2S was burned off at the end of the outlet tube. Samples of reagents or minerals were placed in a 3-in. Alundum boat and inserted within the constant temperature zone of the furnace. Nitrogen gas (6 ppm O2 content) was passed through the tubes at room temperature to remove air in the tube, and continued to pass while the furnace was heating up. When the required temperature was reached, H2S was allowed to enter the tube and the supply of nitrogen was shut off. At the end of the reaction period, the reverse procedure was carried out, the samples being cooled to room temperature in a nitrogen atmosphere. The flow rate of H2S gas was maintained fairly constant in all experiments, although the rate was not determined. SULFIDIZATION WITH SULFUR POWDER - To study the sulfidization of minerals with sulfur (liquid or gaseous), samples of the mineral, together with a measured quantity of sulfur powder, were placed in a stainless steel (SS316) reaction vessel (Fig. 2). The reaction vessel was heated in a vertical resistance coil furnace, the temperature being regulated through a Variac. At the end of the experiment the vessel was cooled to room temperature before the sample was removed. Cooled products were leached with CS2 to remove any free sulfur. SULFIDIZATION OF TIN CONCENTRATES - Fifty gram samples of sulfidized tin ore were prepared by mixing -48 mesh Bolivian tin concentrate with -48 mesh Ottawa sand in various proportions and sulfi-dizing these samples with H2S in the horizontal tube

Jan 1, 1969

By S. C. Miller, R. D. Carter, H. G. Riley

A merhod based on theoretical considerations is presented for estimating srabilized performance of gas wells from data obtained from short-term flow tests. Calculated results are presented for II gas wells evaluated with four-to five-hour flow tests. For seven wells on which comparison war available, satisfactory agreement between measured and calculated performance was obfained. Although the final results are favorable, diflerences in flow capacity k,h obtained from flow tests conducted at different rates were observed. This inconsistency did not seriously affect the calculated results in the test wells. INTRODUCTION In order to negotiate the sale of gas from newly drilled gas wells and to plan otherwise for the orderly exploitation of gas reserves, it is necessary to predict the long-term withdrawal rates physically possible for each well. Since low permeability gas wells do not stabilize except after a lengthy period of flow, the operator is unable, without wasting large amounts of gas, to measure directly the stabilized performance prior to contracting for sale of the gas and connecting the well to a pipeline. This problem is not new, but is one which has recently received considerable attention because of the rising importance of natural gas in the economy and new emphasis on the development of gas reserves. Thus, a method is needed for estimating stabilized performance from information obtained from short-flow tests in which stabilized flow is not attained. Such short tests would be suitable for wells not connected to a pipeline since the tests require only small amounts of time, labor and production of gas. Other requirements of such a method are short calculation time and use of a model of the gas well and its reservoir drainage area which is as consistent as possible with present theoretical, experimental and field knowledge. The method presented herein was developed in an effort to meet these requirements. It has been named the "two-Row method". In addition to the section on the evelopment and description of the method, a description of field application of the method and a discussion of well test results are in- cluded. Finally, there is a section on an additional consideration prompted by some of the test results. The Appendix contains certain details of the derivation of the two-flow method and a sample calculation for Test Well 1. MATHEMATICAL DEVELOPMENT AND DESCRIPTION OF THE METHOD The method is based upon the following assumptions: (1) The system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial. (2) Flow obeys Darcy's law except very near the well where a quadratic flow law is assumed. (3) The system is also characterized by a laminar flow skin factor. The use of a quadratic or turbulent flow term is more fully justified in Refs. 1 and 2. When such a system is opened at the wellbore to constant rate flow from an initial condition of uniform pressure distribution, the radius or drainage quickly passes through the non-Darcy flow and "skin effect" region imrnediately surrounding the well. After the drainage radius or "transient" has passed beyond this region and before it reaches the external boundary, drawdown can be expressed by'." After the drainage radius reaches the external boundary, stabilized flow exists and drawdown is given by'," Eqs. 1 and 2 form the basis for the method. The validity of these equations will now be examined in the light of some well test data previously available. It has been found by actual measurements on gas wells that the relationship between drawdown and flow rate, that is the performance curve, is closely approximated by' Eq. 3 is a straight line on log-log graph paper. Eq. 2 and, for any lixed value of time, Eq. 1 plot as a curve which is slightly concave upward on log-log graph

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967

By R. V. Higgins, A. J. Leighton

A fast method is needed for calculating thoroughly the performance of two-phase flobr~ in reservoir rock with complex geotrietry. The authors present such a method and .slzow its accuracy by romparison with laboratory performrrnce of a five-spot pattern. In making the performance c.alculations, the octant of the five-spot pattern was divided into four channels. The computer tinle required for the culculation was one rllinrcte. This is rorighly about minute per channel. INTRODUCTION By using an electronic computer it is now possible to make reservoir performance calculations in which fewer assumptions and more significant variables are included than was possible with a desk calculator. These performance calculations include such variables as the properties of the reservoir rock and its contained fluids, properties of the rock-fluid system and even different well-spacing patterns. Usually, the more variables that are included the higher the computer cost—and time on the computer is expensive. As a result, a short procedure is needed so the many variables will be included and accurate performance calculations will result. If a reservoir is composed of alternate layers of permeable and impervioi~s rock, and if each conducting layer has a different relative permeability, the combined computer cost for calculating each layer by the present long computer programs (which includes different well-spacing patterns and other variables) would be expensive. The performances of complex spacing patterns in a field often are studied by a potentiometric model, but in these model studies the analogy is to the flow of only a single phase. This paper presents a procedure to calculate the performance of a five-spot pattern with two phases flowing using a potentiometric model as a guide for the dimensions of channels and other related data. The time required to calculate the complete performance is about one minute on the IBM 7090 or four minutes on the IBM 704. The procedure can be extended to calculate the performance of more complex well spacings or an entire field with little modification. The procedure has been tested on the performance of the laboratory water floods of a five-spot pattern reported by Douglas, et alThe range of viscosity ratios in the laboratory experiments is from 0.083 to 754, which is extensive. The accuracy of the performance calculations is excellent as shown by the plotted points in Fig. I. The recovery as influenced by the geometry of the five-spot, taking into account the deviation from true radial character, has been predicted by Aronofsky and Ramey." The authors used a potentiometric model to determine the sweep efficiency to breakthrough of the flooding agent. Dyes, Caudle and Erickson,' by means of a laboratory model, show the influence of mobility ratio on recovery after breakthrough of the injected material. Dyes, et ai, used the X-ray shadowgraph technique and miscible phases in their studies. Craig, Geffen and Morse from the analysis of the performance of their laboratory model correlate many variables, including average saturation of a flooding agent determined by permeability relationship. Craig, et al, used X-ray shadowgraphs to determine the sweep efficiencies in sandstone models. In this report the sweep efficiency is not needed, and the important effect of continuously changing saturation gradients before and after breakthrough

By F. Brons, W. C. Miller

Pressure information for use in material-balance calculations is obtained, where possible, from pressure build-up surveys in shut-in wells. Using proper extrapolation methods, static pressures are obtained which, by averaging, give the field static pressure. Often, particularly in fields with a large number of wells, only spot pressure readings are available. A graphical method is presented here which enables the determination of the corrections which must be applied to the spot readings to give an estinzate of the static pressure. Although the static pressures thus obtained are not so reliable as those from build-up surveys, they are more nearly correct than the spot readings themselves for use in material-balance calculations. INTRODUCTION The average static pressure in a bounded reservoir is determined from a number of individual-well pressure measurements. Where pressure build-up curves are available, extrapolation methods such as described in Ref. 3 may be used to determine the static pressure in the drainage area of each individual well. The possibility of using pressure build-up curves presupposes that each well being surveyed is closed in during a certain period and that, during this time, a pressure recorder is left suspended in the hole. For surveys in fields with a large number of wells, this is usually not feasible for economic and other reasons. Either a large number of pressure recorders would have to be used simultaneously or, using a limited number of recorders, the time interval needed to complete a field-wide survey would be too long. The data obtained for this incomplete coverage are insufficient to permit reliable averaging of the few calculated static pressures over the entire reservoir. In many such cases another type of pressure information is available, i.e., spot pressure readings. During this type of pressure survey all wells are closed-in simultaneously and, after a reasonable shut-in time (e.g., three days), pressure readings are taken by running the recorder in and out of each successive well. These surveys, which are less expensive, are often taken at regular intervals, say once or twice a year. In fact, such surveys are required by the state regulatory bodies in many portions of the united States. These spot readings are normally lower than the static (fully built-up) pressure due to the relatively short shut-in time. No direct means of extrapolating such a reading to the static pressure is available because no pressure build-up record is obtainable. Consequently, these data are inferior to pressure build-up information as a basis for estimating the stabilized reservoir pressure in the drainage area of a specific well. For the common case where abundant spot readings but only a few pressure build-up surveys are available, our studies suggest that it may be possible to design a program which will combine the two types of information and provide an estimate of the average reservoir pressure that is more reliable than that obtainable solely from either set of data. in time) gives a point marked p. The actual static pressure at the time the well was shut-in is marked by p.. The illustration shows that, even if a build-up curve is available and can be extrapolated, a correction

By J. M. McDowell, E. C. Doty, F. Morgan

Improvements in the X-ray method of measuring liquid saturation and saturation distribution are presented. Two identical direct current amplifiers have been added to measure continuously the intensities of the X-ray tube output and the beam which has passed through the core under test. The ratio of these intensities is automatically recorded on a modified Brown Electronik Recorder. Arrangements have been made to drive the core automatically and in synchronism with the recorder chart so as to give a continuous curve showing the variation of liquid saturation along the length of the core. By determining the ratio of the two X-ray beam intensities simultaneously the measurements are made essentially independent of small fluctuation in voltage and current supplied to the tube. This improved technique has been used for investigating the fluid distribution along the length of cores and core assemblies when one fluid is displaced by another. Typical results of these studies are shown. INTRODUCTION In a previous publication' a method was described for measuring the oil saturation in cores which was based on the absorption of X-rays by the fluid content of the core. While that technique is capable of yielding accurate results, the original apparatus was not especially convenient for determining the fluid distribution along the length of a core. This was primarily due to the slow response of the amplifier and galvanometer and the necessity of making a point by point survey of the material under investigation. It was, therefore, not only possible for sections of the core having a variable or unusual saturation distribution to escape observation entirely, but changes in output of the X-ray generator between readings on the unknown and the monitor introduced errors in the measurements that could be avoided only by further stabilization of the source or by taking a large number of readings. This paper describes a method and apparatus for recording continuously and automatically the intensity of the X-ray beam that is transmitted by a core. The variation of saturation with time at a given section of the core may thus be found, or by scanning the core by driving it in synchronism with the recorder the saturation distribution may readily be determined. As indicated in the previous paper, the X-ray method possesses an advantage over the weighing technique in relative permeability experimentation in that the core need not be removed from the holder for the purpose of making a saturation measurement. Furthermore, the actual distribution of fluid and the effects upon that distribution of inflow and outflow heads may be investigated. And errors due to gas bubbles or loss of liquid that are frequently sources of error in volumetric methods of saturation determination are avoided. PRINCIPLE OF THE METHOD The energy in an X-ray beam that has passed through an absorbing material is expressed by the well known formula: E = E0, e-cd (1) where E, is the energy in the incident beam, µ the mass absorption coefficient, c the density or concentration of absorbing material, and d is the length of path of the beam in the absorber. The coefficient µ varies with the wavelength of the radiation as well as with the absorbing material. The output of an X-ray tube varies in an approximately linear manner with the tube current, but to the second or higher power of the anode voltage. Moreover, the quality or wavelength distribution of the radiation varies with the voltage.

Jan 1, 1950

By R. Schuhmann

Distribution of material in the fine sizes of a comminution product generally is well represented by the empirical equation' y = 100 (x/k)a [1] in which y — cumulative percent finer, x = particle size, a -= distribution modulus, and k = size modulus. Charles3 found that the energy consumption in comminution is usefully expressed by another empirical relation, E = Ak(1-D) [2] in which E = energy input per unit volume of material, A = a constant, k = size modulus based on Eq. 1, and n = a constant; (1-n) is the slope of a plot of log E vs log k. Holmes3 has presented energy equations similar to Eq. 2. The constants a and n in Eqs. 1 and 2 have been shown to depend both on the nature of the material and on the comminuting device. Moreover, Charles showed that within experimental error a and n are a — n+1 = 0 [3] Combining Eqs. 2 and 3, E = A k-a [4] In the first sections of this article it is shown that the energy equation, Eq. 4, can be derived directly from the size distribution equation for the fine sizes, Eq. 1. The derivations are made without assuming any of the specific relationships between energy and particle size which have been common in previous literature. For comminution processes in which Eqs. 1 and 4 adequately represent the experimental data, the constant A in Eq. 4 is found to be a simple and useful inverse measure of grindability. That is, A is the energy consumption per unit volume of comminution product finer than unit size as determined from the straight line portion of the log-log plot of the size distribution. These considerations all lead to a unifying hypothesis of comminution mechanism from which both Eq. 1 and Eq. 4 can be derived. Finally, it is pointed out that this hypothesis raises serious questions as to the significance of the Rittinger hypothesis, the Kick hypothesis, and other theories in which energy numbers are systematically assigned to various size fractions of comminution products in order to calculate theoretical energy consumptions. Derivation of the Energy Equation from the Size Distribution Equation: For simplicity, consider the comminution of 100 volumes of a feed material of relatively uniform particle size. The comminution process may be considered as the summation of many individual and independent comminution events. The extent of comminution is most easily expressed as the number of comminution events, z. In the first derivation, the key assumption is that the characteristics of the comminution events in a given crushing or grinding process are substantially constant and do not vary with the progress of the comminution process. Accordingly the characteristics of an average comminution event may be defined. In one such event, a quantity of energy $E is applied to a single particle of size f and volume $v. The crushing of this particle produces fine particles with a size distribution similar to that given by Eq. 1: yw=100(x/ka)a0 [51 In this equation yo, a0, and k0 are used to characterize the product of an individual comminution event rather than the product of the comminution process as a whole. In using Eq. 5, we will not be concerned with values of x close to the feed size f and will therefore assume only that the equation is applicable to the finest sizes of the material. In 100 volumes of total product, the actual volume of product finer than x from a single comminution event, or dy, is given by The total volume of material below size x, resulting from z events, is then given by y = z(dy) =z(dv) (x/ka)a0 =" Eq. 7 reduces to Eq. 1 when we let a, = a and [8a] z =100/dv (ka/k )1 or k = ka (zdv/100)-1/a [8b] Eq. 8a shows that the distribution modulus of the comminution product is the same as for the product of an individual comminution event. Eq. 8b shows how the size modulus of the comminution product k varies with the extent of comminution as measured by the number of events, z, or as measured by the fraction of the feed actually subjected to com- minution1 zdv/100. The energy input to 100 volumes of total feed, or 100E, is the sum of the energy inputs for all the comminution events: 100E =z (dE) [9]

Jan 1, 1961

By F. J. Fayers, S. A. Hovanessian

The one-dimensional displacement equation for a homogeneous porous medium, including the effects of gravity and capillaty forces, has been solved by a numerical method. A finite-difference scheme is developed for obtaining saturation, pressure and fractional flow profiles in waterflood recovety problems. From the numerical examples given, it is concluded that the gravitational forces have a pronounced effect on the saturation profiles and the pressure distribution curves of the system. INTRODUCTION Within the past 20 years, a number of papers have appeared in the literature dealing with the quantitative treatment of waterflood recovery problems. In their celebrated paper, Buckley and Leverett1 described a method for calculating saturation profiles when the effects of capillary pressure and gravity are excluded. Terwilliger, et al,2 included the effect of gravity in their theoretical and experimental investigation of oil recovery problems and obtained close correlation between experiment and theory. welge3 described a simplified method for obtaining the average saturation and the oil recovery from an oil reservoir. The effect of gravity can be included in this calculation. More recently (1958) Douglas, et al,4 presented a method for calculating saturation profiles which includes the effect of capillary pressure. The authors start with the one-dimensional displacement equation (which is nonlinear in the derivative) and, by a change of variable, transform this equation to a semi-linear partial differential equation. This equation is then solved by a finite-difference method on a high-speed digital computer. Fayers and Sheldon5 obtained the solution of the one-dimensional displacement equation by directly replacing the differential equation by a finite-difference equation. The effects of capillary pressure and gravity were in- cluded in their solution. Here, since the movement of the foot of the saturation front is governed by a separate equation, the elapsed time in attaining a particular saturation profile cannot be obtained. The results of these two papers indicate that the inclusion of capillary forces eliminates the triple-valued Buckley-Leverett saturation profiles. In the present paper, the one-dimensional displacement equation for a homogeneous permeable medium, including the effects of capillary and gravity forces, is solved. The method of solution resembles that described in the paper by Douglas, et al, since it involves a similar transformation of variable. However, it extends their work in that the effect of gravity is included and, in addition, pressure profiles may be calculated. Also, the functions of saturation [kro(S), krw(S) and Pc (S) ] required are entered in tabular form rather than as low-order polynomials. This simplifies data preparation and permits the use of a greater range of function types. Unlike the method of Fayers and Sheldon, the corresponding elapsed time required for the development of each saturation profile is calculated and, also, saturation profiles may be calculated after breakthrough. THEORETICAL CONSIDERATIONS The dimensionless form of the one-dimensional displacement equation for a homogeneous porous medium as obtained from Darcy's law and material balance relations can be written as follows.

By Henry Kreisinger

IN evaluating coals as steaming fuel on the basis of chemical and physical properties of the coals, the only properties that can be used are those for which the methods of determination are well established and standardized, and those that affect all methods of burning the coals in the same way. The discussion is written from the standpoint of the coal consumer. The production of steam is the largest field for the use of coal. It includes stations for generating electric power, central heating plants, industrial power plants where steam is used for power generation, factory processes and heating. The field also includes steam locomotives. Small house-heating boilers are not included in this field because such small boilers generally burn specially prepared domestic fuel. The plants using the steam coal are equipped to burn the coal when fired by one of the three general methods: namely, hand, stoker, and pulverized coal. The hand firing is now restricted to small plants and to locomotives. Central stations and larger industrial plants use either stoker or pulverized-coal firing. Some locomotives also use stokers. Pulverized-coal firing has been tried in locomotives but so far its adapta-tion to this use has been handicapped by practical difficulties. The coals used for steaming purposes are generally of the bituminous and semibituminous classes, with volatile matter ranging from 15 to 35 per cent. Coals belonging to these classes are produced in many parts of the United States and form by far the largest part of the total coal production of the country. The value of any coal as fuel for steaming purposes depends on the following factors

Jan 1, 1935

By R. G. Grieves, G. Thodos

A method is presented for the accurate calculation of the cricondentherm and cricondenbar pressures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain six and quite possibly any number of components including paraffins, isoparaffins, olefins, acetylenes, naphthenes and aromatics. The approach is similur to that used for calculating critical pressures and cricondentherm and cricondenbar temperatures. The critical pressure, the normal boiling point and the approximate vapor pressure behavior of each component are all that art required. A stepwise culculation procedure is necessary for mixtures containing more than two components. From an analysis of 22 binary systems and 118 mixtures, the average deviation of calculated cricondentherm pressures from reported values is 2.4 per cent. For nine multicomponent mixtures the average deviation is 1 per cent. Considering 19 binary systems and 108 mixtures, the average deviation of calculated cricondenbar pressures from reported values is 1.7 per cent. For 15 multi-component mixtures, the average deviation is 2.2 per cent. INTRODUCTION A knowledge of the phase behavior in the critical region of multicomponent hydrocarbon mixtures is of value both in industrial processing operations and for the optimum operation of gas condensate reservoirs. Accurate methods of calculating the critical temperatures and critical pressures4,6,7,8,9,23 and cricondentherm and cricondenbar temperatures6p10 of multicomponent hydrocarbon mixtures are available in the literature. If the cricondentherm and cricondenbar pressures could be calculated with equal accuracy, the entire phase diagram of a multicomponent hydrocarbon mixture could be well approximated. The work of Etter and Kay6 is limited to systems containing the normal paraffins and has not been tested on systems containing a heavier component than heptane. In addition, the development of their multicomponent equations is based upon a limited number of mixtures containing from three to six components. Silverman and Thodos have considered systems containing both paraffinic and non-paraffinic hydrocarbons, but their correlation is limited to binary systems and is highly inaccurate for methane systems. Eilerts has done extensive work on the cricondenbar pressure; he has produced an excellent correlation for binary systems. However, his procedure for multicomponent mixtures is chiefly useful for highly complex mixtures requiring a knowledge of the vapor-liquid equilibrium behavior of the mixtures; he has not considered the cricondentherm pressure. The objective of this study was the development of an accurate and rapid method for the calculation of the cricondentherm and cricondenbar pressures of multicomponent mixtures containing all types of hydrocarbons and having a wide volatility range. The approach that was adopted is similar to that used by Grieves and Thodos for cricondentherm and cricondenbar temperaturesl0 and for critical pressures.9 However, methane systems had to be considered separately and a modified stepwise calculation procedure was utilized for the cricondertherm pressure. The correlations were developed in a manner similar to those for critical pressures and cricondentherm and cricondenbar temperatures. Based upon binary data reported in the literature it was observed that the ratios of cricondentherm and cricondenbar pressure to the pseudocritical pressure (molar average), pt/ppc and pp/ppc, respectively, in two-component systems depended upon the mole fraction of the low-boiling component and upon the diversity in properties of the two components. A dimensionless boiling-point parameter Tb,/Tb was chosen to represent the diversity in properties of the components. For a binary system, Ti,'is the molar average of the normal boiling points of the two components involved. For a multicomponent system it is either the molar average of the normal boiling points of all of the components involved or the molar average of the normal boiling point of the pure low-boiling component and of the atmospheric boiling point of the low-boiling-free mixture of

Jan 1, 1965

By D. L. Murphy

As the world's demand for iron ore increases and available sources are depleted, increasing attention will be centered on deposits in the Quebec and Labrador portions of the Labrador Trough. Improvements in the beneficiation of low grade iron ore stimulate exploration for deposits amenable to concentration. This study is centered on the Mt. Wright-Lake Carheil area, lying between Wabush Lake and Mt. Reed in Quebec. Described in detail are the general lithology, structural and economic geology, and mineralogy of the iron-formation. The iron ore province of Quebec and Baffin Island extends from the southern shore of Baffin Island southeasterly through the Labrador Trough to Matonipi Lake, approximately 185 miles northwest of Seven Islands. (See Fig. 1.) Future geologic field studies will very likely extend the limits of the province to include the ferruginous and associated rocks of the Lake Mistassini area. At present, only the deposits near Knob Lake and Lake Jeannine (Mt. Reed) have been brought into production although in the near future several other properties including those at Wabush Lake will commence operations. Within the next decade or so, as the world's demand for iron ore increases and available sources are depleted, this region will undoubtedly become a major mining district. Initial reports of iron-rich rocks in Labrador date back to the investigations undertaken by A. P. Low near the turn of the century. During the ensuing years the country was prospected chiefly for precious and nonferrous base metals because of its inaccessibility and the low tenor of the iron ore. In 1945 an evaluation of the direct shipping ores near Knob Lake indicated that sufficient reserves were present to justify the enormous capital expenditure required to bring the deposits into production. Initial open-pit mining began in 1954 following the construction of a townsite, a 360-mile long railway to Seven Islands, and an elaborate loading dock. Recent improvements in the beneficiation of low-grade iron ore have stimulated widespread exploration for deposits amenable to concentration, particularly along the southwest extension of the Trough. The Mt. Wright-Lake Carheil area, bounded by longitudes 66O45' and 67 30' west and latitudes 52 30' and 52 45' north, is a nearly rectangular tract of slightly less than 500 sq miles. It is in Saguenay County, 174 miles north-northwest of Seven Islands (see Fig. 1). Except for the periods of the spring breakup and the fall freeze up, the area is readily accessible by air.

Jan 1, 1962

By J. H Lightfoot, Stuart E. Buckley

Deep drilling has led to the development of numerous pools which yield products consisting predominantly of natural gas accompanied by high-gravity water-white or slightly straw-colored liquid. Yields of this "condensate" or "distillate," as it is commonly called, vary from traces to over a hundred barrels per million cubic feet of gas. The fact that this liquid exists in the reservoir as an integral part of the gas phase and is formed by condensation attending simultaneous reduction of pressure and temperature after its entry into the well was apparently first recognized about 1932.l The nature of the physical phenomena involved in this condensation was clarified for the petroleum industry in 1933 by a report by Sage, Schaafsma, and Lacey2 of the results of an experimental investigation of the behavior of mixtures of methane and propane. Search of the literature has revealed that the condensation of liquid from gaseous mixtures by isothermal pressure reduction is not a new discovery. Kuenen was led by theoretical considerations to expect such a phenomenon, and in 1892 reported3 the results of experimental verification of his deductions. This was followed by a considerable amount of work on the phase relations of mixtures by Kuenen 4,5,6 and by 0thers,7,8 mostly in European universities. Since 1930 investigators in this country have reported experimental data on the properties and behavior of a number of hydrocarbon systems."' It is now generally known that in the presence of natural gas the solubility in the gas phase of the heavier and normally liquid hydrocarbons increases rapidly with increasing pressure above 500 to 1000 Ib. per sq. in. Thus, at pressures of 2500 Ib. per sq. in. or higher a gas may contain appreciable quantities of hydrocarbons with normal boiling points as high as 6o0°F., which in turn may be condensed from the gas phase by the reduction of their solubility therein attending pressure reduc-tion. The combined effects of pressure and temperature on condensation and vaporiza-tioi of mixtures are most readily visualized by the study of phase diagrams. Since the construction, application, and interpreta-tion of such diagrams have been described fully in readily available literature, 18-20 this paper will not review the theoretical aspects of this subject. While the behavior of gas-distillate sys tems is now clearly understood, at least in a qualitative way, most of the published work deals with the results of laboratory studies of binary systems or of synthetic mixtures of hydrocarbons, and although several systematic investigations of the effects of separator pressure and temperature on the yields of liquid product from gas-distillate wells have been made21,12z it appears that in all cases the limitations have precluded the construction of even approximately complete phase diagrams for naturally occurring systems. There has 2

Jan 1, 1941

By V. V. Valleroy, A. J. Cornelius, B. T. Willman, G. W. Runberg, L. W. Powers

This paper reports the results of an investigation into the use of steam as a recovery agent. High oil recoveries by steam, as much as 100 per cent greater than by water flood, were demonstrated in laboratory linear cores. The principal mechanisms responsible for the additional oil have been identified as (1) thermal expansion of the oil, (2) viscosity reduction and (3) steam distillation. Although the first two also enhance recovery by hot water flood, they are effective with steam injection at lower volumes of produced fluids. Steam-zone residual oil saturation is independent of initial saturation. Recovery by steam is greater for lighter oils because they contain a greater fraction of steam-distillable components. Equations for the steam-zone advance in a radial homogeneous system are presented; these show the effect of reservoir parameters. The displacement in the hot liquid region moving ahead of the steam front can be estimated by a calculative procedure which is suggested. The method which takes account of the temperature gradient is applicable also to a regullar hot water flood. INTRODUCTION Considerable industry interest in the application of heat to petroleum reservoirs as a means of increasing oil recovery has been evident in the past few years. Three thermal methods which have attracted interest are (1) the injection of saturated steam, (2) the injection of hot water and (3) the injection of air into a formation to support combustion within the reservoir. Much has been published in recent years on underground burning, and there have been several papers on hot water injection. Although there are several early reports1,2 on steam injection, no recent experimental studies have been published on this recovery method. In our early thermal methods research, a marked improvement was found in oil recovery by steam injection compared with recovery by water injection in the same system. An investigation subsequently was undertaken to obtain quantitative information from which to evaluate steam injection as an oil recovery method. In addition, it was desired to investigate the recovery mechanism of hot water as compared with that of "cold" water and steam injection. An understanding of the mechanisms by which these processes operate would generally define favorable field conditions for their successful application. This paper reports the results of the laboratory study and suggests calculative approaches to estimate the field behavior of the steam-injection process. EXPERIMENTAL METHOD EQUIPMENT A schematic diagram of the equipment used for the displacement tests is shown in Fig. 1. Essentially, it consisted of a core with suitable means for recording temperatures, a feed system for injecting steam or cold or hot water, and a receiving train for condensing, cooling and measuring produced fluids. For the steam experiments, air-free saturated steam was piped from the boiler to the core face in insulated lines that were provided with steam traps and bleeds to remove any condensate formed in the lines due to heat losses. A 12-in. section of the line immediately ahead of

By J. H. Lightfoot, Stuart E. Buckley

Deep drilling has led to the development of numerous pools which yield products consisting predominantly of natural gas accompanied by high-gravity water-white or slightly straw-colored liquid. Yields of this "condensate" or "distillate," as it is commonly called, vary from traces to over a hundred barrels per million cubic feet of gas. The fact that this liquid exists in the reservoir as an integral part of the gas phase and is formed by condensation attending simultaneous reduction of pressure and temperature after its entry into the well was apparently first recognized about 1932.l The nature of the physical phenomena involved in this condensation was clarified for the petroleum industry in 1933 by a report by Sage, Schaafsma, and Lacey2 of the results of an experimental investigation of the behavior of mixtures of methane and propane. Search of the literature has revealed that the condensation of liquid from gaseous mixtures by isothermal pressure reduction is not a new discovery. Kuenen was led by theoretical considerations to expect such a phenomenon, and in 1892 reported3 the results of experimental verification of his deductions. This was followed by a considerable amount of work on the phase relations of mixtures by Kuenen 4,5,6 and by 0thers,7,8 mostly in European universities. Since 1930 investigators in this country have reported experimental data on the properties and behavior of a number of hydrocarbon systems."' It is now generally known that in the presence of natural gas the solubility in the gas phase of the heavier and normally liquid hydrocarbons increases rapidly with increasing pressure above 500 to 1000 Ib. per sq. in. Thus, at pressures of 2500 Ib. per sq. in. or higher a gas may contain appreciable quantities of hydrocarbons with normal boiling points as high as 6o0°F., which in turn may be condensed from the gas phase by the reduction of their solubility therein attending pressure reduc-tion. The combined effects of pressure and temperature on condensation and vaporiza-tioi of mixtures are most readily visualized by the study of phase diagrams. Since the construction, application, and interpreta-tion of such diagrams have been described fully in readily available literature, 18-20 this paper will not review the theoretical aspects of this subject. While the behavior of gas-distillate sys tems is now clearly understood, at least in a qualitative way, most of the published work deals with the results of laboratory studies of binary systems or of synthetic mixtures of hydrocarbons, and although several systematic investigations of the effects of separator pressure and temperature on the yields of liquid product from gas-distillate wells have been made21,12z it appears that in all cases the limitations have precluded the construction of even approximately complete phase diagrams for naturally occurring systems. There has 2

Jan 1, 1941