Search Documents

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge's method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey's theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey's method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge's6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

By L. H. Reiss, J. Groult, L. Montadert

Many fields, where the reservoir is composed of sandy layers, show great complexity because of the lack of continuity which results from a particular type of seditnentation. This complexity may be a factor of prime importance in studying reservoir drainage, especially in the displacement of oil by water or gas. This paper describes work done to study this heterogeneity and its influence in the exploitation of the fields of the lower Devonian sandy shale in the Fort Polignac Basin in Algeria. Two lines of research were followed. 1. On the border of the hasin, a detailed sedimentological study established rules for the distribution of the sandy zones, which furnished hypotheses necessary for devising a geologic model of the reservoir. A certain number of sandy zones were recognized and characterized by their sedimentation. 2. In the interior of the basin, production logs deterrnined the productivity profile. Study of these point data led to a description of the distribution of oil and gas phases in the reservoir so that the continuity of different zones between the wells could be deduced. These techniques were applied to the study of the lower Devonian sandy shale reservoir of Zarzaitine. The sedimentological study guiding the correlation work, with the aid of the conventional static data (logs and cores), led to the establishment of a useful geologic model of the reservoir. Production logs, thanks to the contribution of dynamic elements, permitted corrobora-tion of the model. In particular, the logs confirmed separation, as proposed by the sedimenfological study, of two comparable sandy groups by revealing the continuity of a thin impermeable layer, invisible on conventional logs. INTRODUCTION The heterogeneity of sandstone reservoirs considerably influences the distribution of fluids, as well as their displacement in the reservoir. Rational development of a field and maximum recovery of hydrocarbons are closely tied to an exact and precise knowledge of the intern architecture. We will try to show how it is possible I refine the knowledge of the reservoir by seeking qui different and varied sources of information. Productic geology, the group of techniques which leads to idei concerning reservoir architecture and its geometric descris tion, furnishes the framework of this undertaking. Of these techniques. those used so far are qui varied. Descriptive as well as quantitative, the technique partake as much of classical geology as of reservo science. They include (1) geologic studies of the enti basin, especially a stratigraphic and sedimentologic study of the reservoir member at its outcrop, detail1 observations of cores taken from wells and, finall oriented cores; (2) interpretation of geophysical dat (3) interpretation of electric logs, nuclear and sonic lo and dipmeters: (4) petrophysical study of the reserve on samples taken from wells; (5) reservoir dynamics (stu of the production history and, in particular, interpretati of production logs); and (6) analysis of the fluids a thermodynamic studies. This incomplete list could also include (1) study the migration of the fluids, (2) geochemistry (chemi characterization of the oils), (3) regional hydrodynamii and (4) injection and detection of tracers. If most of these techniques are classic, attention i nevertheless be called to applying two of the techniql in conjunction: the study of the reservoir member its outcrop and the study of dynamic distribution fluid phases in the reservoir through use of product logs, used in the study of the sandy shale reservoirs the lower Devonian of the Fort Polignac Basin, Alge in particular, their application to the study of the 2 zaitine field where they have led to a new and satisfact model of the reservoir will be described. THE ZARZAITlNE FIELD Situated in the interior of the Fort Polignac 8 about 150 kilometers (km) north of the lower Devol outcrop, the Zarzaitine Devonian oil field was discovc On Ian. 28, 19'' (Figs. 1 and l8). It is in the fern a monocline truncated by two perpendicular fault: great throw, with a downthrown section to the north\: Erosion at the top of the structure has caused one of it to disappear. It is a reservoir with a gas cap an m coefficient on the order of 1:5, open to the aql On the northeast, with the oil zone covering an Of 106 sq km and with initial oil reserves of 222 mi

Jan 1, 1967

By E. R. Landry

This paper reviews the current state of electric furnace technology, with emphasis on some of the more recent developments. A large section of the paper is devoted to the use of increasingly high power inputs by arc furnace operators. Operating practices are then covered, including some of the newer techniques which have been tried.

Jan 1, 1967

By A. H. Larson, L. G. Twidwell

The influence which Au, Ag, Sb, Bi, Sn, and Cu have, both individually and collectively, on the activity coefficient of sulfur in liquid lead at 600"C zuas studied by circulating a H2S-Hz gas wlixture over a specific lead alloy until equilibrium was attained. Subsequently, the H2S concentration in the equilibrium gas mixture and sulfur concentration in the condensed phase were deterruined. The elements gold, silver, and antinzony (above 8 at. pct) increased the activity coefficient of sulfur. Bismuth had no apparent effect. Tin (above 3 at. pct) and copper decreased the coefficient. The influence of an individual element, i, on sulfur is best reported as the interaction parameter, riS, which is defined as The values o these first-order interaction zus are: ESzu = —55.0. These interaction parameters are used to predict the activity coefficient of sulfur in six fouv-component alloys and one seven-component alloy. Comparisons are made with direct experimental determinations. INTERACTIONS in dilute solution have been studied by many investigators. Most of the experimental work has been confined to solute-solvent interactions in simple binary systems and solute-solute interactions in ternary systems. Dealy and pehlke"~ have summarized the available literature on activity coefficients at infinite dilution in nonferrous binary alloys and have calculated from published data the values for interaction parameters in dilute nonferrous alloys. Interaction parameters are a convenient means of summarizing the effect of one solute species on another in a given solvent. Only a few investigators have studied interactions of the nonmetallic element sulfur in a metallic solvent. They are as follows: Rosenqvist,~ sulfur in silver; Rosenqvist and Cox,4 sulfur in steel; chipman, sulfur in alloy steels; Alcock and Richardson,% ulfur in copper alloys; Cheng and Alcock,' sulfur in iron, cobalt, and nickel; Cheng and ~lcock,' sulfur in lead and tin. The only reported work on the Pb-S system in the dilute-solution region is that of Cheng and Alcock.' Their investigation involved a study of the solubility of sulfur in liquid lead over the temperature range 500" to 680°C. The results may be summarized by the following relationship: S (dissolved in lead) + Pb(1) = PbS(s) log at. %S = -3388/T + 3.511 Experimentally, it was found that Henry's law was valid up to the solubility limit of sulfur in lead, i.e., at 600°C up to 0.43 pct. Their investigation did not include the study of sulfur in lead alloys. More accurate calculations could be made in smelting and refining systems if activity coefficients of solute species could be accurately predicted in complex solutions. One of the objectives of this study was to compare the experimental data with the values calculated from the equations derived from models for dilute solutions proposed by wagner9 and Alcock and Richardson. A temperature of 600°C was chosen as the experimental temperature to attain reasonable reaction rates and to minimize volatilization of the condensed phase. EXPERIMENTAL Materials. The Pb, Au, Ag, Sb, Bi, Sn, and Cu used for preparation of the alloys were American Smelting and Refining Co. research-grade materials. All were 99.999+ pct purity except the antimony and tin which were 99.99+ pct. The initial alloys prepared for this study consisted of twenty-one binary alloys, eleven ternary alloys, and one six-component alloy. The constituent elements were mixed for each desired alloy and were placed in a crucible machined from spectrographically pure graphite. The crucible was placed in a vycor tube which was evacuated with a vacuum pump and gettered by titanium sponge at 800°C for 8 to 12 hr. After the gettering was completed, the chamber containing the titanium was sealed and removed. The remaining sample chamber was placed in a tube furnace at 800°C for 2 hr and quenched in cold water. The final operation consisted of homogenization of the alloy for 1 to 2 weeks at a temperature just below the solidus for the individual system. The resulting master alloys were sectioned into small pieces and a random choice made for individual equilibrations. Cobalt sulfide (Cogs8) used to control the gas atmosphere in the circulation system was prepared by passing dried HzS for 24 hr over a Co-S mixture heated to 700°C in a tube furnace. This material was then mixed with cobalt metal to give a two-phase mixture which, when heated in hydrogen to a particular temperature, produced a desired H2S/H2 gas atmosphere in the circulation system. A Cu2S-Cu mixture also used in this study was prepared in a comparable manner. Apparatus for Equilibrium Measurements. The experimental technique of this study required apparatus

Jan 1, 1967

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By D. E. Wolfson

The Bureau of Mines carbonized 35 coal samples from Logan and Mingo Counties) W. Va., at 900° C using the Bureau of Mines-American Gas Association method, and determined yields of products and physical properties of cokes. Twenty samples represented the No.5 Block, Stockton, Chilton, Cedar Grove, Lower Cedar Grove, Alma, No. 2 Gas, and Eagle beds in Logan County, and 15 samples represented the Coalburg, Buffalo Creek, Winifrede, Cedar Grove, Lower Cedar Grove, Alma, and No.2 Gas beds in Mingo County, W. Va, Chemical analyses were determined for all coals. All coals in both counties rank as high-volatile A bituminous. Most of the coals are chemically suitable for metallurgical coals. Coals with accept-able ash and sulfur contents can be used as major constituents in coking blends with higher rank coal to yield coke with satisfactory properties for blast furnace use.

Jan 1, 1967

By H. O. Jahns

This paper describes the application of regression analysis for obtaining a two-dimensional areal description of heterogeneous reservoirs from short-term pressure-time data such as that obtained in interference tests. The method replaces the time-consuming trial-and-error procedure commonly used to match field data on an electric analyzer or digital computer with a systematic search which is programmed for a computer. The computer program adjusts the properties of a reservoir model automatically until a least-squares fit is obtained between observed and calculated pressure data. The reservoir is simulated by a single-phase, compressible, two-dimensional model. It is divided into a number of homogeneous blocks whose transmissibility (kh/F) and storage (Fch) values are varied to obtain the least-squares fit. The reliability of these values is determined from their standard deviations and correlation coefficients. Although the method is rigorously applicable to single-phase flow only, multiphase flow can be handled provided saturation changes are small during the test. Possibly the method can also be used to obtain a reservoir description from pressure-production history, but this application is outside the scope of this work. The paper includes, in addition to a description of the numerical procedure, a discussion of some of the problems associated with the method. Rules are given to help in selecting the number of homogeneous blocks and deciding upon their arrangement. The uniqueness of a reservoir description is considered. Finally, the use of the method is illustrated by the interpretation of field data from two interference tests. INTRODUCTION Pressure data from short-term transient tests, such as single-well and interference tests, are widely used to obtain reservoir properties. These tests are usually analyzed by assuming a simple reservoir model; very often, a homogeneous one is used. As a result, analysis of the transient data from each well frequently gives different values for reservoir properties. The problem then arises to combine all these differing results into a more detailed picture of the reservoir. One technique is to simulate the reservoir with a digital computer or with an electrical analyzer and to adjust the reservoir parameters by trial and error until the simulated pressure data are in reasonable agreement with the observed pressure data for all wells. Although this method has been used for both transient tests and pressure-history data, it is time-consuming and subjective. A second technique uses regression analysis to replace the trial-and-error procedure with a systematic search that can be programmed for a digital computer. Use of regression analysis in reservoir description was proposed recently by Jacquard and Jain.1 They divided the reservoir into a number of homogeneous blocks whose properties are varied until a least-squares fit is obtained between observed and calculated pressures. However, they did not consider their technique to be operational, mainly because of ".. .(I) the lack of experience in using the method .. . notably for the improvement of convergence; and (2) limitations imposed by the insufficiency of available computers". 1 While the analysis presented in this paper applies the same general principle used by Jacquard and Jain, the specific method is significantly different. Some differences are (1) the regression problem is solved in a different way which requires less computer time in most cases; (2) a stepwise solution, in which the detail in the reservoir description is increased from step to step, is used to improve convergence; and (3) the reliability of the estimated reservoir .properties.. as measured by their standard deviation and correlation coefficient, is estimated'

Jan 1, 1967

By V. O. Struebing, K. A. Johnson, F. H. Ellinger, C. C. Land

The system Pu-Ce-(Pu,Ce)Co2 has been investigated by thermal, micrographic, and X-ray diffraction methods. The principal characteristics of this system were found to be: 1) a "eutectic" valley extending across the phase diagram from the 6 Pu/Pu3 Co eutectic at 408°C to the 7 Ce/Ce24Co11 eutectic at 42Z°C, the maximum temperature in the valley being 443°C; 2) three ternary inter mediate phases, each having an appreciable mutual solid solubility for plutonium and cerium; 3) eleven four-phase equilibria comprising one eutectic, four peritectics, and six Class II reactions; 4) a pseudobinary section between the ternary phase A and the y Ce solid solution, and another such section between the 6 Pu solid solution and the 77 phase (Pu3 CO) containing cerium in solid solution. ALLOYS of plutonium that are completely molten at temperatures in the vicinity of 500°C are of interest for use in the liquid state as fast reactor fuels. Since exploratory work1,2 on Pu-Ce-Co alloys had indicated that a "eutectic" valley extended across the ternary system between the Pu-co[6 Pu/n (Pu3co)] and the Ce-CO[? Ce/? (ce24co11)] binary eutectics, thus allowing alloys of any plutonium concentration to be selected within the temperature limitation, the present investigation was begun for the purpose of accurately defining the "eutectic" valley. However, when the existence of three ternary intermediate phases and the consequent complexity of the liquid-solid equilibria were discovered, the investigation was expanded to cover the complete system, except for the corner between the isomorphous series PuCo2-CeCo2 and pure cobalt. 1) EXPERIMENTAL The purities of the alloy components are given in Table I. Alloy buttons, each weighing about 5 g, were prepared by arc melting the accurately weighed components (±0.1 mg) under a purified He-Ar atmosphere. Weight loss as a result of melting was usually less than 1 part in 1000; therefore, the buttons were not chemically analyzed. One half of each button was re-melted by induction heating in vacuo in a magnesia crucible and furnace-cooled, so as to produce a somewhat coarser microstructure than that of the rapidly cooled arc-melted button. About 20 g of each alloy to be examined by thermal analysis were remelted in vacuo in a tantalum crucible of the inverted thermocouple well type. The thermal, X-ray diffraction, heat treating, and micrographic methods used in this investigation were essentially the same as have been described in earlier reports.3,4 Two etchants for these alloys were found to be satisfactory: A) 8 (parts by volume) orthophosphoric acid, 5 ethanol, and 5 glycol; B) 7 tetraphosphoric acid, 36 water, and 57 2-ethoxyethanol. Etchant A was used at 5 to 6 v for 15 to 20 sec, and etchant B was used at 10 v for 10 sec. 2) THE LIMITING BINARY SYSTEMS 2.1) Plutonium-Cerium. The phase relationships of this system have been investigated by Ellinger, Land, and Cramer,5 and by Selle and Etter.6 The over-all agreement of their phase diagrams is good, but there are some discrepancies in the positions of some phase

Jan 1, 1967

By David S. Cahn

For many years the green and dry compressive strengths' of pellets have been, considered as possible criteria of pellet performance during induration. Although widely quoted throughout the industry, little attention has been given to standardizing the unfired pellet compressive strength test or to determining the characteristics of the test. The purpose-of this study was to compare the constant deformation rate test with that using a constant loading rate in addition, the dependence of these tests on such pellet properties as size, weight, and method of balling was to be investigated.

Jan 1, 1967

By Jack Smith

New turbine blades designed to resist ash erosion were tested by the Bureau of Mines in a pilot-scale coal-fired gas turbine. After a total of 1,963 hours, the rotor blades were only slightly eroded; and they had an estimated useful life of 20,000 to 30,000 hours. Changes in aerodynamic design of the blades successfully concentrated the ash at the outer casing, but stepped sidewalls for the rotor drum and wear strips at the bases of the rotor blades are still needed to minimize wear at these points. Stepped sidewalls and wear strips did not adequately protect the stator blades. Ash slightly eroded their leading edges, above the inserts, and cut notches into their trailing edges near the bases. Stator blade life was estimated at 5,000 to 7,500 hours. Blade erosion in these tests was less severe than with previous designs, but more efficient combustors and ash separators are required to reduce the amount and size of ash entering the turbine, to achieve acceptable blade life for commercial operation.

Jan 1, 1967

By G. L. Gates

The Bureau of Mines in cooperation with the City of Long Beach) Calif.) studied some physical factors effecting the deposition of barium sulfate scale in oil-well rotary pumps. Freshly formed barium sulfate was obtained by combining sulfate in sea water with barium from a barium chloride solution. The rate of scale deposition was measured quantitatively by weighing the pump impellers and diffusers. The results of the tests showed that coating the impellers with plastic significantly reduced the rate of barium sulfate deposition.

Jan 1, 1967

By Arnulf Muan, Klaus Schwerdtfeger, L. S. Darken

The equilibrium ratios CO/CO, of. a gas pharse coexislitrg with selected oxide phase assemblages of the system "FeO"-MnO-SiO2 and metallic iron have been determined at 1150°C. The data obtained are combined wilh the thermodynamic data on fayalite (Fe2SiOi) and "FeO"-MnO solid solutions to derive the actiuity-coiHpQsiUon relations for the various oxide components of the olivine (Fe2SiO4-,lIn2Si04) and pyroxenoid (standard free energies of formation of ferrosilile (FeSiO,), tephroite , and rhodonite from their oxides have been evaluated. In the Appendicc~s by Darken and Schivcerdtjeger the ideal solution law for solid deri oed. THE quaternary system Fe-Mn-Si-0 is of interest in many important phases of iron- and steelmaking. For instance, reduction in the blast furnace of ores containing Fe-Mn silicates, deoxidation of liquid steel with Mn-Si alloys, formation and subsequent alterations of oxide inclusions in steels, and the oxidation of manganese- and silicon-containing steels during processing are all governed largely by the thermodynamic relations existing in the system Fe-Mn-Si-O. Although the equilibria between liquid iron and Fe-Mn silicates have been studied in detail previously,' experimental data are lacking for equilibria among the solid phases. The present paper provides such data for a selected temperature, 1150°C. I) PREVIOUS WORK The general feature of subsolidus phase relations in the system iron oxide-manganese oxide-silica under strongly reducing conditions is known.' Fig. 1 represents an isothermal section at 1150°C. The phases present, in addition to silica, are three solid-solution series, viz.. the oxide solid-solution series "Fe0"-MnO (manganowiistite), the orthosilicate solid-solution series Fe2SiO4-Mn2SiO4 (olivine), and the metasilicate solid-solution series FeSi03-MnSi03 (pyroxenoid). The "FeO"-MnO solid solution is continuous between the two end members wustite and manganositeV3 as is the orthosilicate series between the two end members Fe2SiO4 and Mn2 SiO4. 4 The metasilicate solid solution, on the other hand, extends only part way from the MnSi0.l end member toward FeSiO because the latter end member is unstable relative to the phase assemblage fayalite + silica. The letters a through m in Fig. 1 identify end members of the individual solid phases dealt with [for instance, c designates fayalite (FeZSiO4) in equilibrium with wustite, e designates fayalite in equilibrium with silica, and so forth]. In addition to the phase relations there are some thermodynamic data available for the end members of the silicate solid solutions and for the "FeO" -MnO solid solution. The free energy of formation of fayalite (Fe2Si04) from its oxide components has been studied by several investigators. Unfortunately, the data obtained show considerable discrepancies. According to recent tabulations by Kelley, the change in standard free energy Affn for the reaction

Jan 1, 1967

By L. M. Kaas

Production planning techniques, particularly a production scheduling model are proposed for use in taconite mining in Minnesota. Applications of the scheduling model, on a two-year mine plan, which evaluates the equipment involved, are detailed. Data and tables are included. INTRODUCTION Fifty years ago Walter Swart, a mining engineer investigating the Minnesota taconites for the Jackling interests, made this report: "There is no doubt in my mind that this taconite mining is a game worth playing and that it can be played to a successful finish. It seems now only a question of men, money, methods and details."' The prosperous industry envisioned by Swart and other taconite pioneers is today a reality. Its growth has been spurred by technological advances in mining and beneficiation, depletion of natural ore reserves, and improvement of the economic climate. However, the massive problems that were solved to make this development possible are only a foreshadowing of the questions of "men, money, methods and details" that will be encountered in the future. The increasing demands for higher quality taconite pellets have already forced producers to examine better methods of quality control for their operations. THE NEED FOR PRODUCTION SCHEDULING The maintenance of fairly constant crude ore characteristics is of major importance to mill quality control. Through careful mixing of ore from storage bins, the concentrator staff can partially reduce the deviations of the feed analysis from the analysis desired. However, the limited effectiveness of bin mixing systems has caused some mining companies to utilize large capacity blending piles to reduce grade variations. Although the analysis of material being reclaimed from such a pile is relatively uniform, raw materials still must be scheduled from the mine into the pile in such a manner that the average analysis of all material is close to the feed analysis prescribed by the mill. PRODUCTION SCHEDULING MODEL Minnesota taconite mining companies have used three analyses to estimate how a certain ore will react to concentration. These factors are the crude magnetic iron content of the ore, the iron content of the plant concentrate produced from the ore at a standard grind, and the grindability index of the ore. The actual laboratory tests used to determine these factors have been described by Apuli.2 These data are obtained from analysis of diamond drill cores, and are the established criteria used by mine planning engineers to formulate taconite production schedules. Variations in the grade of mine output should be controlled by any production scheduling system. However, to be truly practical, production scheduling methods must also consider many of the less tangible variables that affect mine operation. These include pit layout, available equipment, and ownership of ore leases. The number and scope of the variables inherent in this type of mine planning suggested that a computer simulation could be used to solve the problem. The preparation of information for use with a computer simulation program was the initial phase of this project. Meaningful data was developed from logs of diamond drill cores. Analyses were related to their respective areas of influence by means of mine maps, Fig. 1. The three assay variables and tonnages for each area were punched on cards. In addition to these values, a location code number was punched in the first eight columns of the cards. This code is explained as follows: Column 1 Bench number Column 2 Quadrant location as determined by mine reference lines Columns 3-5 East-west coordinate Columns 6-8 North-south coordinate The heart of any simulation techniques is a well engineered mathematical model. To develop a model for taconite production scheduling, the relationship between the three assay variables had to be established. Because changes in each analysis factor affect the mill circuits in different manners, a weighting system was used to determine the relative im-

Jan 1, 1967

By Marius S. Darrow, William B. White, Rustum Roy

The PbS-PbTe systen has been studied by quench-ing and D.T.A. techniques f?om 400' to 1150°C. Runs were made in evacuated silica tubes so that all equilibria are at the vapor pressure of the system. Lattice parameters of the quenched salnples , measured by X-ray diffraction, show a complete crystalline-solution series existing over a narrow temperature range between approximately 805" and 871°C. An exsolution dome extends from a maximum of about 805"C (approximately 30 mole pct PbTe) to 1 and 96.5 pet PbTe at 400°C. A narrow melting region, deternined by D.T.A., extends form 918c (mp PbTe), The shapes of the liquides and solidus curves imply the existence of a minimum at 871°C at approximately 65 pct PbTe. THe exact composition of the minimum could not be established due to the very narrow two-phase region. At compositions containing less than 50 pet PbTe, liquidus temperatures begin to increase, while the solidus remains almost flat to about 15 mole pet PbTe before beginning to vise toward the mp of PbS (1075 C). LEAD sulfide and lead telluride are isostructural (NaC1 type) semiconductors whose electrical and optical properties have been extensively studied and used in recent years. If appreciable crystalline solution exists between these compounds, the variation of physical properties with composition could be of interest. The purpose of this investigation was to determine the extent, if any. of crystalline solution, and to obtain the phase diagram for the system. To the knowledge of the authors, only three studies of the system PbS-PbTe have been reported, and, in chronological order, each investigation found an increasing amount of crystalline solution. In 1956, Yamamoto reported finding no evidence of crystalline solution between the compounds. Sindeyeva and Godov-ikov,' in 1959, found very limited crystalline solution. but only under conditions of excess tellurium concentration. Finally Melevski s3 investigation in 1963 indicated that one solid phase exists in the region from PbS to 7 pct PbTe and from 82 pct PbTe to PbTe at 886'C, with an eutectic at 55 pct PbTe at that temperature. Detailed data on the solvus boundary were not given. EXPERIMENTAL EQUIPMENT AND MATERIALS Commercially produced PbTe and PbS powders were used as starting materials. Batches of specific mole percent composition were accurately weighed and mixed in a plastic bottle, in a shaker mill. An analy- sis of impurity content is given in Table I for pure PbS and PbTe and for two randomly selected batches after the powders were mixed. Individual samples, ranging in weight from 0.2 to 0.5 g, were sealed in evacuated silica tubes which had been thoroughly washed and rinsed with acetone and distilled water. Thus all data taken were at the pressure of the system. Subsolidus relations were studied down to 400°C by heating the samples in a vertical tube furnace for 24 hr. The sealed tubes were quenched in water with quench time from the hot zone not exceeding 1 sec. Temperatures were measured by a chromel-alumel thermocouple and controlled to 53°C for most runs. The number and composition of phases present were determined from powder X-ray diffraction patterns taken at room temperature on a Norelco diffractome-ter, using silicon as an external standard. Above 850°C quenching techniques were, in general, found to be unsatisfactory, and differential thermal analysis (D.T.A.) was used to determine melting relations. The evacuated tubes were recessed about 1 cm at one end to accommodate the differential thermocouple. Al203 was used as the reference material in a similar tube containing the other side of the differential couple. For temperature measurements, a separate thermocouple was placed in the recess of the tube containing the sample to be measured, thus providing an opportunity to obtain thermal, as well as differential, analysis. All thermocouples for these measurements were Pt-Pt 10 pct Rh. Temperature and differential curves were recorded separately on synchronized strip-chart recorders. Thermocouples and recording equipment were calibrated using NaCl and gold standards, using the melting points 801" and 1063 C, respectively, which span most of the temperature range of interest. Heating and cooling rates generally were from 4 to 7°C per min. It was found, in fact. that rates ranging from 1.5 to 25°C per min did not significantly change the data obtained.

Jan 1, 1967

By H. F. Haller

Large volume sand plugs are placed in tunnels to stem underground explosions at the Atomic Energy Commision's Nevada Test Site. Placement through long, narrow accessways is like stemming coyote holes and blast chambers in mining and construction. Experience with several other methods showed that pipeline conveying could be preferable. Small scale tests justified predicting a placement rate of 20 y. p. h. at 1000 feet, with a modified gunite placer. The initial production trial placed over 2000 c. y. and was estimated to effect a small cost saving after paying for new equipment. Initial problems included duct replacement and supply delays, dust, and other personnel hazards. Operations analysis and SOP'S reduced delays during the second field trial, minimized hazards, doubled production. Re cent re search and operations analyses produced data on abrasion resistant ducts, identified minor equipment modifications which can improve production, developed techniques for impact compaction to 90% ?maximum? and for extension of the practical range to 2000 feet. Twenty-two illustrations include an air consumption-distance chart.

Jan 1, 1967

By R. F. Hehemann, S. T. Zegler

The w phase in zirconium alloys containing more than G pct Nb can form in a difjUsionless manner during quenching or with composition change during aging at temperatures below 550°C. The latter treatment establishes a metastable equilibrium between an w phase containing 5 pct Nb and a ß phase containing 45 pct Nb. The superconducting transition temperat~ire of quenched Z?-10 pct Nb and Zr-15 pct Nb alloys is increased by aging at 375°C and both alloys exhibit the same Tc after prolonged aging. The superconducting transition both in the quenched and in the aged alloys is relatively sharp (occurring ovev a temperature range of <0.5°K), and is interpreted in terms of the microstructure that characterizes the ß-w) mixture and proximity effects postulated for thin-film superconductors. At elevated temperatures alloys of the transition elements from groups 4 to 6 in the periodic table exhibit wide ranges of mutual solubility in the bcc structure. Except for compositions that are high in a group 4 element? the bcc structure can be retained at room temperature by quenching and several systematic studies of the influence of composition on superconductivity in bcc alloys have been reported.&apos;&apos;2 The high-temperature bcc form (ß) of the group 4 elements—titanium, zirconium, hafnium—is not retained upon quenching but instead transforms marten-sitically to the hcp structure (a). The Ms temperature can be lowered by alloying the group 4 metals with a solute from a higher group and with a critical amount of solute the martensitic transformation can be suppressed. The 3 phase, however, is not retained without modification. Rather, a metastable phase, w (in particular in titanium- and zirconium-base alloys), occurs over a fairly wide range of compositions above that required for suppression of the martensitic transformation.3,4 As in the formation of a, the transformation of ß to w occurs in a diffusionless manner during quenching and the temperature for its initiation (ws) is lowered by increasing solute content. The ws temperature can be well below room temperature and reversibility of the transformation allows w that forms on cooling to low temperatures to disappear again when the sample is returned to room temperature.5 In Zr-Nb alloys the w phase can also be formed during aging at temperatures below 550°c.5 In this treatment, a metastable equilibrium is established between a niobium-rich ß phase and a zirconium-rich w phase. The present study of Zr-Nb alloys was undertaken to examine the influence of microstructure: as controlled by composition and heat treatment, on the superconducting transition temperature, Tc. MATERIALS AND PROCEDURE The alloys were prepared from zirconium crystal bar that was spectrochemically pure except for 30 ppm of Al, 40 ppm of Cu, 100 ppm of Fe, and 1 ppm of each B and Pb. The niobium used was 99.5+ pct pure containing according to the suppliers&apos; (Stauffer Chemical Co.) analysis a maximum of 50 ppmO, 60 ppm of N, 30 ppm of C, 5 ppm of H, 200 ppm of Ta, 50 ppm of Mo, and 15 ppm each of Fe, Co, and Ni. The alloys were prepared in the form of buttons weighing approximately 4 g by arc melting on a water-cooled copper hearth in an inert-gas atmosphere. All alloys were heat-treated in the ß-phase region for 3 days at 1000°C and water-quenched. For heat treatment, the specimens, 0.5 cm in diam by 1.0 cm long, were wrapped in zirconium foil and sealed in an inert-gas atmosphere inside quartz capsules. Quenching was done by breaking the capsules under water. Subsequent aging treatments of Zr-10 pct Nb and Zr-15 pct Nb alloys were done by placing the specimens in a lead bath at 375°C and quenching in chilled brine. The phases present in the quenched and aged alloys were identified by optical and X-ray metallography. X-ray powder patterns were taken at room temperature with a 114.6-mm-diam Debye-Scherrer camera and CuK-a radiation. Superconducting transition temperatures were determined from magnetic-permeability measurements6 in a 10 oe field. The superconducting transition temperature was determined from a plot of galvanometric deflection as a function of temperature by extrapolation of the linear portion of the curve to zero deflection. In Figs. 1 and 3 the arrows indicate the temperature range of the transitions. RESULTS AND DISCUSSION Superconductivity in Quenched Alloys. In Fig. 1 Tc for alloys quenched from 1000°C is plotted against nominal composition. Indicated also in Fig. 1 are results of metallographic studies. Alloys containing less than 5 pct Nb transform martensitically to a during quenching from 1000°C and Fig. 1 demonstrates that Tc increases significantly as the niobium concentration in a is increased. On the basis of the Bardeen-Cooper-Schreiffer theory of superconductivity7 it seems most likely that this increase in Tc results from an increase in the density of electronic states as zirconium is alloyed with niobium. This increased density of states has been revealed clearly in low-temperature specific-heat measurements.8 The martensitic transformation to a is suppressed when the alloy content exceeds 5 pct Nb. In the range from approximately 6 to 30 pct Nb, quenched alloys consist of a mixture of B and w phases and the amount

Jan 1, 1967

By Adam Chrzanowski, Eugene E. Derenyi, Peter Wilson

"The rapid growth of the mining industry in Canada necessitates a reevaluation of existing mining surveying techniques and the development of new methods. In this connection, the Department of Surveying Engineering at the University of New Brunswick initiated a research program, and the major points of investigation conducted to date are as follows: 1. Based on practical experience and theoretical studies, all factors influencing classical shaft-plumbing procedures were analyzed. Refined methods were developed to increase accuracy. 2. The performance of electronic and optical distance-measuring instruments in underground conditions were investigated. In this connection, test measurements were made in various mines, employing four types of Tellurometers, two models of Geodimeters and a proto-type optical range-finder. 3. Tests were conducted and methods developed for the mapping of large underground chambers by photo-grammetric means. The aforementioned points are discussed at some length in the paper. The authors have also indicated areas where further research is needed and suggested ways and means to improve and modernize mining surveying as a whole throughout Canada."

Jan 1, 1967

By J. N. Maycock, P. W. M. Jacobs

Previous measurements of the thermoelectric power of ionic crystals are reviewed briefly. It is concluded that, while extensive measurerements are available on systems in which the electrode M has a cation in cornon with the sall MX, there is as yet little modern work on systems in which the metal electrode M&apos; does not haue a calion in common with the salt. We haue now measured the thermoelectric power of- single cryslals of KC1 and of KC1 + SrCl2 between plati.nu?n electrodes over a much wider range of temperature than before. individual ionic heats of transport, which, while probably not very accurate, are of a reasonable order of magnitude, are obtained from this data on the basis of an assumption that charge carriers do not cross the electrode/salt interface. RELATIVELY little work has been done on thermal diffusion in ionic crystals. Reinhold and coworkers1-I measured the thermoelectric power and Soret effect in silver and cuprous salts and their mixed crystals. One electrode was kept at constant temperature and electromotive force was measured as a function of the temperature gradient. This is not a satisfactory technique if the results are to be compared with theory; rather should the electromotive force be measured for small temperature differences as a function of the mean temperature of the two electrodes. They used electrodes having the same cation as the salt. Patrick and Lawsons measured the thermoelectric power of AgBr and of AgBr containing 0.25 mole pct CdBr2 using silver electrodes and temperature gradients of 10°C in the temperature range 100° to 400°C. The results were analyzed using the Teltow9 values of the ratio of the mobility of the two types of charge carrier (cation vacancies and interstitials) found from conductance measurements and the wagnerlo estimate of the heterogeneous contribution to the thermopower, Ohet They obtained values of qi = 0.017 ev and q$ = -0.385 ev for the heats of transport of interstitial Ag ions and cation vacancies, respectively. Their data has been reanalyzed by Haga using the corrections to the Wagner formula for ehet given by Howard andLidiard; Haga found ev. Christy et a1.14 have measured the thermopower of AgBr containing 0.01 to 1.0 mole pct CdBr,. They analyzed their results in terms of the theory of Howard and Lidiard and found that excellent agreement with the theory was obtained if q5 + : + hF was set equal to 0.70 ev at 175°C and equal to 0.58 ev at 250°C. Since Teltow&apos;s value for hF, the enthalpy of formation of a Frenkel defect pair, is 1.27 ev, this means that qt + q: increases in magnitude from -0.57 ev at 175°C to -0.69 ev at 250°C. This result is in poor agreement with Patrick and Lawson&apos;s value of -0.37 ev for pure AgBr (for which the corrections to the Wagner estimate for ehet are very small). christy1&apos; has also reported thermopower measurements for pure silver halides and for the systems he finds that qr + qt for AgBr increases in magnitude from about -0.55 ev at 150°C to about -1.1 ev at 400°C. The temperature dependence of the sum of the heats of transport in AgCl is very similar in magnitude. The sum qi + q cannot be decomposed into the individual interstitial and vacancy heats of transport without a knowledge of the entropy contribution from the change in the lattice vibrational frequencies around the interstitial and vacancy. Despite this it seems an inescapable conclusion from Christy&apos;s analysis that and -6 = Q Ag are quite strongly temperature-de pendent and that -q: is much greater in magnitude than the corresponding activation energy for vacancy migration, +0.37 ev. It remains a mystery as to why the data for pure AgBr are fitted so well by the assumption of temperature-independent heats of transport if association of impurity ions and vacancies is neglected.l The values for q$ and q: so obtained are also much more in line with those expected from a simple kinetic interpretation of the heats of transport due to wirtz.17 Christy and Hsueh have also measured the thermopower of CuCl between copper electrodes; 0 for pure CuCl is comparable with that for the silver halides but, owing to the limited solubility of divalent cation impurities it was not found possible to deduce values for the heats of transport. In contrast to these extensive measurements on cells of the type M MX 1 M, there have been relatively few attempts to measure the thermoelectric power of cells of the form M&apos; 1 MX I M&apos;. (We specifically exclude from discussion cells of this type in which the metal electrodes M&apos; are in contact with gas containing a controlled partial pressure of either M or X,, since in these cells the electrochemical reactions are analogous to those of the first type.) The thermoelectric power of Pb12, PbC12, TlC1, TlBr, TlI, NaNO3, and KNO 3 between platinum electrodes was measured by Thie1e.l3 In discussing these measurements Holtan20 asserts that the thermopotential initially fluctuated wildly and even changed sign before approaching a steady value. A high degree of irreproducibility is also evident in the work of Nikitinskaya and Murin21 who first studied the thermopower of NaCl and KC1 between platinum electrodes. Much better reproduci-bility was achieved in our earlier measurements on KC1;22 these were, however, confined to a relatively narrow temperature range and one of the objects of the work reported in this paper was to extend the temperature range of the measurements of the thermoelectric power 8 of KC1. The only other published

Jan 1, 1967

By Seymour G. Epstein

Using a modified capillary-reservoir technique, electromigration of dilute solutions of silver in liquid bismuth has been measured with relatively good precision. The effects of experiment duration, temperature. current density, and silver concentration have been investigated. Chemical diffusion coefficients of silver in liquid bismuth were also measured at temperatures ranging from 300°to 600°C. The effective valence of silver in bismuth was calculated as a function of temperature from a comparison of electric mobilities and diffusion coefficients. In an alloy containing 1 wt pct Ag, the electric mobility of silver increased from 1.65 x X sq cm per v-sec at 300°C to 2.38 x 10-4at 600°C; the diffusion coefficient in-creased from 3.0 to 17.5 x 10-5 sq cm per sec over the same temperature range. Calculated values of effective valence decreased from 0.28 at 300°C to 0.10 at 600°c. ALTHOUGH there has been considerable progress in the theoretical treatment of electron-transport properties in liquid metals, the mechanisms of mass-transport phenomena in liquid metals are still not clearly understood. This is especially true for the case of electro migration, alternatively termed elec-trotransport or electrodiffusion, which is mass transport induced by an applied electric field. A thorough knowledge of the mechanisms governing electromigration will considerably aid in understanding the electronic state of an impurity atom in a liquid metal. Electromigration has been long observed and studied in liquid metals and alloys. However, no satisfactory theory to explain the observed effects has been developed. The disparity in the present theoretical treatments, summarized by Verhoeven,1 has been attributed to the difficulty in describing the momentum transfer between electrons and ions and to the lack of precise measurements in liquid metals, especially at elevated temperatures. Prior to the last decade most of the electromigration studies were qualitative, in that only the direction of migration of alloy constituents was determined with any degree of reliability. These results have been summarized by schwarz,2 who was able to measure the electric mobility of several elements in mercury. More recently electromigration in liquid metals has been measured by a variety of techniques, with a varying degree of success. Mangelsdorf3,4 made precise measurements on dilute amalgams at room temperature, using a technique relating changes in composition to changes in alloy resistivity, but was unable to make measurements at elevated temperatures. Using a modification of this technique, Verhoeven and Hucke5 measured electromigration and resistivities in liquid Bi-Sn alloys. However, the resistivity technique is not sensitive to small composition changes and is impractical for very dilute alloys. Verhoeven and Hucke6 also measured electromigration of several solutes in liquid bismuth by adapting the capillary-reservoir method for determining diffusion coefficients in liquids. The reported precision of these measurements was poor, however, primarily due to the large uncertainties in the chemical analyses to determine concentration changes. Belashchenko7 has summarized the recent Russian investigations of electromigration, many of which, however, are of questionable accuracy. The present study was initiated to develop an experimental technique which would be both versatile and sufficiently sensitive to yield precise measurements of electromigration of solutes in dilute metal alloys at elevated temperatures. Using this technique, a comprehensive study of silver in liquid bismuth was made to determine the effects of experimental variables on the apparent electric mobilities. Chemical diffusion coefficients of silver in liquid bismuth were also measured over the same temperature range for comparison with the electromigration data. EXPERIMENTAL TECHNIQUE Electromigration Measurements. Electromigration of silver in liquid bismuth was measured by the capillary-reservoir technique, similar to that used by Verhoeven and Hucke,6 but modified to yield more

Jan 1, 1967

By Earl G. Alexander, W. R. Runyan

The nature of silicon growths which form on amorphous silicon dioxide during epitaxial deposition in unmasked areas was investigated. Octahedral silicon crystals nucleate and grow on the oxide at pinholes, adjacent to foreign matter, and probably at oxygen-deficient sites. Generation of these oxygen-deficient sites on the oxide surface is shown to be thermody-namically favorable under usual deposition conditions. The conditions which affect saturation (SiCl4 concentration, flow rate, temperature) determine the number of sites on which silicon will nucleate and grow. Reaction product HCl is likely responsible for suppressing silicon growths on the oxide near exposed areas of silicon. Evidence is given that the oxide-silicon interaction occurs in hydrogen only and involves water as a transport agent. SINCE shortly after the first epitaxial layers of silicon and germanium were grown over large areas, efforts have been underway to grow selectively on specified portions of the surface. Such selective growth should be particularly applicable to monolithic integrated circuits. Usually the approach has been to use a mask to delineate the desired area,&apos;-7 but the VLS methoda can also be used, and direct programmed atomic beam deposition has been suggested.9 The masking techniques may embody two forms: in one the mask is a material which can later be removed along with any growths attached to it; in the other, attempts are made to prevent growth on the mask, so that over a whole slice the only deposition that will occur is in the unmasked regions. The present work is concerned with only the latter case, and specifically with silicon. Joyce and Baldrey1 suggested the use of thermally grown silicon dioxide as a mask material and showed that for the proper choice of operating conditions (total flow, SiCl4 concentration, total deposition time, and temperature) spurious growth could be minimized and that if the deposition times were not too large, and the criteria not too stringent, a usable process occurred. Later Schnable et al 4 observed that after extended deposition times, during which most of the silicon dioxide became covered with growth, there was always a band of oxide around the exposed silicon area that was quite free of spurious silicon crystallites. This latter observation is particularly intriguing since it extends the hope that with proper precautions masking can be much more versatile than it presently is. EXPERIMENTAL PROCEDURE Silicon substrates oriented to within 1.5 deg of a (111) plane were prepared by etch polishing lapped

Jan 1, 1967