Search Documents

By L. A. Rapoport, W. J. Leas

As a preliminary, consicleration is given to the conventional definition of relative permeability and to the conditions governing the simultaneous flow of oil and gas through porous media. For the conditions of flow prevailing throughout most of a gas drive reservoir, the oil and gas can reasonably be supposed to be in capillary equilibrium with each other. Under these conditions, and these conditions only, the relative permeability to liquid can be expressed as a function of saturation. The relative permeability to liquid in that case is dependent upon the distribution of fluids which itself is shown to be related to the capillary pressure, and, in turn, to the saturation. As a consequence, relative permeability to liquid can be expressed in terms of the volume and surface area of a network of liquid channels bounded by the rock and the gas phase. While the volume of this network can be evaluated accurately, the surface area cannot. However, for any such volume, maximum and minimum values of the corresponding surface area can be calculated from capillary pressure data. It is then possible to establish for any saturation the limits within which the value of the relative permeability to liquid must lie. As a consequence of the theoretical development, the validity of an experimental method for measuring relative permeability to liquid which utilizes a stationary gas phase is demonstrated. In this method capillary barriers are cemented to the ends of the core sample to permit the maintenance of capillary equilibrium between the two phases. At the same time, this procedure eliminates undesirable secondary phenomena such as end effects, fissure effects, etc., the presence of which adversely affect the results of other laboratory methods. The results obtained by theoretical calculations, and experimentally, are discussed. In view of the overall precision that can presently be obtained in reservoir calculations, the agreement between the calculated and measured relative permeability to liquid data can be considered satisfactory. In conclusion, for reasons of economy and simplicity, the procedure of calculating limiting relative permeability to liquid curves from capillary pressure data is indicated for general engineering purposes. It is shown that the above procedure can easily be extended to the cases where connate water is present. Its use for reservoir studies is particularly recommended in conjunction with the method for measuring relative permeability to gas' which simultaneously yields the capillary pressure data necessary for the calculations. THEORETICAL Definition of Relative Permeabilities — Basic Equations for Heterogeneous Flow The equations by which the relative permeability concept is defined and upon which the formulation of all of the gas-oil flow problems rests at the present time are expressed as: V, = — Grad PL = — Grad PL ....(la) PL µL, k k KG VG = —K - Grad PG k- Grad Pc .... (lb) Mo where ,. and G refer to liquid and gas; V is the volumetric rate of flow per unit gross area. µ the viscosity, Grad P the potential gradient. and k the specific permeability of the porous medium.* (For horizontal flow, Grad P becomes the pressure gradient; otherwise, gravity must be included.) According to these expressions, each of the constituent phases is considered similar to a homogeneous system where the volumetric rate of flow is proportional to the pressure gradient, and for each of which the constants of proportionality, k, and kG, are termed effective permeabilities. by analogy to the specific permeability as defined by Darcy's law in its original form. In order to obtain a convenient basis of comparison, the effective permeabilities are referred to the specific permeability, k, of the considered porous medium, with the help of the relations: k1. = KI. k kQ = Kn k..........(2) K,. and Kr. are defined as the relative permeabilities to the liquid and to the gas phase respectively, and frequently expressed in per cent of specific permeability. It may be seen that Equations (1); which appear to be a direct generalization of Darcy's law, correspond to the assignment at any given time of a set of "local" permeabilities to each point of the porous medium, and represent in a differential form the two fluid flow system as a simple superposition of the individual single flow systems. The above interpretation implies that the effective or relative permeabilities are independent of pressure or rate of flow,

Jan 1, 1951

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey' have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks' have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.' Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method

By Nicholas J. Grant, Michio Yamazaki

Cast copper alloys containing 10, 20, and 30 pct Ni and 0.75 to 0.80 pct Al were machine-milled into chips, then comminuted in a rod mill to fine flake powder utilizing a number of processing variables. The powders here internally oxidized, mostly at 800°C, in a low-pressure oxygen atmosphere. The consolidated powders were hot-extruded into bar stock. Room-tenmperature tension tests, stress-rupture tests mostly at 650°C, but also at 450° and 850°C, and hardness measurements after various annealing temperature treatments to study alloy stability were perfomted. Excellent room-temperature strength, high rupture strength at 650°C, and resistance to recrystallization at 1050°C were obtained. Problems in optimizing conditions for internal oxidation of Cu-Ni base alloys are discussed. THE interesting high-temperature properties of SAP' have stimulated considerable effort in the study of more refractory alloy systems where the potential for high-strength alloys at high temperature is great.2-13 A number of methods have been utilized to produce the desired fine, hard particle dispersions, of which internal oxidation2,9,7 of dilute solid-solution systems offers considerable promise by virtue of the potential for producing ultrafine, well-dispersed oxides. While most of the published works are concerned with pure metal matrices, a number of investigators have studied the effects of solid-solution strengthening.10,19 Use of more complex alloy matrices (for example, aging systems) has been unsuccessful because overaging still occurs at high temperatures in the oxide-containing alloys.14.15 Solid-solution strengthening is, however, effective at very high temperatures9,10 and might be expected to contribute importantly to the strength of oxide-dispersion strengthened alloys. For this study, internal oxidation of solid-solution alloys of copper and nickel, containing small amounts of aluminum, was chosen as the method of alloy preparation. PREPARATION OF ALLOYS Three copper alloys containing about 10, 20, and 30 pct Ni and each containing 0.75 to 0.80 pct A1 (enough to yield about 3.5 vol pct alumina) were prepared as air-cast ingots measuring 2.5 in. diameter by 6 in. high (see Table I for the analyses). Processing steps for all the alloys were as follows (also see Table 11): 1) Homogenization of the ingot at 982°C (1800°F) for 45 hr in an argon atmosphere. 2) Machine milling of ingots into fine chips. Average thickness was about 0.1 to 0.2 mm. 3) Hydrogen reduction of chips at 593°C (1100° F) for 1 hr to reduce copper and nickel oxides. 4) Rod milling of chips to finer powders. 5) Hydrogen treatment of powders as in step 3. 6) Internal oxidation of the powders. 7) Hydrogen treatment of oxidized powders as in step 3. 8) Hydrostatic compression of evacuated powders. 9) Sintering of compacts in hydrogen. 10) Hot extrusion. Variations in processing among the alloys were made in steps 4, 5, and 10 (see Table 11). In the past, two methods were utilized to internally oxidize alloy powders. Preston and Grant3 surface-oxidized dilute Cu-Al powders to obtain the necessary amount of oxygen to oxidize the solute metal (aluminum and silicon), and then permitted the formed copper oxide to diffuse and react with the solute in an argon atmosphere. Bonis and Grant4 exposed Ni-A1 and other nickel alloys to an oxygen pressure derived from the decomposition of nickel oxide at a preselected temperature, in an argon atmosphere. Both methods are applicable and can be modified to generate a range of oxygen pressures for oxidation of the solute but not the solvent metals. Procedure I: Surface Oxidation of Alloy A3, Cu-10Ni-0.76A1. Powders of -20 to +28 mesh were surface-oxidized at 500°C (932°F) to obtain the desired amount of oxygen for oxidation of the aluminum to alumina; the powder was then sealed in Vycor and heated at 900°C (1652°F) for various times up to

Jan 1, 1965

By Craig S. Hartley

Expressions are obtained for the energy changes associated with the reaction of (a& (111) slip dislocations on intersecting (110)planes in anisotropic bcc metals. An energy criterion for assessing the likelihood of dissociation of the products of such reactions is also presented. It is found that the "burrier reactions" which form a(100) dislocations at the intersection of two active {110) slip planes are more energetically favorable in metals which exhibit a high value of Zener's anisotropy factor, A, than those which have a low value. The results are presented in a form which permits the stacking fault energy to be obtained from a measurement of the separation between par-tials in a dissociated configuration. However, until accurate calculations or measurements of the stacking fault energies involved are available, it is not possible to assess the physical importance of dissociated dislocations. In a recent paper,' the energy changes associated with several types of reactions between two slip dislocations, (a/2)(111){110), in bcc structures were calculated.* Isotropic elasticity and the approxima- tion v = -3- were employed. The purpose of this work is to present calculations of the energy changes for many of the same reactions using anisotropic elasticity. The problem of dissociation of a(100) and a(110) dislocations is also considered, and maximum fault energies for which dissociation will be energetically favorable are calculated for several bcc metals. Two general types of reactions are considered; those for which the reactant (a/2)(111) dislocations have long-range attractive forces and those for which the reverse is true. An example of the former is: (a/2)[lll] + (a/2)[lll]-a[l00] while the latter are typified by: (a/2)[lll] + (a/2)[111] -a[011] Only reactants lying in different slip planes are considered; therefore, the products must lie along (111) or (100) directions, which are the intersection of two {llO} planes. It will be assumed that the reactants and products are infinitely long parallel dislocations, since in this case the energy change associated with the reactions is a maximum.' THEORY The self-energy per unit length of a straight mixed dislocation in an anisotropic medium can be written? where b is the Burgers vector, K is an appropriate combination of the single-crystal elastic constants, and R and ro are, respectively, outer and inner cut-off radii of the elastic solution. The energy given by Eq. [I] does not account for any variation of the core energy with orientation. This could be manifested by an orientation dependence of the core radius or, equivalently, the Peierls width, of the dislocation. However, the energy contribution due to this source is expected to be small, and current models of the dislocation core are not sufficiently accurate to justify such a refinement. It has already been shown that for the isotropic case the energy contributions due to nonzero tractions across the cores of the reactants and products exactly cancel one another in the reaction.' Accordingly, it will be assumed that this contribution to the total energy change in the anisotropic case is small. In the subsequent discussion it is also assumed that the core radii of the reactant and product dislocation are the same and that, where stacking faults are formed, the faulted region is bounded by the centers of the partials. Consequently only changes in elastic energy due to the reactions will be considered. When the dislocation is parallel to either the (111) or the (100) directions, K may be written:375 K = (Ke sin2 a + Ks cos2 a) [2] where K, and Ks are the combination of elastic constants corresponding to an edge and screw dislocation lying along the same direction as the mixed dislocation, and a is the angle between the direction tangent to the dislocation line and the Burgers vector. Eq. [2] should not be confused with the isotropic approximation to the variation in energy with line Orientation.6 It should be noted that the essentially isotropic expression for K is a result of the characteristic symmetry of the (111) and (100) directions and is not, in general, valid for other dislocation directions in anisotropic cubic metals. The energy* change for a reaction in which the re- actant and product dislocations are parallel perfect dislocations can be written: where Ep and E, refer to the self-energies of the products and reactants, respectively. For dislocations parallel to (100) and (111) directions, Eq. [3] becomes:

Jan 1, 1969

By B. G. Matson, M. A. Whitfield, G. R. Dysart

This paper describes the development of u computer program to calculate changes that occur in the length of tubular goods due to temperature and pressure changes during stimulation operations. Due to the numerous variables involved and the uncertainty of all static and dynamic conditions that could exist, it becomes a staggering task for individuals charged with completions to perform the necessary mathematical calculations. The computer program permits advance calculations for several sets of conditions. INTRODUCTION In the Delaware basin of West Texas alone, 50 wells were contracted or drilled to 15,000 ft or deeper in 1965. Deep well activity is continuing in this and other areas on an expanding scale. Many of these deep wells require extensive stimulation for successful commercial production, and during these operations, pressures and temperatures are encountered that have a pronounced effect on the length of tubular goods. This length change during a large-volume, high-pressure stimulation treatment utilizing fluids considerably cooler than bottom-hole temperature can be of such a magnitude that permanent damage to casing and tubing will result unless mechanical design, pressures and fluid temperatures are evaluated and controlled. These pressure and temperature effects can be calculated. However, the process lends itself well to computer solutions because of the mathematical nature of the problem and the calculating hours involved in arriving at an answer. The engineering-hour demand becomes more severe as tapered strings are involved. On initial treatments on a given well, surface pressure and injection rate conditions are unknown, and offset well conditions have not proven to be a reliable method for making predictions. For these reasons, it has become rather standard procedure for operators to compensate for these uncertainties by placing unnecessary pressure and fluid temperature restrictions on stimulation design. On a number of occasions treating fluids have been preheated to as much as 160F as a means of compensating for thermal contmction resulting from pumping cool fluids. The maintenance of packer seals has been treated by Lubinski, Althouse and Logan',' and the problem of therma1 effects on pipe has been explored by Ramey." These works were expanded and the results made applicable to everyday oilfield terminology before submitting them to computer programming. The pressure and temperature effects on tubing movement previously mentioned occur simultaneously as fluid moves through the pipe. The pressure changes, for purposes of explanation, are categorized here as to the various effects these pressures have on a tubing string. These divisions are (1) the piston-like results of forces acting on horizontal surfaces exposed to pressure, (2) swelling or ballooning of the tubing along its entire length due to the forces of pressure acting against the tubing walls, (3) the elongation of tubing due to frictional drag and (4) corkscrewing of the pipe due to internal pressure. Thermal changes are also of great importance, as their results may be more significant than any of the pressure effects. Steel is an excellent conductor of heat and the earth is a relatively poor conductor. It has been calculated that pipe temperatures at depths of more than 20,000 ft approach within as little as 25" the temperfature of the surface fluid after pumping for 2 hours, or a drop in temperature in some treatments of more than 220F. The equations presented in this paper were developed for computer programming and simplicity of input information; therefore, numerical constants such as Young's modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson's ratio for steel (0.3) are included with unit conversion factors. The moment of inertia of tubing cross-sectional area with respect to its diameter was changed to a constant times (D' — d') where D is outer diameter and d is inner diameter. Units in the equations are length in feet, diameter in inches, density in pounds per gallon, pressure in psi, rate in barrels per minute and time in hours. PISTON-LIKE REACTIONS A change in tubing internal dimensions and the exposure of other horizontal surfaces to different pressures on the inside and outside of the tubing result in a reaction much like a piston under pressure. Such is the case when the internal diameter changes in a combination string of pipe, when seals of a slick joint assembly are subject to pressure and in the end effects of a tubing string. The change in tubing length due to the piston effects of a slick joint packer is affected by the various diameters involved, the tubing pressure Ap,, the casing pressure ,Ap,, length of pipe L, densities of fluid in the tubing before and during pump-

By K. A. Slagle

Hydraulic analysis of the wellbore has become increasingly inzportant for designing cementing operations and selecting equipment, materials and techniques to complenzent modern well-c-ompletion practices. Non-Newtonian fluid technology has advanced beyond the point where former empirical methods of analysis adequately define the hydraulic system and fluid properties. In view of these factors, this paper describes a series of rheological calculations which have been found practical, through field usage, for assistance in selecting a cementing program. A relatively simple laboratory method using standard viscometric equipment is suggested for determination of the rheological properties of slurries, and clrrta are presented on some of the more common cementitrg conzposition.A. A criterion for divergence from laminar-flow characteristics has been proposed. Usefulness of the calculations is indicated by examples of cementing operations where they have been used. INTRODUCTION With the changing aspects of well-completion practices during the past few years, it has been increasingly important to have a relatively simple method of analyzing the flow conditions existing in the well during cementing operations. This is particularly true in view of the improved economics toward which most of the changes have been directed. Rheological characteristics of slurries used for cementing should be a major consideration in the trend toward smaller casing sizes, either single or multiple strings. Receiving increased attention is the practice advocated in 1948 by Howard and Clark' of attaining turbulent flow with the fluids circulated during a primary cementing operation. While there may still be a difference of opinion concerning this technique, most available information indicates that superior primary-cementing results are generally obtained when high displacement rates are employed. Fluid properties of the slurry to be used must be available, as well as calculation methods, to determine what flow rates should be attained and the probable consequences in terms of frictional pressure and horsepower utilization. It would certainly be inappropriate to attempt high displacement velocities if sufficient pressure might be developed to create lost circulation. Since cementing slurries are non-Newtonian fluids, it is not possible to define their rheological or fluid properties by the single factor of viscosity and then make calculations for the quantities just described. Because the shear stress-shear rate ratio is not constant: it becomes necessary to establish at least two parameters for adequate fluid-flow calculations. It is not the purpose of this paper to delve into the mathematical development of non-Newtonian technology, nor to discuss the arbitrary classification system under which a single fluid may resemble two or three different classes depending upon experimental conditions. Rather, the intention is to present a useful series of calculations based on a concept applicable to both Newtonian fluids and to the preponderance of non-Newtonian fluids encountered in the oil-producing industry. Development of this approach was begun some 32 years ago,' and has most recently been brought to fruition by Metzner and his co-workers at the U. of Deleware. Some non-Newtonian fluids encountered in the petroleum industry, other than cementing slurries, have also had the benefit of this method of analysis."' The two parameters required to define the fluid are usually denoted by the symbols n' and K' and, for the purposes of this discussion, are called "flow behavior index" and "consistency index", respectively. These two slurry properties permit calculation of the Reynolds' number and the "critical" velocity, or the velocity at which departure from laminar flow begins. EXPERIMENTAL DETERMINATIONS The two principal instruments used for rheological studies are the pipeline (capillary-tube) viscometer and the rotational viscometer. When conveniently possible, a capillary-tube viscometer (where the pressure drop and flow rate of the material can be measured) is the better method for rigorous determination of the flow behavior index and consistency index for non-Newtonian fluids. With pressure-drop data at various flow rates, it is then possible to prepare a logarithmic plot of shear rate as the abscissa-shear stress as the ordinate. For fluids which do not exhibit time-dependency, these data will usually produce a straight line. The flow behavior index n' represents the slope of this line, while the consistency index K' becomes the intercept of this line at unity shear rate in accordance with the mathematical derivations associated with this concept of rheology. Due to the difficulties anticipated in maintaining a uniform, pumpable cement slurry for the time interval required to obtain measurements from the pipe viscometer, the n' and K' data reported herein were obtained using a direct-indicating rotational viscometer (Fig. 2). The

By Harry C. Swanson

TRANSPORTING concrete from mixer to forms has always been a problem. Twenty-five years ago this task was generally accomplished by means of wheelbarrow or concrete buggy. On large dam jobs, as the number of these projects increased, the gantry crane or highline came into use. Today several methods of handling concrete are employed on smaller surface construction jobs, for example, transit-mix trucks or dumpcrete trucks, which have crawler cranes with buckets for placing concrete into forms. In 1944, during early stages of developing Mather mine A shaft, several large underground concrete jobs were necessary. At this time the Cleveland-Cliffs Iron Co, purchased the first pump-crete machine, introduced by the Chain Belt Co. of Milwaukee. The machine was used to pour approximately 200 cu yd of concrete for a dam, or bulkhead, located 400 ft from the shaft. Concrete was mixed on surface, lowered down the shaft 1000 ft in a 2-cu yd bucket hung under one skip, spouted into the bowl of the pumpcrete machine from the bucket, and pumped directly into the forms. Since the day of the first pipeline concrete in 1944 to the present time, other equipment and other methods have been developed to permit transportation of concrete by pipeline through vertical and horizontal distances totaling 1 mile from mixer to forms. Much of the efficiency in present handling of underground concrete can be credited to the Bethlehem Cornwall mines, where concrete was transported through 6-in. pipe for great distances down an inclined shaft and along levels into forms.' During initial development of Mather mine B shaft, with concrete work under way on two or more levels at one time, the pneumatic concrete placer, Fig. 1, was selected as best adapted for underground concrete transportation. The 3/4-cu yd pneumatic placer is a small machine readily moved from one location in the mine to another. It can be equipped with two sets of mine car wheels, which will permit moving on regular mine tracks. It is therefore possible to send concrete through the pipe at great velocity; the pipeline is clean after each shot except for the film of cement adhering to the inside. With the proper slump in the concrete, it is possible to shoot concrete 2000 ft with this machine, using the mine supply of compressed air at 95 psi. This equipment was first used at Mather mine B shaft to concrete slusher drifts, Figs. 2 and 3, and finger raises located about 2000 ft from the shaft. In several instances there were bends into crosscuts and up vertical distances into the forms. For the first pours two placers were used. The first was located near the shaft where the concrete could be spouted into it from a 2-cu yd concrete bucket on the cage. The second was set on the side of the drift at a point approximately 1500 ft from the shaft. The concrete was shot directly into the second placer from the first unit and from the second machine directly into the forms. After completion of several pours with the two machines, a trial pour with only one placer located at the shaft proved that the second placer could be eliminated. Since then all pours have been successfully completed with only one placer underground. As production of iron ore from the mine increased and the development program expanded, use of the cage for handling mine supplies and concrete became a major problem. This brought about the first attempt at shooting concrete vertically down the shaft for 2600 ft. A 6-in. pipeline with victaulic couplings installed during shaft sinking was used for the trial. One placer was set on surface 250 ft from the collar of the shaft so concrete could be spouted directly into it from the mixer. This machine shot the concrete 250 ft horizontally on surface to the shaft, 2600 ft vertically down the shaft, and 100 ft horizontally into the second placer located near the rib of the shaft station or plat. The second machine shot the batch into the forms, about 2000 ft. Total distance horizontally and vertically was 4800 ft. The entire time cycle for a ¾-cu yd batch of concrete from the mixer on surface to the forms underground totaled about 5 min. During the past two years the two-placer method from the mixer on surface to the forms underground has proved a very efficient means of transporting underground concrete. Advantages of using pipeline concrete are as follows: 1—Interference with normal mining operation is eliminated. When the cage, skips, mine cars, or mine openings are used for transporting concrete and materials used for making concrete, mine operation suffers in one way or another.

Jan 1, 1955

By W. G. Pfann

In zone-melting, a small molten zone or zones traverse a long charge of alloy or impure metal. Consequences of this manner of freezing are examined with impurerespect to solute distribution in the ingot, with particular reference to purification and to prevention of segregation. Results are expressed in terms of the number, size, and direction of travel of the zones, the initial intermsofsolute distribution, and the distribution coefficient. IF a charge of binary solid-solution alloy is melted and then frozen slowly from one end, as for example in the Bridgman method of making single crystals,' coring usually occurs, with a resulting end-to-end variation in concentration. Such coring, or normal segregation, is undesirable where uniformity is an object. On the other hand, for certain systems, it can be utilized to refine a material by concentrating impurities at one end of the ingot.'. ' In the present paper a different manner of freezing will be examined with respect to the distribution of solute in the ingot. A number of procedures will be indicated which have in common the traversal of a relatively long charge of solid alloy by a small molten zone. Such methods will be denoted by the general term zone-,melting, while the process described in the preceding paragraph will be called normal freezing. It will be shown that, in contrast to normal freezing, zone-melting affords wide latitude in possible distributions of solute. Segregation can either be almost entirely eliminated or it can be enhanced so as to provide a high degree of sttparation of solute and solvent. A number of simplifying assumptions will be invoked which, while not entirely realizable in practice, nevertheless provide a suitable point of departure for more refined treatments. Moreover, our own experience with zone-melting has shown that, for certain systems at least, the analysis holds quite well. The present paper will be confined to a discussion of principles and a general description of procedures. Comparison with experiment is planned for later publication. Normal Freezing Before considering zone-melting, segregation during normal freezing will be reviewed briefly. If a cylinder of molten binary alloy is made to freeze from one end as in Fig. 1, there usually will be a segregating action which will concentrate the solute in one or the other end of the ingot. If the constitutional diagram for the system is like that of Fig. 2, then the distribution coefficient k, defined as the ratio of the concentration in the solid to that in the liquid at equilibrium, will be less than one and the solute will be concentrated in the last regions to freeze. If the solute raises the freezing point, then k will be greater than one and the solute will be concentrated in the first regions to freeze. The concentration in the solid as a function of g, the fraction which has solidified, can be expressed by the relation: C = kC0 (1-g)k-1 [I] where C, is the initial solute concentration in the melt. Eq 1 is based on the following assumptions: 1—Diffusion in the solid is negligible. 2—Diffusion in the liquid is complete (i.e., concentration in the liquid is uniform). 3—k is constant. Concentration curves representing eq 1 for k's from 0.01 to 5.0 are plotted in Fig. 3. This equation, in one form or another, has been treated by Gulliver,³ Scheuer,4 Hayes and Chipman5 for alloys and by McFee2 for NaCl crystals. It is derived in Appendix I. It should be pointed out that the k which is calculated from the phase diagram will be valid only in the ideal case for which the stated assumptions are correct. In all actual cases, the effective k will be larger than this value for solutes which lower the melting point, smaller for solutes which raise the melting point, and will probably vary during the beginning of the freezing process. For simplification it will be assumed that the ideal k is valid. Zone-Leveling Processes The processes of this part are designed to produce a uniform, or level, distribution of solute in the ingot. Single Pass: Consider a rod or charge of alloy whose cross-section is constant and whose composition, C2, is constant, although permissibly varying on a microscopic scale." Such a charge might be a rapidly frozen casting or a mixture of crushed or powdered constituents. Cause a molten zone of

Jan 1, 1953

By J. H. Swisher, E. O. Fuchs

Rate equations were derived to describe the kinetics of internal oxidation of cylinders and spheres. The derived equations for cylinders were checked experimentally by means of sub scale thickness and electrical conductivity measurements on Cu-Cr alloy wires. The properties of the internally oxidized samples were examined with conductivity applications in mind. It was possible to produce uniform dispersions of Cr2O3 in copper with an initial chromium content as high as 3 wt pct. While electrical conductivities only a few pct less than that of OFHC copper were obtained, the Cr2Os particle size and spacing were too large for effective dispersion hardening. T.HE process of internal oxidation has been used widely in basic studies of the permeability of gases in metals. In a review article, Rapp1 has discussed the principles of internal oxidation in considerable detail. From a technological standpoint, internal oxidation is often considered undesirable, since it is a means by which inclusions can be introduced into an otherwise clean material. Another important aspect of internal oxidation is its use as a means of dispersion hardening a material. Broutman and Krock2 discuss this and other methods for making dispersion hardened alloys. The only internally oxidized material known to the authors which is commercially available is a Cu-BeO alloy.3'4 This alloy is made from Cu-Be alloy powder, using a so-called Rhines pack. It has a tensile strength of 80,000 psi and retains its strength at relatively high temperatures. The objectives of the present study were to derive rate equations for the internal oxidation of cylinders and spheres, to check the derived equations for cylinders experimentally, and to examine the structure and properties of internally oxidized Cu-Cr alloys. The Cu-Cr system was chosen for this study because uniform dispersions are obtainable at high alloy contents, which is a desirable characteristic in dispersion hardened materials. RATE EQUATIONS FOR VARIOUS GEOMETRIES A number of authors5--9 have derived equations to describe internal oxidation kinetics. These derivations differ somewhat in mathematical assumptions and approximations, and all except one of the derivations deal exclusively with the internal oxidation of plates. The exception is a brief treatment of cylindrical and spherical geometries given by Meijering and Druy-vesteyn9 as a part of a comprehensive paper on the general subject of internal oxidation. These authors did not obtain rate data to check their derivations, although they did show that the hardness profile across an internally oxidized sample is directly related to the rate of interface movement. For cylindrical and spherical geometries, a quasi-steady-state approximation is needed to circumvent mathematical complications in obtaining a solution to the basic differential equations. In using this approximation, we consider the concentration gradient of dissolved oxygen in the internally oxidized zone or sub-scale to be the same as the gradient which would be present if there were no movement of the subscale interface. The steady-state approximation introduces an error of about 1 pct in computing the rate of internal oxidation of an Fe-1.0 pct Mn alloy plate, if the present method is compared to the more exact method of Wagner.7'10 The details of the derivations of the rate equations for cylinders and spheres are given in the Appendix, and only the results of these derivations are given below. The final equations obtained by Meijering and Druyvesteyn9 can be shown to be equivalent to our Eqs. [1] and [2], although the two approaches are somewhat different. Cylindrical Geometry. [2] where r1 is the outer radius of the cylinder or sphere, cm, r2 is the radius of the unreacted core, cm, see Fig. l(a), D is the diffusion coefficient of oxygen in copper, cm2 per sec, %O is the concentration of dissolved oxygen at the surface of the specimen, wt pct, %Cr is the initial chromium concentration in the alloy, wt pct, and t is the reaction time, sec. Plate Geometry. The analogous rate equation for a plate has been derived previously for internal oxidation of Fe-Al alloys.8'11 For Cu-Cr alloys, we may write the same equation as follows: [3] where r1 is the half-thickness of the plate, cm, and r2 is the distance from the mid-plane to the subscale intherate is An analysis of Eqs. [1], [2], and [3] shows that for a plate the rate is completely parabolic. The initial

Jan 1, 1970

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By B. L. Landrum, J. Simmons, J. M. Pinson, P. B. Crawford

Patentiometric model data have been obtained to estimate the effect of vertical fractures on the areas swept after breakthrough in water flooding and miscible displacement programs such as gas cycling where the mobility is near one, The data are presented for the case of the fire-spot pattern in which the cemer well is fractured various lengths and orientations, the data indicate that for 10-acre spacing, fractures extetidirrg over 1300 ft in either directior1 from the fractured well may re.srrlt in reductions in sweep efficiencics from 72 to approximately 34 per cent. However. the area swept after break through may be quite largr and only 10 or 12 per cent 1ess than would be obtained if the reservoir were trot fractured. For the specific case when the volume of fluid injected is equivalent to 100 per cent of the pattern vol-unie, the swent area may vary from 80 to 88 per cent, depending on the lenght of the fracture. The former value is that which occurs when the break through or sweep efficiency was orrly 34 per cent and the latter figrrre of 88 per cent is that which is obtained if the reservoir were unfrac-ttm'd. It is pointed out that although the sweep efficiency may he very low in vertically fractured five-spot patterrz.s, the area swept at low water-oil ratios may be only 5 to 10 per cent less than those achieved if the reservoir were unfractured. INTRODUCTION Since the initiation of commercial reservoir fracturing techniques it has been desirable to determine the effect of fractures on the areas swept after breakthrough. Most water flooding or gas cycling projects are continued for substantial periods after the brcakthrough of the injected fluid. Although the sweep efficiency serves as one criterion for rating various flooding patterns. the area swept after breakthrough for various water-oil ratios or percentage wet gas, if cycling. is of perhaps more importance than the sweep efficiency alone. Sweep efficiency data on the vertically fractured five-spot have been presented3. Previous work on the line-drive pattern has shown the effect of vertical fractures on the area swept after breakthrough for the case in which the distance between injection and producing wells divided by the distance between adjacent input wells was equivalent to 1.5 (see lief. 2). The data indicated that for the line-drive pattern it may be desirable to flood or cycle substantially perpendicular to the fractures in order to achieve the greatest recovery for the smallest volume of fluid injected. For this study the center well of a five-spot is assumed as the fractured well. All fractures were assumed to originate at this well and extend into the reservoir for various distances and orientations. All the fractures are straight and are of large permeability compared to the matrix proper. These data are presented to aid the engineer in estimating fractured five-spot pattern performance. ANALOGY The potentiometric model was used in making this study. The model used was 20 20 in. by approximately 1-in. deep. For certain portions of the study one corner of this model was considered to be an injection well and the opposite corner a production well. To simulate vertical fractures a copper sheet was soldered to the wire well and made to conform to the desired length and orientation. In other studies the same model was used except that the four corners of the model might be considered as the corner wells of a five-spot pattern and a fifth well was placed in the center of the model. The well placed in the center of the model was fractured. The total fracture length is L and the well spacing. d. The complimentary fracture angles will be obvious from Figs. 3 and 4. The data obtained on the potentio-metric model assumes the pay to be uniform and homogeneous, the mobility ratio is one, steady-state conditions exist and gravity effects arc neglected. The permeability of the fractures is very great compared to that of the matrix proper. The po-tentiometric model has been used widely both in water flooding and gas cycling projects, and may be used for miscible displacement; how-ever. it is believed that the poten-tiometric model data are more properly applicable to gas cycling than water flooding because the model as-

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.' ' Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.'" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By J. B. Conway

Stress-vuptuve data for several of- the refractory metals are frequently found to yield a linear relationship between the Larson-Miller parameter and the logarithm of the applied stress. In such cases linear stress-rupture isotherms result with slopes bearing a definite relationship to the temperature. It also follows that the stress to produce rupture in a certain period of time will be linear in temperature. Data for several refractory metals have been reviewed and excellent linearity is shown in this type of isochronal plot. In addition, the af ore - mentioned lineavity leads to a linear relation between the log of the stress to produce rupture in a certain time and the homologous temperature. This has been illustrated for the Group VI-A metals, tungsten and molybdenum. EXTENSIVE use has been made of the Larson-Miller' parameter in the interpolation and extrapolation of stress-rupture and creep data. In those cases where this particular parametric approach is applicable a convenient and fairly straightforward procedure is made available for the correlation of experimental stress-rupture data. It is quite common to employ this parameter in the form of a master rupture plot in which the parameter, T(C + log tr), is expressed as a function of log stress. In many cases this functional relationship in log stress is linear within acceptable accuracy and hence the following relation results: where P is the parameter, C is the Larson-Miller constant, T is the absolute temperature, t~ is the rupture time, a is the stress, and a and b are constants. Examples of such a relationship are contained in the work of Green, Smith, and 01son2 dealing with high-temperature rupture behavior of molybdenum and in the work of Green' dealing with the high-temperature behavior of tungsten. In addition, pugh4 has shown a similar linearity for some fairly low-temperature data for molybdenum. It can be shown that when the relationship in Eq. [I] is exhibited certain generalizations can be made concerning the form of the stress-rupture isotherms. For example, rearranging yields: For a given material (constant C) at a given temperature the first term on the right-hand side of Eq. [2] is a constant and hence this equation defines a straight line when log stress is plotted as a function of log-rupture time. This is recognized as the standard form usually employed in this type of data presentation. Such linearity then suggests the linear form of the Larson-Miller parameter. Or, in other words, the linear parametric relationship in Eq. [2] results only when the stress-rupture data are linear on a log-log plot of stress vs rupture time. Another interesting observation can be made in regard to Eq. [2]. It can be noted that the slope of the stress-rupture isotherms is given by - T/b and hence a direct calculation of the constant b is available. It also follows that since the value of b is the same for all temperatures the slopes of the various isotherms on the log-log stress-rupture plot cannot be the same. Indeed, the existence of the relationship in Eq. [2] precludes a system of parallel lines on this common stress-rupture plot. As a matter of fact it further specifies that in addition to being nonparallel the slope of these isotherms must decrease (i.e., become more negative) with increasing temperature. Such a condition is indeed found to exist in the case of the stress-rupture data reported for molybdenum.' As a corollary to the above, it may be stated that stress-rupture data which do not lead to a linear log-log stress-rupture plot or whose isotherms do not exhibit a decrease in slope as the temperature increases will not yield the linear relationship of Eq. [I]. Applying Eq. [2] to two different temperatures and solving for C yields: Eq. [3] affords a simple and rapid method for calculating the Larson-Miller constant from the log-log stress-rupture plot. The slope of a given linear isotherm is measured and the value of "b" calculated based on Eq. [2] as: slope = - -Tb Then at an abscissa value of 1.0 hr (making log tr in Eq. [3] equal to zero) read the stress corresponding to rupture for two different temperatures. Substitution in [3] yields: A value of the Larson-Miller constant can thus be calculated from a few simple mathematical procedures employing data read directly from the log-log plot of the stress-rupture data. Of course, it is not to be overlooked that the above reasoning has been based on the linear relationship of Eq. [I] being applicable. However, if as mentioned above the log-log plot is

Jan 1, 1967

By S. Storchheim

MORE and more attention is being paid to the bonding of metals in their solid states. For a better understanding of this technique for joining metals and how it is affected by changes in temperature, pressure, and time at temperature and pressure, a detailed report concerning nickel to aluminum bonding has been published.' In order to broaden the knowledge accrued, some additional work concerning solid state joining of aluminum to copper and aluminum to zirconium was performed. The investigation of the Al-Cu system was considerably more extensive than the investigation of the Al-Zr system. For the A1-Cu system, not only were tensile sudies made but intermetallic penetration rate investigations also were carried out. The effect of temperature on intermetallic penetration rate for the A1-Cu system was determined at 11 tsi pressure, held 2 min. Procedure Apparatus: The hot pressing technique was the means of solid state reaction used and required the equipment depicted in Fig. 1. The following procedure was involved: The two metals to be reacted were placed in an aquadag-lubricated 18-4-1 tool steel die, 16 in. high by 1.440 in. ID, between punches of 1.366 in. diam made of the same material. A thermocouple well was located in the die body 3½ in. down from the top of the die, while another well was located centrally in the bottom punch 8½ in. from the bottom of the die. This die assembly was located in three cylindrical ceramic heating furnaces placed in tandem. Each furnace was controlled individually by a Variac power transformer. In turn, the die and furnaces were placed in a water-cooled stainless steel pot which could be evacuated. A cover, which contained a centrally located Wilson seal with an 18-4-1 1 in. diam ram running through it, was bolted on the pot. After sealing, the pot was evacuated by a roughing pump to 200 microns pressure, after which a diffusion pump was used to bring the pressure down to 5 to 15 u. At this pressure, the furnaces were turned on. As soon as they started to heat, out-gassing of the entire unit raised the pressure to 30 to 400 p. By the time the specimens were at temperature ready to be pressed, approximately 4/2 hr, the vacuum pumps had re-established the 5 to 15 u pressure. Once the desired temperature was reached, the required pressure was applied for a predetermined length of time to the 1 in. ram, through to the top punch, and to the specimen. When the time for keeping the specimen under pressure had elapsed, the pressure was released, the energizing coil current turned off, and the assembly allowed to cool. After cooling, the die was removed from the pot and the specimen was ejected. Specimen Preparation: Two different types of specimens were made for this investigation. One was for subsequent tensile testing, while the other was for determination of intermetallic alloy zone penetration into the parent metals. Tensile Bars—Commercially rolled copper pieces in. thick or zirconium sheet pieces 1/32 in. thick and 1.366 in. diam were placed between commercial 2-S aluminum rod 1 in. thick and 1.366 in. diam. This sandwich in turn was slipped into a 2-S aluminum sleeve 1.438 in. OD and 1.370 in. ID. This sleeve lined the couple up and prevented the aquadag lubricant from getting in between either the A1-Cu or Al-Zr interfaces. Immediately prior to the specimen assembly, the copper or zirconium was abraded on the flat surfaces with 320 grit silicon carbide paper, producing clean smooth surfaces. The aluminum was chemically cleaned just before assembly by: l—degrease in acetone, 2—distilled water rinse, 3—immersion for 3 min in 5 pct NaOH at 70" to 80°C, 4—distilled water rinse, 5—immersion for 2 min in 50 pct HNO3 solution at room temperature, 6—distilled water rinse, and 7—drying in a blast of gas. After the A1-Cu sandwiches were hot pressed and ejected, the specimens were machined such that the aluminum sleeve was removed, and the remaining aluminum was then threaded; the bar so produced was tested later for tensile strength. In all the instances where Al-Cu couples were tested, the specimens broke during the test at the Cu-A1 interface and never within the aluminum or copper. The ultimate tensile strength values at times showed considerable scatter for a set of given reaction conditions. Because of this, as many as three to five specimens were made for a particular set of conditions. The trend of the average tensile strengths obtained was not as conclusive as was the trend of the maximum tensile strengths, the latter values being obtained under optimum reaction conditions. Therefore, the values of ultimate tensile strength given herein are maximums.

Jan 1, 1956

By H. C. Brown

SINCE 1934 the steel industry has been utilizing the spectrograph for supplementing wet chemical analysis in the production control of electric and open hearth furnaces. This means of control made great strides during the war years because of the general acceptance of the spectrograph and the increased emphasis that was placed on rapid control methods. However, in the post war era, with the demand still on increased production, it became apparent that a still more rapid and economical means of production control was needed. Since the spectrograph had been used mostly in the analysis of low alloying and residual elements, it also became apparent that equipment was needed to extend the spectrographic technique to the analysis of the high alloying elements in stainless steel. For these reasons, companies manufacturing spectrographic equipment were prompted to start development work on direct reading instruments. In June 1949, the Applied Research Laboratories of Glendale, Calif., announced that a direct method of spectrochemical analysis for stainless type steels had been developed. This paper will describe the use of the Applied Research Laboratories Production Control Quantometer in the quantitative control of stainless, silicon, and plain carbon steels being made at the Butler Pennsylvania Div. of the Armco Steel Corp. The Armco Butler Div. has one 70-ton electric furnace and six 150-ton open hearth furnaces. The electric furnace is employed in the making of all types of stainless steel and the open hearth furnaces are used for the production of silicon, wheel, and plain carbon steels. The ARL quantometer was purchased primarily for the purpose of controlling the steelmaking in the electric furnace, but its use has been extended for the analysis of final tests (ladle tests) on a number of different types of stainless, silicon, and plain carbon steels. Because of this additional work by the quantometer, substantial savings in manpower and time have been realized by the laboratory. In the analysis of a set of preliminary tests from the stainless steel furnace, approximately 40 min in laboratory time are saved due to quantometric analyses. Despite the fact that more specialty grades of stainless steel are being made in the electric furnace, the average tons per hour have been increased since the quantometer was put into operation. Specialty grades require more furnace time than regular commodity grades of stainless steel. The installation of the ARL production control quantometer was completed on March 13, 1952. By May 1, 1952, the instrument was calibrated for nickel, chromium, manganese, silicon, and molybdenum, which are the elements necessary for the production control of the stainless steel furnace. Within the following month, training of personnel on the quantometer was achieved and a study of the accuracy of the instrument showed that the results obtained were sufficiently accurate for control purposes. Therefore, on June 11, 1952, the quantometer was placed on production control for all types of stainless steels. Starting September 11, 1952, the instrument was gradually placed on ladle analysis (final tests) as the analytical curves were refined and additional curves were drawn. The quantometer has been relatively free of breakdowns since placing it on production control. The samples from only one stainless steel heat have had to be analyzed by wet chemistry because of instrument trouble. The previously existing heat-time record was also bettered by 15 min on a commodity grade of 18-8 stainless steel. Scope of Control In general, the quantometer determines all elements necessary for the production control of all types of stainless steel heats and for the ladle analysis of various types of stainless steel heats. It is also used in reporting final results for silicon, manganese, chromium, nickel, molybdenum, tin, copper, and aluminum on all silicon steel grades and manganese, chromium, nickel, molybdenum, tin, and copper on several plain carbon steel grades. Table I shows the elements and the concentration ranges of these elements in the various types of stainless, silicon, and plain carbon steel that are determined on the quantometer. A study of the results obtained on ladle test samples of stainless steel types 410, 430, 430 Ti, 446, 301, 302, 304, 304L, 305, and 17-7 PH will be discussed. Also included in the study are the results obtained on ladle test samples of a number of silicon steels. Apparatus In order to take full advantage of the potentials of the production control quantometer, the unit has been placed in an air-conditioned room with relative humidity control. The temperature is maintained at 73'22°F and the humidity at 45&5 pct. The air conditioning serves as a precaution to minimize the amount of adjustment and calibration needed during operation. It also reduces contaminating fumes and dust and thereby lessens the necessity for maintenance on the equipment. The quantometer is composed of three units: the high precision multisource unit, the 1.5 meter vertical spectrometer, and the console. The source unit supplies excitation conditions varying from spark-like discharges to arc-like discharges. The voltage to the source unit is supplied by a motor-generator

Jan 1, 1955

By Benny Langston, Frank M. Stephens

Although little information is available concerning small-scale equipment for dust collection in laboratories, it is possible to adapt standard equipment for laboratory use. Dust from laboratory processes may be collected by cyclone separators, filters, electrostatic separators, scrubbers, and settling chambers. IN recent years much attention has been given to recovery, treatment, and disposal of dusts discharged into the atmosphere from operations of industry. considerable data has been accumulated on both operation and design of dust-collector equipment for commercial installations. On the other hand, there is almost no published data on design and construction of small-scale equipment to handle dust problems that arise in the ore-dressing laboratory. Dust-collection equipment such as multiclones, single-cyclones, scrubbers, chemical and mechanical filters, settling chambers, and electrostatic separators has proved its efficiency for collecting dust in industry. In the laboratory, however, the engineer is faced with the problem of collecting small quantities of dust, inexpensively, without diverting the major effort from the metallurgical problem to the problem of collecting dust produced by the process. For most applications standard dust-collection equipment is too large for use in the laboratory; however, for control of dust in large working areas it is often possible to use a standard dust collector, such as an air filter, with branch ducts to eliminate a health hazard. For example, the well-furnished sample-preparation room containing small jaw crushers, rolls, and pulverizers, in addition to the riffles and screens necessary for preparation of samples, presents a perennial source of dust. The authors' experience has shown that a combination system consisting of overhead branch ducts to the individual equipment and floor ducts with grills, where applicable, connected to a central dust collector effectively removes dust generated in preparation of samples. Fig. 1 is a sketch of a downdraft dust-collector for table installation. Similar systems can be built with floor grids. For portable equipment such as laboratory vibrating screens this type of installation with a steel grill to support the heavy load is reasonably efficient. Overhead branch ducts to individual crushing and grinding equipment, although efficient, must be carefully controlled by dampers to prevent excess loss or a change in the composition of the sample. Change in sample composition can result from excess velocity, which causes selective removal of constituents of low specific gravity. Fig. 2' shows the theoretical effect of terminal velocity on spherical particles of different specific gravities in air and water under action of gravity. Fig. 3 shows the effect of air velocity on composition of CaCO, coal mixtures. Although the entrainment of dust particles in a moving air stream is the basic mechanism by which all dust-collection equipment functions, usually intake velocity of the dust-collection system must be controlled to prevent loss of part of the sample. As an example of what may happen when excess velocities are used, a mixture of 50 pct coal and 50 pct limestone was crushed to —10 mesh and fed to a pulverizer equipped with an overhead dust-collection system. Fig. 4 shows the overhead dust-collection equipment used in this test. The pulverizer was set to give a product 95 pct —100 mesh in two stages. Velocity of air passing over the lip of the pulverizer was measured with an anemometer. After grinding, the finished product was analyzed to show the amount of calcium carbonate present. Fig. 3 shows graphically the increase in calcium carbonate as velocity through the dust-collection duct was increased. These data show that at a velocity of 1 ft per sec little if any of the coal was entrained by the overhead draft. At the maximum velocity, about 6.5 ft per sec, approximately 7 pct more coal was entrained than calcium carbonate. From an operating standpoint, this problem can be remedied easily by dampering the line to control velocity. The lowest velocity commensurate with satisfactory dust control should be used to prevent excess loss and, in some cases, selective dust loss. Collection of Dust in Laboratory Processes As in industry, the engineer desires to collect efficiently the dust produced by processes being investigated on a laboratory scale. However, in the collection of laboratory dusts he is faced with two additional problems: 1—The volumes of gas and the quantity of dust that must be recovered are small when compared with the capacity of standard dust-collector equipment, which must be scaled down in design except for collection of dust from large pilot-plant operations. 2—In addition, because of the variety of problems studied in the process laboratory, the engineer cannot design today a dust collector that will meet the conditions imposed by the processes of tomorrow. Sometimes, therefore, he must compromise collection efficiency to minimize the cost of fabrication and the amount of time diverted from the metallurgical to the dust-control problem. For collection of dust from laboratory processes a cyclone separator, filters, electrostatic separators, scrubbers, and settling chambers can usually be adapted for small-scale operations. The following

Jan 1, 1955

By D. G. Westlake

Polycrystalline tensile specimens of various Zr-0-H alloys have been tested at 298°, 178°, and 77°K. Solute oxygen and hydride precipitates in quenched alloys made individual contributions to the yield strength at 0.2 pct strain which combined to produce a resultant strength increment, a,., Ductility changes which were ohserved can he interpreted in terms of the various oxygen and hydrogen concentrations, testing tem -peratures, and dispositions of the hydride. ADDITIONS of oxygen in solid solution were known to increase the yield and tensile strengths of polycrystalline zirconium as early as 1951.' More recently, the critical resolved shear stress (CRSS) for prism slip in zirconium single crystals was also shown to be affected by the solute oxygen impurity.' This latter work also demonstrated that large increments of strength could be contributed by the finely dispersed zirconium hydride precipitates that are present in quenched Zr-H alloys.3 It was concluded that the combined strengthening due to alloying could be expressed by where to is the increase in the CRSS due to solute oxygen alone and TH is the increase due to finely dispersed hydride precipitates. Eq. [I] is analogous to one used to express the combined strengthening effects of work hardening and neutron radiation damage.4 Eq. [1] was verified only indirectly and for only small amounts of the impurities—up to 0.14 at. pct 0 and 0.63 at. pct H. The present investigation was undertaken to obtain a more direct verification of the validity of the form of Eq. [1] for this system and also to determine the combined effects of oxygen and finely dispersed hydride precipitates on the tensile strength and ductility of polycrystalline zirconium. EXPERIMENTAL PROCEDURE Tensile specimens were machined from the same rolled billet of Kroll zirconium used in the earlier study.' These measured 38 by 4.7 by 0.5 mm and had 10-mm gage lengths which were 2.8 by 0.5 mm. Each specimen was ß-annealed in vacuo at 1173°K for 15.5 hr and a-annealed at 1073°K for 4 hr to D. G. WESTLAKE, Member AIME, is Associate Metal l ur-gist, Metallurgy Division, Argonne National Laboratory, Argonne, III. Manuscript submitted July 17, 1964. IMD______________ give an equiaxed structure with grain diameters averaging 0.06 mm. Oxygen was added by allowing the metal to react with a known quantity of oxygen during the 0 anneal and known quantities of hydrogen were added during the a anneal. Each alloy was encapsulated in Pyrex under vacuum, annealed at 873°K for 4 hr, quenched into ice water, and polished by immersion in a solution of 46.75 vol pct H2O, 46.75 vol pct concentrated HNO3, and 6.5 vol pct HF (49 pct) at 298°K. Special heat treatments given to a few specimens are described in the results below. Tensile tests were done on an Instron machine and were begun within 20 min after the quench, except where specified otherwise. Tests at 298°K were in air, at 178°K in acetone, and at 77°K in liquid nitrogen. All tests were at a strain rate of 8x sec-1. RESULTS AND DISCUSSION Yield Stress at 298°K. The compositions of alloys and the corresponding yield stresses (0.2 pct strain) are given in Table I. A plot of the yield stresses of the oxygen alloys, A, B, C, and D, indicates that varies linearly with CO1/2, where Co is the oxygen concentration, Fig. 1. This is in accord with Fleischer's6 theory for solution strengthening if the oxygen atoms do not cluster, or the cluster size remains constant with increasing oxygen concentration. In Fig. 1, it appears that if one could prepare some oxygen-free zirconium its yield stress would be very low. Therefore, we shall assume that for the oxygen alloys is equivalent to O0, the strength increment contributed by the presence of oxygen. The relationship between0.2and Co is expressed by 0.2 = 31.3 CO1/2, when the yield stress is in kg per sq mm and the concentration is in at. pct. Each of the hydrogen alloys, Al, A2, A3, and A4, contained 0.081 at. pct 0 as an impurity. In Fig. 1, it appears that this small amount of oxygen makes a significant contribution to the strength which cannot be ignored when we evaluate the contribution of the finely dispersed hydride. Let us assume the validity of the following equation: a0.2 = (a2o+a2R)1/2 [2] which is analogous to Eq. [I] for single crystals, and calculate values of UH for the hydrogen alloys by using the experimental values of 0.2 and o (0.081 at. pct) = 8.9 kg per sq mm. For 0.36 at. pct H, oH = 6.47; for 0.72 at. pct H, OH = 11.30; for 2.16 at. pct H, OH = 19.4; and for 3.60 at. pct H,

Jan 1, 1965

By G. H. Miller

With the recent advances which have been made in science and technology and the increased use of mathematics in this area, the question of the best mathematics courses for a mining engineer to take is of major importance. The question becomes even more difficult to answer due to the recent increase in the number of different mathematics courses in the last two decades offered by the mathematics departments. Therefore, the National Study of Mathematics Requirements for Scientists and Engineers (NSMRSE) was designed to provide some answers to these questions. Approximately 10,000 scientists and engineers were selected for the Study, These individuals were placed in two categories: (1) The Awards Group, recipients of national honors or awards and those recommended by the members of the Board of Advisors as having national and international reputations in their areas of specialization and (2) The Abstracts Group, persons exceptionally productive in their research, based on the number of journal articles listed in the last five years in the Engineering Index, Scientific and Technological Aerospace Reports, Chemical Abstracts, Biological Abstracts, and the Physics Abstracts. The NSMRSE Course Recommendation Form and the Instruction and Course Content Sheet were constructed with the aid of the Board of Advisors and other consultants. For the Study, 40 courses were selected by the mathematical consultants. In order to make sure that the basic content of the mathematics courses was the same for all respondents, a brief resume of each of the 40 courses was given. The NSMRSE Course Recommendation Form consisted of seven sections. These sections were as follows: Section 1, 38 different specializations; Section 2, orientation of work (applied through theoretical); Section 3, highest degree obtained; Section 4, place of employment (academic, nonacademic); Section 5, administrative capacity (administrative or nonadministrative); Section 6, age groups (five-year intervals). Section 7 contained the 40 courses which were to be marked according to five categories: (1) Course Length, 3 to 36 weeks; (2) Applied-Theoretical Orientation, a five-point scale; (3) Course Level, freshman through graduate; (4) Knowledge of Course; and (5) Use of Course Content in Work. The Analysis The report of the data is based on the replies received from 44 mining engineers. This group was part of the Awards and Abstracts Group for all engineers. The resume of the recommended courses is reported in quintiles (upper fifth to lower fifth), since recommendations of this kind are not precise. The results of the Study are based on recommendations of the most active research men in engineering in the U.S. today; therefore, the reader should realize that these course recommendations represent an upper bound of mathematics requirements for the Ph.D. in both undergraduate and graduate work. Conclusions and Recommendations 1) Mining engineering students who plan to be active research specialists should take all those courses which are "very highly recommended" (80-10070) and "highly recommended" (60-79.9%). Those courses in the upper two quintiles and recommended by most mining engineers are: first-year calculus, third-semester calculus, elementary differential equations, applied statistics, and machine computation. Courses of "moderate recommendation" (40-59.9%) are: vectors, intermediate ordinary differential equations, the first course in partial differential equations, elementary probability, and linear programming. 2) The great majority of mining engineers indicated that they prefer a course which is concerned primarily with applications. Only the standard courses such as first-year college mathematics, calculus, differential equations, and advanced calculus received a recommendation for 50% theory and 50% practice. Therefore, all mathematics courses given to mining engineers should contain many applications and little theory. Engineers in both the applied and the combination (ap-plied-theoretical) groups indicated a definite need for applications in all courses. 3) In general, recommendations were for mathematics courses to be given for short intervals of time such as 3, 6, or 12 weeks. Only the standard courses mentioned previously received the usual one-semester or one-year recommendation. Therefore, it is of value to combine several related courses into a one or two-semester course so that the mining engineering student could acquire important mathematical knowledge at an early date in order to prepare him for his research. 4) There was little use for the newer courses in modern mathematics such as the functional analysis sequence, the modern algebra sequence, and the group theory sequence. In addition, there were uniformly very low recommendations (0-19.9%) for multilinear algebra, complex variables, mathematical logic, special functions, integral equations, approximation theory, analytic mechanics, integral transforms, and geometric algebra. Therefore, these courses should be given a low priority. 5a) Comparisons among mining engineers with different backgrounds showed that the combination ap-plied-theoretical group recommended more mathematics than the applied group. 5b) There was little difference in recommendations between the administrative group and the nonadminis-trative group. 5c) Analysis of age groups showed that those in the lower age groups gave significantly higher recommendations to courses such as the first course in partial dif-

Jan 1, 1971

By J Price, M. R. Geer, H. F. Yancey

COAL mine operators recognize coal as an abrasive material, because the wear of drilling, cutting, and conveying equipment is reflected as a cost item for replacement of parts. Similarly, industrial consumers of coal experience abrasive wear on all coal-handling equipment. Operators of pulverized fuel plants are doubtless most keenly aware of the abrasiveness of coal, because under the high contact pressures developed between coal and metal in pulverizers, abrasive wear is increased many fold. Moreover, experience in operating pulverized fuel plants has demonstrated that some coals are much more abrasive than others. Hardgrove' stated that maintenance costs entailed by the wear of grinding elements is often a more important variable than the cost of the power required to pulverize different coals. Craig2 also reports that one coal may cause pulverizer parts to wear several times faster than another. It is apparent, therefore, that those concerned with pulverizing coal could profitably employ a method for estimating the abrasiveness of different coals, just as they utilize standard tests for thermal value, grindability, and ash-fusion temperature to assist in selecting the most suitable and economical coal to use in a particular plant. The objective of this investigation was to develop a test procedure that would be suitable for general use in estimating the abrasiveness of coals. However, few, if any, of the standard tests now used for evaluating the properties of coal are the product of a single investigation or the result of a single investigator's efforts. Rather, in each case, a testing procedure was devised by one investigator, used by others on a wider variety of coals, and finally refined completely as the result of the joint efforts of a number of interested people. Thus, the test procedure for estimating abrasiveness developed in the course of this work may not be refined sufficiently in its present form for general use, but it may serve as the starting point from which an acceptable test procedure can be developed. The method has been used thus far on only about a dozen coals, and there has been no opportunity to attempt a correlation between experimental results and actual plant experience. Only wider use of the procedure by other investigators and correlation with plant experience can determine to what extent the method will have to be modified to render it suitable for general application. Test Method Although the literature contains no record of an attempt to devise a method for estimating the abrasiveness of coal that could be used industrially, several investigators have tested properties of coal that are closely related to its abrasiveness. The abrasiveness of a material generally is considered to be related to its hardness, and hardness tests for coal have been employed by Heywood,' O'Neill," and Mathes. Also, the resistance of coal to abrasion, a property that presumably is related to the abrasiveness of coal, was measured by Heywooda and by Simek, Pulkrabek, and Coufalik.2 11 these investigators tested only individual pieces of coal. Since coal is a heterogeneous material having components of varying properties, tests of this type can yield results having little more than academic interest. Only a test method that utilizes a representative sample of coal can give results that are useful industrially. The abrasion tests used for various other materials have been considered for adaptation to testing the abrasiveness of coal. The tests used for metals,7-9 paving and flooring,'" and rubber," cannot be used because coal is not sufficiently abrasive.~ The present experimental work was begun before World War II and was conducted by three research fellows"'" working under a joint agreement between the University of Washington and the Bureau of Mines. After a great deal of preliminary work with a variety of apparatus and materials, a test procedure was developed which consisted of rotating a test disk 2Yz in. diam in a steel mortar containing the coal sample. The shaft carrying the test disk at the lower end and a 100-lb load on the upper end was free to move vertically. The bed of coal in the mortar was kept fluid by low-pressure air admitted through a port near the bottom of the mortar. Measurable wear on an Armco iron disk could be obtained in this test procedure, but, despite extensive efforts to eliminate them, several major disadvantages remained in this test method. First, with most coals the amount of wear on the iron disk did not exceed a few milligrams. Second, a single type of disk was not applicable for all coals. A smooth iron disk gave satisfactory results with both bituminous and sub-bituminous coals, but hardly any wear with anthracite or coke. A disk having studs or projections gave more satisfactory abrasion losses with anthracite and coke and presented no operating difficulties with free-burning bituminous and sub-bituminous coals. It could not, however, be used with caking coals because these coals formed a

Jan 1, 1952