Search Documents

DEVELOPMENT of a chemical process for the extraction of pure metals from mill concentrates or metal scrap has progressed beyond the pilot plant stage and may prove an important adjunct to present smelting and refining methods if initial commercial operations prove economically successful. Developed by the Chemical Construction Corp., the new process involves the treatment of oxide and sulphide ore concentrates by chemical methods, instead of the usual smelting and refining techniques. Several of the many applications are now scheduled for commercial use. Refiners using the new process will prepare ore concentrates by standard flotation methods, introduce the concentrate as a slurry into an autoclave along with water and an acid or ammonia. From the resulting leach solution, recovery of individual metals is made by use of suitable reducing agents. By varying conditions during treatment, different metals in the ore are produced separately as pure powders, which may be pressed into forms ready to market, or in the case of copper, extruded as rods or pipe. The reagents are generally recovered. By manipulating the variables during reduction, selective separation of nickel, and/or cobalt, and/or copper can be made simultaneously. The separation is a continuous operation. Low metal prices are not immediately in view because of the tremendous demand at this time. However, reduced metallurgical treatment costs will allow economical mining of orebodies with lower metal content, permitting increased production to meet demand. Although the practical horizon for use of the new process appears to be limitless, one factor may enter into its employment. Just how low grade a concentrate can be economically used is unknown. In addition, for present practical considerations, only metal below zinc in the electromotive series may be processed. Each commercial application requires specific technique adaptation and pilot-plant data for engineering design. A nickel-cobalt-copper process has been researched and piloted in collaboration with Sherritt Gordon Mines, Ltd., for the firm's Lynn Lake properties. In addition, processes were tailored for cobalt concentrates of Howe Sound Mining Co., and National Lead Co., in view of the urgent need for this specific metal. Major General William N. Porter, president of Chemical Construction, says of the new process, "piloting experience has shown that production cost, from ore concentrates to pure metals should be considerably below current costs." Other savings may be realized by cutting transportation and personnel costs, and by reducing the time lag between mining and pure metal from months to a matter of hours. First commercial use of the process will start sometime this summer. Chemico is slated to finish building a $2.5 million cobalt refinery for Howe Sound Mining near Salt Lake City. The plant is expected to raise world cobalt production by 40 pct. The plant will process 35 tons of 20 pct cobalt concentrates from Howe Sound's Blackbird mine near Cobalt, Idaho, daily. Yearly production is expected to reach 2000 tons of pure metal, about one half of U. S. consumption in 1950. Steps in the Howe Sound application of the process are: 1. Acid oxidation leach. 2. Filtration (residues thrown away are insoluble compounds-gangue, iron, arsenic). 3. Cementation (for Cu removal because Cu content is too low for recovery). 4. Reduction from ammoniacal solution. 5. Separation of Co and Ni as mixed metals (95 pct Co and 5 pct Ni). 6. Recovery of ammonium sulphate. Under construction at the Fredericktown, Mo., mine of National Lead is a $5 million refinery sched- [ ]

Jan 1, 1952

By R. T. Hukki

This paper presents a preliminary analysis of the fundamental relationship between the net energy used and the respective product size throughout the entire range of sizes covered by crushing and grinding, and an attempt to find a sensible correlation between the existing theories. Walker, Lewis, McAdams, and Gilliland' have given the following differential equation of a general form for comminution: dE = -C§ [1] where E is the net energy required per unit weight in a certain process of comminution; x is the factor indicating the fineness of the product; n is the exponent indicating the order of the process, and C is a constant related with the material, units chosen, etc. H exponent n in the above equation is replaced by numerical figures 2, 1, and 1 1/2, the integrated form of the general equation leads to the well known fundamental theories represented by the law of von Rittinger,' law of Kick, and the third theory of comminution by Bond, respectively. The total net energies (E) from infinite feed size to a product of size x are as follows: The net energy required in a certain process of comminution is proportional to the new surface developed according to the law of von Rittinger, to the weight or size of the bodies treated according to the law of Kick, and to the length of the new cracks formed which initiate breakage according to the theory and explanation by Bond. On a logarithmic paper, where particle size is presented on the abscissa and energy consumption on the ordinate (see Fig. 21, all three relationships are represented by straight lines. The slope m of the line according to the law of von Rittinger is equal to -1.0; of Kick, 0; and of Bond, -0.5. Experimental evidence in favor of the law of Rittinger has been presented, e.g., by Gross and zimmerley5 on quartz crushed in a drop weight crusher and evaluated for surface by the method of dissolution, by Dean on magnetite crushed by a similar method and evaluated for surface by the determination of coersive force, by piret7 and co-workers on a group of minerals crushed again by a similar method as well as by compression and evaluated for surface by permeability and gas adsorption methods, and by schellinger8 on a group of minerals ground in a calorimetric ball mill and evaluated for surface by gas adsorption. Experimental evidence in favor of the law of Kick seems to be scant in the field of comminution. On the other hand, in the field of mechanical engineering Kick's law seems to be of fundamental nature in processes such as cutting, pressing, shaping, and rolling of metallic substances. Experimental evidence in favor of the third theory has been provided by Bond." To a large extent, data are based on the vast amount of grindability tests performed in the laboratories of Allis-Chalmers Manufacturing Co. In addition to the devoted proponents of one or the other of the basic theories listed above, certain investigators have indicated that one of the theories might be applied for a certain range of sizes, while another theory might be used for other sizes. In a discussion of a paper by Bond," Dobie" presented a statement at the International Mineral Dressing Congress in London (1952) indicating that 1) for large particles, the law of Kick was approximately correct; 2) for finer particles, von Rittinger's suggestion was nearer to the truth; and 3) Bond's new

Jan 1, 1961

By D. M. James, J. C. Martin

The results are presented of a study of the application of analytical methods to the solution of two-phase flow into single wells. Approximate analytical expressions for the pressure distribution in two-phase flow are found ;;; a number of conditions. The results obtained from the analytical solutions me found to be in good agreement with results obtained by finite difference techniques using a high speed digital computer. Mathematical solutions for four sets of boundary conditions are presented. All of these solutions are composed of a short-term transient plus a steady or quasi-steady state. The rates of decay of the short-lived transients are analyzed. It is found that the durations of the short-term transients may be characterized by a parameter defined as the time constant which can be determined from simple relations. It is shown also that if the outer radius is much greater than the radius of the well. the short term transients decay at rates which are proportional to the square of the exterior radius, and the rates of decay are only slightly dependent upon the radius ratio. Numerical solutions based on finite difference techniques are presented for a number of conditions. The numerical solutions are in good agreement with the predictions based on the theoretical analysis for small and moderate drawdowns. Examples involving large drawdowns indicate that the nonlinearities in the equations of flow do not appreciably alter the longevity of the short-term transients. In all cases the time required for the short-tern transients to disappear is predicted satisfactorily. INTRODUCTION The mechanism by which oil and gas flow into a single well is of vital interest to the petroleum industry. The fundamental equations of two-phase flow which describe this mechanism are nonlinear partial differential equations. Numerical solutions of these equations describing pressure transients have been obtained with the aid of electronic computers. Although solutions obtained in this manner take into account a large number of effects, the reduction of this information to useful generalities is difficult. One method of obtaining generalities is the use of linearized approximations of the nonlinear equations. Since it is possible to obtain explicit solutions of the linearized equation, general properties of the role of pressure in the flow mechanism may be ascertained. Results obtained from this approach are limited to some extent by the linearizing assumptions. The severity of these limitations may be evaluated by comparing solutions of the linear equation with numerical solutions of the more exact nonlinear equations of two-phase flow. In the past considerable amount of work has been devoted to studying pressure build-up using the single-phase flow theory. Unfortunate1y, most prep sure build-up tests involve multiphase flow. A small amount of work has been done studying pressure build-up where the flow is two-phase. The encouraging results of these studies suggest that useful results may be found from additional studies of not only pressure build-up, but also the rapid transients associated with placing a well on production. This paper presents the basic theory of the pressure transients associated with placing a well on production and with closing it in. The paper is concerned chiefly with two-phase compressible flow; however, the results also apply to single-phase flow. The results are based on analytical solutions of the flow equations, and they are verified by numerical solutions using finite difference techniques. Much of the previous work on compressible flow into wells has been confined to single-phase flow. Some work has been done on compressible two-phase steady state flow, and solutions of the equations of flow have been found by finite difference techniques using high-speed computers. Muskat1 presents some rather general solutions to the equations of single-phase compressible flow into wells. Much work has been done on pressure build-up in wells (see, for example, Refs. 2-6). Almost all of the work on pressure build-up concerns single-phase flow with the exception of Ref. 5 and part of Ref. 2. Some work has been done on pressure fall-off in injection wells.718 Muskat9 presents the solution of the equations for radial steady state two-phase compressible flow. Handyl0 extends Muskat's solution to obtain

By S. Kobayashi, A. E. Armenákas, C. Mylonas

Earlier work1-4 has shown that commercial mild steels under static loading at the lowest natural operating temperatures fracture in a brittle manner only when damaged by a suitable history of straining. Notched and then compressed plates have fractured in subsequent tension at loads as low as 10 pct of the limit load and precompressed smooth bars at a subsequent extensional strain as low as 0.01. The comparison of the average net fracture stress with the flow limit stress was shown to be an excellent criterion of brittle or ductile behavior of mild steel structures when only loads and general stress levels are known. The present investigation shows that the deformation at fracture is a far more sensitive measure of brittleness than the fracture stress. Both criteria, deformation and stress, are used to determine the amount of uniform precompression of two mild steels, ABS-B and Project E-steel, resulting in brittle fracture under the strong constraint of a subsequently machined severe circumferential groove. The fracture stress equaled or exceeded the theoretical flow limit Ól = 2.680., ('based on the elevated yield stress at each prestrain) up to prestrains of about 0.20 and then fell off to about 1 at prestrains of 0.60. Yet a prestrain of only 0.05 reduced the elongation at the shoulders by a factor of about 4. The total plastic elongation of a region surrounding a sharp notch in prestrained steel deternines whether or not fracture will be initiated in large structures, hence is a direct and realistic measure of the remaining ductility and provides an excellent test of the material's resistance to embrittlement. ALMOST all engineering structures of mild steel which failed in service in a so-called "brittle" manner exhibited appreciable local plastic deformation (more than 1 pct) and a mixture of "cleavage" and "shear" fracture surfaces. Furthermore, the peak local stress in an advancing crack is always high. Such simple criteria as pure cleavage, absence of plasticity, or fracture stress, so useful in other fields, are insufficient to characterize brittleness in engineering structures at ordnary temperatures. The "brittle)' behavior of engineering structures is caused by the limited ductility at cracks or notch regions, as shown by recent work at Brown University summarized and extended in Refs. 1 to 4. It was also shown that the local ductility of mild steel may be catastrophically reduced by a suitable strain and temperature history similar to those which may occur in real structures. Localized yielding begins at the notch roots at low loads but is contained by surrounding elastic regions. The plastic strains are hence small. They increase slowly with the load up to the flow limit or limit load for an ideally plastic material. Unrestricted plastic flow then occurs. At such strains the real material locally st rain-hardens and fractures. With work-hardening materials no flow limit exists, and the transition from low to high plastic strains is gradual but increases more rapidly at loads close to the flow limit of an equivalent perfectly plastic material. When the available ductility is sufficient for reaching the limit load, the behavior is normal or "ductile". With insufficient ductility, fracture will occur at a low load and will be defined as "brittle". Accordingly the sufficiency or not of the ductility at a notch is reflected in the magnitude of the fracture load or average net stress as compared with the limit load or flow limit. The application of this criterion gave surprising results. Service failures under semistatic loading occurred mostly below the limit load, hence were "brittle". On the contrary all static tests with undamaged commercial mild steels reached the limit load in spite of the deepest notches and temperatures below Charpy impact transition (say about -30°F or higher). Notched bars showed no brittleness even under fully dynamic axial loading up to 1.5 x107 psi per sec at -23°F and about 2.5 x l06 psi per sec at -ll0°F, but broke in a brittle manner at -200°F and l03 psi per sec.' To propagate a fracture at low load, experimenters had to trigger it by superimposing a strong dynamic impact at a notch cooled by liquid nitrogen as in the Robertson5 or Esso tests.' Once started the cracks would run through warmer regions of smaller stress. Obviously the ductility of the apparently undamaged commercial steels was "sufficient" in the laboratory tests but not in the service fracture, where it must have been reduced by some embrittling procedure of fabrication or service. These conclusions were confirmed by the achievement of low static stress (brittle) fracture initiation in unwelded steel after a local reduction of ductility by simple prestraining. Symmetrically notched plates of undamaged mild steel cooled below the Charpy V-notch transition range were tested in central static tension: Undamaged plates withstood loads of limit intensity, but when subjected to prior in-plane com-pressive prestraining perpendicular to the notch axis and to accelerated aging, the plates fractured completely or partially at very low static loads, as low as one-tenth of the flow limit.2-4 However, damage could not be easily related to the strongly variable prestrain around a sharp notch. To obtain easier strain measurements, it was decided to use axially precom-pressed bars and bent bars. Final testing was done by

Jan 1, 1970

By W. R. Ingalls

THE following chronology presents the history of lead-mining in the United States in a brief form and is a useful reference in connection with the statistics of production 1621. Lead was mined and smelted near Falling Creek, Va., to supply the local demand for bullets and shot. This was the first mining and smelting of lead in what is now the United States. 1632. In a report made on the minerals of New England, lead-ore is mentioned. (Bishop, I., 470.) 1650. Supposed beginning of mining by Jesuits in Pima county, Ariz. 1651. Grant of lead-mine at Middletown, Conn., to Gov. John Winthrop.

Jan 9, 1907

By L. D. Mullins, E. B. Elfrink, W. L. Wahl

In the past few years several articles and papers presenting results of solution gas-drive depletion calculations have appeared in the lit-erature. Such calculations are of interest to the oil industry, for investment decisions must often be made before much is known about a reservoir. At other times, an estimate of the possible benefits to be realized from alternate production methods is desirable, and theoretical depletion calculations can serve as a floor or reference level from which to work. In any case, an estimate of ultimate oil recovery based upon engineering data is commonly required. An engineer confronted with the problem of obtaining, for a specific reservoir system, an estimate of ultimate oil recovery by solution gas-drive depletion usually will be forced to perform the calculations himself. This is despite the quantity of data in the literature. Rarely will either experience or the literature provide results from a reservoir system similar in all important respects to the one under consideration, and calculated results are not so plentiful that satisfactory interpolation procedures can be devised. Performing the calculations, however, is a tedious, time-consuming task unless an electronic computer is available, and, in practice, time and manpower are not always available for this purpose. A quick, simple, consistent method was needed for reducing the uncer- tainty in estimated oil recovery from solution gas-drive reservoirs when only minimum information about the reservoir system is available. PROCEDURE Method of Calculation The usual requisite assumptions were made so that the material balance equation could be used to calculate data for the charts. The following assumptions were made: (1) the reservoir is homogeneous and isotropic; (2) oil recovery is due entirely to solution gas drive and neither a gas cap nor a water drive nor gravity drainage is present; (3) the initial reservoir pressure is the bubble-point pressure of the reservoir fluid; (4) initial total liquid saturation is 100 per cent of pore space; (5) interstitial water saturation remains at the initial value as the reservoir pressure declines from the bubble-point pressure to atmospheric pressure; (6) equilibrium gas saturation is 5 per cent of pore space; and (7) oil and gas saturations are uniformly distributed throughout the reservoir at all times. There are no saturation gradients due to a wellbore, nor is the geometry of the reservoir system considered. The material balance equation was written in the form of a differential equation' which was integrated to determine the change in oil saturation for an assigned pressure drop. Formal integration was not possible, so recourse was made to the Runge-Kutta method' of numerical integration. All computations were performed on IBM equipment. Numerical integration yielded the change in oil saturation within the reservoir as the pressure. declined from the bubble-point pressure to atmospheric pressure. The initial oil saturation minus the change in oil saturation yielded the oil saturation at atmospheric pressure. The oil originally in place was obtained by dividing the initial oil saturation by the initial formation volume factor (differential liberation) while the oil in place at atmospheric pressure was obtained by dividing the final oil saturation by the formation volume factor at atmospheric pressure. U1timate oil recovery, expressed as a percentage of the initial oil in place, was obtained by dividing the difference between the oil initially in place and the oil in place at atmospheric pressure by the oil initially in place and multiplying by 100. PVT Data Charts were based upon 135 solutions to the material balance equation. PVT properties of the reservoir fluids were variables in this equation and had to be known as functions of pressure. The required PVT data might have been obtained from either actual reservoir fluid systems or correlated data. However, correlated PVT data were developed and employed in the calculations for the following reasons: (1) it is doubtful that 135 sets of PVT data could have been obtained for the values of variables investigated in this study, and (2) results of recovery calculations from randomly obtained PVT data could not be correlated as well as rssults from selected PVT data. The PVT data used to develop the

By A. R. Gregory

A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy. S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King, 1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. 5- and P-wave velocity measurements on Bandera sandstone and Solenholen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2 and by Gregory 3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fate sand by Gregory, 3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King 1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8 A new shear wave velocity laboratory apparatus has been developed for testing rock samples under elevated pressures. This apparatus employs liquid-solid interfaces as mode converters. The technique for generating pure shear waves at liquid-solid interfaces was developed from studies of critical angle methods.6 Two identical transducers are placed in contact with opposite ends of the rock sample, one transducer acting as an ultrasonic energy emitter or generator and the other acting as an ultrasonic energy receiver or motor. The P-wave energy emitted from the piezoelectric crystal in the generator transducer is converted into pure S-wave energy at the Drakeol-aluminum

By John N. Dew, William L. Martin, `

This paper describes the results of a laboratory investigation conducted to obtain data for an evaluation of the in situ combustion process as a method of producing crude oil from reservoirs. Air and fuel requirements, rates of advance, com-bustion temperatures, and coke and fluid distributions are presented. The mechanism of oil recovery by in stiu combustion is discussed. Five crude oils ranging in gravity from 10.9 to 34.2"' API were prodrrced from a semiadiabatic, uncon-.solidnted sand pack by in situ combustion. Experimental conditions were varied over a wide range in order to determine the inter-relation-.ships of process variables. The mini-mum air flux requirement for self-sustained combustion was found to he less than 10 scf/hr-ft'l. The rate of advance of a self-sustained combustion zone was found to be nearly proportional to the air flux at the combustion front. The effects of presrure and injected air flux were studied in a series of experiments using a 21.2" API crude. A minimum air require-ment was observed at an air .flux of 20 scf/hr-ft'. The oil saturation con-sumed as fuel averaged 5.5 per cent of pore volume. The effect of air pressure was, found to be .small for experiments having high combustion eficiencies. This study should promote a better understanding of the problems and mechanisms involved in labora-tory investigations and field applications of the in situ combustion pro-cess. The data presented will be use- ful in the interpretation of results of field tests. When tempered with volu-metric sweep efficiencies, the data can be used in making preliminary economic appraisals of the process as applied to reservoirs containing high porosity unconsolidated media. INTRODUCTION The purpose of this work was to obtain laboratory data for an evaluation of the in situ combustion process as a method of producing crude oil from reservoirs. In situ combustion basically consists of (1) injecting air into a reservoir through selected input wells to create an air sweep through the reservoir, (2) igniting the crude at the injection well, and (3) propagating the combustion front through the reservoir by continued air injection. By this means, oil is swept toward producing wells in the area. The fuel for combustion is supplied by heavy residual material (coke) which has been deposited on the sand grains during distillation and cracking of the crude oil ahead of the combustion front. Recovery of petroleum by a combustion or heat wave process is not a new idea. F. A. Howard was granted a patent in 1923 on a process in which air and a combustible gas were pumped into an injection well and ignited.' Russian engineers reported on field experiments with a crude oil gasification process in 1935'. Other known early field tests include those conducted near Bar-tlesville and Ardmore, Okla., in 1942." More recently completed field tests include an inverted seven-spot by Sinclair Oil Co. in the Delaware-Childers field, Nowata County, Okla.'; a three-well and an inverted five-spot test by Magnolia Petroleum Co. in Jefferson County, Okla."."; and a test by California Research Corp. in the Irvine-Furnace field in Kentucky. The Worthington Corp. has completed a test in cooperation with the Forest Oil Corp. in Clark County, 111.l Several field tests are now in progress. Sinclair has acquired a 600-acre lease in the Humboldt-Chanute field, Allen County, Kans., and is under way with a large-scale operation. Three field tests are under way in California. The General Petroleum Corp. is conducting an inverted five-spot pattern test in the South Bel-ridge field in Kern County under a cooperative agreement with 11 other companies, including Continental Oil Co. California Research Corp. is conducting a four-well pattern test in Midway-Sunset field near Maricopa, and Richfield Oil Corp. is under way with a test in the Ojai field in Ventura County. Although the published information is valuable concerning results of field tests, little laboratory data other than that reported by Kuhn and Koch5 are available to aid in appraising the tests from a technical or economic standpoint. Engineering data needed to evaluate the process include the amount of air required per barrel of oil recovered, minimum injection rates which will support combustion, rates of advance, air required per unit volume of reservoir cleaned, amount of oil recovered, and amount of oil consumed as fuel. The present laboratory investigation was undertaken to provide these data for that portion of the reservoir swept by the combustion zone. The first phase of this investigation was exploratory in nature and was conducted at injection pressures less than 100 psig. Data were ob-

By A. W. Adamson

A model for transport processes in liquid mixtures is discussed which supposes that the elementary act involves a position exchange between two species and that the exchange is so confined by the solvent cage as to occur nearly isosterically. The rate-determining step, thus, is likened to a bi-molecular reaction and is so treated, using absolute rate theory. The cage model has been applied to diffusion, thermal diffusion, sedimentation and viscosity, but only the first and last of these phenomena are emphasized in the present paper. The model leads to semi-empirical relationships between the absolute value for a digusion coefficient and the activation energy for diffusion, between mutual and self-diffusion coefficients and for the variation of the viscosity of a binary mixture with composition. These are discussed in relation to experimental data for various systems, including hydrocarbon mixtures. It is shown that the proposed viscosity equation and seven other commonly used ones all may be regarded as special cases of a single general relationship; a brief critical analysis is made of the basis of selection of one or the other for data fitting or interpolation. INTRODUCTION AND GENERAL THEORY The present paper covers a brief discussion of a cage model for transport processes in liquid mixtures and how this model may be useful in treating the diffusional behavior and the viscosity of such systems. Since diffusion requires the more detailed treatment, it will be taken up first, and the model then applied to viscosity. There are two types of diffusion coefficients that may be measured experimentally, apart from thermal diffusion quantities. The first is the mutual or binary diffusion coefficient, D which may be defined in terms of Fick's first law. This states that the permeation, or flux P, is proportional to the concentration gradient. In the usual experiment, P is measured relative to a frame of reference fixed with respect to the medium (e.g., the diaphragm in a diffusion cell); as a consequence, the same value of D is obtained regardless of whether P and C refer to Component 1 or to Component 2; i.e., there is only one independent mutual diffusion coefficient for a binary system. In addition to D there will be various self-diffusion coefficients. defined in terms of the gradient in labelled species i and its permeation in an otherwise uniform medium. The thermodynamic approach to mutual diffusion supposes that the actual driving force is the gradient of the chemical potential, i.e., that In the case of a dilute solution of solute, Eqs. 1 and 3 lead to the Einstein equation, If the solution is ideal and the friction coefficient is taken to be then the familiar Stokes- Einstein equation results. Mutual and self-diffusion coefficients can not be related on general thermodynamic grounds; it is necessary to invoke some additional assumptions, i.e., a model; several such have been proposed. Hartley and Crank' supposed the existence of separate, intrinsic diffusion coefficients (Dl and D2) for each component, essentially corresponding to the two self-diffusion coefficients. The two flows can not be independent, however, but must be coupled through the usual restriction that there be no net volume flow. For an ideal solution. one then obtains' Glasstone, et al' treated diffusion in terms of absolute rate theory, but their approach otherwise resembled the previously mentioned one in that each species was considered to move with respect to the general medium in a manner determined by its individual jump distance and specific rate constant. For other than dilute solutions, a coupling of flows leading to an equation such as Eq. 6 would again be present. However, as required by Eq. 6, one does expect that the self-diffusion coefficient for the solute and the mutual-diffusion coefficient for the system become identical at infinite dilution. Lamm4 recognized that there should be three distinctive interactions in a two-component system-1-1, 1-2 and 2-2 — and, therefore, proposed three rather than two fundamental friction coefficients. Mutual diffusion resulted from 1-2 interactions only, and self-diffusion resulted from 1-2 plus either 1-1 or 2-2 interactions. Again, a collective coupling between all motions was imposed to meet the condition of no net volume flow. Laity' has shown how to convert the Onsager equations to a form very similar to Lamm's. Cage Model For Diffusion Work in this laboratory on diffusion in aqueous sucrose solutions made it apparent that three, rather than two, interactions were indeed needed," but considera-

By Carle R. Hayward

The growing importance of metallography has caused a corresponding interest in the improvement of apparatus for preparing specimens of metals and alloys for microscopic examination. The purpose of this paper is to describe an electric heating-furnace, a grinding- and polishing-machine, and a device for mounting specimens, which are used in the metallographical work at the Massachusetts Institute of Technology. These three pieces of apparatus were designed and made in the laboratory of the Institute, and each possesses some original features which may be of interest. Electric Heating-Furnace. The accurate control of heat necessary for metallographic work is best obtained by the electric-resistance furnace. Uufortunately, however, even the platinum-wound furnaces cleteriorate with use and ultimately burn out, while the cheaper resistance-materials are often short-lived. Since the resistance-coil must be replaced from time to time, it is of advantage to be able to make this change with as little trouble as possible, and this is an important feature in the furnace shown in Fig. 1. A is a cylindrical galvanized-iron can, with two handles, B. A porcelain tube, C, passes through the central hole in the cover of the can and rests upon the asbestos disk, D; the bottom of the tube is held in place by the asbestos ring, 3. The space, F, around the tube is filled with powdered magnesia. G is an ordinary assay-crucible, with the cover inverted so as to present a smooth top for supporting a second porcelain tube, H. The latter is wound with " Excello " resistance-tape, 0.014

Jan 1, 1912

By Leon Y. Sadler, E. K. Landis

Froth flotation has been described by several authors" as being analogous to a first-order rate process. Although a few investigators"," have found orders other than one fit their data best, the rate is described most often by the first-order equation, [ ] where r is rate of disappearance of floatable mineral by flotation, C(t) is concentration of floatable mineral in the cell, t is time, and k is first-order flotation rate constant. One of the difficulties in using this rate equation to develop a mathematical description of an actual flotation process is the determination of a suitable value of the constant k. The value of k, which reflects the floatability of the floatable mineral, is dependent on the hydrodynamic conditions in the cell, type and concentration of flotation reagents used, pulp density, and pH, method of froth removal, aeration rate, particle characteristics, and many other factors. Therefore, it is reasonable to expect that each particle of floatable material has a distinct value of k and that a unique distribution of k values exists for given feed and cell conditions. This has been observed and discussed by several investigators, and mathematical descriptions of the batch cell and steady-state continuous cell with distributed flotation rate constant have been developed." The mathematical description in the form of a transfer function of the continuous or flow cell under unsteady-state conditions with a distributed flotation rate constant has not been presented in the literature. The purpose of this note is to make such a presentation. The basic features of the perfectly mixed continuous

Jan 1, 1974

West Virginia University, Engineering Experiment Station, Morgantown, W Va C E Lawall, Director. For copies of Bulletins or other information address the Director Two classes of the Bulletins published by the Engineering Experiment Station are of particular interest in this directory- Research Bulletins, and Extension Bulletins, or bulletins published for use in extension work They are. Research Bulletin 1, The conservation of natural gas in West Virginia (1918), 2, Acceleration stresses in hoisting ropes (1927), 3, Classification of gas wells for operation (1927), 4, Mining methods in West Virginia (1929), 5, The whiteware industry of West Virginia, Extension Bulletin 1, Electricity in mines, series

Jan 1, 1933

(Under this heading will be published notes sent to the Secretary of the Institute by members or other persons introduced by members. ) Member, age 33, single, desires position as null superintendent or charge of mine and mill. Twelve years' milling, cyaniding and flotation in the southwest and Mexico. Will go anywhere. No. 458. Member, mining engineer and geologist, graduated, from- mining college, age 27, single, exempt, with practice in metal mining, surveying of any kind, field geology and ore testing, desires to .engage as geologist or mining engineer with a reputable concern in the west or east of. U. S. Hard worker and reliable. No: 459.

Jan 6, 1918

By W. T. Weller

In Part I, a study of pressure build-up curves calculated for conditions under which both oil and gas flow led to the conclusion that the presence of a dispersed free gas phase in an oil reservoir must be taken into consideration to estimate accurately average reservoir pressure and permeability from build-up curves. Familiar methods based on the assumption of no free gas can be extended to the two-phase case by using total compressibility and mobility in place of single-fluid compressibility and mobility. These methods give correct values for average pressure and permeability when gas saturation is small. Errors become larger as the gas saturation increases. However, for the use to which the results will be put, the methods are satisfactory for reservoir engineering purposes. An improved method of calculating the performance of depletion-type reservoirs is presented in Part 2. Because the mathematical relationships describing simultaneous flow of oil and gas are quite involved. simplifying assumptions are made to provide means of obtaining approximate solutions of reasonable accuracy. One such approximate method now in use is the constant-GOR solution. It involves the assumption that at any instant, the ratio of total gas flow rate (both free and disolved) to oil flow rate is the same at all points in the reservoir. The approach is not applicable unless the free gar saturation in the reservoir is everywhere greater than the critical gas saturation. This paper presents a modification which, by avoiding the constant-GOR assumption, makes the method applicable to all reservoir conditions, and so far appears to be more accurate than the constant-GOR solution and to he comparable in required compuctation time. PART I— BUILD-UP CURVE ANALYSIS INTRODUCTION Pressure build-up characteristics of shut-in wells have been used for many years by engineers to evaluate average reservoir pressure, effective permeability thickness of the pay section, effectiveness of well completion (skin effect) and reservoir size. A number of methods of analysis have appeared from time to time."" Without exception, these methods are based on the assumption that the reservoir contains but one fluid of constant small compressibility and constant mobility. It has been suggesteda" hat these single-fluid methods may be applied to data from reservoirs containing both oil and gas by substitution of some effective total properties of the multiphase system in place of the corresponding single-phase properties. The present investigation was undertaken primarily to evaluate this approach. METHOD A number of theoretical build-up curves were calculated for conditions of two-phase flow, under the assumption of certain reservoir and fluid properties, and were analyzed by single-fluid methods with appropriate total compressibility and total mobility values for the corresponding single-fluid properties. Results of the analyses were compared with the assumed conditions. The theoretical build-up curves were completed by procedures similar to those of West, Garvin and Sheldon." Since these calculations require a considerable amount of computer time, an attempt was made to derive an approximate calculation method. The attempt was unsuccessful for calculating build-up curves, but the effort did result in a new approximate method of calculating the performance of solution gas drive reservoirs, which appears to be an improvement over the constant-GOR method" used previously (see Part 2). The West, Garvin and Sheldon calculations involved the following assumptions: (1) the reservoir is circular and completely bounded, with a completely penetrating well at its center; (2) the porous medium is uniform and iso-tropic, with a constant water saturation at all points; (3) gravity effects can be neglected; (4) compressibility of rock and water can be neglected; (5) the composition and equilibrium are constant for oil and gas; (6) the same pressure exists in both the oil phase and the gas phase; and (7) no afterproduction occurs, i.e., the well is shut in at the sand face for build-up. These assumptions make it possible to describe two-phase flow of oil and gas by the partial differential equations:

Jan 1, 1967

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

There were present eight Directors, the Secretary and the Assistant Secretary of the Institute, and seventeen guests. Messrs. John V. N. Dorr and George D. Barron were appointed members of the Founder Societies Finance Committee on Joint Activities. A request from the New York members of seventeen societies that the Institute send a delegate to Engineering Congress for the purpose of formulating a code for professional conduct was referred to the Engineering Council. There were elected thirty-four members, twenty associates, and thirty-one junior associates. The status of one member was changed from associate to member; and of one junior associate to associate:

Jan 7, 1919

Edwin Ludlow, the forty-first President of the A. I. M. E., died in Muskogee, Okla., on Feb. 10, 1924, after a brief illness of influenza followed by pneumonia. He was born in Oakdale, N. Y. (on Long Island), March 12, 1858, a son of William Handy and Louise (Nicoll) Ludlow. He attended Flushing (L. I.) Institute 1868-75, and Columbia School of Mines (N. Y.) where he was graduated with the degree of Engineer of Mines with the Class of 1879. After graduation, his first employment was as assistant engineer in river and harbor work on the Delaware River near Philadelphia. In 1913

Jan 1, 1924

By S. W. Hazen

The use of operations research and computer techniques continued to expand during 1967. A number of mining companies either obtained computers or were experimenting with their use through service bureaus. This trend was also evident in the Canadian mining industry. Although the annual Symposium and Short Course on Computers and Operations Research in Mineral Industries was not held in 1967, technical periodicals like MINING ENGINEERING, Mining Congress Journal and Engineering and Mining Journal stepped up their coverage on the subject. The annual AIME meeting had eight papers in two sessions and one session on Operations Research was held at the SME Fall Meeting in Las Vegas.

Jan 2, 1968

By S. S. Nosseir

INTRODUCTION General Quebec Cartier Mining Co. (QCM) , Canadian subsidiary of U . S. Steel Corporation, operates two iron ore concentrators at the North Shore of Quebec, Canada. In its own operation at Mount-Wright, over 45 million tons per year of crude ore are ground, fully autogenous in six 32' x 11-1/2 Dominion Engineering mills. At Fire Lake/Lac Jeannine operation, property of Sidbec Normines Incorporated (SNI), over 11 million tons per year of crude ore are now ground, fully autogenous, in twelve 18' x 5' Hardinge mills. The Lac Jeannine concentrator was originally constructed by Q.C.M. in 1960 and. was since operated, in beneficiating the Lac Jeannine ore deposit until fully depleted in 1977, at an annual capacity of 18 million tons of crude per year. In 1976, the new company, Sidbec Normines, was incorporated by Sidbec Corporation, British Steel Corporation and Quebec Cartier Mining Co., in proportional holding of approximately 51 to 41 to 8 percent respectively. Since the acquisition of the Fire Lake ore deposit and Jeannine concentrator by the new company, QCM has continued managing the mine and concentrator operaions for S idbec Normines Inc . Flow-Sheet Design Jeannine concentrator comprises twelve identical mill lines grinding the ROM from 8" top size to 14 mesh nominal size range. Driven by 600 H.P. motor at 13.4 RPM corresponding to 73% critical speed, each mill is operated in a closed circuit with two 5' x 12' vibrating screens dressed with 2.5 x 25 mm slot opening cloth. The ground product or screen undersize at an average grade of 33.4% Fe is beneficiated by gravity separation to produce 66.5% Fe concentrate grade at 42.7% weight recovery. At the rougher stage, each mill line con-tains 128 Humphrey 13- 1/2" pitch spirals which produce a final tailing, middling and concentrate products. The rougher concentrate - so produced - is classified by 64 Humphrey pocket sizers into an underflow, which is a final concentrate product and an overflow which is further upgraded by 64 cleaner spirals. The cleaner concentrate together with the pocket sizer 's underflow are dewatered on a 13'[6] horizontal filter. The rougher middling together with the cleaner tailing are recycled in a closed circuit with the spirals. The rougher tailing after dewatering by hydrocyclones is disposed of by a

Jan 1, 1982

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing