Search Documents

By Donald L. Katz

The compressing of gas into a partially depleted gas drive oil reservoir to bring the contents to a single phase miscible with gas is proposed as a process worthy of serious study. The compressed gas and vaporized oil would be recovered by cycling. This paper explores the Possibility of recovering oil from reservoirs of 6,000 ft or more by this method. INTRODUCTION Many new methods for increasing the recovery of oil from reservoirs are under consideration in the various research laboratories. One method which should be considered for recovering oil From partially depleted gas drive fields is to compress gas into the reservoir until the contents become a single phase miscible with more gas. Natural gas then may be cycled through the reservoir to recover the oil. This procedure ideally will eliminate the capillary effects of the liquid ill the porous media and should permit recoveries similar to those obtained in cycling operations, probably of the order of 80 per cent. This paper explores the information now available for the prediction of the performance of such a project and discusses the problems likely to be encountered in recovering oil by cycling. For many years dry natural gas has been passed through oil reservoirs to recover natural gasoline which would vaporize from the oil. Recently a process of recovering more of the high boiling hydrocarbons has been described by employing high pressures at which the heavier hydrocarbons have a higher volatility. Nevertheless. this reported method does not eliminate the differential nature of the vaporization process and capillary forces are still present in the reservoir. The process under discussion in this paper would minimize differential vaporization and would remove capillary forces which are responsible for low oil recoveries. Higher gas pressures are required to accomplish this removal of the miniscus, but it does not follow that more gas is required for oil recovery. The process is limited to reservoirs deep enough to permit the use of the required pressure — Probabably 6,000 it and deeper. DESCRIPTION OF THE PROCESS Briefly, the process is to compress natural gas into an oil reservoir which is partially depleted so that it has an oil saturation in the range of 40 to 60 per cent of the hydrocarbon pore space. The injected gas will pass through the channels formerly carrying gas to the wells to build up the pressure and dissolve in the oil phase. This injection will continue without production from the field until the reservoir reaches a single phase. If one were to select unique conditions such that the gas composition, the reservoir oil composition and the Per cent saturation would bring the mixture in the reservoir to the critical temperature and pressure simultaneously the phase behavior accompanying injection is shown by Fig. 1. The oil-gas mixture in the depleted reservoir has a critical temperature at C,. Addition of natural gas to the mixture will give a new critical temperature for each addition. Mixture No. 3 is shown to have its critical temperature at the reservoir temperature and reservoir pressure at P: the critical pressure. Under ideal behavior, the entire reservoir would be at iti critical temperature and pressure and the entire hydrocarbon content would be at a single phase. If the reservoir content were a binary system of methane and a heavy constituent. it would be miscible with more methane. However, for complex mixtures the critical locus does not confine the two phase region, and the dew point curves of Fig. 1 occur above the critical locus. If fluid which was injected into the rezervoir at the critical point were of the

Jan 1, 1952

By W. F. Hosford

A series of aluminum single crystals was subjected to axially symmetric flow by drawing through dies. The orientation dependence of the drawing stress, which indicates the resistance to this type of flow, compares reasonably with predictions of Taylor and of Bishop and Hill. Measured lattice rotations likewise are generally consistent with theory. The most successful analyses of the plastic behavior of polycrystalline metal and its relationship to single crystal behavior are those of Taylor1-3 and of Bishop and Hill.4-6 The purpose of this paper is to provide a critical test of some of the basic assumptions involved in these theories. Both treatments apply to fcc metals and assume the shape change of each grain in a polycrystal to be the same as that of the polycrystal as a whole. Thus for a randomly oriented polycrystal extended under uni-axial tension each grain is assumed to undergo axially symmetric flow. dey = dez = -1/2dex, deyZ = dezx = dexy = 0 [l] where x, y, and z are a set of orthogonal axes, x being the tensile axis. Although [011.] oriented grains in bcc and fcc metals may depart markedly from this ideal behavior,7,8 the assumption is probably quite reasonable for most grains. Under a uniaxial tensile stress, ox, the work dW per unit volume expended in producing an increment of axial strain dex is, dW=oxdex [21 Since this work is expended in slip, dw =STi dyi 131 where Ti and yi are the shear stress and shear strain on the active slip systems. In both analyses Ti is taken to be the same on all slip systems, and dependent only on the total prior shear strain on all systems, so dw = Tdy [41 where t is the shear stress required for slip and y is the sum of the shear strain yi on all slip systems. Combining 121 and [4] and defining M = oX/T The ratio M is a generalized Schmid factor. It relates the tensile stress required for axially symmetric flow, ox, to the basic shear stress, 7, and it also expresses the ratio of crystallographic strain, y, to the tensile strain, Ex It depends only on the orientation of the crystal. Taylor evaluated M for various orientations by finding the least value of dy/dEx which could satisfy the shape change. Bishop and Hill, on the other hand, found M by considering the relative work (1/T)dW/dEx expended by stress states which could activate enough slip systems to produce the shape changes. The results of both methods appear to be equivalent, although certain omissions were made in Taylor's original work. The orientation dependence of M, calculated by the second method,9 is shown in Fig. 1. Taylor used the average 2= 3.06 for all possible orientations to predict the tensile stress-strain curve (oxvs Ex) for a random polycrystal from the shear stress-strain curve (t vs y) of a single crystal. Changes in M from lattice rotation were neglected so points on the ox vs ex curve were simply taken at ox = 3.06t and e = y/3.06 from corresponding points on the t vs y curve. Recent extension of the analyses of Taylor and of Bishop and Hill to crystallographically textured metals has been used to relate plastic anisotropy to textures, and to predict the potential texture hardening that could result from texture control.9 Unfortunately there have been relatively few experimental tests of the underlying theory. Taylor found that the tensile stress-strain curve of polycrystalline aluminum did agree reasonably well with that predicted from a single-crystal T-y curve. Had the single crystal undergone easy glide, how-

Jan 1, 1965

By Philip D. Anderson, Ralph Hultgren

Heat contents of extremely pure iron were measured over the range 300"to 1433"K, using a diphenyl ether calorimeter. Results from three samples containing widely differing impurities agreed with one another and with Previously reported results at higher temperatures (1184 ° to 1809°K) by Olette and Ferrier. From these data a table of heat contents of iron has been Prepared which is believed to be more accurate than those previously available. CP values have been derived which are consistent with these heat contents, with the best experimental Cp data, and with the thermodynamic conditions of equilibrium between bcc and fcc iron. The opinion is expressed that low-temperature extrapolations of presently known Cp values for y-Fe are not reliable. It would seem that the thermodynamic properties of iron, which is by far our most important metal, should have been well established long ago. This is not the case, however. As discussed by Olette and Ferrier,' there is objection to every earlier heat content measurement. While the best true Cp determinations agree well at temperatures below the Curie point, there is wide divergence at higher temperatures. Several investigators have examined available data and attempted to calculate thermodynamic functions for iron that are consistent with the two equilibrium temperatures of bcc and fcc iron. The two primary efforts were those of Austin' in 1930 and Darken and smith3 in 1948. The recommended values differ considerably; integration of Austin's Cp values leads to a heat content at 1100°K which is 300 cal per g-atom higher than the value recommended by Darken and Smith. The compilation of Fisher4 in 1949 was based on the selection of Austin. In 1956, Weiss and Tauer obtained yet another set of values by resolving Cp values into components due to dilation, lattice vibration, electronic excitation, and magnetism. The magnitude of contribution by each component and its dependence on temperature was estimated by semiempirical means. Their values are not self-consistent in that they do not form smooth curves when plotted vs temperature; moreover, their calculations cannot easily be repeated since they do not give the values of some of the parameters used. Clearly, experimental values for iron were in a most unsatisfactory state. The presently described measurements of the heat content of iron were therefore undertaken, covering temperatures as high as 1433°K. While the work was in progress, the high-temperature (1193" to 1789"K) heat content measurements of Olette and Ferrier' were published; in the region of temperature overlap the present results agree with Olette and Ferrier well within experimental accuracy. It is believed that the present results, combined with those of Olette and Ferrier, provide for the first time a reliable set of heat contents of iron from room temperature to the melting point. From these data, guided by available Cp measurements, it proved possible to derive Cp, entropy, and free energy function values for @-iron from room temperature to the melting point, and for y- iron in the temperature region of its stability. Exhaustive study led to the conclusion that extrapolations of y-iron functions to room temperature were not reliable enough to be useful; hence these have been omitted from the tables. EXPERIMENTAL Materials. Three samples of extremely high purity iron from different sources, containing widely different impurities, were measured. The first of these, designated "Chipman", was kindly supplied by Dr. K. K. Kelley of the U. S. Bureau of Mines from a lot which was made many years ago by Dr. John Chipman of the Massachusetts Institute of Technology. The second, designated6'Teddington", was supplied by the National Physical Laboratory, Teddington, England, from a lot zone-refined by them. The third sample, designated "Irsid", was supplied by the Institute de Recherches de la

Jan 1, 1962

By Peter A. Allen

INTRODUCTION Lac Minerals Ltd., a significant North American producer of gold and other metals for over three decades, recently completed the financing for its third discovery in an eight year period. From 1975 to 1982, three important gold deposits in Quebec and Ontario: namely Bousquet, Doyon and Hemlo, were located and tested by company geologists. When current development plans are complete, Lac will have injected $400 million into these deposits. A discussion of these projects and their financing reveals some extraordinary challenges faced by top mining management. Responses to these challenges developed by Lac management may interest executives planning to proceed with new North American gold or base metal mines. MINE DE BOUSQUET In 1976, we explored a property containing a large low grade surface gold deposit. The geologial setting was exciting and we tested recent rethinking about the volcanigenic genesis of gold and base metal occurrences. These ideas helped pus to persist in drilling a nearby conductor and 600,000 ton, 0.22 ounce/ton orebody was outlined. With gold at $Cdn.150 per ounce the feasibility was approved and work started February 1978. Availability of a company custom mill 35 miles away reduced our immediate capital needs to $8.5 million. First of all, understanding the geological significance of results at Bousquet was critical. Mine and surface exploration geologists at Lac are continually and closely working together to better understand the nature of gold deposits in the Shield. This management emphasis on both pure and mine geology meant that a much larger ore reserve could confidently be predicted. (A proposal to quickly develop and ship the ore to a nearby company mine was discarded.) Instead of designing just for the initial plan of 500 tons per day we therefore included headframe, shaft and hoisting capacity for 1,700 tons and provided a compartment for future shaft deepening. These prebuilt additions were welcomed within two years. The type of collaboration whereby surface and mine geologists make study and field projects of each other's problems sets the foundation for this predictive ability. At Lac, the traditional responsibility for and supervision of the mine geologist by the mine manager is modified by this collaboration. This adds interest and scope to the geologist's work and lets him think beyond production requirements. Besides such management considerations, dedicated, professional geologists with integrity and lots of practical experience are needed to plan successful developments like Bousquet. Second, projecting future bullion revenue was critical and our experience selling production for future delivery gave us a way to ensure that loan repayment at term would not be jeopardized by falling metal prices. Output from a single mine is often reduced or curtailed for numerous reasons, so "selling forward" is done best and with least risk by companies like Lac with more than one or two producing locations for that metal. At first, a forward marketing program was considered unwise, speculative and unorthodox by analysts and other producers. Declining markets in the early 1980s prolonged gold project capital recoveries however, and critics were better able to appreciate the usefulness of such a tool. The notion that gold miners shouldn't sell forward is rooted in the belief that long-run gold prices will rise as modern social theory implementation in democracies inevitably erodes currency values. Selling against this trend goes against this industry belief or "ethic" and is at odds with the company shareholder who expects the gold price to rise not fall. Nonetheless, an efficient tool is available to the gold producer (not to most other metal producers) to quickly disengage from adverse price trends or to lock in spot prices instantaneously for two years or more of profits. Gold had recently retreated from U.S.$200 to $104 and now had reversed to $150. Lac had been forced to close several of its producers during the late 1960s because central bank-controlled

Jan 1, 1985

By T. M Morris

THE rate of flotation of solid particles determines the percentage recovery of these particles which can be obtained during a given time interval. It is an established fact that the recovery is greatest in the fine sand size range and falls off on either side of this size range. It is also a fact that the history of the development of concentrating devices has been marked by a desire to recover the small particles: tables were developed to treat material too fine to recover in jigs; flotation recovers material too fine to be recovered on tables. As our ores become more complex and require finer grinding, a greater proportion is in the slow-floating size range. It is therefore of economic importance that the effect of size of particle upon rate of flotation be studied. Many studies have been made along this line. The present investigation attempts to quantify the rate of flotation by means of a rate equation, and then compares the effect of size upon the rate of flotation. It was found that the flotation rate constant increases with increase in diameter of particle in the size range smaller than the optimum size particle. The flotation rate constant decreases with an increase in diameter of particle in the size range coarser than the optimum size particle. This suggested two different actions taking place, and these are discussed. Suggestions for further work in order to improve the flotation rate of both fine and coarse particles are made. The usual method of expressing rate is to plot cumulative percent recovery vs time. If we knew the equation of this kind of curve we could replot it as a straight line, and the slope of the line would be an index of the magnitude of the flotation rate. It is easier to compare straight lines than curves. The basis used for obtaining the equation of the recovery vs time curve is the law of mass action. This law states that the "rate of a chemical reaction is proportional to the 'active masses' of the reacting substances present at any given time." This, of course, implies a constant temperature. We can modify this for our use by using the term flotation reaction instead of chemical reaction and the term floatable particle instead of active masses. dC Thus----— kc [1] where C denotes concentration of particles, t denotes time, k denotes a proportionality constant. If the exponent of C is 1, then we have a first order rate reaction. If the exponent of C is 2, then we have a second order rate reaction. For a first order rate reaction, we obtain, after integration, where k is again a proportionality constant and has the dimension of the reciprocal of time. C, is the original concentration or weight or percentage of floatable particles. C, is the concentration or weight or percentage of floatable particles remaining in the flotation cell after time t. If all particles of the mineral to be floated were capable of floating we could take C, as 100 and C, would be 100 - R, where R is the cumulative percentage recovery of mineral at time t. But it is not usual for all of the particles to be capable of floating, because of surface oxidation, etc. Therefore, we must introduce a factor x which is the percentage of un-floatable mineral. Hence the first order equation is It has been proposed by various workers1-' that this equation describes the flotation rate. More recently N. Arbiter4 states that a second order rate equation fits the recovery vs time curve more closely than does a first order rate equation. Morris- as discussed the evidence presented for this statement. It has been concluded that the first

Jan 1, 1953

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 L. W. Pommier, F. B. Brien

The Rheological properties of pulps are non-Newtonian in character and more than one viscosity parameter is necessary to describe their behavior, therefore, the single term 'apparent viscosity' which has been used in the past may be misleading. The study reported in this paper uses a dimensional analysis based on the Ostwald-de Waele shear stress-rate of shear model for non-Newtonian liquids to develop velocity and resistance equations for falling spheres. An electronic monitoring apparatus was assembled for determining the velocity of metallic bodies as they fall through opaque liquids and suspensions. By means of this apparatus the behavior of the viscosity parameters of the Ostwald-de Waele model and the applicability of the proposed equations were examined for a number of non-Newtonian liquids. The validity of the equations for these liquids was confirmed. Preliminary data collected using solid-liquid suspensions established the validity of the form of the basic equations for these systems. The differential movement of solid particles through solid-liquid pulps and suspensions is very important LO many branches of engineering, such as minerals processing, soil mechanics, ceramic slip processing, oil well drilling-mud technology, materials-handling in pipelines, and many others. There is a dearth of fundamental knowledge on the flow resistance encountered by immersed bodies moving through solid-liquid suspensions. The resistance offered by Newtonian liquids to the movement of free settling spheres is well described for laminar and turbulent conditions by Stokes' and Newton's Laws, respectively. However, problems are encountered when attempts are made to extrapolate these simple relationships into solid liquid suspensions. The rheological properties of pulps are such that they are non-Newtonian in character;'" that is, the applied shearing stress and resulting rate-of-shear are related by non-linear functions, and thus, in the simplest case two viscosity parameters are necessary to describe their rheological properties. For this reason the single term 'apparent viscosity' which has been commonly used to explain the flow resistance of suspensions and non-Newtonian liquids may be misleading. It is realized that a rigorous analysis of the rnechanisms of momentum transport of liquid-solid systems is complex indeed. There are many variables that affect the rheological properties of the systems such as density, size, shape, and surface energy of the suspended solid and density and viscosity of suspending liquid, as well as the solid to liquid ratio. However, it was felt that a useful contribution to the knowledge of the behavior of suspensions could be made by examining the systems in terms of the basic fundamentals of nowNewtonian flow. This has led to the development, by means of dimensional analysis, of velocity and resistance equations, using two viscosity parameters to describe the velocity of penetrating bodies and the flow resistance encountered by these bodies, as they pass through non-Newtonian systems. THEORY In order to examine the effects of the rheological properties of solid-liquid suspensions on the movement of immersed bodies it is appropriate to start with non-Newtonian liquids and develop equations, analogous to those of Stokes and Newton, for laminar and turbulent flow respectively, from which quantitative information may be obtained. Further investigation would then be directed at the complex suspension systems, and the effects of the many variables on the velocity and resistance equations. Liquids are classified on a basis of viscosity which may be defined as the relationship between applied shear stress and resulting rate of shear?-'* Pure single phase liquids composed of small light molecules behave as Newtonian liquids, having a single viscosity parameter since the shearing stress and rate-of-shear are related by a linear function. Bingham plastics are similar but have a yield stress. Of particular importance in our investigation are dilatant (shear thickening) and pseudoplastic (shear thinning) non-Newtonian liquids, with and without yield values.

Jan 1, 1968

WELL known member of the Institute, Henry Krumb, survived the Japanese earthquake and has written a most interesting description of his personal experience to a friend in New York, an extended excerpt from which we publish, by permission, below. Not having been written for publication it has increased value as a human document in its vivid description of personal impressions and sensations in the midst of such a catastrophe. The earthquake was the most terrible experience I have ever passed through. We arrived at Miyanoshita, Friday, Aug. 31, about 6 p.m., having left Yokohama about 2:30 that afternoon by automobile. It was a most beautiful ride, although rather warm, but not unusually so, as it has been uncomfortably hot ever since we arrived in Japan. The first part of this automobile trip was along the old highway from Kyoto to Tokyo, which is lined with fine old pine trees. There are 52 villages along this highway and we saw real Japanese country life: naked children, both boys and girls, and men wearing nothing but gee-strings. The country on either side was intensively cultivated, the green rice fields and lotus ponds making a most beautiful sight. The occasional hills were all re-forested with young pines. In places the road runs close to the ocean beaches. The last part of the trip was up a steep mountain road, cut out of the side of the mountain. On one side of the road, in the gulch below, was a roaring river. The road crosses this river several times. Miyanoshita is beautifully located on the side of a mountain about 400 ft. above this river at an altitude of 1400 ft. The grounds of the hotel are laid out in Japanese style, small lakes with fish, waterfalls and dwarf trees making a very pretty sight. On the morning of Sept. 1, after walking around the garden and the hotel, we went to our room and I started to write up a few notes. Mrs. Krumb was reading a guide book. A few minutes before noon she came over to the table where I was writing and said that it was about lunch time. Just at that moment hell broke loose. It was the first earthquake I have ever experienced and I did not immediately realize what was happening. The building was shaking as a mine office might shake when heavy blasting was going on near by. Mrs. Krumb, having experienced many earthquakes, stepped to the open door between our room and the bathroom and called to me to get out from under the chande-lier. The plaster was commencing to fall and a large wardrobe fell over. I then knew what was happening and ran to the win dow. Before I reached it a large chunk of plaster from the ceiling hit me on the head and cut my forehead and nose. Mrs. Krumb ran for the other window just as the marble mantelpiece and the dresser fell. Without any hesitation, she jumped out of the win-dow and I followed suit out of the other. Luckily neither of us sprained an ankle-we had jumped 8 ft. We started to runaway from the building but had not gone more than 25 ft. when a second shock or rather a continuation of the first became so violent and the ground under our feet was trembling so that it was almost impossible to stand up without holding on to a tree. The ground seemed to have become liquid and rocked violently, it seemed to me in a vertical direction. In some later shocks the movement was more horizontal. About this time Mrs. Krumb saw the blood on my face and commenced to break down, but only for a moment: As soon as she saw it was not serious she was all right and during the following days stood the experiences and hardships as well as the best of them:

Jan 11, 1923

By T. M. Morris

THE rate of flotation of solid particles determines the percentage recovery of these particles which can be obtained during a given time interval. It is an established fact that the recovery is greatest in the fine sand size range and falls off on either side of this size range. It is also a fact that the history of the development of concentrating devices has been marked by a desire to recover the small particles: tables were developed to treat material too fine to recover in jigs; flotation recovers material too fine to be recovered on tables. As our ores become more complex and require finer grinding, a greater proportion is in the slow-floating size range. It is therefore of economic importance that the effect of size of particle upon rate of flotation be studied. Many studies have been made along this line. The present investigation attempts to quantify the rate of flotation by means of a rate equation, and then compares the effect of size upon the rate of flotation. It was found that the flotation rate constant increases with increase in diameter of particle in the size range smaller than the optimum size particle. The flotation rate constant decreases with an increase in diameter of particle in the size range coarser than the optimum size particle. This suggested two different actions taking place, and these are discussed. Suggestions for further work in order to improve the flotation rate of both fine and coarse particles are made. The usual method of expressing rate is to plot cumulative percent recovery vs time. If we knew the equation of this kind of curve we could replot it as a straight line, and the slope of the line would be an index of the magnitude of the flotation rate. It is easier to compare straight lines than curves. The basis used for obtaining the equation of the recovery vs time curve is the law of mass action. This law states that the "rate of a chemical reaction is proportional to the `active masses' of the reacting substances present at any given time." This, of course, implies a constant temperature. We can modify this for our use by using the term flotation reaction instead of chemical reaction and the term floatable particle instead of active masses. [ ] where C denotes concentration of particles, t denotes time, k denotes a proportionality constant. If the exponent of C is 1, then we have a first order rate reaction. If the exponent of C is 2, then we have a second order rate reaction. For a first order rate reaction, we obtain, after integration, [ ] where k is again a proportionality constant and has the dimension of the reciprocal of time. C, is the original concentration or weight or percentage of floatable particles. C, is the concentration or weight or percentage of floatable particles remaining in the flotation cell after time t. If all particles of the mineral to be floated were capable of floating we could take C, as 100 and C, would be 100 - R, where R is the cumulative percentage recovery of mineral at time t. But it is not usual for all of the particles to be capable of floating, because of surface oxidation, etc. Therefore, we must introduce a factor x which is the percentage of unfloatable mineral. Hence the first order equation is [ ] It has been proposed by various workers1-3 that this equation describes the flotation rate. More recently N. Arbiter4 states that a second order rate equation fits the recovery vs time curve more closely than does a first order rate equation. Morris' has discussed the evidence presented for this statement. It has been concluded that the first

Jan 1, 1952

By T. M. Morris

THE rate of flotation of solid particles determines the percentage recovery of these particles which can be obtained during a given time interval. It is an established fact that the recovery is greatest in the fine sand size range and falls off on either side of this size range. It is also a fact that the history of the development of concentrating devices has been marked by a desire to recover the small particles: tables were developed to treat material too fine to recover in jigs; flotation recovers material too fine to be recovered on tables. As our ores become more complex and require finer grinding, a greater proportion is in the slow-floating size range. It is therefore of economic importance that the effect of size of particle upon rate of flotation be studied. Many studies have been made along this line. The present investigation attempts to quantify the rate of flotation by means of a rate equation, and then compares the effect of size upon the rate of flotation. It was found that the flotation rate constant increases with increase in diameter of particle in the size range smaller than the optimum size particle. The flotation rate constant decreases with an increase in diameter of particle in the size range coarser than the optimum size particle. This suggested two different actions taking place, and these are discussed. Suggestions for further work in order to improve the flotation rate of both fine and coarse particles are made. The usual method of expressing rate is to plot cumulative percent recovery vs time. If we knew the equation of this kind of curve we could replot it as a straight line, and the slope of the line would be an index of the magnitude of the flotation rate. It is easier to compare straight lines than curves. The basis used for obtaining the equation of the recovery vs time curve is the law of mass action. This law states that the "rate of a chemical reaction is proportional to the 'active masses' of the reacting substances present at any given time." This, of course, implies a constant temperature. We can modify this for our use by using the term flotation reaction instead of chemical reaction and the term floatable particle instead of active masses. dC Thus----— kc [1] where C denotes concentration of particles, t denotes time, k denotes a proportionality constant. If the exponent of C is 1, then we have a first order rate reaction. If the exponent of C is 2, then we have a second order rate reaction. For a first order rate reaction, we obtain, after integration, where k is again a proportionality constant and has the dimension of the reciprocal of time. C, is the original concentration or weight or percentage of floatable particles. C, is the concentration or weight or percentage of floatable particles remaining in the flotation cell after time t. If all particles of the mineral to be floated were capable of floating we could take C, as 100 and C, would be 100 - R, where R is the cumulative percentage recovery of mineral at time t. But it is not usual for all of the particles to be capable of floating, because of surface oxidation, etc. Therefore, we must introduce a factor x which is the percentage of un-floatable mineral. Hence the first order equation is It has been proposed by various workers1-' that this equation describes the flotation rate. More recently N. Arbiter4 states that a second order rate equation fits the recovery vs time curve more closely than does a first order rate equation. Morris- as discussed the evidence presented for this statement. It has been concluded that the first

Jan 1, 1953

By R. A. Morse, P. M. Bridges, L. E. Wilsey, Howard N. Hall, P. L. Terwilliger

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By G. W. White

Over the past several years, the steady progress of gas lift from a mysterious art into a developing science has been traceable through the many fine papers that have appeared in Journal of Petroleum Technology. Consequently, this magazine has rightly become a standard for the serious students of gas lift all over the world. The following discussion is offered in a sincere attempt to accurately present certain topics covered in the subject paper in the light of acknowledged state of the art procedures. INSTALLATION DESIGN The first point to be questioned is the author's installation design, beginning with Eq. 6. The use of 1,200 psi in the spacing equation instead of the 1,275 psi noted as available seems, at first glance, to be a simple arithmetic error. Upon further examination, it appears evident that both this 75 psi plus the 100-psi safety factor must be subtracted from available kickoff pressure because of conditions imposed by the characteristics of the ASC valve (see Fig. B-1). Pressure-charged valve design is incorrectly penalized by spacing the first valve using only 1,100 psi. The use of such a high back-pressure makes it appear that a typical ASC installation would require no more valves than a comparable conventional valve installation. Using a value of 100-psi tree pressure, the well could be valved to 11,000 ft very conservatively with seven pressure-operated valves, but would still require 11 ASC valves. The use of 400-B/D spacing gradient for a 140-B/D rate is not an accepted practice using precharged valves. The author would have done well to note the advantages to be gained by reducing the 500-psi tree pressure, which is clearly explained by Beadle, Harlan and Brown.' In Step 1 of Part III of the continuous flow discussion, the author incorrectly applies the Poettmann and Carpenter curves.2 Item 2 is an incorrect statement. An excellent reference on this point is the paper by Baxen-dell and Thomas.' VALVE MECHANICS The next topic relates to the mechanics of precharged valves. Eq. 2 for the ASC valve can be used to describe the condition of equilibrium, wherein the valve is closed with opening pending. When this equation is used to determine relative values of casing, dome and tubing pressure with conditions approaching the transition value, some strange results are obtained. The following uses the author's symbols and values given on page 689. Fig. B-1 shows the condition that exists at the moment an ASC valve has been uncovered during the unloading process (Condition I). The obvious result is a locked condition. The basis for the curve is shown on the figure which includes 1,000 psig casing pressure to uncover the valve. If injection gas is not bled off the casing until the pressure is decreased by an amount exceeding the stabilizer discharging valve pressure, upward adjustment of the casing pressure is mandatory for the valve to open and unloading operations to continue. In Condition I the dome pressure exceeds the casing pressure by 20 psi, and the tubing pressure required to open the valve is 1,212 psig. This tubing pressure would not exist under the assumed conditions. In Condition I1 the casing pressure is increased to the dome pressure, and the tubing pressure required to open the valve is 1,102 psig. It is apparent that the casing pressure must be increased to a value exceeding the dome pressure before unloading operations could continue. In Condition I11 a maximum casing pressure of 100 psi (stabilizer charging valve set pressure) above the dome pressure is reached, and the tubing pressure required to open the valve is 549 psig. The author's statement after Eq. 5 that "the minimum tubing pressure required to open an ASC valve is equal to the casing pressure at depth less the transition pressure of the valve" is incorrect except for Condition I11 only. Otherwise, the equation is invalid.

Jan 1, 1965

By L. E. Wilsey, R. A. Morse, P. L. Terwilliger, Howard N. Hall, P. M. Bridges

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By H. A. Koch, R. L. Slobod

This paper discusses a new oil recovery process called the "miscible slug process." This process involves the injection of propane or LPG into the reservoir prior to gas injection. The operating conditions are maintained so that the slug is miscible with both the reservoir oil and the injected gas phase. Thus, a miscible phase displacement is achieved, and high recoveries are obtained from the area of the reservoir contacted. The most striking discovery of the laboratory study which included displacements from cores up to 123 ft in length was that relatively small slugs of LPG are effective over reservoir distances in the miscible displacement of oil. This makes possible the commercial application of the process. The important factors controlling the size of the slug are: (I) reservoir length, (2) reservoir fluid composition, and (3) reservoir pressure at the displacement front. Of lesser importance are the effects of injection rate and porous medium type. Questions which have not been completely answered deal with the effects of gravity and reservoir inhomogeneities on the process. INTRODUCTION One of the major problems facing the oil industry today is recovery of the large fraction of the discovered oil which will remain unrecovered unless some new processes are found. Water flooding has been one successful method used to increase recovery over natural depletion. But even successful water floods may leave 25 to 40 per cent of the pore space filled with residual oil. Recently a process which involves the use of thermal energy has been reported1,'*' which may prove useful in increasing recovery from shallow reservoirs containing highly viscous crudes. This paper deals with a new oil recovery process involving use of miscible phases to displace the oil from the reservoir. Miscible phase displacement oil has been an intriguing idea because the elimination of capillary effects in the reservoir leads to 100 per cent recovery in the areas contacted by the miscible displacing phase. The big deterrent in the past to the practice of miscible phase displacement has been the high cost involved in using large volumes of solvent in the process. The use of high pressure gas to achieve a miscible phase displacement of reservoir oil has been reported4,5 as one economically feasible method. However, this process calls for the maintenance of pressures above 3,000 psi at the dis-placement front. The process described in this paper is not restricted to high reservoir pressures; the only qualification being that the slug of material which is injected into the reservoir must be miscible with both the displaced oil 2nd the displacing phase. Two such cases are the use of an alcohol slug followed by water to achieve a mis-cable water flood, and the use of a low boiling hydrocarbon slug such as propane followed by gas to achieve a miscible gas drive. This paper will deal specifically with the propane dug process. This process will obtain a miscible phase displacement at pressures as low as about 1,200 psi at a reservoir temperature 01 150°F. The process is very simple. The propane displaces oil completely at one front because it is miscible with the oil, and gas displaces the propane completely at a second or trailing front because it is miscible with the propane. A small band of propane should thus remain wedged between the gas and oil phases. This idea has been called the "miscible slug process."6 Unfortunately, the propane slug will become diluted as it flows through the reservoir. This results from a mixing action since the fluids flow through a complex reservoir porous medium rather than simple capillaries. Fig. 1 diagrams a miscible slug displacement. The effects of some of

Jan 1, 1958

By M. Prats

The effect of a sand-filled vertical fracture of limited radial extent and finite capacity (fracture capacity is the product of the permeability and width of the fracture) on the flow behavior of a cylindrical reservoir producing an incompressible fluid through a centrally located well has been investigated mathematically. The shape of the lines of equal pressure near the fracture is essentially independent of the size of the reservoir, provided that the field radius is of the order of the fracture length or larger. For reasonable values of production rates and of fluid, reservoir and fracture properties, the total pressure drop between the end of the fracture and the well is gelrerally negligible compared with the pressure drop in the reservoir. With regard to production response, the effect of vertical fractures can be represented by the production response of an equivalent or effective well radius. For a high-capacity fracture, the effective well radius is a quarter of the total fracture length, decreasing with the fracture capacity. When invasion effects are simulated by decreasing the width of the damaged zone with distance from the well, the effect of formation damage around a fracture on the production response is not so seriolls as indicated by the literature. This suggests that frac fluids with a conventional filter-loss response are better than high-spurt-loss frac fluids, provided the effective permeability of the damaged zone is the same. INTRODUCTION This paper considers the effect of the fracture capacity, as well as the formation damage which can result from fracture treatments, on the productive capacity of vertically fractured wells. Other publications, notably those of van Poollen,l,2 consider these same effects. In addition to providing more general results for vertical fractures than are available from the liter- ature, the present paper gives the equivalent well radius of fractures having different lengths and capacities and, also, includes pressure distributions in and around the fractures. The effect of a damaged zone around a fracture on the production response was not found to be so great as that reported by van Poollen. 2 This difference probably stems from the fact that we consider a damaged zone which is widest (but is still small) near the well and thins out toward the extremities of the fracture, whereas van Poollen considers a damaged zone having a uniform width for the entire fracture length. Simple, but adequate, equations which describe the effect of these variables on production response are presented (in Appendixes A and B). Thus, results can easily be extended to values of the variables not specifically considered here. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The liquid is incompressible. It is assumed that, after the well is fractured, the cylindrical outer "boundary" is at a uniform potential, provided that it is not too near the fracture. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below the producing layer (reservoir).3 Fig. 1 indicates the general three-dimensional geometry of the fractured reservoir described above. If gravity effects are neglected, the flaw behavior in the reservoir will be independent of the vertical position ill the oil sand. This means that the flow behavior in the fractured reservoir can be described by the two-dimensional flow behavior shown in the horizontal cross section of the reservoir, Fig. 2. PARAMETERS The effect of a fracture on the pressure distribution in the fracture itself and in an otherwise undamaged (no skin) reservoir indicates that the pressure distribution is a function of three parameters.

By R. D. Begley, A. W. Khair

Introduction The technique being investigated here uses a physical mine model which consists of a plexiglass room-and-pillar model 11.4 x 11.4 cm (4.5 x 4.5 in.). It is placed at the bottom of a larger plexiglass box 25.4 x 30.5cm (10 x 12 in.). The box is then filled with a model material representing geologic overburden. Sand has been selected because of its deformability under gravitational forces. The mine model is capable of creating a cavity and induces movement on the surface of the model material. This surface disturbance can be measured using a photographic technique called holographic interferometry. This paper discusses results of an investigation of the possibility of generating subsidence criteria using physical mine models and laser holographic interferometry. Problem Statement Subsidence appears to be widespread in two states - Pennsylvania and Illinois - although serious subsidence is found in localized areas of numerous other states (HRB Singer Inc., 1977). The US Bureau of Mines estimates that more than 1618 km2 (400,000 acres) of urban areas are threatened by subsidence. Pennsylvania and West Virginia contain 57% of these threatened areas, with the remainder being spread out over 16 states (Johnson and Miller, 1979). Events of subsidence have been reported in urban areas, including Pennsylvania (Gray et al., 1977) and Illinois (Mahar and Marino, 1981). One location has observed three episodes of subsidence beginning in 1950 (Bruhn et al., 1981). These subsidence events all occurred over abandoned room-and-pillar mines. Some mines have been abandoned more than 50 years (Bauer and Hunt,1981) . Available research does not offer much on the problem related to time dependent subsidence, especially those associated with the failure of remnant coal pillars. The physical modeling technique developed at West Virginia University is capable of measuring both instantaneous and time delayed subsidence. Laser Holographic Interferometry Laser holographic interferometry, or holometry, uses the interference properties of laser light to measure small movements of objects. The object requires only optical access and no surface instrumentation. This technique is closely related to photoelasticity. Both techniques rely on the uniqueness of the light reflected off or through an object in any given state of stress. The entire measurement process is done remotely without any surface instrumentation and without using a camera. The interferogram is constructed and visual observation will yield fringes on all areas that yielded or moved. The displacement represented by one fringe was determined by Wilson et al. (1977) and is given by the relationship: [d ='1(1) 2 cos 8 where d = displacement of one fringe 6 = incident angle of laser light A = wavelength of light (632.8 x 10-9 meters for Helium-Neon lasers)] In summary, holometry is capable of detecting very small movements of objects under study. It is a very sensitive process that requires laser light and vibration free experiments. Experimental Program Lab Facilities The primary components of the facilities were the laser and the optical table. The light source was 50 milli-watt, Helium-Neon laser. The optical table was a granite slab weighing nearly one ton. All experiments were performed on the optical table. It consists of a large table sitting on air cells similar to car innertubes. Six air cells were placed on top of two layers of wood and Styrofoam. Figure 1 is a schematic diagram of the vibration free table. Even with these precautions, the experiments still had to be performed at night to reduce potential traffic related vibrations. Room-and-pillar mine models consisting of three different pillar sizes were constructed out of plexiglass. The models were designed to represent three different extraction ratios. Figure 2 shows the models. The extraction ratios were determined by the tributary loading concept (Peng, 1978). These models were the primary component of a seam extraction simulator. A schematic of this simulation device is shown in Fig. 3.

Jan 1, 1986

By J. S. Browning

The U.S. Bureau of Mines has been experimenting with flotation processes to separate the spodumene-beryl ores mined at Kings Mountain, N.C. The success to date as well as the present status of the process is discussed. Further experimentation is underway and studies will be made on the milling costs involved in the flotation process. The pegmatites of the Kings Mountain-Lincoln ton, N.C., area constitute the largest known domestic reserve of beryl and spodumene. The reserve is estimated to contain 90 million tons of pegmatic material with 1,280,000 tons of recoverable Li,O(lithia) as spodurnene.' The pegmatites also contain 0.4 to 0.5 pct beryl disseminated throughout the orebodies. The pegmatites may contain a potential reserve of 240,-000 tons beryl, equivalent to approximately 34,000 tons of BeO. Different flotation methods for separating beryl from feldspar or quartz, or both, have been developed by various investigators.2-7 On the other hand, published information is limited on concentration of spodumene-beryl ores. As the response of beryl and spodumene to flotation is essentially the same, a successful separation of the two minerals depends upon use of a selective depressant or a selective collector for one of the minerals. Several years ago, the U.S. Bureau of Mines, recognizing the need for developing a large domestic supply of strategic beryl, undertook studies to develop improved methods of recovering this mineral. Because of the low beryl content of the pegmatites, economic production of the beryl entails its recovery as a byproduct of spodumene flotation. Thus, the development of procedures that would provide for maximum recoveries of the spodumene, as well as the beryl, was a primary objective of the studies. Tests were made on run-of-mine ore and spodumene flotation tailing from the Foote Mineral Co.'s spodumene concentrator at Kings Mountain (Fig. 1). Pet-rographic analyses of the ore and flotation tailing are given in Table I. Detailed analyses and examinations of the ore and tailings revealed that the beryllium was present in the form of a low-alkali beryl. The crystals were small, rarely larger than 14-mesh diam, and were clear and colorless. About 10 pct of the Li20 content of the ore was present in the mica and feldspar components of the pegmatite. In some weathered ores the associated clays also contained small amounts of lithia. This is typical of pegmatites in the Kings Mountain area. LABORATORY TESTS In the past, a vigorous chemical treatment of the spodumene surfaces with acid or caustic soda has been considered essential for satisfactory selective flotation of the spodumene from the other minerals in pegmatites. Many tests were made during this investigation to determine if chemical treatment of the spodumene could be eliminated in favor of a reagent combination that would depress the other minerals while selectively floating the spodumene. A simple reagent combination and process was developed that gave better spodumene recovery and grade of concentrate than the more elaborate acid or caustic treatment methods previously used. The finely ground ore pulps were conditioned with either an ammonium, alkali, or alkaline earth lignin sulfonate and sodium fluoride, and the spodumene was floated and cleaned using oleic acid as the collector. The spodumene tailing was deslimed, conditioned with sulfuric acid and coco amine acetate, and floated to reject mica. The mica tailing was thickened and conditioned with hydrofluoric acid, then washed to remove the acid. The washed pulp was conditioned with sodium hydroxide for pH control, and the beryl was floated and cleaned with oleic acid as the collector. Table II summarizes the results of the laboratory test work. About 77 pct of the spodumene and 75 pct of the beryl were recovered from the pegmatitic material. The spodumene concentrate assayed 6.1 pct Li20 and 0.01 pct BeO, and the beryl product assayed 3.0 pct Li20 and 1.57 pct BeO. Similar results were ob-

Jan 1, 1961

By A. N. Crownover, H. T. Kennedy, E. P. Miesch, C. H. Bowman

The paper presents correlations of (l) molar volume of gaseous hydrocarbon mixtures with pressure, temperature, composition and properties of the C,-plus fraction; (2) shrinkage of oils during flash and differential liberation of gas, including the calculation of formation volume factor under various conditions; and (3) bubble-point pressure with temperature, composition and characteristics of C, plus. The data on which the correlations are based comprise 1,615 measurements on 900 hydrocarbon systems, including numerous systems containing nitrogen, hydrogen sul-fide and carbon dioxide. In each correlation, the number of data points covered and the accuracy is substantially greater than in previously available work. Thus, the equation yielding molar volumes of gases has an average deviation of 2.04 per cent, applied to mixtures having temperatures up to 313F and pressures up to 9,800 psis, compared to 2.37 per cent for the Benedict-Webb-Rubin equation applied to the same data, and 4.53 per cent for the method based on the law of corresponding states. he equations resented are all explicit in the dependent variable, and require no iteration on the digital cow puter. INTRODUCTION The ease and accuracy of determining the composition of hydrocarbon mixtures, compared to the difficulty of measuring their properties under reservoir conditions, makes it desirable to utilize composition as the key to physical behavior to the greatest possible extent. As a result, there are available correlations between composition, or easily measured characteristics dependent on composition, and practically every important engineering property of reservoir fluids. The task confronting us is one of finding more exact relationships between important variables rather than extending correlations to new properties. This paper describes new correlations of molar volumes of gases, formation volume factors, and bubble-point pressures with composition, temperature and where possible, pressure. Each correlation is obtained by employing a sufficiently large amount of data so the calculated properties are probably as least as accurate as the measurement on which they are based. MOLAR VOLUME OF GASES Although many equations of state have been proposed for pure gases, only a few methods are applicable to hydrocarbons at conditions comparable to those in petroleum reservoirs. Still fewer arc useful in describing the behavior of mixtures, with which the petroleum engineer is largely concerned. The correlation presented here involves procedures similar to those of Alani and Kennedy The van der Waals equation is modified to make a and b functions of temperature instead of being constants for each material.This is a cubic equation, which may have either one or three real roots. The lowest root corresponds to liquid volume, while the highest applies to gas. When the procedure was applied to 703 pressure-temperature points of 164 gases of known composition and volume, the average deviation was 12.08 per cent, and the standard deviation 8.15. The above calculations were made using the constants derived for liquids by Alani and Kennedy and the mixture equation developed by them. A closer approximation to measured molar volumes is obtained by employing different sets of constants for different areas on the pressure-temperature chart, and by changing the relationship between a,,, and b,,, for mixtures and the a, and b, for individual hydrocarbons. The various sets are shown in Table 1, and the areas for which they are recommended are plotted in Figs. 1 through 7. In these figures, the area designated by zero is in the critical region for pure materials, and the values obtained for them may be unreliable. For normal mixtures, however, the area may be added to any of the bounding areas without substantial error. The new relationships employed are

Jan 1, 1966

By E. A. Gulbransen, K. F. Andrew

Weight loss measurements were made using a sensitive microbalance operating in a high vacuum system. The Langmuir equation was used to Calculate the vapor pressures of the several metals and alloys. Aluminum lowered the vapor pressure of i9,on in a 5.4 Al-94.6 Fe alloy, and of iron and chromium in a 4.8 Al-21.5 Cr-73.7 Fe alloy. The influence of oxide films formed on the alloys on the effective vapor pressures of the alloys was also studied. These studies were used to aid in the inte9,pretation of the high temperature oxidation of alloys based on iron, chromium, and aluminum. RECENT studies on the oxidation of chromium1 and the heat resistant alloy 5 Al-22 Cr-73 Fe showed transitions in the rate of oxidation at 900o and 1050oC respectively. The transition for chromium occurred at a temperature where the rate of evaporation of chromium from oxide-free surfaces equalled the rate of chromium atoms reacting with oxygen. This paper presents new vapor pressure studies on iron, chromium, and the specially prepared alloys 5.4 Al-94.6 Fe, 21.9 Cr-78.1 Fe, and 4.8 Al-21.5 Cr-73.7 Fe. The ternary alloy is similar in basic composition to the patented heat-resistant commercial alloy known as Kanthal. The binary alloys were studied to show the individual effects of aluminum on the vapor pressure of iron and of iron on the vapor pressure of chromium. The purpose of the proposed studies was to test the influence of metal volatility on the oxidation behavior of heat-resistant alloys based on iron, chromium, and aluminum. Since oxide films formed on the metal may limit metal transfer to the surface, the influence of oxide films on metal volatility was also studied. LITERATURE A) Vapor Pressure. The vapor pressure of iron has been studied by Jones, Langmuir, and Mackay,' Dornte and Nrton. Edwards. Johnston. and Ditmar. and McCabe, Hudson, and Paton. Stull and SinkeT have calculated a heat of sublimation at 298° K of 99.83 kcalper g atom. The vapor pressure of chromium has been reported by Bauer and Brunner,' Speiser, Johnston, and Blakburn, Gulbransen and Andrew,'' Vintaikin,o McCabe, Hudson, and Paton, and Kubaschewski and Heymer.12 Good agreement was obtained except for the early work of Bauer and Brunner.' Stull and Sinke7 calculate a heat of sublimation at 298°K of 95.0 kcal per g atom. Vapor pressure measurements on aluminum have been made by Brewer and Searcy,13 Bauer and Brunner,aand Farkas.I4 Stull and Sinke,7 giving Brewer and Seary's' data the most weight, derive a heat of sublimation of 77.5 kcal per g atom at 298 oK. B) Method. Two methods15 are used for measuring the vapor pressure of metals: 1) The Langhuir free evaporation method, and 2) the Knudsen effusion method, In the Langmuir method, used in the present work, the vapor pressure of the metal P is given by the equation: Here M is the molecular weight of the vapor species. T is the absolute temperature. (dw/dt) is the rate of sublimation in g per sq cm per sec and a is the condensation coefficient. a is a measure of the efficiency of condensation of molecules striking the metal surface from the vapor. If condensation results from each collision a is unity. This coefficient has been found to be unity,a, le for metals. These results establish the general validity of the Langmuir free evaporation method as applied to metals. EXPERIMENTAL A) Vacuum Microbalance Method. A modification of the Langmuir free evaporation method was used. Strip specimens were suspended from a sensitive microbalance operating inside of a high vacuum system. For alloys, microbalance methods17 must be used on large area samples at relatively low temperatures to minimize composition changes. The

Jan 1, 1962

By P. F. Gnirk, J. B. Cheatham

Single bit-tooth penetration experiments under static load were conducted on six rocks at confining pressures of O to 5,000 psi using sharp wedge-shaped teeth with included angles ranging from 30 to 120° In general, the force-displacement curves for all rocks exhibit an increasingly nonlinear and discontinuous behavior with decreasing confining pressure. The confining pressure at which a rock exhibits a macroscopic transition from predominantly ductile to predominantly brittle behavior during penetration vmies from about 500 to 1,000 psi for the limestones to greater than 5,000 psi for dolomite. The correlation between calculated values of force per unit penetration based on plasticity theory and experimental values is quite encouraging, even at confining pressures as low as 1,000 psi A qualitative correlation between volume of fragmented rock per unit energy input for a single bit-tooth and drilling rate for microbits appears to exist over a confining pressure range of 0 to 5,000 psi. INTRODUCTION Laboratory experiments utilizing a small-scale drilling apparatus have demonstrated that penetration rates are reduced considerably as a result of increasing the confining pressure from atmospheric to a few thousand psi?-3 This undesirable situation can, in general, be attributed to a combination of decreased efficiency of chip removal at the bottom of the borehole, increased rock-failure strength, and a possible change in the mechanism of chip generation and rock fragmentation with increasing confining pressure. To more fully understand the principles underlying the last circumstance, it is the purpose of this investigation to experimentally study the mechanism of single bit-tooth penetration into dry rock at low confining pressures and, in particular, to establish the confining pressure at which the penetration mechanism may undergo a brittle to ductile transition for various rock types commonly encountered in drilling. Confining pressure as used here refers to the differential pressure between the borehole fluid pressure and the formation-pore fluid pressure. EXPERIMENTAL PROCEDURE Using an experimental apparatus previously described: a single, sharp wedge-shaped tool was forced under a "statically" applied load into an effectively semi-in finite dry rock sample subjected to a prescribed confining pressure. To prevent the invasion of the con fining-pressure fluid into the pores of the rock sample during penetration, the exposed surface of the rock was jacketed with a layer of silicon putty.* Electrical instrumentation incorporated into the apparatus yielded a graphical plot of force on the tool as a function of penetration or displacement of the tool into the rock during an experiment. During the course of the experimentation the following conditions were maintained constant: (1) pore pressure — atmospheric (i.e., the rock was dry); (2) temperature — 75F; (3) rate of loading — essentially static (approximately 0.002 in./sec); (4) bit tooth — a sharp wedge-shaped tool loaded normal to the rock surface; (5) rock surface — smooth and flat; (6) drilling fluid — hydraulic oil; and (7) maximum depth of penetration — approximateiy 0.1 in. In addition, each experiment was performed on a different rock sample so the rock surface is free of a layer of cuttings and of any previous indentation craters. The influence of the corners of a borehole was neglected, since each rock sample was cemented inta a section of aLuminum tubing to simulate a semi-infinite body. Indentation experiments were made on a sandstone, two limestones, a marble, a dolomite and a schist with a 60' bit-tooth over a range of confining pressures from atmospheric to 5,000 psi and at a constant confining pressure of 1,000 psi for a variety of bit-teeth with angles ranging from 30 to

Jan 1, 1966