Search Documents

By Arnold H. Kackman

At the Kaiser Cement and Gypsum plant near San Jose, Calif., one rod-ball mill has taken over the entire raw grinding function for the largest single cement operation in the West. Installed as part of a major overhaul performed on the plant by Kaiser Engineers, the new mill does the work of seven previous mills in feeding the plant's six kilns and keeping pace with its annual capacity of 8.5 million bbl of cement. Installation of the mill, which is the largest in the world, was started in December of 1966 and completed in April of 1967. The over-all modernization program surrounding the new mill has included the modification of existing crushing facilities; the building of new raw materials handling and blending facilities, new raw milling facilities, and new slurry handling facilities; and the installation of a digital-analog control system.

Jan 7, 1967

By B. H. Caudle, R. L. Slobod, E. R. Brownscombe

An earlier publication has discussed three methods lor obtaining relative permeability data on small core sample-: and the apparatus and technique for the capillary pressure displacement method. This paper describes the apparatus and technique for the solution gas and the dynamic displacement methods and presents a routine procedure for obtaining oil-gas and water-oil relative permeability data. Theoretical and experimental considerations are presented to show that the end effect commonly associated with the dynamic flow mechanism is extremely small where constant rates are employed in the flowing phase. An effect of flow rate on relative permeabilities obtained by the dynamic system is found only when gas is one of the flowing phases and this effect is ascribed to a form of channeling in the capillary system. The apparatus and procedure* used to obtain relative permeability data with three phases flowing are described and some preliminary results of the use of this method are shown. INTRODUCTION In a previous publication from this laboratory there appeared a brief discussion of the concepts behind three basic methods for obtaining relative permeability data.' These three methods were called the capillary pressure displacement method. the solution gas displacement method, and the dy. namic displacement method — the names being suggestive of the type of process used for obtaining the desired saturation prior to making the permeability measurements. In the same publication, the apparatus and technique for the capillary pressure displacement method were described and some typical results obtained by this method were presented. This paper will present: 1. The routine procedure used to obtain complete relative permeabilitv data (In small core samples with two phasec flowing; 2. The experimental techniques for the solution gas and dynamic displacement methods for obtaining relative permeability data; 3. The results of some studies on the mechanism of fluid flow through consolidated porous media; and, 4. The preliminary results on the determination of relative permeability for a system in which three phases are flowing. TECHNIQUE AND APPARATUS FOR THE DETERMINATION OF RELATIVE PERMEABILITY WITH TWO PHASES FLOWING In order to obtain relative permeability data on a small core sample which will describe the various phases of reservoir depletion. it is necessary that the methods used duplicate, as closely as possible, the various mechanisms involved in rescbrvoir depletion. It is important that. in any type of relative permedbility determination. the core be originally at virgin reservoir conditions with respect to the fluids within the pore *paces. (for example. oil used to displace water from core until connate water saturation estahlished). and that the subsequent saturation changes in the core be analagous to tile reservoir depletion mechanism for which the data are to he used. PROCEDURE OF DETERMINING RELATIVE PERMEABILITY DATA WITH TWO PHASES FLOWING A standardized procedure has been set up for obtaining relative permeability data and other information on small core samples. This procedure essentially consists of 11 steps as shown below: 1. Determination of tile capillary pressure curve.'! This information is used to help select samples for relative permeability determinations and to charactrrize. in general. the pore size distribution of the sample. 2. Determination of the specific permeabilities to air, oil. and water. These values are those on which the relative pern~eahility calculations are based. 3. Reduction of the core to connate water by an oil drive using the high speed centrifuge. 1. Determination of the critical gas saturation value and the associated relative permeability to oil by the solution gas drive method. 5. Determination of the relative permeabilities to oil and gab in the region where both phases are flowing by the capillary presqure method. 6. Reduction of the sample to residual oil in the centrifuge by an air drive. and the determination of the relative permeability to gas at this point.

Jan 1, 1951

By F. W. Jessen, R. W. Durie

The development in the laboratory of mall salt cavities illustrated that the rate at which salt dissolved was affected by the rate at which fresh water was injected into the cavity. Experiments carried out with small salt cores indicated, however, that the rate of salt removal with forced convection was not significantly different from that with natural convection, provided the flow was in the laminar range. The rate of salt removal was apparently controlled by the water salinity, which was indirectly determined by the rate at which fresh water was injected. With turbulent flow at the salt surface, the rate of salt removal could be increased several-fold with the same salinity in the cavity. A quantitative evaluation of the rate of salt removal under natural convection was attempted by using techniques developed for the analogous heat convection system. The values of the rate of salt removal obtained experimentally were considerably greater than the calculated values due to surface irregularities. To determine the effect of these irregularities, the rate of salt removal was found for a smooth surface by experimentally determining the rate of salt removal as an irregular surface developed and then extrapolating this trend back to initial time. Excellent agreement between the experimental and calculated values then were obtained. INTRODUCTION Underground cavities leached from massive salt formations have been used extensively for many years for storing LPG products. The development techniques for storing methane gases in the liquid state at a low temperature leads to speculation that underground salt cavities of controlled geometry may be quite important in this field. In any case, the increasing number of salt cavities used for storage indicates an impending need for a better knowledge of the washing process to efficiently develop maximum storage volumes within specified regions. It is the purpose of this paper to discuss some aspects of -the mechanism of the dissolution of salt in the formation of cavities of controlled geometry. Although no attempt is made to resolve the problems that may be associated with the storage of methane in the liquid state in salt formations, the fact that salt cavities can be formed at reasonable cost would appear to justify investigation in this direction. Further, underground cavities have the advantage of allowing storage at higher pressures (and hence higher temperatures) than surface installations. Although the thermal conductivity of salt is somewhat higher than frozen earths, 0.017 cal/(sec) (cm)(c) for salt at 32F vs 0.005 cal/(sec) (cm) () for ice at this temperature, this disadvantage is offset by the higher storage temperature in underground cavities. For example, at atmospheric pressure, the temperature of liquid methane would be —258F. In a storage cavity at 2,000 ft, the pressure would be approximately 1,000 psia, which would allow a storage temperature of -9OF. The formation temperature would be roughly 110F at this depth, hence the total temperature differential would be 200F. Under surface storage conditions the corresponding temperature difference would be some 333 F (75 F + 258 F). The results of this difference can be compared more readily by substituting the appropriate values of thermal conductivity, temperature difference, heat capacity and density into the solution of the diffusion equation for the unsteady linear conduction of heat. For the underground salt storage, the heat flux is 15cal/(sec) (sq cm), where t is the time in seconds from the initiation of the cooling. For the surface system, and using the thermal properties of ice at the average temperature, the corresponding figure is 21cal/(sec) (sq cm). Hence, the heat losses through ice under surface conditions would be some 30 per cent higher than the underground salt storage. This ratio would vary, considerablv. Of course, for other conditions, but it indicates that an advantage over surface facilities may be obtained by storage in underground salt cavities.

Jan 1, 1965

By Margaret Dienes, Karl L. Fetters

IT is recognized that ingot structure is important in determining the quality of finished steel. Such elements of ingot structure as the size, shape and distribution of primary crystals; the size, shape and distribution of shrinkage cavities; the degree and kind of segregation (both macro and micro); the distribution of nonmetallic impurities and gas cavities are each the result of solidification phenomena. These structural features vary with variations in ingot practice: Their control requires a knowledge of their correlation with the variables of ingot practice. The structures of full-sized commercial steel ingots have been studied extensively by the Steel Ingots Committee of the Iron and Steel Institute.1 This work has contributed much to our knowledge of ingot structure and the factors that control it. These investigations of necessity have been limited in scope because of their costliness and the great amount of time required. Studies of the structure of small steel ingots have not been very satisfactory because of the rapid rate of freezing of the ingots. Their structures are thus difficult to resolve and to interpret. Medium sizes of steel ingots in the ranges of 150 to 500 lb. do offer opportunity for study, but their cost also becomes greater and their use less expeditious. Experimental models for the study of steel-ingot structure have been made of nonferrous metals. Gathmann2 used lead-tin alloys and reported that they were satisfactory for determination of the location of the major shrinkage cavities and the "planes of weakness" resulting from the pattern of crystallization. Certain features of mold design have been illustrated by these results. Bezdenezhnihk3 has used copper and aluminum ingots similarly and has shown the effect of mold-wall thickness on the depth of columnar grains. These studies of nonferrous models have been useful, but their results have been limited because the analogy with steel ingots obviously is not a close one. Model wax ingots have also been useful in pointing out important features of mold design. Gathmann4 watched the solidification of wax ingots, through a glass-walled mold, and has reported the effects of mold type on the location of the shrinkage cavities. Segregation was simulated by using a mixture of paraffin and stearine and adding a red dye that on freezing segregated preferentially to the stearine. Brearley5 also used wax models to determine .the effects of some of the principal features of mold design. Carlsson and Hultgren6 used stearic acid to show that the columnar crystals of ingots tend to slope upward progressively from the mold wall toward the center of the ingot. As with the nonferrous ingots, the wax ingots have

Jan 1, 1943

By F. F. Schmidt, H. R. Ogden, E. S. Bartlett

Continuing tantalum alloy development studies have been concerned with a more detailed investigation of promising binary, ternary, and more complex tantalum alloys containing Groups IV-A, V-A, VI -A. VII-A, and VII-A metal additions. A primary objective of these studies was to improve the elevated temperature strength of tantalum with minimum degradation of fabricability and low-temperature ductility by alloying. Several alloys were found to exhibit hot tensile strengths in the excess of 25, 15, and 10 ksi at 2700°, 3000°, and 3500°F, respectively, combined with excellent cryogenic toughness and ductility. Particularly outstanding in this respect are alloys in the system Ta-Mo-W. Ternary parameter curves describing hot-strength, low-temperature ductility, recrystallization temperature, weld ductility, and fabricability are presented for the Ta-Cb-V, Ta-Hf-W, Ta-Mo-Hf, Ta-Mo-V, Ta-Mo-W, and Ta-V-W systems. At present, only four metals qualify for large scale structural applications at temperatures above about 1800°F. These are tungsten, tantalum, molybdenum, and columbium. Research and development work on these refractory metals has progressed to the point where efforts are now primarily directed toward the improvement of high- and low-temperature mechanical properties through alloying. High melting point combined with excellent low-temperature ductility, fabricability, toughness, weldability, and high solubility and tolerance for both substitutional and interstitial solutes are important attributes of tantalum. Conversely, its high density and relative scarcity compared with other refractory metals are the major deterrents to the use of tantalum. Early studies in the Battelle alloy development program established base line data on the mechanical and metallurgical properties of high-purity tantalum.' These studies were extended to include a preliminary investigation of the effects of various alloying elements on the established base properties.2,3 Recent studies involved a more detailed investigation of promising binary alloys containing Cb, Hf, Mo, Os, Re, Ru, and W, and ternary Ta-Cb-V, Ta-Hf-W, Ta-Mo-Hf, Ta-Mo-V, Ta-Mo-W, and Ta-V-W alloys. Results of major importance obtained are given in this article. EXPERIMENTAL PROCEDURES Electron-beam-melted tantalum (Wah Chang) and high-purity alloying elements were consolidated as 150-g ingots by multiple nonconsumable arc melting under a partial atmosphere of helium. Homogeneity of the button ingots was checked radiographically. Alloy buttons measured about 1-1/2 to 1-3/4 in. in diam by about 0.35 and 0.50 in. thick. Nominal alloy compositions were checked on each alloy button by weighing after arc melting. In most instances, weight losses were within 3 pct of the total weight charged. Subsequent chemical analyses on a number of selected alloys showed interstitial contents to range as shown below: Range of Interstitial Content, ppm 0 C N H 25-75 50-100 10-50 1-5 Metallic alloying additions were found to agree well with nominal compositions based on weight change data. Fabrication directly from the cast button to sheet was conducted at low (75° to 900°F), intermediate (1800° to 2000°F), or high (3000°F) temperatures, depending upon alloy content and hardness. Evacuated stainless steel or molybdenum packs were used for fabrication at intermediate and high temperatures, respectively. The temperature required to complete recrystal-lization in one hour was determined by hardness measurements and metallographic observations (longitudinal test sections) after vacuum annealing at 180°F intervals in the range 2010° to 2910°F. All mechanical property evaluations were conducted on recrystallized sheet, nominally 0.04 in. thick. Weldability of the various alloys was screened using manual TIG welding techniques. Bend ductility (bend axis perpendicular to weld axis) was the criterion used for evaluation. Tensile specimens employed 0.200 by 1 in. gage lengths, and bend test specimens measured about 1 by 1/4 in. Low-temperature tensile tests were conducted using a crosshead speed of 0.02 in. per min up to the point of yielding, and 0.05 in. per min to fracture. High-temperature tensile tests were conducted in vacuum using a crosshead speed of 0.01 in. per min

Jan 1, 1963

By M. C. Flemings, P. J. Ahearn

Dendvite morphology of unidirectionally solidified Sn-12 pct Bi alloy was determined with the aid of X-ray and macroetch techniques. Columnar growth direction (growth direction of primary arms) in this alloy is [110]; secondary arms grow in the [lii], [111), and [112] directions. Interstices between [112] dendrite arms fill in during early stages of solidification so that well-defined (110) planes are seen in the final structure (parallel to the heat-flozo direction). Similarly, interslices betweex the [111] secondary arrlls fill in to form (112) planes, and interstices between [1111 ] arms fill in to form (112) planes. DENDRITE morphology in alloys is conveniently studied in unidirectionally solidified castings. In this type of solidification, heat is extracted from a single chill surface of the casting while other surfaces are insulated. Solidification progresses from the chill and results in cast grains that are columnar with the long axis of the grain perpendicular to the chill. Dendrite morphology is conveniently studied by examination of etched surfaces parallel and perpendicular to the chill, as previously described.'" Experimental studies on the growth direction for tin are in general agreement. The major growth direction is [110]:-' although this is altered somewhat by supercooling. In one study, tin dendrites pulled from a supercooled melt had a growth orientation 12 deg away from the [110] direction toward the [002] direction.' This orientation could be changed back to the [110] or close to the [110] direction by changing the tempera- ture gradient in the melt. In a study using autoradio-graphic techniques, (111) directions were shown to be directions of preferred growth in addition to the [110) direction.7 PROCEDURE The samples employed in this study were cut from 12 pct Bi-88 pct Sn castings, 1.5 in. in diam and 3.5 in. tall. The casting procedure consisted of heating the charge to 390°C and then pouring the molten alloy directly into a gypsum sleeve mounted on a massive copper chill. Specimens were cut from the ingot, mounted, polished, and etched with Taffs reagent [FeC13 .6H2O, 2 g; HC1, 5 cu cm;H20, 30 cu cm; CzHaOH (95 pct), 60 cu cm].' The dendritic pattern was clearly visible after etching. Photomicrographs at about 8 and 35 times were taken on surfaces parallel and perpendicular to the chill. Laue back-reflection patterns were obtained from three different columnar grains in a singIe polished face 2 in. from the chill. The cut face was polished and positioned in the beam of a Laue camera, so that diffraction occurred from a single dendritic grain. The resulting diffraction pattern was interpreted to determine crystal log raphic orientation in relation to the heat-flow direction and in relation to microstruc-ture of the given grain. RESULTS Fig. 1 shows a horizontal and vertical section of the columnar structure, at low magnification. Light regions are rich in tin, and darker regions tin-poor; the primary phase is tin-rich a phase. The two sections shown were taken orthogonally from the same metallographic specimen. It is readily seen from Fig. 1 that the dendrites have

Jan 1, 1968

By P. Rama Rao, D. H. Sastry, T. R. Anantharaman

Deformation (a) and twin (a') stacking-fault densities in cold - wor ked filings of fcc silver alloys containing 2, 4, and 6 wt pct Sb have been estimated from Debye-Scherrer peak shifts and asymmetries, respectively. Evidence has been obtained for a possible strain-induced transformation in a two-phase alloy containing 9.0 zut pct Sb. Using Fourier analysis of X-ray line profiles, the effective particle sizes and root mean square strains have been determined for the three fcc alloys in the cold-worked condition. The values of the compound fault probability 0.5a + a') evaluated from peak shifts and asymmetries are in reasonable apeement with those derived from peak broadening. COLD work of fcc metals may produce both deformation and twin (also called growth) stacking faults, which can be detected and quantitatively estimated by X-ray diffraction techniques.' In addition to peak broadening, deformation faults produce a shift in peak positions and twin faults an asymmetry of the Debye-Scherrer reflections. Besides stacking faults, a reduction in the size of coherently diffracting domains and distortion within these also broaden the powder reflections of deformed fcc metals and alloys. Fourier analysis' of peak shapes provides quantitative infor- mation about the effective domain sizes, lattice strains, and also the compound fault probability (1.5a + a'), where a and a' represent deformation-and growth-fault density, respectively. The Paterson method2 for determination of a from peak shifts, the warren3 or wagner4 method for evaluation of a' from peak asymmetries, and the Fourier analysis for quantitative estimates of all the contributory causes of peak broadening have, in the course of the last decade, been extensively applied to obtain data for a number of deformed fcc metals and alloys. One of the useful results of earlier measurements of a in several solid solutions is a correlation between a and the solute valence as well as electron concentration per atom (p) in alloy systems based on a particular solvent. The work of Smallman and westmacott5 on Cu-Zn, Cu-A1, and Cu-Sn, of Foley, Cahn, and ~aynor' on Cu-Si and Cu-Ge, and of Rama Rao and Anantharaman7 on Cu-Si and Cu-Sb has revealed a common trend in the increase of a with p and also a with valency of the solute. Cahn and Davies' and Ad-ler and wagnerg have observed similar trends for silver alloys containing cadmium, indium, and tin. Data on a for silver alloys with a pentavalent solute element are not, however, available. It is primarily to fill this need that the present investigation on Ag-Sb alloys was undertaken. Utilizing the Patersonand Warren methods, values of a and a', respectively, have been determined in three silver alloys containing 2, 4, and 6 wt pct Sb. By Fourier analysis of the broadening of Debye-Scherrer peaks, effective particle sizes, lattice strains, and the compound fault probability (1.5a +a') have also been evaluated.

Jan 1, 1967

By J. T. Norton, A. L. Mowry

The monocarbides of the A subgroup elements in the fourth and fifth group of the periodic table in addition to being hard and refractory are of special interest in that they are isomorphous in crystalline structure. They are cubic with a sodium chloride type structure in which the metal atoms are essentially close packed in a face-centered cubic arrangement with the carbon atoms placed in the interstices between. Interstitial structures of this close packed type were first investigated systematically by Egg1 and he gave the rule for their formation, stating that the radius ratio of the nonmetal to the metal atom should not exceed the value of 0.59. The carbides of interest are those of titanium and zirconium of the fourth group and vanadium, columbium and tantalum of the fifth group. Table 1 shows the radius ratio using the Goldschmidt radii for 12 coordination for the metal atoms and the diamond radius for the carbon atom. It will be noted that while there is considerable variation in the size of the metal atom, in all cases the ratio is smaller than the limit of 0.59 placed by Hägg. It has been known for some time that these cubic carbides are soluble in one another, at least to some extent or, in other words, the metal atoms can be replaced, one by another without destroying the stability of the structure. Since the stability of these close packed interstitial substances appears to depend more upon geometry than upon the exact chemical nature of the atoms involved, it is of interest to examine the possibilities of replacement in these carbides in some detail. Hume-Rothery2 has pointed out the importance of the difference in size of solute and solvent atom as a factor in limiting the solubility in simple binary solid solutions. Largely on an empirical basis, he states that if the difference in size between solvent and solute atom is more than 14-15 pct of the solvent atom, the range of solubility is very restricted. The atom size was based on the distance of closest approach in the elements involved. While there is some question as to how one should calculate the size of the metal atom in the carbide structures, reference to Table 1 will show that zirconium is the largest and vanadium the smallest of the group and that the difference is about 15 pct. The Ti-Zr difference is about 9 pct and the others are smaller. Thus one would predict that if the size factor controls the solubility, all of the pairs except VC-ZrC would have wide or complete solubility whereas this latter pair is on the border line and might have restricted solubility. The purpose of the present investi- gation was to examine the solubility of the several pairs of carbides by heating them together until equilibrium was established and then examining the product by X rays. Previous Work Agte3 and his associates prepared various transition metal carbides and determined the melting points of binary mixtures. He concluded from the shapes of the melting point curves that there was extensive solubility in the case of the cubic carbides. Umanskii and his colleagues made an investigation of a number of pairs of the cubic carbides, using X rays and plotted lattice parameter vs. composition curves for the systems TaC-Tic, CbC-Tic, TaC-ZrC and CbC-ZrC. All pairs showed a continuous series of solid solutions. The first two pairs gave a linear relation while the latter two showed a negative deviation from Vegard's law. Kiefer and Nowotny, in a paper which became available after the present work was well advanced, investigated the binary pairs of the five cubic carbides by means of X rays. Relatively few points were obtained and results indicated that in some cases, at least, equilibrium was not reached at the temperatures used. The results indicated that solubility in the VC-ZrC system was not complete. All of the results of previous investigations indicated the desirability of a more detailed study. Materials The raw materials used were mono-carbides of titanium, zirconium, vanadium, columbium and tantalum and were the purest which could readily be obtained commercially. Spectrographic qualitative analysis showed that the CbC and TaC contained less than 1 pct

Jan 1, 1950

By R. L. Slobod, B. H. Caudle, E. R. Brownscombe

An earlier publication has discussed three methods lor obtaining relative permeability data on small core sample-: and the apparatus and technique for the capillary pressure displacement method. This paper describes the apparatus and technique for the solution gas and the dynamic displacement methods and presents a routine procedure for obtaining oil-gas and water-oil relative permeability data. Theoretical and experimental considerations are presented to show that the end effect commonly associated with the dynamic flow mechanism is extremely small where constant rates are employed in the flowing phase. An effect of flow rate on relative permeabilities obtained by the dynamic system is found only when gas is one of the flowing phases and this effect is ascribed to a form of channeling in the capillary system. The apparatus and procedure* used to obtain relative permeability data with three phases flowing are described and some preliminary results of the use of this method are shown. INTRODUCTION In a previous publication from this laboratory there appeared a brief discussion of the concepts behind three basic methods for obtaining relative permeability data.' These three methods were called the capillary pressure displacement method. the solution gas displacement method, and the dy. namic displacement method — the names being suggestive of the type of process used for obtaining the desired saturation prior to making the permeability measurements. In the same publication, the apparatus and technique for the capillary pressure displacement method were described and some typical results obtained by this method were presented. This paper will present: 1. The routine procedure used to obtain complete relative permeabilitv data (In small core samples with two phasec flowing; 2. The experimental techniques for the solution gas and dynamic displacement methods for obtaining relative permeability data; 3. The results of some studies on the mechanism of fluid flow through consolidated porous media; and, 4. The preliminary results on the determination of relative permeability for a system in which three phases are flowing. TECHNIQUE AND APPARATUS FOR THE DETERMINATION OF RELATIVE PERMEABILITY WITH TWO PHASES FLOWING In order to obtain relative permeability data on a small core sample which will describe the various phases of reservoir depletion. it is necessary that the methods used duplicate, as closely as possible, the various mechanisms involved in rescbrvoir depletion. It is important that. in any type of relative permedbility determination. the core be originally at virgin reservoir conditions with respect to the fluids within the pore *paces. (for example. oil used to displace water from core until connate water saturation estahlished). and that the subsequent saturation changes in the core be analagous to tile reservoir depletion mechanism for which the data are to he used. PROCEDURE OF DETERMINING RELATIVE PERMEABILITY DATA WITH TWO PHASES FLOWING A standardized procedure has been set up for obtaining relative permeability data and other information on small core samples. This procedure essentially consists of 11 steps as shown below: 1. Determination of tile capillary pressure curve.'! This information is used to help select samples for relative permeability determinations and to charactrrize. in general. the pore size distribution of the sample. 2. Determination of the specific permeabilities to air, oil. and water. These values are those on which the relative pern~eahility calculations are based. 3. Reduction of the core to connate water by an oil drive using the high speed centrifuge. 1. Determination of the critical gas saturation value and the associated relative permeability to oil by the solution gas drive method. 5. Determination of the relative permeabilities to oil and gab in the region where both phases are flowing by the capillary presqure method. 6. Reduction of the sample to residual oil in the centrifuge by an air drive. and the determination of the relative permeability to gas at this point.

Jan 1, 1951

By L. J. Erck, T. E. Ban

A discussion of laboratory procedures used to determine pellet quality and to simulate handling and firing conditions. Strength-temperature relationships in pelletizing; effect of chemical additives on agglomeration; and crushing, abrasion, and impact resistance of the pelletized products are discussed with respect to five different iron ore concentrates. IN most routine experimental work on products and processes, method:; of operation and product evaluation are established from accurately controlled laboratory tests. With particular emphasis on the pelletization of fine iron ore concentrates, standardized procedures of laboratory production and product testing have been established at Cleveland-Cliffs Iron Co. Research Laboratory at Ishpeming, Mich. A program of this nature was used for investigating the agglomeration characteristics of various iron ore concentrates. The factors which control the process were studied and subsequently controlled for these investigations. Pelletizing as applied to iron ore agglomeration involves the preparation of spheroidal masses from a fine moist concentrate, followed by high-tempera-ture firing to produce a final product. The preparation of moist spherorids is generally designated as balling. This produces the so-called green pellets or balls. After, these have been fired they are called pellets. Balling is accomplished by rolling moist concentrates in a revolving drum until small spheroids form by a snow-balling action. These balls are then heated in a shaft furnace, or a modification of this device, to a temperature high enough to cause agglomeration within the individual balls while they descend countercurrent to hot oxidizing gases. By this process the lour strength balls are converted through various chemical and physical-chemical processes into hard pellets. With correct balling and firing, the pellets are strong, resist abrasive action and handling treatments, and provide a suitable blast furnace and open hearth feed. In 1950 Tigerschiold and Ilmoni developed a working hypothesis to explain the formation and stability of green pellets.' Their work was presented on the basis of geologists' discoveries concerning the soil-water system and C. V. Firth's early postulates of green pellet formation.' Fundamental reactions oc-curing within fired pellets have been described by Cooke and Ban, who showed the microstructure changes caused by different firing temperatures." All these investigators have provided scientific clarifications of agglomeration by the pelletizing process and have opened avenues for greater development. Mineral grain growth and ceramic bond formation are the two basic phenomena imparting coherency to particles within a fired pellet of iron ore. Each of these is established by the transfer of thermal energy to the mineral particles. Grain growth mechanisms are chiefly physical reactions wherein smaller particles consolidate into larger ones, with resulting decrease in surface energy conforming with the second law of thermodynamics. Ceramic bond formation is a chemical reaction usually between basic and acid metallic oxides of the gangue and ore minerals. CaO, MgO, A1203, and FeO of Fe3O4, react in varying proportions with SiO, or Fe2O3, to form a great series of compounds or amorphous solid solutions. These product bond iron mineral particles together by a bridging action or by direct reaction with the surfaces of the particles. The result of the grain growth-ceramic bond formations is a consolidated, coherent compact of gangue and iron minerals formerly held together by much weaker bonds of interstitial water and forces common to dry, unconsolidated particles.4 With respect to iron ore agglomeration, these high temperature bonding reactions have been described by Schwartz, who referred to the bonds in sinters," and by Cooke and Ban, who referred to the bonds in fired pellets." The extent of the bonding phenomena and the consequent strength of fired pellets is governed by

Jan 1, 1954

By L. L. Handy

Two equations are compared for calculating imbibi-tion rates of water into porous media. A linear relation was found to apply when the square of the volume of water imbibed was plotted against time for water imbibing in natural rock materials. The effect of adsorbed organic materials on surface wettability was determined from measurements of rate of imbibition, water saturation and water permeability. INTRODUCTION An understanding of the mechanisms by which water displaces oil from porous media requires an understanding of the role of capillary forces in the displacement process. The capillary forces are determined by properties of the water-oil-solid surfaces. Consequently, some measure of the surface properties of reservoir rocks is necessary if displacement behavior in oil reservoirs is to be predicted from laboratory measurements. Several experimental procedures have been proposed to assign quantitative wettability indexes to reservoir rock surfaces. The most recent of these proposals are those of Bobek, Mattox and Denekas' and Amott.' These indexes are designed to show a continuous variation from the preferential oil-wet to the preferential water-wet systems. They require measuring some property of the rock which is a function of surface wettability. The quantities are measured on unaltered core material and compared with values obtained for known oil-wet and water-wet extremes on the same material. These methods are useful but are semi-empirical in nature. They have the disadvantage that the measured quantities may be functions of other variables in addition to surface wettability. Some understanding of the behavior of water-oil-rock systems may be obtained from studies of the simpler water-air-rock systems. Surfaces of a given wettability are more easily maintained for the latter. Also, certain approximations can be made for the flow of water and air which are not justified for water-oil flow. The theoretical and experimental work discussed in this paper is part of a study of the role of capillarity in the dynamic displacement of air by water. The expec- tation is that observations for the water-air system will lead to generalizations which are applicable to water-oil systems. THEORY Imbibition is defined as that process by which one fluid displaces another from a porous medium as a result of capillary forces only. Differential equations have been derived to describe the displacement of oil by water during imbibition, but these equations have not been solved.".' With certain assumptions, however, solutions are possible for water-air imbibition.

By K. A. Gschneidner, O. D. McMasters

DgIel-ential thermal, nretallographic, and X-ray paramzetric methods were used to establish the Yb-Pb phase dingram. The terminal solid solubilities in the system are less than 0.2 at. pct. Lead additions lower the 792°C transformation tewzperature of pure ytterbiun~ at 787°C by means of an inverted peritectic veaction. Ytterbium additions raise the 327oC melting point of lead to 329°C by means of a perilectic reaction. Three intermetallic compounds Yb2Pb, YbPb, and YbPb, melt congruerttly at 1246", 1116o, and 740°C. respectively. YbsPb3 forms by a peritectic reacliorz at 1250°C. YbzPb crystallizes in the PbCh (C23) type structure and the orthorhombic latliceo constants are a, = 7.478. bo = 5.225, aud C0 = 9.549A. YbsPbs crystallizes in the hfnsSi3 (D8a) type structure. The hexa~nal lattice constants for this phase are a, = 9.325 and co = 6.92d. YbPbl crystallizes in the cubic AuCU3 (Ll2) type structure with a: = 4.8628 + 0.0008A. YbPb is diwrorphic with the structure change occurring al 507°C upon healing. The crystal structures of the two forttzs were not determined. Eulectic reactions occur at 5.1 at. pct Pb and 769oC, 41.5 at. pct Pb and 1086oC, and G8.0 at. pct Ph and 717°C: RELIABLE information concerning the alloying behavior of most of the rare-earth metals is scarce. A systematic study of the rare earth-lead alloy systems now in progress is an effort to supply such information. Determination of the phase diagram of the Yb-Pb alloys is the first step of this study. Variations in the valence of normally divalent ytterbium may be studied by comparing the characteristics of the Yb-Pb alloys to those of the trivalent rare earth-lead and of the divalent alkaline-earth-lead systems. Of particular interest are the melting temperatures of the intermetal-lic compounds since these are known to be high relative to the melting points of the constituent metals. Crys-tallographic and magnetic susceptibility data for the rare earth-lead alloys are to be included in the comparative study. Results from this and subsequent investigations should provide a better understanding of alloy formation. EXPERIMENTAL PROCEDURE Materials. The major impurities in the ytterbium and lead used in this investigation are listed in Table I. The lead was obtained from Cominco Products, Inc., and was specified to be 99.99 pct pure. The ytterbium was prepared in this laboratory by the lanthanum reduction of the oxide followed by a distillation of the ytterbium. Alloy Preparation. The alloys were prepared by

Jan 1, 1968

By T. Seldenrath

In the development of a fluid -operated hammer drill' for accelerated penetration of hard rock formations in oil wells, a research investigation was conducted to evaluate the percussion effects obtained with different design characteristics and to determine the possibilities of percussion in the low frequency range. Tests were conducted on granite blocks and comparable impact forces were measured with a load cell under selected test conditions. These full-scale laboratory tests provided an evaluation of the effectiveness of percussion which was further supported by actual downhole field performance. While much pertinent data have been presented by other investigators the method of evaluation described in this article resulted in good correlation between laboratory and field performance and may be applicable in other low frequency percussion developments. EXPERIMENTAL METHOD For these tests the bit was held stationary and the test block was rotated and forced upward against the bit by a hydraulically-lifted rotary table. "Weight on bit" was calculated from pressure readings on the fluid system of the lifting cylinders—corrected for tare and friction. Percussion was obtained with several tubular hammers of different weights up to 320 lb. These were lifted mechanically and allowed to drop by gravity from various heights onto a slidably-mounted bit sub or "floating anvil" for calibration of the load cell response with known kinetic energy of the hammers. During evaluation of percussion effects, the hammers were hydraulically operated to provide a range of percussion frequencies up to 1,020 blows per minute delivered to the anvil. To obtain comparable test readings, the load cell of the strain gauge type, was attached to the bottom of the anvil or bit sub in place of the bit. The load cell rested on a typical granite test block mounted on the drilling table. The latter was hydraulically supported as in the actual drilling operation. The effect of support resistance was the same for all tools when calibrating. Impact forces so determined werethusdirectly comparable with one another for the purposes of this investigation. When a roller bit was substituted for the load cell for actual drilling, the impact force would probably be reduced due to the resiliency of the bit body, legs. cones and teeth and the penetration of the teeth into the rock. However, since the bits and rock specimens were uniform in these tests, an evaluation of actual tool performance, in terms of parameters previously determined with the load cell at the selected standard conditions of calibration, proved satisfactory. Throughout these tests, 8¾-in. standard tricone, hard rock toothed roller bits and carbide-insert roller bits were used. Rotary speeds were 50 and 100 rpm; static weight on bit was applied up to 40,000 Ib, and percussion frequencies up to 1,020 blows per minute were available. During this investigation it was found that toothed bits on which the teeth had developed flat ends about 3/16 in. wide did not change appreciably for the duration of a test. Hence, all toothed bits used in this study were in this condition. Conventional reference curves obtained in these tests show that increasing weight on hit by 2:1 yielded an increased penetration rate of approximately 2.7:1, while an increase in rotary speed of 2:1 netted an increased penetration rate of roughly 1.7:l. Results are thus of the same order as those reported for hard rock drilling by other investigators." CALIBRATION Under the selected standard test conditions previously described, with the lifting table hydraulically supported, and the anvil interposed between hammer and load cell, the hammers were dropped from known heights. Thus, with the kinetic energy of the hammer known, the corresponding maximum impact force shown by the cell was determined. (It may be of passing interest to note here that special drop tests with the hammers showed that no appreciable reduction in impact force resulted from introduction of the anvil or hit sub.) From these calibrations, it was possible to determine the theoretical kinetic energy available in the hammer from the maximum impact force shown by the load cell while operat-

By C. L. Thornburg

FLANKING the Valley of Mexico on the northeast is a mountain range known as the Sierra de Pachuca. This northwesterly-trending range is about 30 miles long and 5 miles wide, its summit attaining an elevation of more than 10,000 feet above sea level, or 2000 feet above the valley floor. Pachuca, a town of 50,000 inhabitants, lies nestled at the southwest base of the range, 60 miles northeast of Mexico City. Three miles to the east, just over the summit and on the northeasterly slope, is the mountain town of Real del Monte with a population of about 20,000. These are the two principal towns in a 40-sq mile area whose yield has brought the district to its high rank among the world's silver producers. Total production probably exceeds 1.25 billion ounces of silver and 4.5 million ounces of gold. Exploitation by the Spanish was under way by 1530. As the mines were deepened operations became handicapped by the inadequate method of handling water with horse whims and bull skins. To offset this disadvantage a two-mile drainage tunnel was started in 1749 and completed ten years later. John Taylor of London acquired the mines from the third Count of Regla in 1824 and formed the Compania de Real del Monte y Pachuca. Two years later Cornish pumps were brought from England. A Mexican company purchased the British interest in 1848 and within a few years resumed an old project of driving a second, lower three-mile drainage tunnel, completing it in 1857. After American interests acquired control of the property in 1906 operations expanded to an unprecedented scale, despite periods of political strife, unstable prices, labor problems, and a growing burden of taxation. An idea of the scale of operations may be gained from records of the last two decades prior to 1947, the last year for which production information is at hand, during which time production ranged from 1.4 million to 1.1 million tons per year, and the combined development and exploration, exclusive of diamond drilling, amounted to about 18 miles annually. Pachuca differs from Zacatecas and some other Mexican silver vein districts that flourished and declined in earlier centuries in that a substantial part of its total production was made after 1900. While the district's accelerated activity in this century was largely due to the successful application of the cyanide process to the treatment of the ores, improved pumping methods, and new mechanical equipment, it was at the same time stimulated by the opening up of important new orebodies as a result of active and extensive exploration of the vein system. Some of these orebodies were found on veins which terminate upward far below the surface, and though these particular veins do not crop out, they show relationship with a type of alteration which may be their surface expression. The district has been studied by many able geologists, and a large store of information has been accumulated, much of which is recorded in published as well as private papers. Wall Rocks Extrusive Rocks: The Sierra de Pachuca is a thick accumulation of Tertiary eruptive rocks consisting mainly of flows, breccias, and tuffs. Numerous recognizable volcanic vents within the range itself manifest that the range was built up by material ejected from these vents. The maximum thickness of the volcanic pile is not known, but from deep exploration and from relationship with Mesozoic sediments to the north and east it is inferred that the average thickness in the main part of the district may be at least 6000 ft. The extrusive rocks comprise four general types, which in ascending order of age are: andesite, rhyolite, dacite, and basalt. The veins are confined largely to the andesite, the dominant rock of the region. Characteristically this rock is augite andesite, but it includes comparatively acidic as well as basic layers, the acidic strata being

Jan 1, 1952

By Byoung Whie Lee

In order to investigate the change in surface energy as a function of cold work ad environment, the contact angles between a liquid and the surface of elec-tropolislzed aluminum with various degrees of cold work were measured. Specinlens cut from aluminum rod (99.99 pct Al) were annealed at 600°C for 24 hr and then cold compressed parallel to the drawing di-rectzon of the rod in increments of approximately 10 pcl from 0 to 70 pct. The sections, taken transversely to the drawing direction, were mechanically and electro chemically polished. Contact angles between the sessile drop of distilled water, glycerin, and glycol and the specimen surface were measured optically. The differences in interfacial energy were calculated from the contact angle data. The effects of chemical washing and surface roughness of the specimen on the contact angle were also studied. The liquid contact angle was a linear function of the degree of cold working (D) ad changed markedly with differed liquids. THE plastic deformation of polycrystalline metals induces a preferred orientation in crystal structure. A preferred orientation changes the atomic density of a metal surface. This change in surface atomic density brings about the changes in surface interactions of the metals."' Also, the plastic deformation of metals raises the internal energy by increasing the stored energy.' Consequently, the total Helmholtz free energy of the metal increases and, in turn, the surface free energy of the metal. In order to study the effect of plastic deformation on solid surface phenomena, the liquid contact angle of a sessile drop made on cold compressed and elec-trochemically polished aluminum was measured. The preferred orientation in aluminum, due to uniaxial compression, is known to be [110] direction parallel to the compression axis. However, this description is only a first approximation of a rather complex orientation. In actuality, half the crystallites are more than 10 deg from this orientation, regardless of the amount of compression." For randomly oriented fcc crystals, the average surface atomic density can be calculated as follows: The number of atoms in a unit volume is 4/a3, where a is a lattice constant. Hence, the number of atoms per unit surface area is about (4/a3)2/3 = 2.52/a2 However, when the preferred orientation took place along [110] direction, the surface atomic density is as follows: The area of (110) plane of a unit cell, considering ridges and valleys, is BYOUNG WHlE LEE, Member AIME, is Research Metallurgist, Atomic Energy Research Institute of the Republic of Korea, Seoul, Korea. Manuscript submitted May 29, 1968. IMD 2 x J2- a x Ga = 6a2 Since this plane contains 4 atoms, the surface atomic density is 4i6a2 = 2.3/a2 Upon changing from random to preferred orientations, the surface atomic density decreases about 10 pct. The increase in the stored energy and the change in surface atomic density due to plastic deformation would certainly alter the surface tension of a metal. EXPERIMENT High purity aluminum rod of 6-mm diam (99.99 pct Al, obtained from Johnson, Matthey & Co., Ltd, London, England) was cut into sections approximately 1 cm long. The transverse cross section was mechanically polished. Then, the specimens were annealed in air at 600°C for 24 hr. Thus, eight annealed specimens were uniaxially compressed at room temperature from 0 to 70 pct of cold working with an increment of 10 pct. The extent of cold working was determined by the percentage of reduction in axial length to the original length. The specimens were further polished mechanically and, then, electrochemically with the Modern Battelle Type 1 electropolishing solution with a current density of 0.3 to 0.4 amp per sq cm and an electrode potential of 8 to 12 v at the solution temperature of 78" to 82°C for approxirr~ately 5 min. The Modern Whlte Diffusing Screen \t/ --------------- re Goniometer for Mercury Lamp ^-^ Sessile Drop ¦=• Condenser c=> Microscope l_*~_i bJ Mlcromster Filer with Protmcter Fig. 1—Schematic diagram of apparatus set-up for liquid contact angle measurement.

Jan 1, 1970

By M. G. Benz

The superconducting performmce of diffusion-processed Nb3Sn is influenced by its micro structure. High isotropic transverse current density may be achieved in this material by a process which forms a precipitate of ZrO, within the Nb3Sn. FOR an ideal type-I1 superconductor, little or no transport current can be carried in the mixed state; i.e., little or no transport current can be carried above the lower critical field H,,, where the field penetrates abruptly in the form of current vortices or fluxoids, even though full transition to the normal state does not occur until the upper critical field H,,.' Fortunately, nonideal type-I1 superconductors can be readily obtained and these carry large transport currents up to the upper critical field H. Both theoretical and experimental investigations have attributed this current-carrying capability for nonideal type-I1 superconductors to pinning of the fluxoid lattice by heterogeneities in the microstructure of the superconducting material. These heterogeneities may take the form of dislocations or dislocation clusters,2"5 grain boundaries: structural imperfections introduced by phase transformations; radiation damage,8"10 or precipitates.11"15 Nb3Sn formed by diffusion processing is a type-I1 superconductor. Heterogeneities are needed for high superconducting critical currents above H,,. This paper will cover: a) what the microstructure of diffusion-processed NbSn looks like; b) what changes in the microstructure take place when the system is doped with precipitates, and c) how these changes in microstructure influence the superconducting critical currents. EXPERIMENTAL Preparation of Samples. Diffusion processing was used to form the Nb3Sn. The procedure used was as follows: a) coat the surface of a niobium tape with tin; b) heat-treat this tape at a temperature above 930°C to form a layer of Nb3Sn at the Sn-Nb interface. Such a layer of NbsSn is shown in Fig. 1 The thickness of the NbsSn layer formed was controlled by the time and temperature of the heat treatment. The same general procedure was used for preparation of both undoped samples and samples doped with a precipitate. An additional step was included in the preparation of the doped samples which consisted of internal oxidation of zirconium to form ZrOn. The details of the doping process will be reported in a later paper. Sample Testing. The Nb3Sn tape samples were soldered to a copper or brass shunt. Current and voltage leads were then attached to the sample in the usual four-probe resistance measurement configuration. The sample was cooled to 42°K. In some cases it was cooled in the presence of a high magnetic field and in other cases with the field turned off. The results were the same for both cases. The samples were oriented in a configuration with field transverse to current but could be rotated such that the angle between the field vector and the wide side of the tape sample could be changed. Measurements up to 100 kG were done in a superconducting solenoid and measurements above 100 kG in a water-cooled copper magnet at the MIT National Magnet Laboratory. Once the test field was reached, the current in the sample was increased until voltage was detected across the sample. The critical current was taken as the current at which voltage was first detected in excess of background noise. In most cases this was 1 to 2 x 10~6 v for a— in.-wide sample carrying several hundred amperes with a in. separation between voltage leads and with a 10 "-ohm shunt resistance. RESULTS AND DISCUSSION Microstructure. Examination of the microstructure of the undoped Nb3Sn shows rather large-diameter (1 to 2 columnar grains growing outward from the niobium surface toward the tin surface. As the layer is made thicker by longer diffusion times, these grains grow longer. Few new grains are started. Transmission electron microscopy shows little or no second-phase material within the bulk of the Nb3Sn layer. The microstructure of a diffusion-processed NbsSn layer changes quite drastically when the system is doped so as to form a precipitate within the NbsSn layer. Instead of large-diameter columnar grains of NbaSn forming, smaller-diameter (0.5 to 1 ) equiaxed grains of Nb3Sn decorated with the precipitate form. Fig. 2 shows a transmission electron micrograph of a Nb3Sn layer doped with zirconium oxide. This layer has been etched so that one may look between the grains

Jan 1, 1969

[ ]

Jan 1, 1938

By S. G. Gibbs

A new method for predicting the behavior of sucker-rod pumping systems is presented. The pumping system is described by a flexible mathematical model which is solved by means of partial diflerence equations with the aid of computers. Polished rod and intermediate-depth dynamometer cards can be calculated for various bottom-hole pump conditions. The technique permits simulation of a wide variety of operating conditions, both normal and abnormal. The data generated with the new technique are useful in refining the criteria for design and operation of sucker-rod systms. INTRODUCTION Sucker-rod pumping systems are used in approximately 90 per cent of artificially lifted wells. In view of this wide application, it behooves the industry to have a fundamental understanding of the sucker-rod pumping process. Oddly enough, our understanding has been rather superficial. This is evidenced by the semi-empirical formulas which have been used as the basis for design and operation of sucker-rod installations. Thoueh- we have realized the limitations of our methods for many years, it has not been computationally feasible to use more refined techniques. With the advent and widespread use of digital computers, it is now possible to handle the mathematical problems associated with sucker-rod pumping. This paper summarizes a computer-oriented method which can provide greater insight into the sucker-rod pumping process. It is hoped that this technique, and techniques which may evolve from it, will prove to be the tool needed by industry to obtain the most efficient use of rod pumping equipment. THE MATHEMATICAL MODEL Prediction of sucker-rod system behavior involves the solution of a boundary value problem. Such a problem includes a differential equation and a set of boundary conditions. For the sucker-rod problem, the wave equation is used, together with boundary conditions which describe the initial stress and velocity of the sucker rods, the motion of the polished rod and the operation of the down-hole pump. Of these items, the wave equation, the polished rod motion condition and the down-hole pump conditions are of primary importance. Discussion of the mathematical model centers about these factors. ROD STRING SIMULATION WITH THE WAVE EQUATION The one-dimensional wave equation with viscous damp- is used in the sucker-rod boundary value problem to simulate the behavior of the rod string. This equation describes the longitudinal vibrations in a long slender rod and, hence, is ideal for the sucker-rod application. Its use incorporates into the mathematical model the phenomenon of force wave reflection, which is an important characteristic of real systems. The viscous damping effect postulated in Eq. 1 yields good solutions, even though nonviscous effects such as coulcomb friction and hysteresis loss in the rod material are present. Fortunately, the nonviscous effects are relatively small, so the viscous damping approximation used in the wave equation is adequate. The coefficient v is a dimensionless damping factor which is found in field measurements to vary over fairly narrow limits. For mathematical convenience the gravity term is omitted in Eq. 1. The effect of gravity on rod load and stretch can be treated separately, as will be noted later. Since Eq. 1 is linear, the legitimacy of this procedure is easy to demonstrate. POLlSHED ROD MOTION SIMULATION The motion of the polished rod is determined by the geometry of the surface pumping unit and the torque-speed characteristics of its prime mover. By determining the motion of the polished rod, we formulate an important boundary condition. From trigonometrical considerations it can be shown that the position of the polished rod vs crank angle 0 is given by (see Fig. 1) These equations are obtained from the general solution of the "four-bar" linkage problem and can be used to describe the kinematics of any modern beam pumping unit.' If prime mover speed variations are disregarded, the angular velocity of the crank is constant, and Eq. 2 can be used to predict the position of the polished rod vs time. However, the constant-speed condition leading to constant crank angular velocity is only approached in practice: hence, it is better to make provisions for prime mover speed variations in the mathematical model. The speed at which the prime mover runs is determined by its torque-speed characteristics and the torque imposed upon it. The torque that the prime mover "feels" is the net torque arising from the polished rod load and the opposing torque from the counterbalance effect. The

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.

By T. Cantwell, N. C. Rasmussen, H. E. Hawkes

Beryl is a mineral that may be difficult to distinguish from quartz by casual field inspection. The easily recognized green color and hexagonal crystal form of coarse-grained beryl are by no means universal, even in beryl from pegmatitic deposits. If it occurred as a fine-grained accessory mineral in an igneous rock, it would almost certainly escape detection unless samples were submitted for petrographic or chemical analysis. There may be substantial deposits of some beryllium mineral, other than beryl, that has been overlooked because that mineral also closely resembles the common rock-forming minerals. A reliable and simple method of identifying beryllium minerals and determining the beryllium content of a rock would be helpful in exploration. This article describes preliminary experiments in applying nuclear reaction to the qualitative identification of beryl and to the semiquantitative determination of the beryllium content of rock samples. Gaudin,1,2 the first to apply a nuclear reaction in detecting beryllium minerals, developed a method that irradiates the sample with gamma rays, which react with beryllium nuclei to produce neutrons. The neutrons are then measured with standard equipment. The cross section for this reaction is about 1 millibarn. The cross section is a measure of the probability that a reaction will take place, for example, between a beryllium nucleus and an incident gamma ray or alpha particles.3-5 At 1-millibarn cross section for the reaction, satisfactory performance required a source strength of the order of 1 curie (3.7 x 10"' disintegrations per sec, where each disintegration releases one or more gamma rays). The reactions will not take place if the gamma radiation is below a minimum energy, in this case 1.63 mev. The size of the source and the energy of the radiation made heavy shielding necessary for these experiments, both to reduce the background count of the neutron counter and to safeguard personnel. The original discovery of the neutron by Chad-wick in 1932 resulted from experiments with another nuclear reaction, induced by bombarding beryllium with alpha particles in which the products are carbon-12 and neutrons. The equation for this reaction is as follows:' " ,Be" + ,He'? 6C12 + 8,n' [1] re-particle neutron In the above nuclear equation (Eq. 1), the sub- script number indicates the number of protons in the nucleus (the atomic number) and the superscript the total number of neutrons and protons (approximately the atomic mass). For the alpha-neutron reaction the cross section is about 250 milli-barns, or 250 times that of the gamma-neutron reaction used by Gaudin. The positively charged alpha particle is repelled by the positive charge of the beryllium nucleus; it must, therefore, have a certain minimum energy in order to approach close enough to the beryllium nucleus to react. For reaction with the beryllium nucleus, the lower limit of the alpha-particle energy is 3.7 mev. The alpha-neutron reaction, with polonium-210 as an alpha source, was selected for the present experiments. Alpha particles are emitted by polonium-210 at 5.30 mev, which is adequate for the reaction with beryllium. Furthermore, this isotope of polonium emits alpha particles with negligible associated gamma radiation, thus eliminating the necessity of shielding. The half-life of polonium-210 is 138 days. Inasmuch as alpha particles carry a possible charge and are large compared with most nuclear particles, their energy is rapidly dissipated in passing through matter. Their range in standard air is 3.66 cm,3 and they penetrate only a few tens of microns into a mineral sample. The short range in air can be minimized by preparation of a flat sample surface that can be brought very close to the alpha source during analysis. On the other hand, short range of alpha particles in air lessens the radiological health hazard and makes it possible to use this method without shielding. It must be emphasized, however, that the alpha emitters are potentially very dangerous if they enter the human body. Polonium must be handled with extreme caution. The literature has reported experiments on the yield of neutrons from reaction of alpha particles with beryllium nuclei. Feld" reports that in intimate mixtures of polonium and beryllium, 3 x 106 eutrons per sec are produced per curie of polonium. Elsewhere in the same reference it is stated that a sandwich-type source yields about one third as many neutrons as an intimate mixture. A table of neutron yields for full energy polonium alpha-particles on thick targets as reported by Anderson7 is the basis of Table I. From Table I it can be deduced that the elements most likely to interfere, i.e., those that also produce neutrons when bombarded by alpha particles, are boron and fluorine. These data also show that it will probably not be possible to determine very small quantities of beryllium in rocks because of the masking effects of the major elements, sodium, magnesium, and aluminum. The neutrons emitted in the alpha reaction are detected by another nuclear reaction. Either of the

Jan 1, 1960