Search Documents

By C. F. Weinaug

In recent years, several methods have been proposed for correlation of physical properties and behavior of hydrocarbon systems. These correlations were found to be unsatis. factory for predicting differential gas liberation and shrinkage data now obtained experimentally. To achieve the desired accuracy, methods were developed to adjust the correlations with two experimental pieces of data. The required data are the saturation pressure at reservoir temperature and the specific volume or density at these conditions. The saturation pressure is used to adjust the vapor-liquid equilibrium-ratio correlations. The density measurement is used to determine the apparent density of liquid methane. With these adjustments, differential gas liberation and shrinkage are predicted within the experimental error. EQUILIBRIUM RATIOS Equilibrium ratios published by the National Gas Assn. of America' were used as a basis for the calculation. The convergence pressure for the sample was determined by using a hypothetical binary system, with methane as one component and the remainder of the system lumped together to form the other component. For the critical temperature of the heavy constituent, the weight average of the critical temperatures of the ethane and heavier components was used. The correlation of pseudo-critical temperature, density and molecular weight of Matthews, Roland and Katzl was used in this calculation for the hexane and heavier fractions. Points for the critical properties of methane and the hypothetical heavy constituent were placed on a chart of the critical loci of hydrocarbon systems. A critical loci for the hypothetical binary system was drawn, and the critical pressure of this system was determined at the reservoir temperature. This critical pressure became the convergence pressure for which equilibrium ratios were read from the NGAA4 charts. The NGAA equilibrium constants for C, through C, were used for the C,; through C, fractions since they were separated on the bases of normal paraffin boiling points. The correlation of Hoffman, Crump and Hocott5 as used to estimate the equilibrium constant for the C fraction. The NGAA eauilibrium-ratio charts are based primarily on paraffin systems and have to be corrected for the presence of naphthenes and aromatics. The correlations in the literature proved too unsatisfactory for determining these corrections. EXperimental data, the saturation pressure, is used with the correlation of Soloman2 to make this adjustment. The correction factors of Soloman were cross-plotted so that the factors for ethane and heavier components were made a function of the factor for methane. a2 = 0.7927 at + 0.1646 .........(1) where a, = correction coefficient for for methane NGAA equilibrium ratio, and a2t = correction coefficient for ethane and heavier NGAA equilibrium ratios. This equation and the equation for the bubble point are solved simultaneously for the correction coefficients. The bubble-point equation is anKnXn = 1.0 . . . . (2) The equilibrium constant for the saturation pressure at reservoir temperature is used. This procedure not only corrects for the relative volatility between components, but also assures that the results will be directly related to the experimental saturation pressure. METHANE DENSITY After determining a set of equilibrium ratios, the apparent density of condensed methane at 60°F and 14.7 psia was estimated. This was found to be necessary because the correlations of Standing3 ere found to be unsatisfactory. The law of additive volumes was employed to make this calculation.

By C. Feng

HE sensitivity of twinning in molybdenum alloys containing 20 or more at. pct Re was reported by Dickenson and Richardson&apos; and Sims et a1.&apos; The latter3 also studied the effect of twinning on plastic deformation of these alloys. Recently, in this laboratory, the twinning planes in a series of Mo-Re alloys were successfully identified by tracing their boundaries on two known surfaces of the crystals, with the aid of stereographic projection. The nominal alloy composition and the results of the determination are listed in Table I. Metallic powders of molybdenum and rhenium with 99.9 pct purity were mixed, pressed and arc-melted into buttons under purified argon in a water-cooled copper hearth; the buttons were annealed at 2400 °C in vacuum (10-&apos;mm Hg) for 1 hr, remelted under argon and finally reannealed in vacuum at 2400°C for another hour. This treatment yielded a homogeneous alloy. Large-grain specimens devoid of subgrains were obtained with grain diameter varying from 1/8 to 1/2 in. A typical specimen weighed about 15 g. Orientations were determined by taking X-ray back reflection photographs for each grain on a polished surface of the specimens. Two surfaces were cut and polished; the angle between the two surfaces was measured with an accuracy of *2o. The twins were produced by compression in a hydraulic press at room temperature. Usually, several sets of twins could be obtained within each grain. Except in cases where no twins were formed, (Crystals 7, 8, and 9) the twinning planes were identified consistently as the (112) family over a wide range of alloying compositions and surface orientations of the crystals, as indicated in the stereo-graphic triangle in Fig. 1. The twins, however, were not broad enough to accommodate a conventional X-ray beam and no identification could be made of the twinning direction. Nevertheless, the direction must be one that lies in a plane of symmetry perpendicular to the twinning plane. There are only two such possible planes in a bcc lattice, namely the (110) and the (100) planes. Of these, only the (110) planes are perpendicular to (112), and the <111> direction is the intersection of the two perpendicular planes. It can thus be concluded that <1l1> is the twinning direction. The twinning elements of molybdenum have been reported as the (112) planes and the <111> directions,4 which are also the operative elements in other body-centered-cubic metals, such as W, a-Fe, and ß -brass.5 Thus, the present finding of the twinning elements in Mo- Re alloys supports the concept that twinning formation is mainly governed by the lattice structure. It is of interest to note that for specimens containing 12 at. Pct Re or less, no twins were found after being plastically deformed to as much as 30 pct reduction in height by compression, while for the 35 at. pct Re alloy, twins were formed with less than 1 pct deformation. This, of course, does not exclude the possibilities that for Re-Mo alloys containing 12 pct Re or less, twins could be formed under different conditions, such as at liquid air temperature or under impact loading. The present results, however, suggest that under the same conditions twinning is increasingly sensitive with increasing rhenium content. Consequently, the plastic strain at which twinning starts was found inversely proportional to the amount of rhenium in solid solution with molybdenum. This plastic strain necessary to cause twinning may be used as a criterion in the study of twinning characteristics.

Jan 1, 1961

By S. A. Hovanessian

Analytical solutions are obtained for calculating the pressure distribution in rectangular fields due to injection and/or producing wells located anywhere within the field. The field is assumed to be homogeneous with either constant pressure or no-flow boundaries. These solutions are extended to include the calculation of average reservoir pressure and permeability from pressure build-up curves. Numerical examples are given to illustrate the application of equations to practical problems. INTRODUCTION Pressure distribution analyses are of considerable importance in the field of reservoir mechanics. A number of papers dealing with this subject have appeared in the petroleum literature. perrinel in 1956 presented a summary of the methods available for pressure build-up calculations and discussed the applicability of each in some detail. More recently (1958), Hazebroek, et a1,2 discussed a theoretical method for obtaining pressure distribution in a radial field. Also in 1958, Nisle3 gave a theoretical solution of the pressure distribution in a field extending to infinity with a partially penetrating line source at the origin. Mathews, et a1,4 applied the method of images to existing solutions for pressure distribution in radial fields and obtained expressions for calculating pressures in bounded reservoirs. However, a literature survey revealed that a general analytical solution for the pressure distribution in a rectangular field, with injection and/or production wells located anywhere within the field, was not available. Such a solution could be used to approximate the pressure distribution in a field the shape of which is nearly rectangular. Furthermore, it could be applied to pressure calculations pertaining to a single well taken from a system of wells which were drilled in a rectangular pattern. In this case the boundaries of the well drainage area are considered to be the perpendicular bisectors of the lines joining the well and its nearest neighbors. The basic problem considered in this paper consists of a rectangular field with either constant-pressure 01 no-flow boundaries and with either an injection or production well located anywhere within the field. The physical properties of the rock and fluids present are considered to be constant. Analytical expressions are obtained for the pressure distribution within the field during the production or injection periods. These expressions, which are infinite series, can be used to evaluate pressure at any point within the field as a function of time and position. For fields containing more than one well (production and/or injection), the solutions for each well acting independently are superposed. That is, at any point in the field, the pressures due to each injection 01. production well acting alone can be algebraically added to give the pressure resulting from all wells acting simultaneously. Evaluation of the series solutions given in this paper is readily accomplished by means of an electronic digital computer; and, in most cases, sufficient accuracy is obtained by setting the upper limits of summation at 20. These solutions are extended to include the calculation of average reservoir pressure and permeability from build-up data. A number of numerical examples are given to illustrate application of the solutions to actual field problems. ANALYTICAL SOLUTIONS The partial differential equation representing the pressure distribution in a homogeneous rectangular field (see Fig. I.), having a single source or sink of strength (2, can be written*

By J. C. M. Li

The steady state shape of a shrinking cylindrical grain boundary of miform boundary energy is shown to be circular. This is based on the principle of either the minimum rate of entropy production or the minimum thermokinetic potential. The circular shape corresponds also to a state of minimum (not maximum) velocity. The steady state shape of a grain boundary moving between two inclined surfaces is a circular cylindrical arc whose position and curvature will be afjected by the nature of the surfaces. The minimum thermokinetic potential is the only valid criterion in this case. THE preceding paper1 reports an interesting experiment in which it is shown that a grain boundary moving between two inclined surfaces assumes a circular cylindrical arc whose axis nearly coincides with the line of intersection of the two surfaces. The driving force is well defined since it is the energy of the grain boundary itself. The velocity of the grain boundary is found not to vary linearly with the driving force. Several interesting questions arise: 1) What principle determines the shape of a moving grain boundary? 2) What is the effect of the inclined surfaces? 3) How do we understand the nonlinear velocity-driving force relationship? It is attempted in this paper to discuss the first two questions based on irreversible thermodynamics and in a following paper to discuss the third question based on dislocation theory. STEADY STATE SHAPE OF A SHRINKING LOOP (A CYLINDRICAL GRAIN BOUNDARY) In irreversible thermodynamics, there are two criteria for the steady state. One is the minimum rate of entropy production2 which is based on the symmetry of phenomenological coefficients and is applicable to linear force-flux relations with constant coefficients. The other is the minimum thermokinetic potential3 which is based on the integrability of a certain Pfaffian differential equation and is applicable to nonlinear systems with variable coefficients. In the case of a shrinking grain boundary loop (a two dimensional version of a cylindrical grain), the flux of atoms across the grain boundary can be taken as the local boundary velocity v measured perpendicular to the boundary (the number of atoms transferred per unit area per second is equal to the product of atom density and the velocity of the boundary). The local thermodynamic driving force can be taken as y/p with y being the grain boundary energy (invariant with respect to the orientation and curvature of the grain boundary) and p being the radius of curvature. Such driving force can be understood by considering the loop as an elastic string with tension y. As shown in Fig. 1, across a length ds along the loop, the direction of y changes by an angle d@ which is ds/p. Because of this change of direction, the tension exerts a force perpendicular to the curve equal to yd@. The force per unit length is then y/p. According to the experimental observation, the velocity v is a power function of y/p, where M is the velocity at unit driving force and &apos;n is a constant. The rate of entropy production per ds along the loop is vyds/pT, where T is the temperature. For a given length L of the loop, the total rate of entropy production is The thermokinetic ptentia13 has only a differential definition and its existence depends on the integrability of this definition. In this case, Eq. [3] happens to be a total differential and the thermokinetic potential can be defined as which turns out to be the same as the rate of entropy production, Eq. [2], except for a constant factor n + 1. Strictly speaking, since the flux-force relation is not linear, the rate of entropy production may not be a minimum at the steady state. However, because of the proportionality between Eqs. [2] and 141, the state of minimum thermokinetic potential happens to be the same as that of minimum rate of entropy production. The steady state shape thus can be obtained by minimizing either integral [2] or [4] for a given L.

Jan 1, 1970

By R. K. McGeary, B. Lustman

The textures produced in zirconium by cold and hot rolling, and by recrystallization above and below the transformation temperature were determined. Thermal expansivities were measured in the thickness, transverse, and rolling directions of preferentially oriented zirconium and were correlated with the texture scatter in these directions. REVIOUS investigations have indicated that minor differences between hexagonal close-packed metals of similar axial ratio may appear with respect to the textures produced both on cold rolling and on subsequent recrystallization. In the case of magnesium, beryllium, and titanium, metals of axial ratio similar to that of zirconium, the ideal orientations produced by rolling are fundamentally the same, although marked variance is reported in the degree and type of scatter about the mean orientation; in those instances where recrystallization textures were observed, they were reported to be similar to the rolling textures. Measurement of the anisot-ropy of thermal expansion of both rolled and re-crystallized zirconium could not be correlated satisfactorily with the textures reported for the above metals, and therefore a study was made of the preferred orientations produced in zirconium. Reported below are the textures produced in zirconium by cold and hot rolling, and recrystallization above and below the transformation temperature, together with the results of thermal expansion measurements. Determination of Preferred Orientation Two types of zirconium were investigated: 1— "crystal bar" zirconium obtained from the Foote Mineral Co., produced by the thermal decomposition of zirconium tetraiodide, and 2—zirconium ingot obtained from the Bureau of Mines prepared by melting sponge zirconium in a graphite resistor vacuum furnace in a graphite crucible. The major impurities present in the two materials used are listed in Table I. Several of the pole figures were later checked with 0.03 pct hafnium crystal bar material and the results were identical with those to be shown for the 1.5 pct hafnium material. The materials were cold rolled to 0.014 in. in thickness as shown in Table 11. Specimens were cut from the 0.014 in. thick rolled sheets and etched to thicknesses of 0.002 to 0.010 in. Such specimens were used for exposures up to a 50&apos; to 60" angle between the beam and plane of the specimen; for higher angles a wire shape, similar to that described by Bakarian,&apos; was formed on an end of the original 0.014 in. sheet. A fine-bladed abrasive cut-off wheel was used to slot the sheet and to form the cylindrical cross-section. The wire shaped ends were then etched to 0.006 to 0.010 in. in diam. Although absorption of X-rays in the wire-shaped specimens does not vary with angle of rotation, the line width around the diffraction rings was not uniform, because the wire was narrower than the X-ray beam, and this condition caused some uncertainty in the estimation of azimuthal intensities. Furthermore, scanning was not practicable with this type of specimen so that spottiness of the rings due to large grain size was excessive for specimens which had been heated above about 650°C. Nevertheless, satisfactory information could be obtained for high angle exposures from the negatives by the use of both types of specimens. Transmission Laue photograms were taken using unfiltered molybdenum radiation (47.5 kv, 18 ma) and a 0.025 in. pinhole. With the film 8 cm from a 0.005 in. thick specimen exposures of about 30 min were adequate. For specimens with a coarse grain size, a device that scanned about 0.15 sq in. of sheet surface was used. An attempt was made to plot the pole figures by use of an X-ray spectrometer as described by Norton.&apos; However, for the particular technique used, the intensity variations obtained were not considered definite enough to give reliable results, especially for the large grained recrystallized and transformed specimens. This method was therefore abandoned in favor of the standard photographic method. Nine exposures were taken of each specimen: seven exposures with the beam perpendicular to the rolling direction and at 0°, 10°, 20°, 35", 50°, 65", and 80" to the transverse direction, and two exposures with the beam perpendicular to the transverse direction and at 60" and 80" to the rolling direction. Additional exposures were then made where necessary. The intensity variations of the diffraction rings were estimated by eye. It was usually possible to estimate 3 degrees of intensity from the photograms but in some cases 2, 4, or 5 degrees were estimated. Experimental Results The preferred orientation was determined for the following treatments: 1—cold-rolled, 2—low temperature rolled, 3—cold-rolled surface layer, 4— cross-rolled, 5—hot-rolled, 6—recrystallized below the transformation temperature, and 7-—recrystallized above the transformation temperature. I—Cold-Rolled Textures: The slip plane in hexag-

Jan 1, 1952

By R. E. Smallman, I. Milne

The deformation behavior of single crystals of Nb-Mo alloys has been investigated with particular reference to the influence of composition, orientation, and temperature. Strong solid-solution hardening was observed reaching a maximum at the equiatomic cotrlposition and can be attributed to the difference in atomic size between niobium and molybdenutrz. Changes in the form of stress-strain curve, as shown by a high work-hardening rate and restricted elongation to fracture, were observed at a composition of Nb-85 pct Mo and are attributed to the presence of MozC DreciDitate. Conjugate slip was only extensive in dilute alloy samples; at the 50/50 composition deformation rnainly occurred by primary slip, and the onset of conjugate slip gave rise to failure by cleavage on (100). The variation of yield stress of Nb-50 pet Mo with orientation was consistent with slip on (011)(111) slip systems. The temperature deperndence of the yield stress between -196" and 250°C was similar to that of pure bcc metals, but at a much higher stress level; no evidence for twinning %as found. IN recent years the deformation behavior of various pure metals in groups VA and VIA has received considerable attention, but surprisingly little work has been carried out on binary alloys made by mixing metals from the two groups. Such an investigation would be of interest since single crystals of metals of group VA have been shown to deform characteristically with a multistage deformation curve1"3 while a parabolic type of deformation curve has been reported for most of the group VIA metals.4&apos;5 It has been suggested by Law ley and Gaigher~ that the difficulty encountered in obtaining multistage deformation curves for molybdenum in group VIA was possibly because of the presence of a microprecipitate of MozC which they observed even at carbon contents as low as 11 ppm. Recently a multistage deformation curve has been reported for molybdenum ," although the stages are not so definitive as those for group VA metals. The binary alloys of the particular refractory metals which have been investigated in single-crystal form include Ta-w,&apos; Ta- Mo,&apos; and Nb- Na." While a large amount of hardening was observed for alloys of the Ta-W and Ta-Mo systems, associated with room-temperature brittleness for alloys approaching the equiatomic composition, Ta-Nb remained ductile over the complete composition range with little or no solution hardening. Other systems have been investigated by hardness measurements on polycrystalline material and a discussion of the hardening of these alloys has been presented by ~udman." The purpose of the present investigation was to examine the deformation behavior of Nb-Mo alloys in detail, with particular reference to alloy composition and single-crystal orientation. In this way it was hoped to shed some light upon the restricted ductility of these alloy specimens. 1) EXPERIMENTAL PROCEDURE The starting materials were obtained in the form of beam-melted niobium rod and sintered molybdenum rod of suitable dimensions. Since niobium and molybdenum form a complete solid-solution series at all temperatures, alloy single crystals were produced by melting the two constituents together in an electron bombardment furnace (EBM). To produce specimens free from segregation a molten zone was passed over the length of each rod six times in alternate directions at a speed of 10 in. per hr. Typical specimens were analyzed for interstitial impurities by gas analysis and for metallic impurities by spectrographic analysis. The results of this analysis are shown in Table I. Many of the tensile specimens were also analyzed (after testing) by scanning the gage length in an electron beam microanalyzer, from which it was found possible to predict the approximate composition of a specimen from the original proportions of each element in the EBM. The tensile specimens were made with a gage length of 0.5 in. and diameter of 0.075 in., using a Servomet Spark machine. By careful machining on the finest range for the final i hr of this technique, surface cracks could be reduced to the level where they were easily removed by electropolishing in a solution of nitric and hydrofluoric acids. The specimens were strained at a rate of 10 4 sec-&apos; using friction grips designed to prevent accidental straining and maintain a good alignment before straining. The orientations of the individual specimens tested are shown in Fig. 1 and the corresponding compositions listed in Table I1 together with collated experimental data. 2)RESULTS a) General Deformation Behavior. The effect of composition on the room-temperature deformation curves of similarly oriented specimens is shown in Fig. 2. The yield stresses of the pure constituents, while not the lowest reported to date, were at least comparable with existing data. Although the solution hardening was large for alloys at either end of the phase diagram, and comparable with the Ta-W solution-hardening data of Ferris et a1.,8 the low work-hardening rate characteristic of niobium was sustained until a composition of Nb-85 pct MO had been reached. Associated with the peak yield stress ob-

Jan 1, 1969

By Bob J. Rawicki

TYPES OF WET DUST COLLECTORS Basically, there are two types of wettype dust collectors. One is mechanical, incorporating pumps, motors, fans, sprays, filters, or flooded beds. These come in many forms, but their operating principle is basically the same. The screen or flooded bed is wetted by a series of sprays. Polluted air is drawn in through a fan; here the dust particles impinge on the screen and are flushed off into a settlement tank. The air then passes through spray eliminators to the atmosphere. A great deal of time and money has been spent developing this type of collector. In laboratory tests where clean water at constant pressure, air at constant volume, and dust at low but constant mass to volume were applied, these collectors gave high percentage dust collection efficiencies. Theoretically and technically, they are very good. In practice, however, under rugged mining conditions where there are no constants, and dust load, particle size, water pressure, and volume of air vary with time, their efficiency varies dramatically. Some of the most common problems are: 1. clogging of screens 2. blocking of nozzles 3. replacement. of damaged pump stators Any of the above result in loss of air flow and dust collection efficiency. Generally, there are so many mechanical moving parts, something is always breaking down. This type of collector is very expensive in parts and labor to maintain. The second type of dust collector, which I personally developed, with the help of the British National Coal Board of Great Britain, has none of the above mentioned problems and it works on a completely different principle. My collector, the Mark III Precipitaire Wet-Type Dust Collector, (an improved version of previously successful models) has no internal moving parts, no flooded beds, screens, spray nozzles, etc. It is highly efficient - 99.8% efficiency on total collected dust and 97% on respirable dust. The pressure drop through the collector is constant, hence, air flow remains constant. There is no maintenance expenditure, be it labor or parts. The only maintenance that there is, is desludging. This, of course, varies with the dust load. With a low dust load, desludqinq may only be required every 4 weeks, but this is a very simple operation that can be handled by unskilled men and it takes only a very short time. My collector is powered by a ventilation fan that is placed on the clean air side for a longer fan and motor life. Dust is extracted solely by water action. The collector&apos;s operating principle is: Dust laden air is drawn in along ducting into a tapered scrubber section, shaped to create a self-induced curtain of water. This action washes the dust particles from the air and the collected dust settles into the bottom of the tank. Cleaned air is then exhausted via a series of spray eliminators where it is completely dried and on into the atmosphere. These collectors have been sold in Great Britain, to the National Coal Board, and throughout Europe. They have been working very well for many years now. GENERAL APPLICATIONS Dust is very dangerous; silicosis and black lung kill and disable many workers each year. Many of the mine explosions can be attributed to coal dust. By the use and application of dust collectors, these deaths and disasters can be avoided. There really is no excuse. However good and efficient a dust collector may be, it&apos;s positioning and installation arrangement is critical. The collector may remove 99% of dust supplied to it, but it takes the skill and experience of a dust control engineer to design and fit ducting, hoods, etc. to remove polluted air at source and deliver it to the collector. In some cases this is very simple; in the case of longwall machines, it is very difficult. In general, dust collectors should be used to collect dust at every point of dust generation. I have noted a great reliance on dilution as a means of dust control. I personally think this is a dangerous and unreliable method. As mines become deeper and the use of larger and more productive equipment becomes more common, reliance on dilution will be impossible and dust collectors will be essential.

Jan 1, 1982

By H. L. Hartman

The estimation of the magnitude of pressure changes which occur in mine ventilation circuits is of primary importance to the mining engineer in making changes in an existing mine or in projecting the ventilation requirements of a new mine. It is a formidable task in many cases, because no adequate theory is available and empirical data are unreliable or even lacking. The common practice of increasing one&apos;s estimate of the mine static pressure by an arbitrary "safety factor" to compensate for the unknown is a poor substitute for accuracy. Generally little difficulty is encountered in calculating friction losses in mine airways. As long as the friction factor can be reliably estimated, the well-known Atkinson formula gives accurate results. However, it is in the determination of shock losses that the ventilation engineer encounters difficulty and is apt to resort to an arbitrary allowance. Shock losses may constitute as little as five or as much as twenty-five percent or more of the overall pressure drop in the mine, and such allowances are generally completely unsatisfactory. Since other shock losses have been covered by McElroy (1), this paper deals with a long-neglected source of shock and pressure change in ventilation circuits, divergence and convergence of airflows, here referred to by their more common mining terms, splits and junctions. PREVIOUS WORK There is practically no fundamental information available in ventilation literature on the splitting and junction of airflows, either in mining or mechanical engineering periodicals or handbooks. Fluid mechanics textbooks seemingly avoid the subject of divergence and convergence, except for brief empirical treatment. The only investigation reported for mine ventilation systems was by Weeks, et a1 (21), in 1933. Although applicable to the case of line brattices in coal mines, it cannot be applied to the general case of mine splitting and junctions involving change of direction with little or no area change. In dealing briefly with the subject of shock losses at splits and junctions, McElroy (1) suggests only that they be considered comparable to bends, with losses at junctions being increased 50 pct to allow for interference of merging streams. He recognizes the importance of quantity ratio but concludes that exact quantitative data are lacking to compute losses at splits and junctions. A very capable summary and analysis of the state of knowledge in industrial ventilation and fluid mechanics regarding losses due to divergence has recently been reported by Gilman (3). In considering all earlier work dealing with both air and water as fluids, he compares results on a common basis and obtains good empirical agreement. For divergent flow, the controlling variable in each branch is the quantity ratio, although the deflection angle is of secondary importance. Empirical equations are presented for straight and 90-deg branches with varying velocity ratios, assuming branches smaller than the main duct. These formulas are modifications of the Borda equation for abrupt contraction and expansion. The only information available on shock losses in convergent flow are approximate experimental data presented by Alden (4) and the Manual of Recommended Practice for Industrial Ventilation of the American Conference of Governmental Industrial Hygienists. Alden considers the effects of both quantity ratio and deflection angle. The data indicate an effect opposite to that observed in divergence, the shock loss varying directly with quantity ratio. SHOCK LOSS THEORY Because pressure losses due to shock have been found to bear a constant relation to the mean velocity of flow in a given conduit, they are frequently expressed as a dimensionless function of the velocity head, termed the shock loss factor. For airflow in ducts and mine airways, the shock loss may be represented by the formula (1): Hx - XHv, (1)

Jan 1, 1961

By A. L. Benham, W. E. Dowden, W. J. Kunzman

The mechanism involved in creating a miscible fluid displacement of typical reservoir fluids through porous media by light hydrocarbon mixtures usually involves a number of mixing and separation stages in the contact area. This mechanism has been illustrated on ternary diagrams for miscible displacement of reservoir fluid by rich hydrocarbon gases. A relatively simple method has been devised for calculating the approximate conditions for a miscible displacement of reservoir fluid by rich, light hydrocarbon gases or by LPG mixtures. Miscibility is favored by increased pressures, decreased temperatures, light reservoir fluids and light hydrocarbon-displacing fluids rich in the heavier components. The conditions for miscible displacements of five typical reservoir fluids by various light hydrocarbon mixtures were determined for a range of temperatures from 70" to 260°F and pressures from 1,000 to 3,000 psia. A correlation of these conditions was obtained using the variables of temperature, pressure, C,+ molecular weight of the reservoir fluid, C,+ molecular weight of the displacing fluid, and mol per cent methane in the displacing fluid. This correlation has been applied to data appearing in the literature with very good success. A partial check of the correlation has been made using a step-wise procedure with a windowed PVT cell. INTRODUCTION The term, "miscible fluid displacement", may be defined as any oil-recovery displacement process where there is an absence of a phase boundary or interface between the displaced and displacing fluids. A water flood would not be a miscible fluid displacement, although cycling in a condensate reservoir would be. Another example is displacement of oil by gasoline. Not so obvious, however, are the conditions under which fluids (such as propane, mixtures of propane and methane, or other similar systems) will give a miscible fluid displacement of oil. There are four standard miscible-fluid-displacement processes being used for the recovery of oil from an "nderground reservoir—(1) high-pressure dry-gas miscible displacement, (2) rich-gas miscible displacement, (3) miscible slug using a slug of LPG followed by a lean dry gas under miscible conditions, and (4) gas driven-slug process using a slug of LPG followed by a lean dry gas under immiscible conditions. The first problem faced by the engineer designing a miscible displacement process is to determine the conditions for miscibility. A number of papers have included discussions of the process of obtaining miscibility in these processes.1-8 However, methods for calculating the conditions for miscibility have not been published. This paper includes a method for calculating the conditions for a miscible displacement of reservoir fluid by a rich gas (LPG diluted with methane or natural gas). In this paper, a rich-gas miscible displacement process is defined as an oil-recovery process in the reservoir. The cost of a rich-gas miscible displacement will be decreased by diluting the LPG with methane or natural gas as much as possible. Therefore, it is important to be able to predict the maximum dilution for miscibility. MECHANISMS The mechanisms for obtaining miscibility between a displacing mixture of light hydrocarbons and displaced reservoir fluid often involve a number of displacing and mixing stages. As stated previously, the cost of a rich-gas miscible displacement is decreased by diluting available LPG with natural gas. The problem is to predict the maximum dilution allowed while maintaining miscibility. To get a picture of the displacing and mixing stages, it is helpful to consider these complex hydrocarbon mixtures as pseu do-ternary systems. Graphical representation of a complex system as a pseu do-ternary system can be thermodynamically rigorous if enough variables are used to describe the character of each pseudo-component in one of the phases. However, even if the characterization of each pseudo-component is not complete, pseudo-ternary representation does serve to illustrate the mechanism for obtaining miscibility and can be used to make some general conclusions concerning the conditions for obtaining miscibility. A phase diagram for the ternary system, methane-

By R. W. Parsons, P. R. Chaney

Oil recovery by the imbibition mechanism can be important in fractured carbonate reservoirs with a bottom water drive. Laboratory experiments were performed on water-wet carbonate rocks to model this process. The results are interpreted in light of the applicable scaling laws. The imbibition production behavior of a pillar of rock subjected to a slowly rising water table can be synthesized from data on total immersion experiments on small rock samples. If the rate of water table rise is slow enough and if the physical properties are such that the imbibition zone is of a small vertical extent, then only the final oil recovery on the small samples need be known. Applying laboratory imbibition data to a real reservoir should be done with caution because of possible wettability problems and the unknown behavior of two-phase flow in fractures. INTRODUCTION The spontaneous taking up of a preferential wetting fluid into a porous medium with the simultaneous expulsion of the contained fluid is termed imbibition. One of the many physical imbibition processes is the production of oil by the taking up of water into a water-wet reservoir rock. In certain types of heterogeneous reservoirs, e.g., fractured, it is realized that imbibition alone may be the dominant mechanism. Theoretical studies of imbibition oil recovery have followed several lines of attack. The behavior of certain mathematically tractable and conceptually simple models have been analyzed. Birks1 considered a bundle of vertical capillaries model, while Perotti et al.2 and Kelemen3 used a model of a series of vertically oriented fractures. If the oil production behavior from a block of reservoir rock is known, then the cumulative behavior of a vertical stack of these blocks subjected to a rising water table can be determined by a method outlined by Aronofsky et al. 4 Knowing the capillary pressure and relative permeability functions should uniquely set the imbibition behavior for a given sample. Solutions of the partial differential equations describing the process, however, present some formidable problems. Digital computer techniques have been applied to displacement processes in which the capillary pressure has been included.5"8 More recently, similar techniques have been used by Blair 9 to calculate countercurrent imbibition behavior. The fact that imbibition behavior is described by the usual fluid flow equations (Darcy&apos;s Law and continuity equation) implies that the derived scaling laws should apply equally to this situation. Use of model experiments to study imbibition has been employed by Graham and Richardson 10 and Mattax and Kyte,11 who actually tested the scaling law applicability. This study is concerned with scaled model imbibition experiments on water-wet carbonate rocks. Of particular interest is the slowly rising water table experiment. This simulates a highly fractured reservoir with the water table rising uniformly through the vertical fracture system. Experiments were performed on different sized samples and with different rates of rise of the water-oil interface surrounding the rock pillar. Total immersion imbibition experiments on smaller rock samples were performed to test the basic linear scaling laws and to see if these simpler tests could be used to synthesize the rising water table experimental results. THEORY IMBIBITION SCALING LAWS Scaling laws for the flow of two incompressible immiscible fluids in a porous medium have been presented by Rapoport.12 These are derived directly from Darcy&apos;s Law for the individual phases and the continuity equation. Results from a laboratory model experiment will duplicate those from some reservoir prototype if a constant proportion is maintained among the three forces acting on the fluids — the capillary pressure gradient, the gravitational gradient and the flowing pressure gradient. Operation of the model (i.e., fluid injection or withdrawal) must be conducted in accordance with specific equations. If the

Jan 1, 1967

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

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

By G. S. Cole, G. F. Bolling

Recent experitnents on crystal growth have indicated that thermal convection in the liquid influences the conditions of solidification and subsequent crystal characteristics. However, the important columnar-equiaxed transition in castings has not been examined as a function of this free convection. An experiment is presented here to demonstrate a change in the transition as convection is suppressed. Lead- and tin-based alloys were cast in stainless-steel molds either without special procedures or with thin grids placed horizontally across the tnold diameter in portions of the castings. While the grids allowed some diffusion of heat and mass transport to occur during solidification, massive fluid flow was hindered in the vegiotz of the grids. A larger percentage of columnar grains appeared in the grid region than in other parts of the special castings or than in all regions of normal castings. This structural difference can be explained in terms of thermal convection. THERE is little information about the possible effects of natural convection which may occur during the solidification of a casting. Ruddlel-3 has stated that the movement of liquid metal produced by convection might assist in the transfer of heat, but in the absence of confirmatory evidence he suggested that this idea be treated with caution. chalmers4 considers, but does not assess, the importance of fluid motion in casting, while Tiller5 suggests that, in liquid metals which have extremely high thermal diffusivities, the effects of thermal convection may be negligible. In addition, throughout the literature, there are statements indicating the possible presence of thermal convection but without proofs or evaluations of its influence. Natural convection may be divided into two categories: thermal convection, where liquid flow is due to density differences caused by temperature gradients in the liquid (e.g., liquid near a cold wall falling with respect to hotter liquid further away from the wall); and solute convection, where flow is due to density differences as a result of concentration gradients, see Wagner.6 In most instances thermal convection predominates over solute convection. Nevertheless, even thermal convection has not been studied in any detail and is difficult to handle, since it can vary quantitatively from one experiment to another depending upon the exact positions of heat sources and sinks or the attitude of the container with respect to gravity. In spite of these difficulties it is possible to realize that a higher rate of heat transfer must result with convection than in any situation where conduction acts alone. Thus, the degree of increased heat transfer is the important undecided question in considerations of casting. It will be demonstrated here that a reduction in convective flow leads to a suppression of the columnar-equiaxed transition in a casting and that this suppression is related to a decrease in the degree of heat transfer. Fortunately it is possible to do this by comparative rather than quantitative measures. The insertion of horizontal plates or grids into a simple mold does not appreciably influence heat transfer per se, but in reducing the allowable height through which colder liquid may fall freely it reduces natural thermal convection. This leads to a reduced fluid velocity, a different flow pattern, and an altered heat transfer which can be measured. EXPERIMENTS The alloys were prepared from elements of 99.9 pct nominal purity and cast into an uncoated 2-in.-diameter stainless-steel mold insulated at the bottom with a suitable refractory so as to ensure predominantly radial heat flow during freezing. Casting was performed either in empty molds or in molds containing grids or plates in various configurations . The grids were constructed of aluminum screen (1/32 in. thick) which was cut to the desired shape, cleaned in aqua regia, and coated with a mixture of 2 parts aqua dag graphite -1 part Sauereisen cement mixture, thus supplying a hard, protective, and insulating layer approximately 0.006 in. thick. The

Jan 1, 1965

By J. Kostuik, M. J. de Bastiani

The method of mining the uranium orebody at Denison mine in the Elliot Lake district of Ontario is discussed. Details of the development of haulage-ways and rooms are giuen, as well as methods of ore haulage. The milling and beneficiation stages are noted briefly. Located in the Elliot Lake mining district of Canada, the Denison mine lies between Sudbury to the east and Sault Ste. Marie to the west, Fig. 1. This area is overlain by Huronian sediments filling a basin in the Keewatin basement complex. Subsequent folding of these sediments developed a synclinal structure, known locally as the Quirke Lake Trough, which forms the northern portion of the now famous big Z structure, the geological trademark of this region. Plunging to the west and extending nine miles from the granite outcropping in the north to the Keewatin lava contact to the south, the syncline has dips of 10" to 45" on the north limb and 10" to 20" on the south limb, Fig. 2. The Denison orebody, located on the north limb of the syncline, lies between 800 and 3000 ft below the surface. The main ore zone is within 50 ft of the Missisagi-Keewatin contact (see Fig. 3), and it varies in thickness from 6 to 32 ft. This ore zone has a composite dip of 19" to the south and is composed primarily of two pyritized conglomerate beds separated by a band of pebbly quartzite stratum. The upper conglomerate bed of the main zone varies from 2 to 15 ft in thickness, while the lower conglomerate varies from 2 to 12 ft. Over extended areas, mining is carried out over the entire section comprising the two conglomerate layers and intercalated quartzite; in other areas, the conglomerates are selectively mined, depending on the thickness of the quartz-ite. The ore zone contains approximately 6 pct sulfides and 80 pct silica. The predominant uranium-bear ing minerals are brannerite, uraninite, and an unidentified mineral (probably thucolite). Some of the uneconomic metallic minerals identified to date are pyrite, pyrrhotite, and chalcopyrite. MINE LAYOUT Shafts: Underground exploration and development (Fig. 4) were performed through two shafts 1/2 mile apart. The No. 1 shaft, sunk to a depth of 1856 ft with the main station established at the 1600-ft ore level, is a 250-sq ft, 5-compartment shaft capable of handling 4500 tpd. Two 5x5-1/2-ft compartments serviced by a 12-ft diam double-drum hoist, are used for hoisting 190-cu-ft bottom-dump skips in balance. A 30-man cage, offset by a counterbalance, operates within a 5-1/2x11-1/2-ft compartment. The fifth compartment is used as a manway, pipe, and cable-way. This shaft is also provided with an intermediate pumping station used during the initial mining phase as well as loading-pocket and spill-pocket installations and pumping facilities. The 480-sq ft No. 2 shaft, located down-dip from No. 1 shaft, has a 2776-ft depth. Serving as the main hoisting shaft with a maximum capacity of 9000 tpd, it has eight compartments, four of which are used for hoisting ore in 10-ton capacity bottom-dump skips. These four compartments, measuring 5x6-ft in section, are serviced by two 12-ft diam double-drum semi-automatic hoists. A 50-man cage is operated in a 15x8-ft compartment, and it is counterbalanced by a 12-ton weight. Two compartments are used for ventilation, manway, electric cables, and pipes. A station was constructed at the 1325-ft horizon of the No. 2 shaft for relay pumping. The main sta-

Jan 1, 1961

By Hans Allenius, R. T. Hukki

This paper describes in quantitative terms the effect of sharpness of classification on the performance of the closed grinding circuit. The analysis is based on a large number of laboratory experiments designed to simulate the two most common industrial closed grinding systems. The experimental results show quantitatively the degree of improvement achievable by more efficient classification. In an earlier paper1 the senior author has presented a trend showing qualitative analysis of mill and classifier performance in tile closed grinding circuit. According to this analysis, the master key to a major improvement appeared to be effective removal of finished fine material from the classifier sands, or in other words, improved sharpness of classification, leading simultaneously to a substantial reduction of the circulating load. One purpose of this paper is to present the results of a quantitative investigation of the closed grinding circuit. Another is to set forth the essential pertinent variables revealed by the quantitative experiments carried out by the junior author. In the two most common systems of closed-circuit grinding, (1) feed to the circuit is introduced into the grinding unit, normally the ball mill, operating in closed circuit with the classifier; or (2) feed to the grinding circuit is introduced into the primary grinding unit, normally the rod mill, discharging into the classifier in closed circuit with the secondary grinding unit, normally the ball mill. In the following discussion, these two systems are analyzed separately in the indicated order. BALL MILL - CLASSIFIER CIRCUIT The Test Apparatus and Procedure: The test appara- tus included: (1) a F 195 mm x 220 mm laboratory ball mill rotated at a speed 77% of the theoretical critical speed, and a 7-kg batch of F 20 mm-F 50 mm steel balls; (2) a conventional 65-mesh test sieve and a Ro-Tap sieve shaker serving as a classifier; and (3) a Permaran instrument manufactured by Outo-kumpu Co. for surface area determinations by the permeability method. As no continuous closed-circuit experiments could be brought about on a laboratory scale, the process was broken down into alternating grinding and sizing steps in such a way that the new feed plus the returning sands always formed a batch of 1000 gm. For each selected grinding period and for each selected sharpness of classification the basic steps were repeated six times. Steady-state conditions were then reached. Crystalline vein quartz was used as a test material. Feed to the process consisted of -10-mesh fraction of this quartz crushed in rolls. This fraction included 15% of-65-mesh material. Experimental Results: The general flowsheet is shown in Fig. 1. Table I gives the essential data obtained representing steady-state conditions under the indicated set of variables. The table is based on 110 grinding experiments, 440 screen analyses and an equal number of specific surface area determinations. Fig. 2 shows the relationship between the cumulative net energy consumption and the cumulative number of mill rotations as evaluated by separate experiments. Note that this relationship is not linear but is instead represented by a slight curve. The data given in Table I are presented in graphical form as follows: Fig. 3 shows the produced -65-mesh material in grams per minute vs. time of grinding in minutes; the sharpnesses of classification at 65 mesh were 100%, 75% and 50%. It is clearly indicated that the highest-capacity figures call for relatively short grinding times and for the sharpest possible classification. Fig. 4 presents the specific surface area on the final -65-mesh product in square centimeters per gm vs. time of grinding in minutes. In conventional mineral dressing processes, a final product characterized

Jan 1, 1969

By S. T. Yuster

An analysis has been made of the factors responsible for tilted fluid contacts in petroleum reservoirs. The factors are static and dynamic with the former being controlled by those variables responsible for the capillary rise of fluids, and the latter include rate of formation tilt and hydraulic circulation. The explanation that rate of formation tilt may be responsible for tilted contacts appears untenable. The role of some of these factors in the migration, fractionation. accumulation and possibly the prospecting for oil is indicated. INTRODUCTION While it is generally considered that the segregation of gas, oil, and water in a petroleum reservoir results in approximately horizontal planes separating these phases. there are factor; which modify this simplified picture to a very great extent. Some of these factors are static. and the departure of the fluid interfaces from horizontal may represent a static equlibriun. Some of these factors may be dynamic and the interface positions may be non-equilibrium and represent a steady .state condition. Of course combinations of these two situations may also exist. The reasons most commonly given for this behavior are hydraulic circulation of fluids, change in character of the formation, and a rate of tilt of the formations more rapid than the rate at which the fluids may re-establish an equilibrium. Russell&apos;s Principles of Petroleum Geology&apos; contains an extensive discussion of tilted fluid interfaces and their importance in the migration and accumulation of petroleum. Russell cites examples of tilted water tables at Lance Creek, Wyo.; Rock River, Wyo.: Cut Bank. Mont.; Cushing, Okla.: and Graham; Okla. Also given was the Wasson Field of West Texas which has saucer shaped oil-water and gas-oil interfaces with the saucer being deeper in the case of the oil-water contact. The Bradford Field in Pennsylvania has a tilted water table. San Ardo. Coalinga None and some other California field.; have tilled contacts. A recent paper by King Hubbert discuses the effect of hydraulic circulation of the water underlying the oil in a given stratun and it- implications in the accumulation of and drilling for the oil. The foregoing would indicate that tilted fluid interfaces are not uncommon and that it would be of importance to underbtand the forces which are responsible. Further. it would be of interest to examine certain idealized models of situations which may cause tilting with the thought in mind of deriving relatinnships between the important variables. These relationships may only hold exactly for the model* assumed but the results may apply semi-quantitatively to actual formations and at least indicate trends. The discussion will be divided into the effect of static factors and the effect of dynamic factors. STATIC FACTORS It has already been mentioned that one reason given for tilted contacts is a variation in formation properties. Specifically this is a change in the pore size of the formation capilaries. If the internal surface of the reservoir is preferetially wet by water, capillarity will cause the water 10 rise to a greater extent above the free water table in a small capillary than in a large one. Consequently the water-oil contact will be higher in the lower permeability sands than in the higher permeability sand.;. The free water table may be defined as the position the oil-water interface would take if capillarity were absent. Since the phenomenon is the capillary rise of water in an oil environment, a derivation of this shold give a relationship between the variables. Fig. 1 is an idealized model of the situation with a uniform bore circular capillary ubstitttting for the tortuous non.circular formation pore. The upward force.; on the column of water will be equated to the forces downward to give the equilibrium situation. forces upward = force downward where r is the radius of the capillary in ceritimetera. a is the interfacial ten4on oil-water in dynes per centimeter, 0 is the

Jan 1, 1953

By Milton F. Goodrich

The idea of using water-only cyclones to scalp clean coal from the feed to other devices has recently been gaining in popularity. 1,2,6 An indication of this popularity is that water-only cyclone scalping is already employed in a few commercial plants. Water-only cyclones scalp the feed to a heavy media vessel in the Warwick plant,3 to Deister tables in the Rowland plant,4 and to a heavy media cyclone circuit in the Amherst Coal Company&apos;s McGregor plant.5 Despite this trend, we found surprisingly little information in the literature on the performance of water-only cyclone scalping. Therefore, we undertook a study of its technical merits, limiting our investigation to a water-only cyclone scalper in conjunction with heavy-media cyclones. We felt this to be a reasonable combination of equipment, since both processes operate well on the same size coal. First, we developed simple mathematical relationships revealing general characteristics of the system. Next, we compared the clean coal yield and quality obtainable by the system on four different coals with that obtainable by heavymedia cyclone circuit alone and water-only cyclones middlings recirculation circuit alone. This was done by using computer program to simulate the operation of the three systems. Finally, we considered the effect that implementation of the system would have on a typical preparation plant circuit. The System The system chosen for study [(shown schematically in Fig. 1)] comprises two stages: a primary stage (PI) consisting of water-only cyclones, and a secondary stage (P2) consisting of heavy media cyclones. The water-only cyclones separate raw coal at a low specific gravity, producing a relatively small amount of high quality clean coal and a refuse containing much good coal. This refuse is then recleaned in the secondary stage, where the sharper separation of heavy media cyclones, operated at a higher specific gravity, reclaims additional clean coal. Finally, the clean coal from both stages is combined as the system&apos;s clean coal, while the refuse from the secondary stage becomes the system&apos;s refuse output. Theoretical Analysis Two assumptions lie behind the following theoretical analysis. The first assumption is that partition curves can represent not only the performance of each stage but also the performance of the system as a whole. These curves show partition number as a function of specific gravity, where partition number is either the percent or fraction of a raw coal with a given specific gravity reporting to clean coal. The second assumption is that these curves are independent of the raw coal characteristics. At any specific gravity the system&apos;s partition number is solely a function of the partition numbers of each stage. For convenience, we use the following variables in describing this relationship: g = specific gravity, PS(g) = the partition number of the system at specific gravity g, P1(g) = the partition number of the primary stage at specific gravity g, P2(g) = the partition number of the secondary stage at a specific gravity g. For this system, we can show algebraically that the relationship is: [Ps (g) = PI (g) + (1- P1(g)) P2 (9)(1)] With this equation the system&apos;s partition curve can be determined if the partition curves of each stage are known. By definition, the partition numbers P1(g) and P2(g) must each be less than, or equal to, one. Placing these constraints on Eq. (1), we arrive at the relations: [Ps(g) % P1(g) and Ps(g) % P2(g)(2)] Thus, at any specific gravity the system&apos;s partition number is always greater than that of both stages. [Figure 2] shows typical partition curves. Curve I represents the primary stage, which consists of water-only cyclones operating at a separating gravity of 1.3. (The separating gravity or the

Jan 1, 1979

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,&apos; van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method&apos;&apos;, as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By G. R. Countryman, G. H. Thomas, I. Fatt

Procedures are described for scaling up turbulent friction pressure drops of inelastic fluids and laminar friction pressure drops of inelastic/elastic fluids in circular ducts. The laminar scale-up calculation method does not require the use of arbitrarily defined flow models, and can be performed directly from laminar capillary viscometry data, or by means of a numerical solution of the integral equation for Poiseuille flow, inputing true shear rate-shearing stress data from other simple shearing experiments. e.g. Couette and plane Poiseuille flows. The turbulent scale up procedure, based on the Dodge-Metzner correlation for inelastic fluids, requires evaluating the characteristic rheological parameters at the existing wall stress under laminar flow conditions. Since the wall stress is the quantity sought in a pipe flow scale-up problem, a trial-and-err or solution is indicated, and therefore a calculation method has been developed and programmed for a medium-size magnetic-drum-memory digital computer. Laminar flow data from other somces may also be adapted to this turbulent scale-upmetbod for inelastic fluids. Comparisons between predicted and experimentally observed friction losses for generalized Reynolds numbers ranging from 3 — 32.000 indicate an error band width of less than 10 per cent in circular tubes ranging in size from 0.5 to 1.89 in. in diameter. INTRODUCTION Considerable effort has been devoted to the problem of establishing an accurate procedure for predicting the laminar and turbulent flow behavior of non-Newtonian materials in pipes from data obtained with small scale laboratory or pilot plant equipment, as evidenced by the variety of reviews on the subject.1-3 Although there does not appear to be a universally accepted scale-up method for laminar flow, the various methods currently employed can be divided into two principal categories depending upon whether a functional relationship between shear rate and shearing stress is assumed, or whether a generalized correlation is adopted. The methods of Hedstrom4 and of Weltmann are typical of the first category, while the generalized correlation of Metzner and Reed6 illustrates the second category. Under appropriate conditions either approach yields a reasonably satisfactory description of laminar flow behavior. In contrast, the prediction of non-Newtonian transitional and turbulent friction losses is much more complex. For example, certain types of materials have been observed to sustain laminar flow at much higher flow rates than Newtonian liquids. When turbulence does develop, the level of agitation is damped to such an extent that the friction factor is lower than that of the original low viscosity Newtonian medium. Extremely broad transition regions and "diameter effects" have also been noted in the turbulent flow of certain chemically complex materials.&apos; Aside from the latter complications, which are attributed in part to viscoelasticity, the existence of a variety of scale-up procedures for predicting transitional and turbulent friction losses is due principally to the chioce of viscosity parameter used in the definition of the Reynolds number. In contrast, Dodge and Metzner have demonstrated that the friction factorReynolds number relationship for inelastic fluids can be satisfactorily generalized provided the fluid property parameters are evaluated at the existing wall shearing stress.* However, the quantity usually sought in a pipe flow scale-up problem is the wall stress, and hence a trial-and-error procedure is required when dealing with non-power law materials. In this report we discuss two calculation methods: (1) for determining the wall stress, and hence the friction pressure drop, for inelastic fluids under conditions of turbulent flow by means of an iterative procedure, and (2) for predicting laminar friction pressure drops without recourse to arbitrarily defined flow models.

By G. W. Doepel, W. P. Sibley

A three-dimensional, multilayer, mathematical model has been developed for predicting performance by mis-c.ible displacement. Areal and vertical coverage, total gas-oil ratio, propane-oil ratio and oil-recovery factor for the multilayer system are computed as functions of hydrocarbon pore volumes of solvent injected. The technique has been programed for use on the IBM-7090 digital computer and was used in calculating perfortnance for continuous miscible displacement of a five-spot, 50-layer .system. Results reported are based on the use of published laboratory data for a single, uniform five-spot layer. However, the technique is sufficiently generalized to calculate performance of any system of layers of any pattern configuration and mobility ratio, provided data are available to describe perfortilance of a single layer. The technique has also been adapted for predicting oil recovery resulting from various gas-driven slug sizes. The results presented for both continuous and slug-type injection show, that permeability variation has 0 significant effect on oil recovery but that solvent composition, within the range studied, has a lesser effect on oil recovery. INTRODUCTION Experience from many of the miscible-phase displacement projects currently in operation indicates that the injected solvent flows preferentially through the more permeable zones in the reservoir. This preferential solvent flow in the high-permeability streaks results in a reduction of vertical coverage (C,) and an attendant decrease in oil recovery. The multilayer technique presented in this paper uses areal coverage (C.,) data for a single layer from the literature and from laboratory measurements of displacement efficiency (herein referred to as microscopic coverage, C,,) to calculate areal coverage (C,) and vertical coverage (C,) for the multilayer system. Oil recovery for the multilayer system is obtained as the product of (Ca) (Cm) (Cv). It is assumed that C, in a single layer is 100 per cent; therefore, (C,) (Cm) is a measure of oil recovery from that layer. Data presented in the literature for a single layer generally deal only with areal coverage and/or microscopic coverage —but seldom specifically with oil recovery. Nevertheless, this technique is flexible in that either Ca and C,, of an individual layer or a direct measurement of oil recovery for the layer may be used. In developing this technique, the problem arose concerning a method of determining the relative amounts of solvent that will enter each layer during any increment of solvent injected. To solve the problem, a dimensionless equation was derived for expressing solvent injected into each layer as a fraction of total injection. A finite-difference method was applied in integrating the equation to determine oil recovery for the system of layers as a function of pore volumes injected. The technique subsequently was programed for solution on a digital computer. The method was applied to predict performance of a deep California reservoir under continuous miscible displacement and has been extended for calculating oil recovery as a function of slug size for gas-driven slug injection. The technique may be applied to any petroleum reservoir that can be adequately represented as a layered system, by considering permeability variation,1-3 permeability range and average permeability. Each layer is then assigned its appropriate permeability, porosity and hydrocarbon saturation. RESERVOIR VOLUME CONTACTED BY SOLVENT It has been recognized for some time that high microscopic coverage (C,) can be obtained by injecting a miscible agent, such as propane, into an oil reservoir. These high values of C,,, have been demonstrated in the laboratory and are normally between 80 and 100 per cent, depending on solvent and reservoir oil composition as well as on reservoir temperature and pressure. In miscible-phase displacement, because there is no interface between solvent and oil, effective permeability to solvent is the same as that to oil. (This assumes that oil saturation has not changed and that water is immobile.) Consequently, mobility ratio (M) simply reduces to that of viscosity ratio . Perrine4 and Habermann5 have shown that large

By W. O. Philbrook, K. M. Goldman, G. Derge

THIS is the third in a series of papers from the Metals Research Laboratory dealing with the transfer of sulphur across the iron-slag interface in a carbon-saturated system. The first paper&apos; showed that the transfer obeyed first-order kinetics. The second&apos; presented detailed evidence in support of the three-stage mechanism: FeS(Fe) FeS(slag) FeS(slag) + CaO(sing) ? CaS(slag) + FeO(slag) FeO(slag) + C ? Fe + CO. This emphasized the role of an iron-sulphur complex as the transfer agent across the interface and as an essential intermediate in the overall process. The present paper shows the influence on the process of the individual elements, carbon, silicon, manganese, and phosphorus, the alloying elements normally present and of interest in commercial blast furnace operations, aluminum, which has a pronounced accelerating influence, and nickel and copper, which are common minor impurities in the raw materials. The experimental methods and procedures for interpretation of data are the same as described earlier.1,2 Induction heated, graphite crucibles were used, and sulphur-free slag was added to the molten iron containing sulphur and the selected alloy. Slag samples were taken at regular time intervals and analyzed for sulphur and other components of interest. The authors wish to state specifically that these experiments were not intended to simulate industrial blast furnace operations. However, they do represent a study of the kinetics and mechanism of sulphur transfer at the slag-metal interface under specified laboratory control of the major variables in the system. It is believed that the results can be usefully applied to the blast furnace if consideration is given to the differences which are known to exist between the two conditions. A rigorous evaluation of the accuracy of rate studies is difficult. Recognized sources of error include sampling and analysis of both slag and metal, temperature measurement and control, selection of zero time, composition of the gaseous atmosphere over the reaction system, and extraneous sources of sulphur from the graphite crucible and reagent materials. In this study carefully standardized procedures were followed consistently in order to minimize the errors. It is believed that the reproducibil-ity of duplicate runs indicates that errors have been rendered negligible relative to the observed effects from intentional variables. For example, original experimental points for duplicate runs are shown in Figs. 1, 2, 7, etc. In most cases in this paper the drawn curves represent an average of at least two heats. The spread in sulphur values at a given time ranges from 0.02 to 0.08 pct, which is a small frac-tioi of the total values ranging up to 2 pct S in slag. In the systematic study of the influence of alloy additions, each element was added over a range of concentrations. The two slags used in this study were designated as "acid" or 1530 and "basic" or 1545. They had the nominal compositions given in Table I. Rates were measured in the temperature range 1500" to 1650°C. Data will now be presented to illustrate the various features of the study. Silicon The first alloying element examined was silicon, because its presence is unavoidable in any silicate slag system through the reaction y/2SiO2(slag) + xM?y/2Si(liq.Fe) + MxOy where M is any reducing agent such as iron or carbon. It is known that the attainment of equilibrium by this reaction is extremely slow and probably not reached in the blast furnace." The influence of silicon on the final stages of the sulphur reaction has also been studied.&apos; It has been shown5 hat silicon increases the activity of sulphur in iron. For such reasons silicon may also be expected to influence the kinetics of the reaction. The observed effects will now be described. The general procedure was to melt 530 grams of ingot iron in a graphite crucible, 3 in. OD, 2 1/4 in.

Jan 1, 1955