Search Documents

By K. E. Gray, M. S. Seth

Equations of elasticity for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure were formulated. Using an analogy between thermal and porous body stresses, it was shown that the solution for a transversely isotropic porous body may be obtained by incorporating body forces and the stresses due to a boundary load into the corresponding solution for the thermal stress problem. Equations of elasticity ad the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant pressure or a constant rate of flow is maintained at the wellbore. The vertical and tangential displacements are zero from the conditions of the problem. A numerical analysis was made of the solutions obtained by using a digital computer to determine the relative influence of each system variable. Considering rock as a porous body with internal fluid pressure generally gives results significantly different than considering the rock to be nonporous; the directional character of rocks leads to significant differences as compared to results based upon the common assumption of isotropy. Stress gradients are high near the wellbore but die out away from the well. Radial stresses are compressive or neutral, whereas tangential stresses are tensile, neutral or compressive, depending upon the boundary conditions and physical properties of the system Vertical stresses are compressive for an unbounded system. For constant wellbore injection rate, the vertical stress is proportional to the rate of fluid injection and decreases with time, whereas the radial and tangential stresses increase with time. At a given location, the radial displacement generally is very dependent upon time. INTRODUCTION A realistic appraisal of the state of stress in subsurface rock formations would be of considerable interest and use to the petroleum industry. For example, knowing the state of stress in proximity to a wellbore would be of fundamental importance in designing a fracturing operation or, more important, of clearly understanding the conditions necessary to produce rock failure of desired dimensions and geometry. Understanding conditions necessary for rock failure at the wellbore would also be of utility in a preventive sense. For example, borehole stability is an important consideration for many rock formations, and knowledge of the stress state at and near the wellbore under conditions of substantial pressure gradient due to fluid flow would be of great value. When fluid flows through a porous body which is initially at some uniform stress level, the following forces generate stresses at any point in the body. (1) Forces due to nonuniform pressure distribution. With increasing pressure the elements of a body are compressed. Such compression cannot proceed freely in a continuum when the pressure is not uniform throughout, and thus, stresses due to flow of fluid are set up. (2) Pore fluid pressure. This gives rise to normal stresses whose value at any point is the product of areal porosity and the fluid pressure. Although the fluid exerts uniform pressure, the stresses it creates in an anisotropic body may not be the same in all directions since the areal porosity in an anisotropic porous body is a direction-dependent quantity. This consideration leads to the concept of directional porosity. (3) Body Forces. The flow of the moving fluid imparts a force to the rock matrix, the magnitude of which depends upon the pressure

Jan 1, 1969

The income of the Institute is derived mainly from dues, advertising in MINING AND METALLURGY and sale of publications. These sources are fortunately supplemented by the interest from invested funds now amounting to over $500,000, a considerable portion constituting endowments for especial purposes. Aside from the medal funds the principal endowments are the James Douglas Fund for support of the Library, the Rocky Mountain Fund and the Seeley W. Mudd Memorial Fund. The income from the two last named is available for the support of research and a variety of other purposes. It is allocated by the Board of Directors upon recommendation of standing committees in each case. An active effort is being made to increase the endowments of the Institute, it being felt that unusual opportunities exist for research and public and professional service through the organization.

Jan 1, 1939

By Michael A. Toole, O&apos

INTRODUCTION The tremendous capital investment required to produce a low grade ore deposit demands a reliable answer to the question: "How much does it cost to drill a well to produce the uranium that geologists have indicated is there in the ground?" However, you will find the answers given will be many and varied, depending on whether they come from the operating company, the drilling contractor, a geologist, a reservoir engineer, purchasing agent, or whoever. Each generally considers only a limited area of the total operation. The operating company usually depends upon the drilling contractors, service companies, and consultants to supply expertise. Because it is often difficult to see the "whole picture" or even to agree upon what the "whole picture" is, planning a well program and its costs is often done piecemeal. Frequently, costs "saved" in one area of the program are needlessly spent in another, because the effect of one area upon the other was overlooked. Since the drilling fluid is such an influential part of the drilling program, it should be given utmost consideration when planning the overall program. As W. D. Lacabanne wrote in 1954: "The drilling fluid system is understandably called the heart of the rotary oil drilling rig. Any other type of rotary rig.... should benefit by the incorporation of mud fluids in the drilling scheme." And G. R. Gray mentions that "The driller recognizes the drilling fluid as one of the useful tools available to solve drilling problems." However, in minerals drilling, only in recent years has the drilling fluid been considered as more than a tool to get through special problem areas. Although there are many similarities between drilling for oil and gas and drilling for minerals, the differences in the drilling equipment used justifies designing specific fluids for the minerals drilling industry. NL Baroid/NL Industries, Inc. has been a leader in introducing such fluids and in providing technical know-how to the minerals drilling industry. The purpose of this paper is: (1) to discuss the selection of drilling fluids to meet specific drilling conditions during both exploration and production and (2) to show the interrelationship among factors present in the exploration and production phases that influence total well costs. DRILLING FLUIDS FOR DRILLING PROBLEMS The drilling fluid is a tool that can be used to improve drilling performance by improving hole cutting, cleaning, stability, and formation productivity. Properly formulated and maintained drilling fluids enable the drilling operation to be carried out with increased efficiency and lower total (overall) costs. However, it should be noted that not all drilling problems can be solved by even the most carefully prepared and maintained drilling fluids. Of the many possible drilling problems encountered in a well before reaching target depth, this paper will discuss only those most likely to be present in the shallow [61 to 610 m (200 to 2000 ft)] drilling operations in South Texas. LOST CIRCULATION Loss of circulation is the most common problem encountered in drilling. Because the losses occur under varied conditions, it is often difficult to determine the exact causes. "Lost circulation" or "lost returns" means the partial or complete loss of drilling fluid to voids in the formation. "Loss of water" while drilling with water may take place into any permeable section and should be distinguished from "water loss" or filtration of fluid through the filter cake of mud solids laid down on a permeable formation. "Loss of water" can frequently be stopped by the addition of colloidal sized clay particles such as high yield bentonites, whereas "water loss" may be controlled with organic polymers. Subsurface conditions that lead to loss of circulation can be classified as: (1) natural fractures, (2) induced fractures, (3) unconsolidated or highly permeable formations (loose gravel), and (4) cavernous formations (crevices and channels). Loss of circulation may occur whenever the borehole pressure exceeds the formation pressure. The greater the differential pressure, the more likely it is that circulation may be lost. To stop loss of whole mud, voids must be bridged so that a filter cake can be laid down on the permeable section. The plugging material must be of the proper size and shape to offer greater resistance to the fluid flow around it than the flow up the annulus. A plugging composition that satisfies these requirements may not be able to be handled by the small rig pump available.

Jan 1, 1979

By Herman Dykstra, R. L. Parsons

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.' THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy's Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. 'Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If' .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By N. A. Gokcen, J. Chipman

A revised treatment of the authors' published data eliminates the complex relation previously proposed between concentration of silicon and activity coefficient of oxygen in liquid iron. Revised values of the thermodynamic properties of the liquid solution are presented. IN a recent experimental study of the reaction SiO2 (s) = Si+ 2O; Kf, = [% Si] [% O]² [1] the authors' found a substantially constant equilibrium product in liquid iron at 1600°C of 2.8x10-5 They also reported extensive data on the reactions: SiO2 (s) + 2H2 (8) = Si- + 2H2O (g); K'2= [% Si] (H2O/H2 [2] and H, (g) + 0 = H2O (g); K'3 =( H2 O [31 (H2) [%O] From the results on reaction 3 and earlier data of Dastur² on this same reaction in the absence of silicon, they determined the activity coefficient of oxygen, f0, on the basis of the definition K3 = (H2O)/ (H2)f0 [% 0] where K, is the equilibrium constant and f0, is taken as unity in the pure Fe-0 system. Similarly values of fsi were deduced from results on reaction 2. In a more recent study" of analogous reactions in the system Fe-A1-0, it was found impossible to reconcile the results on reaction 3 with Dastur's data; accordingly the latter were ignored and the equilibrium results were extrapolated to find a value of K, at zero concentration of aluminum. This procedure failed to locate the cause of the discrepancy but it did yield reasonable values of activity coefficients. It also avoided introduction of the complex empirical relation between the oxygen activity coefficient and the concentration of the added element. The same type of discrepancy exists for system Fe-Si-0.' In the earlier paper an attempt was made to fit both sets of data by a single curved line (Fig. 6 of ref. l), the form of which is contrary to the theoretical requirement of a finite slope at infinite dilution. In the light of experience on the Fe-A1-0 system the discrepancy must be recognized as one which can be resolved only by more refined measurements. Accordingly Figs. 6 and 10 are retracted. It is pointed out also that until the discrepancy is resolved Figs. 7, 8, and 11 are subject to some uncertainty. Qualitatively the following conclusions still appear valid: 1—The activity coefficient of oxygen is reduced by addition of silicon. 2—In dilute solutions the activity coefficient of silicon increases with its concentration. 3—With respect to equilibrium in reaction 1, the above effects are approximately compensating. The discussion of K'1 in the previous paper requires no revision. It was pointed out that the constancy of the product [% Si] [% 0]² ndicated a compensating effect of the activity coefficients of silicon and oxygen. Therefore, as a very good approximation, K1 = K'1 and the following average values are suggested both for K, and K', at the temperatures 1550°, 1600°, and 1650", respectively, 1.0x10-", 2.8~10-" and 5.5 ~lo-'. Revision of the thermodynamic treatment is necessitated by the recent appearance of new data, based on a combination of combustion and solution calorimetry,' which yields for the heat of formation of low-cristobalite from the elements, the value —209,330 ±250 cal per mol at 25°C. This is about 4000 cal larger than the value previously accepted. The new value for cristobalite is used, together with Kelley's tables of high-temperature heat contents" and entropies and with Korber and Oelsen's' heat of fusion of silicon to obtain the following equation for the standard free energy of cristobalite in the temperature range 1700" to 2000°K: Si (1) + 02 (g) = Si02 (crist.); ?F° = -217,700 + 47.OT [4] The free energy of solution of 0, in liquid iron is:8 O2 (g) = 20 (in Fe); AF° = -55,860 - 1.14T [5] and these two equations are combined to give: Si (1) + 20 = SiO, (crist.); AF° = -161,840 + 48.14T [6] ?F°1873 = -71,700 cal. From the experimental value of K, = 2.8x10-5, Si + 20 = SiO2 (crist.); ?F°1873 = -39,000 cal. [7] The combination of Eqs. 6 and 7 yields the free energy change when liquid silicon dissolves in iron to form the dilute solution of unit activity (1 pct). Si (1) = Si; ?F°1873 = -32,700 cal. [8] The heat effect in this process according to Korber and Oelsen' is an evolution of 28,500 cal per gram

Jan 1, 1954

By Hokuichiro Ohmachi

INTRODUCTION Metallic minerals have been formed only through complex geologic processes which took place at certain stages of the earth's histrory. Their concentration, abundance, and distribution are, therefore, restricted geologically, and very small in global scale. Ore is a mineral deposit or mineral concentration from which metals can be economically extracted by the contemporaneous technology. By this definition, several factors should be considered before any mineral deposit is regarded as an ore. The grade of mineral concentration and the scale of the reserve are most important. Demand for the metal and its price are the other factors. Advancement of technology in mining and extraction are also vital. Copper, for example, is now commonly extracted from the ores containing less than 0.3% of the metal. However, in the 1700's the common minerable grade was about 13%, and in the 1900's between 5 to 2.5%. This decrease in minerable grade is a function of not only the demand for copper, but also technological progress in mining and metal extraction in this case. Availability of ores has recently been subject to political factors, which were not the primary concern in the past. The petroleum of the Middle East is the prominent example. The metallic minerals are classified into three categories. The first one is the minerals free from these factors mentioned above. Iron and manganese belong to this category. The second category include aluminum, copper, nickel, cobalt, titanium, lead and zinc. These metals can be provided by the future progress in technology which enables the use of lower grade deposits. The third one represents the metals whose occurrence is geologically limited, and, thus, subject to the political factors. Niobium and tantalium are the example. In this article, I discuss these minerals in detail to give an outline of the factors which brought their concentration, distribution and availability. SELECTED MINERAL CONCENTRATIONS Iron Major iron ore deposits are the features of the relics of ancient continental crusts. They are found in all the continents where the Precambrian rocks expose. They do not occur in the oldest rocks (3,000my old) nor in the youngest ones (600my old). They are sedimentary rocks, normally, exhibiting alternations of bands of iron oxides (magnetite and hematite) or iron silicates (especially greenalite), and bands of silica, variously described as jasper, quartzite and chert. The best known banded iron ore deposits exist in the Lake Superior region of the North America, and their distribution extends at intervals to Labrador of Canada. Other examples occur in the U.S.S.R. Brazil, Venezela, India, the mainland China, and southern and western Africa (Liberia and Guines) (Table 1). Recently large deposits have been discovered in the Hammersley region of the western Australia, and northern and central region of Brazil. These ore deposits were originally formed in shallow seas where simple but abundant life existed. After the deposition, enrichment processes , related to tropical weathering, have brought about wholesale removal of silica to produce the deposits of best quality. Ores with 50-60% iron represent the best products in the Precambrian fields. Since local reserves of this high grade ores approach to exhaustion, benefication processes have been used effectively to upgrade much leaner ores (for example, taconite with 20-25% iron in Minesota, U.S.A.). There are two other types of iron deposits. Ironstone, a sedimentary rock containing goethite, chamosite and siderite, first appeared in early Paleozoic times in the stratigraphic record, and reached its Zenith in the Jurassic. Similar bedded iron deposits are found in the belt from the Cleveland Hills to Oxford in England, and in the "minette" oolitic ores of Alsace-Lorraine in Luxenburg and France. The iron content of these ironstones seldom exceed 30%, but they usually contain calcite and are self-fluxing. They also have phosphorous as undesirable impurity. Because of the grade and impurities, they require more fuel than the Precambrian ores. The third type of iron ore is associated with igneous activity and consists of magnetite and some hematite, and apatite. This deposit is believed to be a product of magmatic differenciation. It occurs in Kiruna of Arctic Sweden. The iron ore, therefore, present no global shortage problem. They extend to considerable depths. Their concentration is largest among the metals. If low grades were treated, the resource can be stepped up much substantially. Manganese Manganese is a ferro--alloy metal, thus, essential to the manufacture of sound steel. Manganese comes from manganese oxide and silicate ores. Many minerals contain manganese, but only a few oxides, silicates, and, in some places, carbonates (rhodecrosite) are mined as ore. Types of manganese deposits are bedded, massive,

Jan 1, 1982

By G. A. Swanquist, E. M. Morales

INTRODUCTION The idea of water management and control at the Church Rock Mill operations began to take shape in February 1979. At that time, we were already investigating the feasibility of decreasing the fresh water requirements so that the solids would become the limiting factor in tailings impoundment utilization. The area for solution evaporation could be kept at a fraction of the normal requirements under the standard process of full water usage. The Church Rock Mill is an acid leach circuit followed by solids/liquid separation with thickeners in counter current decantation, and solvent extraction. Following the normal design of acid leach circuits, reuse of tailings solution was not incorporated in the original mill process design. INITIAL WATER CONTROL INVESTIGATIONS The investigations to decrease the fresh water requirements centered around modifying the grinding circuit from the present semi-autogenous grinding (SAG) mill in closed circuit with hydrocyclones, to open circuit grinding with a rod mill. The open circuit grinding with the SAG mill and rod mill in series had the potential of decreasing the water requirements for grinding and leach dilution by approximately 50% or 1.4 m3/min (300 gpm). The grinding pulp density would be maintained at 70 to 72% solids, and the leach dilution to 50% solids would be accomplished with acid tailings liquor recycle. In such a grinding circuit arrangement, the SAG mill would provide the primary or coarse grind, and the rod mill would be used for the fine grind. By the SAG mill and rod mill series grinding method of water control and other secondary water controls in various places downstream from the grinding circuit, the required necessary evaporation area was estimated at 120 acres of liquid surface. A second method of water control at grinding was investigated. A two-stage cyclone classification circuit appeared to have a good potential of achieving the same water reduction at a much lower capital and operating cost. However, in retrospect, this would not have been a viable method since a high slime recycle load would have been established hindering classification. The use of reagents to neutralize the acid tailings solution was not considered seriously at that time, since it would have materially increased operating costs, although it would have also allowed more tailings solution recycle and consequently, less fresh water usage. However, with the tailings solution deposition area available at that time, it was not then necessary to incur the high cost of neutralization. The control expected by the series grinding of semiautogenous and rod mills would have been sufficient to maintain a water consumption/evaporation equilibrium well in line with the available land area. IMPLEMENTATION OF NEUTRALIZATION OPERATIONS During the summer of 1979, the UNC Church Rock Mill experienced a tailings dam breach which resulted in a prolonged mill shutdown. Upon resumption of operations at the end of October 1979, tailings deposition was restricted to a small portion of the tailings impoundment area. Figure 1 shows the general tailings area and the limits of the present deposition area in the central part including the borrow pits. These borrow pits had been excavated to provide materials for tailings dam construction. Immediately after resumption of operations, it became evident that it would be necessary to control the quantity of liquid to be evaporated because of the small confined area available for tailings solution deposition and to maximize the deposition time in the tailings area. The water control required had to be exercised on a large scale, and to be in operation as quickly as possible. An obvious solution was to reuse the tailings liquor in mill process. Immediate steps were taken to install the necessary equipment for tailings neutralization on an interim basis. Anhydrous ammonia was selected as the primary neutralization reagent since it was the quickest system that could be placed in operation. Previous laboratory tests indicated fair results with ammonia neutralization. Such a system required a minimum of installed equipment and handling. INITIAL NEUTRALIZATION OPERATIONS Actual neutralization operations began on November 26, 1979. The raffinate solution which normally would have been discarded was pumped to a 3.7 m (12ft) diam by 4.3 m (14ft) tank for reagent contact, see Figure 2. At this tank, anhydrous ammonia was added directly from the tanker trailers and controlled at pH 7.0 nominally. Agitation was provided by air sparging. The neutralized product formed a highly viscous slurry in the grinding circuit which resulted in pumping and cyclone classification problems.

Jan 1, 1982

By J. J. Gilman

The data presented indicate that the critical shear stress and strain-hardening Thedatapresentedrate of a zinc monocrystal depend on the orientation of its slip direction with respect to its external boundaries. The tendency of a crystal to form deformation bands also depends on its shape. THE plastic behavior of pairs of zinc monocrystals in which both members of the respective pairs had the same orientation with respect to the longitudinal axis, but each had different orientations with respect to their rectangular external shapes, were compared in this investigation. The purpose of the investigation was to see what influence the shape or surface of a zinc crystal has on its mechanical properties. In a previous investigation of triangular zinc monocrystals,1 anomalous axial twisting was observed which seemed to be related to the triangular shape of the crystals. Wolff,' in 400°C tensile tests of rectangular rock-salt crystals bounded by cubic cleavage planes, found that, of the four equivalent slip systems, the two with the "shorter" slip directions yielded and produced slip lines at lower stresses than the other two. This observation and the work of Dommerich³ as formulated by Smekal4 as a "new slip condition" for rock-salt: "among two or more slip systems permitted by the shear stress law, with reference to the formation of visible slip lines by large individual glides, that slip system is preferred which has the shortest effective slip direction." More recently, Wu and Smoluchowski5 reported essentially the same effect for ribbon-like (20x2x0.2 mm) aluminum crystals at room temperature. Experimental Chemically pure zinc (99.999 pct Zn), purchased from the New Jersey Zinc Co., was the raw material. Glass envelopes, containing graphite molds and zinc, were evacuated while hot enough to outgas the graphite but not melt the zinc. At a vacuum of about 0.2 micron the envelopes were sealed off and then lowered through a furnace at 1 in. per hr so as to melt and resolidify the zinc and produce mono-crystals. One-half of one of the molds is shown in Fig. la. Each mold consisted of four pieces from a cylindrical graphite rod that was split longitudinally and transversely at its midpoints. Rectangular milled grooves 0.050 in. deep and % in. wide formed the mold cavity when the split halves were assembled with twisted wires. Fig. lb shows the specimen shape obtained when the top and bottom mold-halves were rotated 90" with respect to each other. Good fits prevented leakage and excess zinc was necessary to provide enough liquid head to fill the mold completely. In removing soft crystals from the molds it was impossible to avoid small amounts of bending. However, manipulations were carried out whenever possible with the crystals protected by grooved brass blocks. All specimens were annealed prior to testing. From the top and bottom sections of each crystal, X-ray specimens and tensile specimens 7 to 8 cm long were sawed. The tensile specimens were annealed inside evacuated tubes for 1 hr at 375°C. Next the crystals were cleaned and polished by 2-min dips in a solution of 22 pct chromic acid, 74 pct water, 2.5 pct sulphuric acid, and 1.5 pct glacial acetic acid.' Cleaning was followed by a 10-sec dip in a 10 pct caustic solution, then washed in water and alcohol, and dried. This treatment results in a bright surface covered by an invisible oxide film. The testing grips were a slotted type with set screws and were supported in a V-block during the mounting operations in order to avoid bending the crystals. A schematic diagram of the recording tensile-testing machine is shown in Fig. 2. The machine has been described elsewhere.' The head speed was 0.3 mm per sec for all tests. The crystal orientations were determined by the Greninger X-ray back-reflection method with an estimated accuracy of 1. Description of Crystal Geometry A schematic picture of a rectangular zinc mono-crystal is shown in Fig. 3. ABD designates the front edge of a basal plane (0001) of the crystal, the only active slip plane for zinc at room temperature. Of the three possible (2110) slip directions, the active one is indicated by an arrow. Cartesian coordinates are taken parallel to the specimen edges. The normal, n, to the basal plane (n is parallel to the hexagonal axis) has the direction cosines a, ß and ?. X0 = 90 — y is the angle between the longitudinal axis and

Jan 1, 1954

By R. L. Parsons, Herman Dykstra

A numerical method for solving partial differential equations in steady state fluid flow is described. This method, known as the "relaxation method," has two advantages over analytical methods: (1) practically any problem can be solved, and (2) a solution can be obtained quickly. A disadvautage is that the solution is not general. The method is applied to core analysis and relative permeability measurement to calculate constriction effects and to calculate the true pressure drop measured by a center tap in a Hassler type relative permeability apparatus. Further applications are suggested. INTRODUCTION Many problems in fluid flow cannot be solved analytically because of the nature of the boundary conditions. For many problems, however. an exact answer is not necessary because boundary conditions are not exactly defined or the parameters describing the porous medium are not accurately known. The relaxation method can be used to obtain an approximate answer easily and quickly for the flow of incompressible fluids in porous media. The method can also be used for other types of problems, such as determining the stress in a shaft under load. or the temperature distribution during steady state heat flow. In this discussion only calculations concerned with the flow of fluids in porous media will be considered. The method was introduced by R. V. Southwell in 1935.' THEORY The treatment given here follows that given by Enimons.2 Consider a porous medium to be replaced entirely by a net of tubes of equal length and uniform cross-sectional area as shown in part in Fig. 1. Assume that the net of tubes behaves exactly like the porous medium which it replaces; that is, the net can be made fine enough to reproduce exactly the porous medium. Assume also that Darcy's Law can be used to calculate the flow from one point to another point through these tubes. The flow from point 1 to point 0 is KA . ------ P-P) .......(11 where a is the distance between points: K is the "permeability" of a tube; A is the cross-sectional area of a tube; is the viscosity of the liquid in the porous medium; and (P1 — P0) is the pressure difference between point 1 and point 0. In like manner the flow can be calculated from points 2, 3, and 4 to point 0. The net flow into point 0 is Qo = KA/µa (P1 + P2 + P33 + P4-4P0) . . (2) MB For an incompressible fluid the net flow into point 0 will be zero or, Q. = 0. This says that at point 0 fluid is neither being accumulated nor depleted. 'Therefore. P1 + P2 + P3 + P4 - 4P0 = 0 .... (3) . If. now. with specified boundary conditions. the pressure i.; known at a finite number of points in a given region, as at the points shown in Fig. 1, Equation (3) will be satisfied at every point. If, on the other hand, the pressure is not known, the pressure can be guessed at these points. Then. unless the guess is perfect. Equation (3) will not be satisfied at all of the points. When Equatiol~ (3,) is not satisfietl. let d = P1 + I?, + P, + P, - If' .,....(4) where 6 is an apparent error and is called the residual at point 0. Equation (4) shows how much the pressure guess is in error at point 0 with respect to the surrounding points. A positive residual means that the pressure is too low, and a negative residual means that the pressure is too high. To bring the residual, 6. to zero in order to satisfy Equation (3). it is necessary to make changes in the pressure guesses. Equation (4) shows that a +1 change in Po will change the residual at point 0 by -4. A +1 change in the pressure at any of the four surrounding points will change the residual at point 0 by +l. Thus it can be seen that a change at any point will affect the residual at that point and the four surrounding points. By changing the pressure from point to point, all of the residuals can eventually be brought nearly to zero and the problem will be solved. This procedure is the essence of relaxation methods and is used to relax the residuals so that Equation (3) is satisfied at every point. The procedure can be most easily explained in detail by solving a simple problem. as Southwell says, "To explain every detail of a practical technique is to risk an appearance of complexity and difficulty which may repel the reader. A

Jan 1, 1951

By S. W. Ross

In 1936 a description of the plant (Fig 1) and process employed by the Electrolytic Zinc Co. of Australasia Ltd. for the recovery of zinc from zinc concentrate by the electrolytic process was prepared. † During the twelve years which have elapsed since the preparation of the earlier paper, several major changes in the metallurgy of the process have been introduced. It is the purpose of the present paper to give a general description of these changes and thus to bring up to date the description of the plant and process. Summary The major changes in Risdon practice since 1936 have been: 1. Replacement of two stage roasting by a preliminary roast followed by the flotation of all the leach residue and the roasting of the flotation concentrate. 2. Screening of all calcine fed to the pachucas. 3. Continuous leaching of calcine and improved classification of pachuca discharge. 4. Close control of hydrogen ion concentration during purification for iron removal. 5. Recovery of cobalt as a good grade oxide. 6. Production of part of the zinc output in the form of "four nines" metal (99.99 pct purity). 7. Closer spacing of electrodes thus increasing the potential output of cathode zinc per cell by 50 pct. Changes which are in prospect and for which construction work is proceeding at the present time involve the starting up of: 1. Two suspension roasters. 2. A contact acid plant to produce 150 tons‡ of acid per day and replacing the existing Mills Packard chamber plant. 3. Extra power station capacity to permit greater current flow to existing cell room units. This will increase the output of cathode zinc from about 245 to 290 tons per day. Plans for the future envisage the building of an ammonium sulphate plant, the first unit of which will produce about 50,000 tons per year, and improved treatment of zinc plant residue for the recovery of zinc, lead and other metals. At the end of this paper tables of metallurgical data are presented relating to the year ending June 30, 1948. Details of Changed Practice Output In 1936 the production of cathode zinc amounted to about 200 tons per day. This has since been increased to about 245 tons per day while plant extensions are practically complete which will permit of an output of about 290 tons per day in the near future. Roasting Division ROASTING POLICY A major change has occurred in the roasting policy. Twelve years ago the method in use was to carry out a two stage roast in the first stage of which sulphide sulphur was reduced to about 6 pct. The pre-roast calcine was re-roasted in modified Leggo furnaces using coal as fuel, sulphide sulphur being reduced to about 0.8 pct. The whole procedure was described on pp. 482491 of the earlier paper. Although this roasting procedure had certain advantages it possessed some distinct disadvantages. For instance, it appeared uneconomical to heat up the entire input of pre-roast calcine to roasting temperature by the expenditure of fuel in order to oxidise a few percent of sulphide sulphur. It was argued that if the pre-roast calcine were leached and a process could be developed for the recovery of a zinc sulphide concentrate from the leach residue, this concentrate, small in weight compared with the pre-roost calcine, would probably roast autogenously, thus virtually eliminating the expenditure of fuel as well as greatly increasing the weight of sulphur oxidised per square foot of furnace hearth area. The obvious method of producing a suitable concentrate from leach residue was by flotation. It will be recalled (see p. 495 of the earlier paper) that when two-stage roasting was practised the leach residue was classified, the granular fraction was ground and floated while the slime fraction was thickened, filtered and dried ready for shipment to a lead smelter. This process worked quite successfully. However, when trials were made of leaching a calcine carrying several percent of sulphide sulphur the granular fraction still floated well, but the slime fraction carrying 8-10 pct sulphide sulphur yielded very poor results when subjected to flotation. This fact held up the application of "pre-roast" leaching for many years. However, successful flotation of the slime fraction of leach residue was finally achieved and in August 1940 the slime flotation plant began operation, while the leach-

Jan 1, 1950

By S. D. Baumer, P. M. Hulme

IN the late thirties, the National Carbide Co. cooperated with C. E. Wood, of the U. S. Bureau of Mines, in his investigation of the relative merits of various desulphurizers, including soda ash, caustic soda, and calcium carbide. Laboratory tests showed that carbide, when it could be made to react, is an excellent desulphurizing agent for molten iron. Sulphur content can be driven to lower levels and higher extractions obtained with carbide than with actionsany of the more common reagents. Wood's results1 are shown in Table I. Unfortunately, as the Handbook of Cupola Operation puts it, the chemical fact that carbide is a good desulphurizer was of only academic interest because it was found to be extremely difficult to devise a practical means to make it react with molten iron. Calcium carbide is formed in the electric furnace at 4000°F and above, and its softening point is probably at least 500 °F above the usual working temperatures encountered in iron and steel practice. Consequently, carbide does not form a true slag but floats as a dry powder on top of the metal and only a very small portion of it ever comes in actual contact with the iron. Stirring with a rabble, or pouring the metal over the carbide, increases the efficiency only slightly. Extractions of 20 to 30 pct can be obtained in this manner, but conventional soda slag treatment can do better than this and do it more cheaply. All attempts to lower the melting point of carbide in order to obtain a reactive, liquid slag have so far proved fruitless. Directly under the arc in a metallurgical electric furnace, carbide becomes highly reactive. Excellent sulphur removal can be obtained without any slag other than a thin layer of carbide." imilarly, good results are obtained by adding small amounts of carbide to the finishing slag in double-slag arc furnace practice. To react a liquid with a solid, it is axiomatic that the liquid has to wet the solid before anything can happen. If the solid is heavier than the liquid, the problem is easy, but it becomes more difficult when the solid is much lighter than the liquid, as in the case of carbide and liquid iron. Wood recognized this problem and solved it in a unique fashion. The results shown in Table I were obtained by spinning the carbide beneath the surface of the molten iron by means of a refractory centrifuge. This technique allowed each particle of the finely divided carbide to come into intimate contact with the metal and to be wetted thereby. Wood's centrifuge technique was successful in the laboratory where it achieved excellent and consistent results. Some attempts were made to expand this method to commercial practice, but serious difficulty was encountered in obtaining a refractory centrifuge head that would be economically feasible. About this time the war intervened and the project lay dormant for several years. In 1944, it was revived. It was suggested that the carbide could be blown into the metal with a carrier gas in an attempt to eliminate the necessity for the expensive and brittle centrifuge. The idea was first tried out in a fairly large ladle of iron using natural gas as the carrier. Considerable sulphur was removed, but it was quite obvious that the use of natural gas was not practical. Attempts then were made to blow carbide into molten iron using, in turn, nitrogen, argon, carbon dioxide, air, and oxygen. The latter two gases proved unsatisfactory. Calcium evidently prefers oxygen to sulphur because in the tests calcium oxide and carbon dioxide were produced, the sulphur still being untouched in the iron. Nitrogen, argon, and carbon dioxide gave much better results, although the efficiencies and extractions were erratic, and only a few isolated tests approached the results obtained by Wood. Table II shows typical results obtained with these gases. The sulphur removals were interesting, sometimes even encouraging, but it is evident that such erratic behavior could not be tolerated in commercial practice. A number of different types of equipment, such as sand blasting machines, refractory guns, and the like can used to blow the solid into the metal. All types required relatively large quantities of gas in order to maintain the flow of solid carbide through the system and into the metal. It was observed that the bubbles of gas breaking through the surface of the metal contained quantities of unreacted carbide. The liquid metal never came in contact with these particles and if it cannot wet them it cannot react with them. The initial work had shown that carbide had great possibilities as a desulphurizer. In practice

Jan 1, 1952

By N. J. Grant, C. C. Wang

The Ti-Cr-O ternary system has been studied in detail near the titanium-rich corner within the limits of 10 wt pct 0, and 20 wt pct Cr. Studies were extended, but not in detail, to the region beyond 25 wt pct 0, (50 atomic pct) and 62 wt pct Cr (60 atomic pct). Four isothermal sections at 1400°, 1200°, 1000°, and 800°C are presented as well as two vertical sections at 1 and 2 wt pct 02. DURING the last decade much interest has been shown in the development of high strength titanium alloys for high temperature and corrosion resistant applications. Extensive research is being carried out at present, as the current literature indicates, in order to study the properties of titanium and to develop improved alloys. Two of the important alloying elements in commercial titanium alloys are chromium and oxygen and it would be desirable to know their combined influence upon titanium. For this purpose the present work was carried out to investigate the titanium-rich corner of the ternary system TiICr-0. The binary systems Ti-Cr and Ti-0 have been published recently. The Ti-Cr system was studied by several investigators " and their results are in close agreement. The eutectoid decomposition of the B phase has been shown to be extremely sluggish. TiCr, was the only intermetallic compound found in this binary system and was formed at 1350°C by a transformation from the p phase. TiCr? was established as the cubic C 15 (MgCu,) type of structure with 24 atoms per unit cell and was designated as the y phase. This terminology will be adopted in the present work. There was disagreement about the actual composition of this compound among the several investigators, although it is evident from their data that the compound probably has a solubility range of about 2 to 3 pct and is in the vicinity of 65 pct Cr. It has been indicated recently that a high temperature modification of this y phase (TiCr,) existed at a temperature above 1300°C." ' This high temperature modification was identified as a hexagonal C 14 (MgZn,) type of structure with 12 atoms per unit cell. The exact transformation temperature from the high temperature phase to the low temperature phase has not been established. A considerable hysteresis was observed and, due to the sluggishness of this transformation, the high temperature phase often co-existed with the low temperature phase at temperatures below 1300°C. A preliminary study of several Ti-0 compounds and the Ti-0 system had been carried out by Ehr-1ich."-"' The most complete binary Ti-0 system was the one reported recently by Bumps, Kessler, and Hansen." The first intermediate phase found in the system was the 8 phase which formed by a peritec-toid reaction of the phases a and Ti0 at temperatures below 925 °C. This reaction is extremely sluggish. The structure of this 8 phase was tentatively identified by these authors as being tetragonal and the lattice constants were found as c,, - 6.645A, a,, = 5.333A and c/a = 1.246A. Experimental Procedure The raw materials used for this investigation were TiO,, electrolytic chromium, iodide titanium, and sponge titanium. The TiO, was in the form of powder of chemically pure grade (99.8 pct pure). The chemical analysis of the electrolytic chromium was: 0, 0.50 pct; Fe, 0.07; Cu, N, and C, 0.01; and Pb, 0.001. The oxygen in the chromium was calculated as part of the final oxygen content of the alloys. The alloys were prepared by the cold crucible method using a tungsten arc. The entire system was evacuated and flushed with purified helium three times and then filled with helium. Each alloy was melted, turned over, and remelted at least four times to insure homogeneity. The total melting time was generally from 6 to 10 min. A master alloy of 25 pct 0,-75 pct Ti was prepared to facilitate alloying by melting compacts of TiOl powder with either iodide or sponge titanium, yielding the compound TiO. It was found necessary to bake the TiO, powder compact at about 150°C to remove adsorbed moisture. This was done to prevent the disintegration and spattering of the compact when the arc was struck. TiO, powder dissolved quite readily into the melt and no other trouble was encountered.

Jan 1, 1955

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. N. Poundstone

The mining and ventilating system used in development work in the Pittsburgh Seam in northern West Virginia is discussed. The seam conditions and the nature of the accompanying methane gas are described. The type of equipment and the mining cycle will be discussed, showing how they are well suited for very gaseous development work. Face ventilation in development work is possibly the fastest growing problem of the industry. The coal mines of the future will be faced with the prospect of mining from under increasing depths of cover. Consequently, larger and larger amounts of methane gas probably will be found. The Pittsburgh Seam, in northern West Virginia, is an example of an area already faced with this problem. At the present time, most of the development work being done in this seam lies beneath 500 to 1200 ft of cover. The Pittsburgh Seam in this area has always been very gassy—even near the outcrop—and the recent development work has been accompanied with extremely large volumes of gas. In many cases, a single development section has liberated in excess of 1,000,000 cfm in 24 hr. This problem of heavy gas liberation was the chief concern, several years ago, when continuous mining equipment was first considered at Christopher Coal Co. All of us were apprehensive about the liberation that would accompany rapid extraction in a single working place. However, the experience during the past few years has shown that this ability to mine only one place at a time, is actually the key to working this type of coal. With all of the mining or advancement concentrated in one place, the ventilation can also be concentrated. By this it is meant that continuous mining permits the active working place to be ventilated with a maximum amount of fresh air, taken directly from the intake source without first passing another working place. Continuous mining (and a good ventilating system) also permits a much greater concentration of attention or vigilance to the actual working place. There are two things that are very important to the mining of coal having high rates of liberation. First, adequate volumes of air are necessary. Second, and perhaps more important, a mining and ventilating system must be used that will provide an uninterrupted flow of air to every portion of the working section. Liberations of this magnitude take only a few seconds of interruption for a dangerous accumulation of gas to occur. With adequate volumes of air available, the ability to concentrate ventilation more than offsets the concentration of gas emission that is inherent with continuous mining. The mining system used with continuous mining equipment at Humphrey Mine is similar to the system used at many mines in the area for development work. This system is designed to favor ventilation, realizing that other efficiencies are meaningless if the equipment must stop because of ventilation difficulties. This plan is especially well suited to minimizing ventilation interruptions. Basically, the overall plan of mining is to develop headings into virgin coal and encircle or block out large areas. These blocks are generally at least 2000 ft sq. The purpose of this blocking out is to bleed gas from the area before pillaring. Experience has shown this method to be quite effective, even in the most gassy areas. The gas in this field seems to migrate or flow readily from the solid coal into the outside return headings of the development work. The numerous clay veins and slips that are found in the area are extremely good avenues for gas flow. A block of coal, surrounded with development headings, usually bleeds off readily; and since it is cut off from the virgin coal, it is not subject to gas migration, through the seam, from this source. However, the outside return of the encircling development work, adjacent to the virgin coal, may liberate gas for years. This liberation from the outside ribs in the virgin coal is the reason split ventilation is used in development work. If split ventilation were not used, there would, in many cases, be a serious build-up of gas in the intake before it could reach the working face. Fig. 1 shows a typical development section having seven headings. The two outside places on each side are returns, and the three center headings serve as intakes. This section is equipped with a ripper-type continuous mining machine. An off-track loading machine is used to load from a surge pile on the

Jan 1, 1961

By J. R. Elenbaas, J. A. Vary, D. L. Katz

The development of gas fields, oil fields and aquifers for storing natural gas is treated from two main vieu.-points: (I) the volumetric storage capacity for gas in a given situation and (2) the prediction of the number of wells required for the delivery of gas. Other experiences in the design and operation of storagc fields are incluclerl INTRODUC TION Storage of natural gas in underground reservoirs near the terminus of long distance pipelines has been the prime factor in opening the space heating market to the natural gas industry. Storagc has permitted a major. increase in both the load and the load factor of pipe-lines; some are now operating at steady load throughout the year. Thus, underground storage has been responsible for the rapid increase in demand for natural gas in recent years. Three types of reservoirs have been used for gas storage: natural gas reservoirs, oil reservoirs, and waterbearing sands or aquifers. This paper presents the factors to be considered when developing gas storage reservoirs of these three categories. There are two prime considerations tor any storage reservoir: (1) the volume of gas which a given reservoir will store advantageously and (2) the number 01 wells needed to provide the required peak deliverability. These two problems will be considered for the three types of reservoirs just noted STORAGE IN PARTIALLY DEPLETED GAS FIELDS Early storage operations consisted of replenishing the natural gas in a depleted gas field situated adjacent to the market. Today, newly discovered fields near the market may be considered for storage, and this discussion applies equally to both types of reservoirs. For reservoirs originally containing gas or oil, the question of the impermeability of the cap rock nor-mally does not arise. However, such fields are likely to have many wells drilled either to or through the reservoir under consideration. Positive assurance must be obtained that such wells are or can be made mechanically tight. Corroded casings may need to be lined or permanently plugged. Abandoned wells should bc reopened and properly cemented. The volumetric capacity for gas storage depends upon space available in the porous rock as well as pressure and temperature of the gas in the reservoir. The production-pressure decline data on partially depleted gas reservoirs without water drive permit calculation of the reservoir space for gas. Isopachous maps of sand volume and porosity data for the reservoir rock provide an alternate method of calculating the pore volume for water-drive reservoirs. The pressure range selected for the storage cycle depends upon ()) the safe upper limit of pressure. 2) the flow capacity of wells and (3) compression requirements when injecting gas into the reservoir or delivering to market. Normally, gas and oil fields have pressures at discovery in the range of 0.43 to 0.52 psi/ft of depth. Pressures of around 1.0 to 1.2 psi/ft of depth appear to lift the overburden1-3 and invite uncontrolled movement of fluids in the porous rock. Some top pressure is normally selected for a storage reservoir ranging from below discovery pressure for deeper reservoirs to 0.65 psi/ft of depth for shallower reservoirs. Pressures to 0.66 psi/ft have been experienced without difficulty. The lower pressure limit is set by water intrusion accompanying low pressures, reduced flow capacity for wells at lower pressures and compression requirements. Depletion-type gas reservoirs often encounter water problems in the later stages of gas production. Such water intrusion may be due to movement from the surrounding aquifer. Accordingly, displacement of this water back into the aquifer by gas pressure and subsequent surges of water corresponding to the gas storage pressure cycle must be considered. Storage fields often produce in four months a volume of gas equal to its initial content. Rapid decreases in reservoir pressure occur, such as 20 psi/day. Accordingly, closed-in pressure observation wells which reflect the pressure in the bulk of the reservoir are required for following the operation of the reservoir. It has been found that a plot of observation wellhead pressures against gas content, Fig. 1, is very useful in observing operation of the field, checking the inventory and predicting future behavior. The plot is based on a given quantity of base or cushion gas in place. The injection and withdrawal curves may spread depending upon the homogeneity of the reservoir rock. permeability of the rock, well spacing and flow rates.

By R. W. Regehr, R. A. Burmeister

Epitaxial growth of GaAs 1-xPx on germanium substrates was achieved using an open tube vapor transport system. The compositional range of 0.3 < x < 0.4 was examined. The best results were obtained with (311) orientation of the germanium substrate. The physical and chemical properties of the resulting layers were investigated using several techniques. Spectrographic analyses of the layers indicate substantial incorporation of germanium into the GaAs t-X Px layer. Evidence is presented which indicates that this incorporation occurs via a vapor phase transport process rather than by solid phase dijfu-sion. Electrical measurements suggest that the germanium thus incorporated behaves predominantly as a deep donor in the compositional range of 0.33 < x * 0.40 and has a deleterious effect upon the luminescent properties of GaAs1-x Px. The increasing technological importance of GaAs1-xPx for use in light-emitting devices has led to an evaluation of several aspects of existing growth processes. The method most commonly used to prepare GaAs1-xPx for electroluminescent device applications is vapor phase epitaxial growth on GaAs substrates.&apos;-4 In a typical electroluminescent diode structure the active region of the diode is entirely within the epitaxial layer and thus the electrical properties of the substrate are relatively unimportant since it is effectively a simple series resistance (assuming hetero-junction effects to be negligible). The use of germanium rather than GaAs as the substrate material is of interest for several reasons. First, GaAs of reasonable structural quality has been epitaxially grown on germanium4-2 and it is reasonable to expect that GaAs1-xPx could subsequently be deposited on the GaAs layer. Second, germanium substrates are readily available with both lower dislocation densities and larger areas than GaAs. Finally, single crystals of germanium are more economical than GaAs single crystals. The principal objective of the present investigation was to test the feasibility of growing GaAs1-xPx epi-taxially on germanium substrates, and to evaluate the properties of such layers with regard to electroluminescent device requirements. The approach used was to a) demonstrate epitaxial growth of GaAs1-xPx on germanium, and b) characterize the relevant structural, electrical, and optical properties of the GaAs1-xPx layers. The possibility of germanium incorporation into the grown layers was of special interest since there was some indication of this in previous studies of GaAs growth on germanium.5&apos;11,12 Although a study of the electrical properties of germanium in GaAs1-xPx was not an intent of this investigation, several features of the electrical properties of the layers grown in the present study which appear to be due to germanium are described. EXPERIMENTAL PROCEDURE The open-tube vapor transport system used for the epitaxial growth of GaAs1-xPx is illustrated in Fig. 1. This system utilizes the GaC1-GaC13 transport reaction and is similar in most respects to the larger system described elsewhere.&apos; The germanium substrates were n-type, with a resistivity of 40 ohm-cm (Eagle-Picher Co.). These were cut to the orientations of {100), {111), and (3111, and were mechanically polished and chemically etched in CP-4 (5 min at 0°C) prior to growth. In some cases, a GaAs substrate was employed in addition to the germanium. The orientation of the latter was {loo}, and they were also mechanically polished and chemically etched prior to growth. The initial composition of the deposited layer was pure GaAs. After approximately 10 microns of GaAs was deposited on the germanium substrate, the phosphorus content of the layer was gradually increased over a distance of approximately 15 microns to the desired concentration and maintained at this value throughout the remainder of the growth. Typical operating parameters used during growth are given in Table I. Selenium was used as a n-type dopant in several runs to facilitate comparison of the electrical properties of the layers grown on germanium with those of layers grown on GaAs substrates, which are usually doped with selenium. The concentration of H2Se in the gas phase was adjusted to a value which would normally yield a carrier density of 1 to 5 x 101 7 at room temperature in layers grown on GaAs substrates. The terminal surfaces of the epitaxial layers were examined by optical microscopy for structural characteristics. Laue back-reflection photographs (Cu radi-ation) were also made on the terminal surface to verify the epitaxial nature of the deposit. After these steps

Jan 1, 1970

By Elwood G. Haney, R. N. Parkins

Stress corrosion cracking of two heats of 18 pct Ni maraging steel in rod form immersed in an aqueous solution of 0.6N NaCl at pH 2.2 has been studied on un-notched specimens stressed in a hard tensilf machite. Austenitizing temperature in the range 1830 to 1400 F has been shown to have a marked influence on the propensity to crack, the loulest austenitizing- temperature producing the greatest resistance to failure. In the nzosl susceptible conditions, the cracks followed the original austenile grain boundaries; but when tlze steels zcere heal treated to inproze their resistance to stress corrosion, the cracks becatne appreciably less branched and slzouqed significant tendencies to become trans granular. Electron metallography of the steels indicated the presence of snzall particles, possibly of titanium carbide, along- the prior austenite grain boundaries and these particles u:ere more readily detectable in the structures that were most susceptible to cracking. Crack propagation rates, which appeared to be dependent upon applied stress and structure, were usually in tlze reg-ion of 0.5 mm per hr and may, therefore, be e.xplained on tlze basis of a purely electrochetnical ,nechanism. However, there is some ezliderzce from fractography that crack extension may be assisted by ttlechanical processes. Anodic stit)zulation reduced the tiwe to fracture, although cathodic currents of small magnitudes delayed cracking-; further increase in cathodic current resulted in a sharp drop i,n fracture litne, possibly due to the onset of hydrogen ewbrittlement. THE use of the high strength maraging steels, with their attractive fracture toughness characteristics, is restricted because of their susceptibility to stress corrosion cracking in chloride solutions. Although this limitation has resulted in investigations of the stress corrosion susceptibilities of these steels, there have been few systematic studies aimed at defining the various parameters that determine the level of susceptibility. It is the case that the usual tests have been performed with the object of defining some stress or time limit, on unnotched or precracked specimens, within which failure was not observed,&apos; but while such results may be of some use in design considerations, they are necessarily concerned only with the steels as they currently exist and not with their improvement to render them more resistant to stress corrosion failure. This omission may be considered unfortunate because the indications are that stress corrosion in maraging steels shows dependence on structure in following an intergranular path, and since experience with other systems of intergranular stress corrosion crack- ing is that susceptibility may be varied by modifying heat treatments, a similar effect may be expected with maraging steels. It is sometimes from such observations that a fuller understanding of the mechanism of stress corrosion crack propagation begins to emerge, leading in time to the development of more resistant grades of material. The present work was undertaken to study only one aspect of the influence of heat treatment upon the cracking propensities of the 18 pct Ni maraging steel, namely the effect of austenitizing temperature, although certain ancillary measurements and experiments have been undertaken. EXPERIMENTAL TECHNIQUES Most of the measurements were made on a steel, A, having the analysis shown below, although a few results were obtained on a steel, B, having a slightly different composition. Both steels were supplied in the austenitized condition, A as 3/8-in-diam rod and B as 1/2-in.-diam rod. Cylindrical tensile test pieces were machined from the rods: the overal length was 2 1/2 in., the gage length 1 in. and the diameter 0.128 to 0.136 in. The stress corrosion tests were carried out with the specimens strained in tension in a hard beam testing machine, the necessary total strain being applied to the specimen over a period of about 30 sec, after which the moving crosshead was locked in position and the load allowed to relax as crack propagation proceeded; the load relaxation was recorded. The load was applied after the specimen had been brought into contact with the corrosive solution, the latter being contained in a polyethylene dish having a central hole through which the specimen passed, leakage being prevented by the application of a film of rubber cement. The specimen was in contact with the solution for over half of its gage length and the solution was exposed to the air during testing. The solution was prepared from distilled and deionized water to which NaCl was added, 0.6N, and the pH adjusted to 2.2 by HCl additions. The composition of the solution

Jan 1, 1969

By L. R. Raymond, J. L. Hudson

This paper proposes a comparatively short-term (8 to 10 hours) well test for detecting and characterizing well-bore damage and for measuring mean formation permeability. The proposed test is made by injecting fluid at constant pressure, recording injection rate as a function of injection time. After one to four hours of injection, the well is shut in and fall-off of bottom-hole pressure is obtained as a function of shut-in time. Formation permeability is estimated by an iterative technique. First, a value of formation permeability is assumed. Then, a plot of the recorded injection rate as a function of dimensionless time is made, using the assumed pertneability value. From the slope of the injection-rate curve. a new value of formation permeability is calculated. If the new value agrees with the original assumed value, the assumption was the correct formation permeability. If the values do not agree, the process is repeated using the new permeability value in the calculation. Convergence is rapid, and a reliable permeability value results. Pressure fall-off data are used to check the result. Graphs of pressure and injection rate us functions of time given in the paper show that changes in permeability of the formation in the neighborhood of the wellbore are disclosed by this technique. Thus, the short-term test can he used to detect formation damage. Also, a rough measure of the radial extent of damage can be inferred, which is helpful in designing stimulation treatments. The mathematical model used for this work was a single-zone, horizontal reservoir with a damaged zone in which permeability decreased continuously as radial distance to the wellbore decreased. This model is more realistic than the usual two-zone, discontinuous permeability model used in published works; calculations indicate the realistic model is valid. Vertical variations in horizontal permeability were studied with this model, and results indicate that the permeability measured by the short-term test is the mean horizontal permeability for the vertical interval tested. The proposed short-term test thus should be useful in detecting and characterizing formation damage and in measuring formation permeability needed in calculating reservoir transmissibility. INTRODUCTION To plan the most efficient production or injection schedule for a well and to design or evaluate the optimal stimu- lation treatment, it is necessary to know the properties of the reservoir adjacent to the well, particularly the reservoir transmissibility and characteristics of a damaged zone, if one exists. Several techniques for determining reservoir transmissi-bility from well tests have been presented in the literature. 1,2,3,4 All these techniques rely on conducting constant-rate well tests that often are difficult to execute. A constant-pressure well test is generally easier to carry Out. and this paper contains the first available method for the analysis of constant-pressure well tests. Determination of wellbore damage from transient well tests has been the subject of several papers."" From these studies it is apparent that information necessary for determination of the characteristics of a damaged zone is available shortly after the transient well test is initiated. Consequently, it may not be necessary to carry out an extensive well test (for example, a pressure build-up test) if the primary purpose of the test is to detect the existence of wellbore damage. All previous studies of well testing to determine wellbore damage have been based on the two-zone perrneability model. In this model the damaged zone has a permeability k,, extending to a radius r,,, and the formation permeability k obtains from r, to the drainage radius r,.. Consequently, there is a discontinuity in permeability at r = r,,. This discontinuity can be eliminated by assuming a continuous variation in permeability through the damaged zone. The effect of this assumption on transient well tests is discussed in following sections of this paper. In addition, all formations have within them vertical permeability variations associated with lithology changes throughout the zone of interest. This paper also considers the effect of these variations on transient well tests. ANALYSIS OF CONSTANT-PRESSURE WELL TESTS The mathematical analysis associated with the injection of fluid at constant wellbore pressure into a single-zone, horizontal reservoir completely filled with a fluid of small and constant compressibility and constant viscosity is given in Appendix A. In this analysis it is assumed that the well is located at the center of an undamaged, circular drainage area. From this analysis, the formation permeability can be obtained as follows. 1, Estimate a value for the formation permeability k. 2. Prepare a plot of injection rate q vs

Jan 1, 1967

By W. A. Oates, J. A. Lambert, P. T. GaIIagher

Monte Carlo methods have been used to compute the arrangements of interstitial atoms dissolved in tetrahedral sites in bcc lattices. It is assumed that the presence of an interstitial atom "blocks " a certain number of neighboring sites and prevents their occupancy. Sites "blocked" by more than one filled site are allowed for. The computed values of. the mean occupation number (defined as the ratio of the total number of sites blocked to the number of solute atoms are used to calculate the configurational entropies of the solutions. These entropies are compared with those resulting from previous theoretical studies of this problem and also with available experin~ental data for the p Zr-H, Nb-H, V-H, and Ta-H systems. Evidence is also given that the "blocking" explanation of low limiting compositions in these systems, rather than this being due to initial limitations on the number of sites available, is probably correct. THE ideal partial configurational entropy of mixing of an interstitial solute in a metal is given by: where p is the number of interstitial sites per metal atom and Xi is the atomic fraction of the interstitial. For the bcc lattice. which we shall be concerned with in this paper, the interstitial positions are shown in Fig. 1. It can be seen that for the tetrahedral sites, p=6. whereas for the octahedral sites, p = 3. Different emphasis has been placed on the relative importance of energy and entropy effects in determining deviations from ideality in interstitial solid solutions. In some cases the same system, e.g., Fe-C, has been described by the contradictory regular and athermal solution models indicating that the enthalpy and entropy functions, derived from equilibrium data, are frequently not accu.rate enough to differentiate between these treatments. However, for certain metal-hydrogen solutions the equilibrium data is available over sufficiently wide ranges of temperature and composition to permit a reasonably accurate determination of the compositional variation of the heats and entropies. Hoch&apos; has attempted to interpret the results of interstitial solid solutions in terms of a regular solution model. In the case of the Ta-H system where 13 = 6, this model entails fitting the experimental relative partial entropies of solution, asH, to the equation: where ASgs is the relative partial excess entropy of solution of hydrogen. Hoch found that the results of Mallett and Koeh1 could be fitted to this equation with an approximately constant value of AF up to XH = 0.25. However, it is apparent from the solubility isotherms in this system which become asymptotic to the composition TaH that, since (Xh /6 - ~Xh ) becomes infinite only at TaH6, it is necessary that AS<&apos; tends to infinity at TaH. In other words, the low saturation composition of TaH, instead of the anticipated TaH,, eliminates the possibility of applying regular solution theory to such systems. Rather large negative excess configurational entropies must exist at higher hydrogen concentrations in order to explain the lower saturation values. To account for these low limiting compositions and excess entropies two distinctly different approaches have been followed. Rees and many others1-l2 have assumed that not all interstitial sites are crystallographically equivalent with respect to the interstitial addition; that is, in Eq. [I] p is less than the value anticipated from geometrical considerations. To describe, say, a bcc metal-hydrogen system with a limiting composition of MH by this approach one would consider that p = 1 in the first instance instead of p = 6.&apos;j3 In some cases, nonintegral values of B have been taken in order to improve the fit with the experimental data over limited ranges of composition. The other approach which has been used to explain the low saturation compositions is to assume that, although all sites are available for occupancy, strong repulsive interactions exist between the neighboring interstitial atoms, and hence occupancy of any site excludes or blocks a certain number of neighboring sites from being occupied. Earliest treatments of this concept considered the exclusion of an integral number, of nearest-neighbor sites from being occupied at all concentrations. In this case, the partial configurational entropy is given by: These early treatments failed to allow for the overlap of the blocked sites which will arise at all but the very lowest concentrations. More recently attempts have been made to calculate the effect of this decrease in the number of blocked sites on the configurational entropy. Using the quasichemical treatment of interstitial solid solutions as given by Lacher and assuming that an infinite repulsive interaction energy existed between the solute atoms. atom obtained an approximate configurational entropy applicable to the blocking with overlap case:

Jan 1, 1970

By V. J. Sikora, T. S. Hutchinson

This paper presents a method for making a water-rlrive ana1gsis without prior knowledge of aquifer geometry and uniformity using a standard desk calculator. Although it is necessary to know the initial oil in place to use this method, this is a minor limitation to the scope of (he method since in most reservoirs it is possible to obtain reasonably accurate volumetric oil-in-place estimates. An equation is developed which relates pressure at the water-oil contact to the water influx rate as a function of time by a factor called the resistance function. The resistance function introduces the composite effect of the aquifer geometry and flow resistance disrribution. It is pointed out that the characteristic shape of this function makes it possible to start with on approximation of the function and successively improve the approximation until the correct resistance function curve is obtained. In this fashion water-drive estimates can be made without the limitation of assuming simplifed aquifer shapes and flow distributions. This is the novel feature of this development. Methods are given for extrapolating the final curve to calculate future aquifer behavior. Equations are developed for adjusting the pressures and water influx rates where it appears that possible errors in these quantities make it difficult or impossible to obtain a useuble resistance function curve without this adjustment. Application of the pressure build-up analysis techiiique to estimate. some of the aquifer properties is also presented INTRODUCTION A great many oil pools are the result of oil accumulation in some type of trap in an otherwise large and continuous porous stratum. The void space of this stratum outside of the oil pool itself is filled with water or brine. In analyzing performance of the oil pool surrounded by this water aquifer, it is quite necessary. in most cases, to include behavior of the water. When the pressure at any point in a fluid system is lowered, such as by opening a well in an oil sand, fluids in the immediate neighborhood of this point will begin flowing towards this lower pressure sink. As pressure in this area drops due to flow towards the sink, fluids from farther out will start to flow towards the lower pressure. As more and more fluid is removed from the system, the distance from the sink or well within which flow is occurring will continually increase; that is, thc region of disturbance will grow. If some rigid boundary such as a fault is reached by the disturbance, this area will cease to grow; but if some movable boundary is reached, such as a water-oil contact, the area will grow on out into the water, although rate of growth may and almost always does change. The relation between amount of pressure drop and amount of fluid flowing at any point and at any time in an aquifer depends on such factors as compressibility and viscosity of the fluids, the porosity and permeability of the rock, geometry of the whole system and withdrawal rate or pressure drop. With these factors known, it is theoretically possible to calculate the pressure-flow behavior of the system, However, in practice true solutions to the problem are next to impossible due to complexity of reservoir systems. A number of approximate methods of solution have been developed based on various simplifying assumptions. One frequent assumption is that the water-oil contact can be located and equations defining oil reservoir and aquifer solved separately by assigning values to various parameters so as to match past pressure and production history. Usually the properties of the aquifer are not known, since few, if any, wells are drilled through the aquifer; but the water influx and the pressure at the reservoir boundary are known over some time. If it can now be assumed that the aquifer is circular, pie-shaped, or linear, and that it has uniform properties, it is possible to fit a theoretical dimensionless curve to the past aquifer performance history and therefore calculate future water drive. These theoretical curves are available in the literature,1,2,3 Of course, the assumption of uniform aquifer properties is almost always somewhat in error. Fortunately, however, moderate variations from the average have little effect on behavior of the system; and, hence, the best fit of a theoretical curve is frequently satisfactory. In other cases variations in aquifer properties and geometry are large enough that none of the available theoretical curves will give an acceptable fit of the data. In these cases, methods of fitting the data with various electrical analyzers have been developed.1,2&apos; It is the purpose of this paper to present an approach