Search Documents

By W. S. Sanner

Experiments were conducted to determine the quantity and purity of ultraclean anthracite that could be prepared in the laboratory, using conventional separating techniques. A low, a medium, and a high-volatile anthracite were studied. Specific gravity separation followed by grinding, screening, and refloating at lighter specific gravities yielded small quantities of material purified to 1.0% ash, as follows: from low-volatile anthracite, 0.7%; from medium-volatile coal, 1.0%; and from high-volatile coal, 2.3%. The reduction in ash had no appreciable effect on the volatile matter or sulfur content of the coal. The limited quantities of material containing 1.0% or less ash obtainable indicate that if ultraclean coal is desired, a practical approach should include the processing units as an integral part of a convc.ntiona1 dense-media washing plant. Anthracite preparation practices during much of the industry's history were geared to maximum recovery of the larger sizes, pea through egg, which were used in hand-fired, domestic furnaces. Those sizes were cleaned to meet trade requirements, and the smaller sizes were usually higher in ash content. The impression, thus, was that ash content of anthracite always increased as the sizes got smaller. It was presumed that it was impossible to wash the finer sizes to extremely low ash contents. Following World War 11, the classical domestic market began to decline and the finer sizes assumed an increasingly large share of the market. For example, in 1949l 61% of the production was still pea and larger. Fig. l2 shows that -by 1955 these sizes accounted for only 50%, and in 1965, they were only 39% of the total anthracite produced. With the reduction in the market, increased quality was demanded, especially in export markets, and an increasing number of dense-media magnetite washers were installed in new and existing cleaning plants. In 1965, the coarse coal prepared at 59 breakers in the anthracite fields was sampled by the U.S. Bureau of Mines (USBM); approximately 50% of these plants were equipped with dense-media magnetite washers.3 Since 1950, more attention has been focused on the use of finer sizes for numerous industrial uses, including iron-ore sintering, pelletizing, foundry coke-making, cement manufacture, production of calcium carbide, and of carbon electrodes, etc. Most industrial uses require relatively clean anthracite, such as 10 to 11% ash for sintering and foundry coke, 7% for aluminum melting-pot cathode carbon, and 5% for carbide production. Hence, cleaning methods have been improved to the extent that dense-media washers, hydraulic cyclones, and flotation cells have been introduced as fine-coal-cleaning devices. It is now possible to wash the buckwheat, rice, and barley sizes (9116 X 3132-in.) to ash contents as low as 8 to 10% in dense-media devices, but washing the -3132-in. sizes to those limits presents more difficult problems, due mainly to the fineness of the material. This USBM investigation was conducted to determine the amount and purity of ultraclean anthracite that could be prepared in the laboratory, using conventional separating techniques. The goal was to reduce the ash content to a level acceptable for industrial carbon uses, including electrode manufacture, which generally requires an ash content of 1% or less. While low-ash content is of primary importance in selecting industrial carbons, the amounts of other constituents are also important, including iron, silicon, titanium, nickel, vanadium, and sulfur. They were not considered in this work, since their occurrence in deashed anthracite has been studied elsewhere.4 The results of the present investigation were used to design a flowsheet, suggesting a method for producing ultraclean anthracite as part of conventional preparation practice, using dense-media vessels, crushers, pulverizers,

Jan 1, 1970

By C. M. Sellars, F. Maak

Electvomotie -force measurements hazle been made on ten Au-Ni alloys at temperatures 7754 825O, 900O, and 935°C using galvanic cells with solid electrolyte. Partial and ivtegral thermodynamic functions are derived from the determined activities. which show la9,ge positive deviations jvom Raoult's law. The enthalpies of mixing- are in close accord with calorinletric data and are of comparable accuracy. The results are used to determine the miscibility gap and the spinoidal curL1e, and are qlalitaticely interpreted in terms of the large size difference and electvonic interactions between gold and nickel. Tabulated nalues of the actielities and partial and inteqal free energies, entropies , and enthalpies of mixing are given in the appendix. The Au-Ni binary alloys are of particular interest thermodynamically as despite the large size difference between gold and nickel atoms they have a simple phase diagram.' This shows complete solid solubility above a miscibility gap with a maximum at about 71 at. pct Ni and 810°C. Activity measurements (over a range of temperature) have been made previously on this system using an electrolytic-cell method with a liquid chloride eletrolyte. The accuracy of this data has, however, been queried3'4 as it leads to values for the enthalpy of mixing which are too high when compared with calorimetric measurements.46 Also, the entropy of mixing is too high when compared with the values derived from specific-heat measrements. It was therefore considered worthwhile redetermining the activity data using a galvanic cell with a solid electrolyte (0.85 ZrO, 0.15 CaO). This method is based on a technique introduced by Kiukkola and aner'' and has been used to determine the thermodynamic data for Cu-Ni alloys.' THE ELECTROLYTIC CELL The cell used to determine the activity of nickel in Au-Ni alloys was: The cell reaction, which is discussed for a similar cell by Rapp and Maak, gives a direct measurement of the activity of nickel in the alloys as where aNi(alloy) is the activity of nickel in the alloy, 3 is Faraday's constant, E is the measured electromotive force, R is the gas constant, and T is absolute temperature. For the valid application of this type of cell the arrangement must be ther modynamic ally stable, which means that: 1) The free energies of formation of the oxides fulfill the relation 2) The vapor pressures or dissociation pressures of each component must be so small that during the time of the experiment the amount of any material transported by the gas is negligible. 3) Side reactions which change the alloy composition must be negligible. These reactions are possible both with the surrounding gas atmosphere and with neighboring solid material. The first two conditions are fulfilled by the above cell.' The third condition is approached by care in the experimental technique. The factors which control the attainment of equilibrium in this type of cell have been considered by Rapp and aak, who conclude that metallic diffusion in the alloy is the rate-control ling step. A lower temperature limit for the practical application of this cell is therefore set by the condition that equilibrium must be attained in reasonable times. EXPERIMENTAL MATERIALS AND TECHNIQUE The starting materials were mint gold (99.99 pct Au) and Mond nickel (99.9 pct Ni). Alloys, nominally 5 at. pct Ni and then each 10 at. pct across the system, were produced by vacuum-levitation melting compacts of mixed drillings of the pure metals. The nickel used for the reference tablets in the cells was melted in a similar manner. The alloy ingots, cast in boron nitride molds, were cold-forged and

Jan 1, 1967

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967

By C. C. Martin, M. C. Fuerstenau, R. B. Bhappu

Experiments revealed that quartz could not be floated in conductivity water at any pH with a long-chained sulfonate as collector. Various cations, Fe+++, Al+++, Pb++, Mn++, Mg++, Ca++, are shown to function as activators when the pH is such that the cation hydrolyzes. As the solubility product of the various cation-sulfonates was exceeded at the concentrations involved, precipitated cation-hydroxy-sulfonates, e.g. Fe (RSO3)2OH must be functioning as the collector in these systems. A number of interesting and important phenomena were revealed during a study of the response of beryl to sulfonate flotation.' Perhaps the most important observation is the effect of various cations on flotation response. Experiments showed that in the vicinity of pH 3.0, about three ppm Fe+++ increased recovery of leached beryl from about 60 pct to 100 pct, whereas about four ppm Ca+ + reduced the recovery from 60 pct to 2 pct. With the exception of ferric iron, all of the other cations investigated reduced flotation recoveries drastically at this pH. As postulated in the paper Sulfonate Flotation of Beryl,' even with low concentrations of polyvalent cations and sulfonate (on the order of one ppm), cation-sulfonates are precipitated immediately in solution. That is, blueish-white clouds were noted to form after given additions of salts to known solutions of sulfonate. Furthermore, it was also shown that one-hundred-fold more ferric iron than necessary to form the precipitated cloud had to be added to effect complete flotation. In other words, some form of solid ferric sulfonate was found to function as the collector in these systems. The form of the compound was suggested to be Fe(RCO3)2OH and that bonding to the surface occurs through the hydrogen of the hy-droxyl. As the form of this compound was shown to be dependent on the extent of hydrolysis of ferric iron, it is to be expected that ferric iron will be unique in some pH region. Further, if the presence of hydroxyl in the neutral precipitated compound is of primary importance, then other cations should function similarly but at different values of pH, since they will hydrolyze at different concentrations of hydroxyl ion. To determine whether the premise of hydrolysis is correct, the study was extended to quartz systems in the presence of various cations at various values of pH. EXPERIMENTAL MATERIALS AND METHOD Sodium alkyl aryl sulfonate was chosen as the collector for the experimental work. The reagent, sup plied as a courtesy sample by the Shell Chemical Co., has the following physical properties: Physical form Solid — Finely ground Sulfonate content, pct wt 95 to 97 Molecular Weight 450 to 470 Number of carbon atoms in hydrocarbon chain 25 to 30 Melting Point, OF 250 to 260 Solubility Soluble in water with gel formation at concentrations higher than 25 pct by weight sulfonate This particular reagent was chosen because of its high active sulfonate content, long-chain length, and solid form. With the exception of the frother and sulfonate, all other chemicals were of reagent grade quality. Conductivity water was used in the experimental work. This water, made by passing distilled water through an ion exchange column, had an average measured conductivity of one micromho. Experiments were conducted at room temperature in a small glass flotation cell with the following procedure: 1) A redetermined amount of water and salt solution (e.g. CaC12 . 2H2O) were combined and the pH adjusted to a given value. 2) A given amount of sulfonate was added so that final solution volume was 130 cc. 3) One drop of Dow 250 frother was added. 4) Five g of quartz were added and the system conditioned for five minutes. 5) The pH of the system was measured (termed flotation pH).

Jan 1, 1963

By J. C. Todd

The volume of air needed to move the combustion wave through each acre-foot of the reservoir is a very important quantity for engineering economic analyses. A new method, which involves backflowing the injection well, has been developed for determining the air requirement .for the forward combustion thermal pilot in the Delaware-Childers field. The method depends on measurement of the total volume of air injected and measurement of the pore volume in the zone swept out by the combustion front. The basic premise of the method is that the pore volume of the swept zone contains air, while the formation around this zone contains flue gas, oil and water. The pore volume of the swept zone is determined by blowing down the pilot section of the formation through the injection well after the pilot is shut down. Oxygen balance is used to differentiate between air from the swept zone and flue gas from the surrounding formation. Data indicate that air'is produced from the swept zone by relatively efficient miscible displacement by the flue gas. The pore volume occupied by the air in the swept zone is determined by simple PVT relationships. The average temperature of the swept zone is calculated by heat balance from the total amount of heat generated during the pilot and losses to the surrounding formation. The new method was applied successfully at the conclusion of the pilot in the Delaware-Childers field; the air required for the reservoir conditions involved was calculated to be 19.2 MMscf/acre-ft of reservoir swept. INTRODUCTION One of the most important pieces of information to be found from a forward combustion pilot is the air requirement; that is, the volume of air needed to move the combustion wave through each acre-foot of the reservoir. Compression costs per barrel of oil recovered, which are a direct function of the air requirement, are a large portion of the expenses incurred in the thermal recovery process. Thus, the profitability of thermal recovery depends largely on the air requirement, which is a measure of the amount of fuel deposited and burned in the thermal recovery process. Several attempts have been made to correlate crude oil properties with fuel deposition and, hence, air requirement.l*2 Other attempts have been made by our laboratory to find the air requirement by material balance, gravity survey, pressure transient analysis,3 magnetic survey and the poten-tiometric model. Better methods are needed to predict more accurately the economics of fieldwide application of the process. Another method of finding the air requirement consists of operating a pilot thermal project for several years. This is followed by an extensive coring program to define the volume swept by the combustion wave. The air requirement is found from the total injected air divided by the volume of the swept region found from the coring operation. This method is costly; an adequate coring program usually costs several hundred thousand dollars. A thermal pilot was operated in a watered-out portion of the Delaware-Childers field, Bartlesville sand, from Nov., 1960, until May 3, 1965. One of the primary purposes of this pilot operation was to determine the air requirement for economic calculations which determine the advisability of expanding the project to the entire 128,000 acre-ft field. Barnes4 has described the result of this pilot operation,, including calculation of the air requirement from limited coring data. The value obtained, as calculated, was 15.8 MMscf/acre-ft. One purpose of this report is to describe a new method of finding the air requirement from a thermal pilot operation. A second purpose is to describe the application of this method to the Delaware-Childers thermal pilot. THEORY AND DESCRIPTION OF THE METHOD The combustion zone in the forward combustion process is propagated from injection w producing

Jan 1, 1970

By H. C. Brassfield, R. E. Fryxell, E. C. Duderstadt, E. A. Aitken, P. K. Conn

Permeation rates of hydrogen through are-cast tungsten were measured at temperatures from 1300" to 2600°C with hydrogen pressure differentials of 1 and 0.1 atm across isothermal membranes. Rates were nearly Proportional to square root of pressure, were inversely ProPortional to membrane thickness, and were temperature-dependent according to the equation: p = permeation coefficient in cu cm (stp)-mm per sq cm-min-atmlr2, R = gas constant in cal per mole-OK, T = absolute temperature in OK. Permeation rates of oxygen through tungsten were deternlined in the range 1965" to 2300°C and at calculated oxygen (02) pressures of 3.3 x 10-4 to 2.2 x at?n using sealed capsules containing Kmown quantities of tungsten oxide. Permeation coefficients obtained in this manner had a temperature dependency in the range of 40 to 45 kcal per mole of permeating species. Permeation coefficients for oxygen were higher than those for hydrogen in the temperature range studied. This is probably the result of oxygen solubility being much greater than hydrogen solubility at these temperatures at a gzven pressure. The purpose of this investigation was to determine experimentally the permeabilities of hydrogen and oxygen through polycrystalline tungsten at temperatures ranging from 1300° to 2600°C. Permeation rates of hydrogen through tungsten have been measured at temperatures up to 925°C,but no previous measurements have been reported for the temperature range covered in this investigation. Some diffusion measurements3&apos;&apos; and solubility measurements3 have been made for hydrogen at temperatures below 2200°C. No previous permeability studies of oxygen through tungsten are known to the authors; however, efforts have been made to determine diffusion rates5-&apos; and solubilities PERMEATION THEQRY Permeation has been defined by ortton" as the overall steady-state flow process from the gas phase on one side of a membrane or wall to the gas phase on the other side and includes the steps of adsorption, solution, diffusion, exsolution, and desorption; in a given temperature range, any of these processes may be rate-controlling. For diatomic gases, which dissociate prior to diffusion, the quantity of gas which permeates through a membrane or wall at a given temperature is directly proportional to the surface area of the membrane, the length of time of permeation, and the difference between the square roots of the pressures on each side of the membrane, and it is inversely proportional to the thickness of the membrane when diffusion-controlled, thus: q = total amount of gas permeating a membrane, p = temperature-de pendent permeation coefficient, A = area of membrane exposed, / = time, p i = gas pressure of high-pressure side, p2 = gas pressure of low-pressure side, d = thickness of membrane. When the surface reactions are not rate-controlling steps and the solution of the gas in the metal follows Henry&apos;s law, the permeation coefficient, P, is equal to the product of diffusivity and solubility. The permeation coefficient is also temperature-dependent according to the equation: Po = permeation constant, Q = activation energy, R = gas constant, T = absolute temperature. MATERIALS The tungsten materials used in this investigation were greater than 99.9 pct pure. The interstitial impurities ranged from < 1 to 50 pprn C, 7 to 35 ppm O, less than 15 pprn N, and undetectable quantities of hydrogen. Previous high-temperature studies of are-cast tungsten have shown that thermal treatment similar to that which the hydrogen permeation cells experienced reduces the carbon content from starting levels of 50 to 60 ppm to about 30 ppm. Metallic impurities were in the low (1 to 50) ppm range. The argon gas used in this investigation contained about 8 ppm N. None of the following gases were

Jan 1, 1968

By R. L. Schmidt, J. D. Verhoeven

The diffusion coefficients and electrotralzsport mobilities of aluminum, gallium, and arsenic have been determined in molten germanium with the capillary reservoir technique. The diffusion coefficients are slightly higher than those determined by the solidification technique of Brirton et a1.1 It was found that aluminium and gallium migrate to the cathode in molten germanium) whereas arsenic migrates to the anode under an electric field. A possible technique for producing p-n junctions by means of electric field freezing is discussed in view of the experimental data. In the past decade or so the literature has abounded with articles concerning the solid-state properties of the semiconducting elements germanium and silicon. It is surprising by comparison how little has been published concerning the liquid-state properties of these materials. In the production of doped germanium or silicon the diffusion coefficient of the dopant in liquid germanium or silicon is an important quantity. Yet, to the best of our knowledge, no one has ever made a direct measurement of the diffusion coefficient of solutes in either molten germanium or silicon. A number of authors&apos;-4 have determined some diffusion coefficients by the indirect technique developed by Burton et al.1 This technique appears to give numbers of the correct order of magnitude but it does depend upon a number of assumptions in the analysis and also requires prior knowledge of the viscosity of the melt. In the present work a more direct measure of the diffusion coefficients of arsenic, gallium: and aluminum in molten germanium has been made. It has been shown independently by two groups of authors5, 6 that, when an electric current is passed through the interface of a solidifying metal containing an impurity, the resulting electrotransport in the liquid boundary layer at the interface offers a mechanism for control over the redistribution of the impurity accompanying the solidification. These authors show that the degree of control is proportional to the mobility, U of the impurity atoms in the liquid boundary layer, where the mobility is defined as the drift velocity relative to the solvent per unit electric field. In principle the results apply to the solidification of a doped semiconductor and offer an interesting technique for control of the doping of a semiconductor. In the present work the mobilities of aluminum, gallium, and arsenic were determined in molten germanium in order to evaluate the potential usefulness of electro-transport as a technique for solute control when applied to the solidification of doped germanium. EXPERIMENTAL PROCEDURE Both the diffusion coefficients and the mobilities were measured using high-purity alumina capillaries and a resistivity technique for analysis. A schematic diagram of the electrotransport apparatus is shown in Fig. 1 and the details of the cell design are shown in Fig. 2. For the mobility experiments electrical contact with the bottom of the capillary was made with graphite, as shown in Fig. 2. In order to prevent electrical contact with the melt the graphite was coated with a thin layer of alumina using a plasma spray torch. The cell assembly, a graphite electrode, and another assembly containing two blank capillaries were fastened to the lavite disk shown in Fig. 1. To make a run the system was evacuated and the doped germanium was outgassed at the operating temperature of 1030°C. The cell and blank capillaries were then submerged beneath the melt and a slightly positive pressure of helium was introduced to fill them. After passing current through the cell for a given time the cell and blank capillaries were raised from the melt and the furnace was shut off. The temperature of the cells dropped below the freezing point of germanium within 1 to 2 min after the furnace was shut off.

Jan 1, 1968

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study&apos; that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,&apos; flow in packed columns,&apos; fracture flow capacity&apos; and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By Milton E. Wadsworth, W. Martin Fassell, William H. Dresher

A study of the rate of dissolution of molybdenite (MoS2) in alkaline solution was carried out under carefully controlled conditions. Effects of temperature, oxygen over-pressure, and KOH concentration were evaluated. Studies were made in the temperature range 100°C to 175°C and in the pressure range 0 to 700 psia of oxygen. Under these conditions molybdenite was found to leach according to a linear mechanism. Both oxygen over-pressure and KOH concentration were found to control the rate of leaching. The mechanism has been explained in terms of adsorption of oxygen at the molybdenite surface followed by configurational rearrangement of the adsorbed molecules. The hydroxyl ion dependency is believed to be diffusion-controlled. Laboratory batch studies have shown that molybdenite can be readily dissolved in alkaline solutions under moderate conditions of temperature and pressure. Commercial application of this process to the production of ferro-molybdenum and molybdenum chemicals is promising in view of the ease of dissolution of molybdenite and the relatively noncorrosive conditions involved in the process. HIGH temperature-high pressure techniques have long been used to great advantage in the organic chemical industry, the petroleum industry, and the paper industry. Only recently, however, have these methods been used to extract metals from their ores on a commercial scale. The Chemical Construction Corp., together with interested producers of nickel and cobalt, has done much to develop methods of producing nickel and cobalt powders by the use of high temperature-high pressure techniques.&apos; The Howe Sound Co. is currently operating a plant near Salt Lake City in which the new Chemical Construction Corp. technique is applied to cobalt-arsenic sulfide concentrates. In this method the concentrates are leached in a dilute solution of sulfuric acid under air pressure at a high temperature. Near Edmonton, in Alberta, Canada, the Sherritt Gordon Co. is also using the new method in conjunction with its ammoniacal leach process.&apos; National Lead Co., at Fredricktown, Mo., is operating a process similar to that of Howe Sound. In spite of the rapid growth of high temperature-high pressure hydrometallurgical processes, there has been very little work done on fundamental principles involved in this type of process. Conse- quently the literature lacks any significant data in this field. Oxidation of pyrite,3 phalerite,4 galena,5 and recently molybdenite" in alkaline solutions has been reported to some extent. This investigation was initiated, therefore, to broaden the application of high temperature-high pressure hydrometallurgy and also to contribute some insight into the basic mechanisms of oxidation in an aqueous medium under these conditions. In addition to the immediate consequences of these processes much valuable information remains to be learned about fundamental reactions when they are carried out at temperatures and pressures slightly higher than those normally encountered. Leaching of molybdenite by alkaline solutions under oxygen pressure has recently been reported in the Russian literature by E. S. Usataya,8 who found that leaching rate increased with an increase in the pH of the solution—the strongest effects being observed at a pH of 10—and also with an increase in temperature. Usataya also has reported the formation of a protective film on the surface of the molybdenite when the leaching solution was weakly alkaline, neutral, or acidic. He believes that this oxidation process explains the migration of Mo+ + in the oxidation zones of molybdenite ore deposits and the subsequent impoverishment of the ores. Equipment: The high temperature-high pressure reaction unit used in this work was especially designed for kinetic studies of this nature, emphasis being placed on accuracy of measurement and precision of control. As details of the unit&apos;s construction and operation have been discussed elsewhere,7 Only basic features and principles of operation will be discussed at this point.

Jan 1, 1957

By E. E. Underwood, L. L. Marsh

Aluminum single crystals, alloyed with 0.042 atomic pet Cu and 0.11 and 1.1 atomic pct Mg, were subjected to constant stress creep tests, tensile tests, and hot hardness measurements within a temperature range of 300° to 866OK. Calculations based on Dorn&apos;s temperature-compensated time parameter, 6, gave a value of DH, = 27,000 cal per mol for the activation energy of early creep in aluminum single crystals. Correlations have been obtained for aluminum alloy single crystals with the parameter E for solid solution strengthening, as well as with the parameter F for solid solution hardening, by using a valence of three for aluminum. Limited measurements on tensile specimens show that the slip band density tends to decrease with increasing temperature and with decreasing solute concentration. INCREASING interest is being shown in the mechanisms of plastic deformation in single crystals during tensile and creep testing. The complexity of deformational processes in polycrystalline materials has led to a search for simpler experimental conditions. Numerous creep and tensile investigations have been conducted with pure, metallic single crystals. To a lesser degree, the effects of alloying elements on the tensile properties of single crystals have been determined. However: the literature dealing with the effects of alloying additions in single crystals under creep conditions is vanishingly small. This paper represents an attempt to narrow this gap in the knowledge of the subject. Previous investigations of creep behavior at the Battelle Memorial Institute1-1 have been of great value in the analysis of the present single crystal data. It is equally desirable, however, to ascertain the extent of correspondence between the behavior of single crystal and polycrystalline materials. For this purpose, correlations, similar to those developed for polycrystalline aluminum alloys by Dorn and co-workers, have been made with the single crystal data. The results from this study have tended to confirm and extend those correlations to the case of alloyed single crystals. Materials and Procedures Three dilute binary aluminum alloys were prepared for this investigation from 99.99 + wt pct Al, 99.8 wt pct Mg, and 99.92 wt pct (electrolytic) Cu. The nominal compositions of the alloys were 0.042 atomic pct Cu, and 0.11 and 1.1 atomic pet Mg. Precautions were taken to avoid contamination of the stock during melting and casting of the alloys. After a heat treatment of 8 hr at 925°F the alloys were extruded, then machined into) threaded tensile specimens with a 3 in. reduced section and a 0.505 in. diam. Spectrographic examination showed less than 47 ppm metallic impurities in each alloy. The single crystals were grown by the strain-anneal method, with a critical strain of about 11/4 pet giving the optimum results. In general, 3 in. crystals were obtained with the Al-Cu alloy, but smaller crystals in the magnesium alloys necessitated the use of 2 in. and 1 in. gage lengths with the low and high magnesium alloy specimens, respectively. After an electrolytic polish, the orientation of that portion of the specimen containing the largest single crystal was determined from Laue back-reflection photographs. Tensile tests were conducted at a constant load rate of about 2 lb per min. Creep runs were made in a constant temperature room, under constant stress at the higher creep temperatures, and constant load at the lower temperatures. The eloneation was measured to within ±5 microin. by a specially designed capacitance extensometer. The ex-tensometer arms were attached to the 3 in. specimens at the shoulders of the test piece or, where the crystals were smaller than 3 in., by knife-edge grips.

Jan 1, 1957

By N. M. Levine, W. M. Fassell

In this laboratory study the problem of aerative conditioning to separate chalcopyrite and pyrite from cobaltite was simply effected with a sulfy-drate collector and pH by proper choice of mixing variables. It was shown that for a given impeller type and geometry, selectivity is governed by a relation between the Froude number and a modified power number multiplied by impeller tank-diameter ratio wherein the power factor is expressed as power intensity or horsepower per unit volume of contained pulp. In applying this approach to other systems, other criteria would be used in place of selectivity. The purpose of this paper is to show how solutions to problems of aerative conditioning of pulps may be enhanced by development of mixing parameters in a given system through application of a fundamental approach to agitation. The term aerative conditioning refers to a combination of processes which include mechanical gas dispersion, gas solution in the liquid phase, transport of dissolved gases and other reactants to the reaction zone, physico-chemico interaction within the liquid phase and/or at the solid surface, and, finally transport of products away from the reaction zone. Important elements of these processes are: 1) Air-Liquid Contact Area: This area limits the rate at which a given gas may dissolve in a given liquid. Thus, an aerative conditioning system should be designed for optimum gas dispersion. 2) Diffusional Factors: Such factors as concentration gradient, diffusion film thickness, and diffusivity may govern the rate at which dissolved constituents are transferred to the zone of interaction at a particle surface and the rate at which reaction products are transported therefrom. Thus, such systems would be designed for maximum turbulence to effect good mixing and to minimize diffusion film thickness. 3) Solids-Turbulent Liquid Contact Time: Sufficient time should be provided for interactions to take place at solid-liquid interfaces. The region of maximum turbulence for mechanical systems being in the impeller zone, it is important that the number of passes through this zone be optimized. This means that a proper balance should be struck be- tween flow and turbulent power to achieve optimum conditioning time. Too high a flow rate and insufficient turbulence may be improper if the reactions are diffusion limited. On the other hand, too much turbulent power and too little time in the reaction zone due to insufficient flow power may also be limiting. 4) Geometry and Impeller Design: The design of the conditioning vessel and the type of impeller employed are important factors in determining the relative distribution between flow and turbulent power and the degree of gas dispersion. Maximization of power input would be to no avail if the distribution between flow and turbulent power were improper or gas dispersion insufficient. 5. Physical Chemistry: The reactions between dissolved gases and solid surfaces may be limiting if gas solubility, activation energies and/or free energies are not favorable. J. H. Rushton, E. W. Costich, and H. J. Everett&apos; have summarized the mechanical elements of the combination of processes listed in a fundamental mixing equation. This equation was derived by dimensional analysis and relates the physical variables for a single impeller, centered in a cylindrical, vertical axis, flat-bottomed tank. It is, where T is tank diameter in feet; H is liquid depth in feet; C is height of impeller off tank bottom in feet; S is pitch of impeller; L is length of impeller blades; W is width of impeller blades; R is number of baffles; D is impeller diameter in feet; Np is the power number, pg/dN3D5; P is power in foot-pounds per second; g is the gravitational constant, feet per second; J is width of baffles; B is number of im-

Jan 1, 1961

By G. T. Pruitt, T. L. Irby, P. B. Crawford, H. Ferrell

In water flooding peripheral, center-to-edge, line-drive or water-encroachment patterns the question has arisen, "when should a producing well be converted to a water-injection well?". It is realized that, if the producing well is shut in when water breakthrough occurs, some oil could be left in an unswept area. Yet, operation to high water-oil ratios may require processing much additional water. If a well is shut in, "where should the allowable be transferred?". A review of the literature indicates that very few data are available to the engineer to assist him in answering these questions. It is the purpose of this paper to present information which may be of value in selecting the flooding technique and the injection and producing wells and in scheduling injection and production rates. Certain phases of these problems can be studied in the laboratory. The Potentiometric model is well suited for studying many waterflooding problems. This model assumes that steady-state flow exists, that the mobility ratio is equal to one and that gravitational effects are neglected.&apos;-&apos; The three-phase sand model also has been recognized as an effective means for studying certain reservoir problems. Both models were used in making the study reported here. This paper discusses seven cases utilizing various waterflooding techniques. Data for the first five cases were obtained on the potentiometric model. Data for Cases VI and VII were obtained by constructing and flooding an unconsolidated sand pack (a three-fluid-phase system) and by determining the swept area and cumulative oil recovery for various types of operations. APPARATUS AND PROCEDURE The potentiometric model reservoir was 30 in. long X 10 in. wide by about 1/2 in. deep. Four wells were equally spaced in a line array along the lateral center line. (See Fig. 1.) Well No. 1 initially was an injection well for each of the patterns reported. The cases studied with the potentiometric model and sand model are described herein. The laboratory sand model reservoir was a 30- 10-in. rectangular Lucite vessel packed with sand to a depth of about 0.65 in. Four wells were equally spaced in a line array along the lateral center line. The reservoir material was an unconsolidated sand into which the reservoir fluids were mixed to yield fluid saturations of 55 per cent oil, 25 per cent water and 20 per cent gas. The water was added to the sand and the material mixed, after which the oil was added. Kerosene and water, the fluids used, gave a viscosity ratio of 1.64:1. In packing the reservoir, care was taken to assure constant fluid saturations and constant porosities. A new mixture was prepared for each run. During water flooding, a colored water was injected by a constant rate pump. The injection rate was 1.1 B/D/ft of sand thickness. For this series of tests, each of the producing wells was open to the atmosphere and, as in the potentiometric model study, Well No. 1 initially was an injection well for each of the patterns. KESULTS OF POTENTIOMETKIC MODEL STUDlES Two cases (Cases I and 11) involving equal pressure injections were studied, with the following results. Case Water injection was initiated at Well No. 1; the producing wells (Nos. 2, 3 and 4) were operated at constant and equal bottom-hole pressures. When water broke through at Well No. 2 and achieved the selected water-oil ratio, Well 2 was converted to a water injection well, resulting in two injection and two producing wells. The injection pressures at Wells 1 and 2 were maintained equal and constant; pressures at the producing wells were held equal, but they were lower than those at the injecting wells. A similar process was repeated when water broke through at Well No. 3. The flood was terminated when

By F. H. Deily, Jean M. Quets, Milton E. Wadsworth, John R. 222-000-000-012 Lewis, D. S. Rowley, R. J. Howe

Samples of synthetic magnetite were reduced in hydrogen-water vapor atmospheres in the temperature range 450o to 900oC. The reaction was always surface controlled, indicating the final products of reaction are nonfirotective. Three separate cases have been identified from the kinetic results: Case I, reduction of magnetite to iron below 570°C; Case II, reduction of magnetite to iron above 570°C; and Case 111, reduction of magnetite to wustite above 570 oC. In Cases I and 111 the kinetics have been explained by the formation of oxygen anion vacancies in the magnetite surface followed by associated diffusional processes. In Case 11 the observed kinetics have been related to the concentration of iron cation vacancies in the wustite layer formed, followed by associated steady-state diffusional processes. The concentration of cation vacancies in the wustite is maintained constant by the magnetite-wustite equilibrium. THE reduction of the oxides of iron has been a subject of study and interest for many years. However, few attempts have been made to interpret the results in terms of the mechanism of reduction of an iron oxide by either hydrogen or carbon monoxide. Stalhane and Malmberg&apos; using CO, H2, and CO-Hz mixtures, developed the rate expression: where m = mass of oxide reduced, t = time, k = a constant, s = surface area, P = partial pressure of of CO or H,, and P, = partial pressure of the reducing gas at equilibrium. The rate of reaction per unit area was proportional to the difference between the partial pressure of the reducing gas and its pressure at equilibrium. Hansen, Bitsianes, and Joseph2 studied the reduction of pellets of hematite to magnetite with mixtures of CO-CO,. They concluded that below 450°C the reaction proceeded according to a model in which CO reacts with the hematite surface producing an oxygen ion vacancy. The rate controlling step was then assumed to be oxygen ion diffusion through vacancies resulting in the formation of magnetite. They derived an expression giving the rate of reduction proportional to the cO/CO, ratio. Recently McKewan3 published results for the reduction of magnetite spheres in atmospheres below 570°C. He explained his results by an equation of the form where K, is the equilibrium constant for wustite-iron equilibrium attributed to the formation of wustite as an intermediate in the reduction of magnetite. According to McKewan, Kp is an equilibrium constant for H, -HO adsorption in which oxygen is deposited on the surface from the H,Oand subsequently acts as a poison for reduction by hydrogen. Sample Preparation and Experiments. In this work the authors used the same apparatus and experimental procedure employed in the hydrogen reduction of magnetite. Sintered disk samples weighing approximately 1.5 g and measuring 1.2 cm in diameter and 0.25 cm in thickness were prepared by a method described previously.4 Partial pressures of hydrogen and water vapor were controlled by bubbling hydrogen through a series of flasks containing distilled and deionized water maintained at constant temperature in an oil bath. The total flow rate of the reacting gases—hydrogen and water vapor—was approximately 0.5 liter per min. The entire external system was maintained at a temperature above the oil bath to prevent condensation of water vapor. This was accomplished by wrapping electrical resistance elements over all surfaces. The glass column containing the weighing system-a McBain balance suspended from a gold chain-was enclosed in an aluminum tube heated by nichrome resistance wire enclosed in an insulating material. The temperature of the glass column was held constant by means of a series of thermocouples connected in parallel to a Leeds and Northrup temperature controller. This was necessary since variations in temperature resulted in expansion and contraction of the glass column and could be readily observed with the cathetometer. The spring extensions were measured through a window with a gaertner cathe-

Jan 1, 1962

By J. C. Griffiths, D. W. Bennion

A mathematical model, which does not assume a priori that stratification exists, but was designed to test for the stratification was developed. The model segmented the reservoir horizontally into areas of common variance, then divided it vertically into strata (if strata were present). Next, trend surface techniques were used to determine the lateral extent and variation of each stratum. The model was tested on two reservoirs (sandstone and limestone) using porosity as the reservoir rock property. The sandstone reservoir contained approximately 60,000 samples from 2,000 cored wells, while the limestone reservoir had approximately 24,000 samples from 430 cored wells. Within the areas of common variance tested, the model was able to distinguish four separate lithological units or zones in the sandstone reservoir, four in the Marly section and seven in the Vuggy section of the limestone reservoir. INTRODUCTION To accurately predict the movement and production of fluids from a reservoir rock, it is necessary that the reservoir engineer have a knowledge of the rock properties and their distribution throughout the reservoir. For this study, a reservoir rock is defined as a solid containing interconnecting holes or voids which occur relatively frequently and are dispersed within the rock in a regular or random manner. The fraction of bulk volume that these interconnecting voids occupy is called porosity. The rock must also possess sufficient permeability to allow fluids to move through it and it may be composed of one or more strata or lithological units. A stratum or lithological unit as defined in this study is a body or volume of reservoir rock whose properties are so distributed that its lateral and vertical extent can be traced throughout or through a portion of the reservoir. Operationally, the existence of strata are question- able unless the variance within strata is less than among strata. Within the past few years computers have been widely used as a tool for calculating movement and production of fluids from reservoirs. Most prediction models require some type of a mathematical description of certain reservoir rock properties. In the past, most of these models have been relatively simple. Some have assumed vertical variation but no lateral variation, while others have assumed lateral variation but no vertical variation. In general, most of these models have not used functional relationships to predict reservoir rock properties as functions of position. The object of this study was to develop a model which would predict lateral and vertical variations of reservoir rock properties. Since samples are usually available from only a small portion of the total reservoir rock, it seemed logical that if measurements from these samples were to be used to infer the properties of the actual reservoir, the data should be treated statistically as a sample from the total population (reservoir). Therefore, the model developed in this study to predict reservoir rock properties was a stochastic model. The model was tested using porosity (a macroscopic reservoir rock property) measurements from two reservoirs — the first, a sandstone reservoir from which 60,000 samples were obtained; and the second, a limestone reservoir with 24,000 samples. Several investigators have proposed methods to determine the vertical and areal variation of a reservoir rock property. Stiles: Dykstra and Parsons 2 and Suder and Calhoun3 have all developed waterflooding prediction techniques which assume vertical variation but no lateral variation. Stiles developed a method to segment a reservoir arbitrarily using frequency distribution. Law4 suggested that porosity has a normal frequency distribution and that permeability has a log-normal frequency distribution. Elkins and Skov 5 proposed that when a reservoir is segmented vertically, the

Jan 1, 1967

By A. N. Chirico

This paper reviews the three major distillation processes: multiple effect (LTV) evaporation, multi-stage flash distillation and vapor-compression forced circulation evaporation. Scale preventative measures are discussed for all saline water applications. Operational data is presented on the Freeport, Texas, demonstration plant, which is the first government facility to produce one million gallons per day of fresh water for a municipality. There has been much publicity concerning the water crisis that confronts many peoples today. There is little doubt that the arid regions - which include the Middle East, parts of Africa and the Caribbean — have suffered to no slight extent for the lack of this most precious commodity. But what about our own country? Government experts have extrapolated demands and flatly predict that by 1980 we will be faced with a shortage of 80 billion gal per day. The validity of this prediction has been disputed by others; however, it appears certain that there will always be a natural abundance in some areas, a shortage of supply in other areas. Distributing our total natural water resource equitably would be impossible. Transmission of water by means of aqueducts can be an expensive undertaking. The controversial Feather River Project in Calif. is an example. The estimated cost for this venture is a staggering three billion dollars. As our population continues to increase and our standard of living rises, agricultural and industrial requirements will be proportionate. What can be done to avoid a crisis? We can conserve and use our water for industrial and irrigation purposes more judiciously: we can attempt to eliminate pollution of our natural resources from sewage and industrial wastes; we can develop, while there is time, saline water conversion processes for the massive production of low cost fresh water. The scope of this paper will be limited to only the distillation type processes that are commercially feasible. The first step in making potable water is to find a source of raw material. An inexhaustible supply is provided by oceans in coastal regions, but the supply of water is more critical in the inland communities of our country. The source of raw material in these. regions is brackish well water and this must be utilized. There is a significant process difference when handling brackish water compared to seawater and this brings out an important fact: No single conversion method can be universally applied to solve all the water problems. Composition of the raw water and degree of distillate purity are factors which must be considered before a particular process is selected. Local conditions at the plant site, fuel costs, labor costs, disposal facilities, etc. are other equally important factors which must be evaluated. WATER CONVERSION PROCESSES Conversion processes commercially feasible today are: 1) multiple effect evaporation, 2) multiple flash distillation, 3) vapor compression distillation and 4) electrodialysis. The first three are distillation processes which take potable water out of the raw water, while the fourth is a membrane process which takes out the salt. In addition to the above there are a good number of other promising methods for production of potable water, many of them in the basic research stage. Freezing is another promising process with commercial potential. This process hinges upon the formation of ice crystals which are free of salt occulsions. Separation and washing of the crystals are problems that must be solved before this process can be considered commercially feasible. Freezing processes have several advantages mostly attributable to the low temperature of operation. This low temperature reduces the scaling and corrosion problems encountered in higher temperature operations. Several different freezing processes are being considered and will be tested in the pilot plant at Wrightsville Beach. One process is a flash freeze process, a second uses a secondary refrigerant such as butane for freezing, and a third process employs a secondary refrigerant particularly controlled to produce large ice crystals. Of the above processes, only the distillation processes are reliable when water of high purity is re-

Jan 1, 1963

By J. E. Dorn, L. A. Shepard

LOW and Gensamer&apos; demonstrated a number of years ago that the yield point phenomenon in mild steels was associated with the presence of fer-rite soluble carbon or nitrogen. More recently the yield phenomenon in body-centered-cubic metals containing interstitials was rationalized by Cottrell&apos; in terms of a simple dislocation model. Interstitial atoms interact with dislocations in two ways; they cause not only local expansions but induce local tetragonality in the lattice. Consequently, interstitial~ interact with the hydrostatic tension and shear components of the stress about dislocations. They tend to migrate toward the expanded regions of edge dislocations and to assume sites that relieve the shear stresses of screw dislocations. Thus, a dislocation saturated with solute atoms constitutes a lower free energy state than that obtained when the dislocation threads through the average composition regions of the matrix. A greater stress will be required to separate the dislocation from its atmosphere than to move the dislocation through the matrix. This factor gives rise to the upper yield stress, which is required to unleash a series of dislocations in a localized region. This local yielding is propagated across the specimen to form a thin band of plastically deformed material known as a Lueder&apos;s band, making an angle of about 45" to the stress axis. Once the band has formed, deformation continues at the lower yield stress by the spreading of the Lueder&apos;s band in the direction of the applied stress. Undoubtedly the spreading of Lueder&apos;s bands at the lower yield stress is accomplished by the high stress concentrations at the band fronts, which serve to induce continued unlocking of new dislocations in advance of the migrating band fronts. Cottrell and Bilby have shown that the dependence of the yield point on temperature can be deduced by assuming that thermal fluctuations aid the stress in unlocking small dislocation loops from their solute atmospheres. Once a loop that exceeds a critical size has been nucleated, the entire locked diqlocation is released and can migrate. Fisher4 simplified this analysis by assuming that the locking forces were short range, so that if the dislocation loop were displaced only one Burgers vector from its atmosphere, it would be unlocked. Applying his model to the special case of delayed yielding under a constant stress of the order of the upper yield strength, he demonstrated that the delay time 7 for yielding should depend on stress and temperature according to where A and B are constants, G is the shear modulus, and u and T are the resolved shear stress and absolute temperature, respectively. Cottrell and Bilby, Fisher, and Fisher and Rogers&apos; have shown that the above deductions are at least in qualitative harmony with the experimental facts. A number of investigators5-0 have shown that the yield point phenomenon can also be induced in sub-stitutional alloys of face-centered-cubic metals. In general such yield points are not as pronounced as those encountered in body-centered-cubic metals containing interstitials. The yield point phenomenon in these materials is usually enhanced by prestrain-ing at low temperatures and aging at intermediate temperatures. Undoubtedly the yield point phenomenon induced by strain aging substitutional alloys also results from locking of dislocations. But the locking of dislocations in substitutional alloys of face-centered-cubic metals differs somewhat from interstitial locking of dislocations in body-centered-cubic metals. Substitutional elements in face-centered-cubic lattices cause only radial displacements of the adjacent lattice points. Consequently, only the edge components of dislocations can be locked by the mechanism suggested by Cottrell. Additional locking, however, can be obtained by the Suzuki mcchanism.10 In face-centered-cubic metals, dislocations exhibit lower energies when they are present in the

Jan 1, 1957

By H. I. Kaplan, C. Laird

Pulsating contpresslon experiments have been carried out on coppev single crystals in order to test the adequacy of mechanisms which have been suggested for stage I cvack grouth when tension-compression loading is the mode of- applying stress. Although dejorrrzation continues slowly for many millions 01 cycles under Pulsating compression stresses both within slip bands formed in the first cycle and by the creation of new bands, the morphology of these bands is not typical of fatigue slip bands in general and neither initiation of cracks nor propagation of artificially induced cracks was observed. It is therefore concluded that 1) initiation and growth of cvacks in fatigue is a consequence of repeatedly reviersed, microplastic deformation rather than intermittent, unidirectional deforma-tion, 2) gvowth of- cracks in polycrystalline specimens FAILURE of a plain specimen subjected to cyclic tension-compression stresses is known to take place by two stages of cracking.1,2 The first stage is characterized by propagation along slip planes roughly at 45 deg to the stress axis, in the outermost layers of the specimen, and at a rate of the order of angstroms per cycle. The average macro orientation of a stage I crack is normally at right angles to the stress axis,3 however, except in a few cases, e.g., single crystals in suitable orientation, when both the micro and macro directions of propagation are at 45 deg to the axis. After a stage I crack has penetrated to such a depth that the stress on the remaining section of metal is significantly increased, stage II crack propagation succeeds stage I, at right angles to the stress axis both micro and macroscopically. Direct evidence of the processes occurring during a subjected to pulsating compression stresses is caused by the setting-up of local tensile stresses due to interactions between grains, and 3) both unslipping mechanisms and operation of the "plastic blunting process " within a single slip band are adequate descriptions 01- stage I growth, being fundamentally similar. The slip mode of a material, however, is considered to control whether unslipping or the plastic blunting process takes place. Thus materials of planar slip mode fauor unslipping while those hacing a wavy mode, and capable of cross slip, can fail in stage I by the plastic blunting process. Similarly, and perhaps more importantly, reversed torsion and tension-compression modes of applying stress may be expected to cause the unslipping mechanism and the plastic blunting process, respectively. single stress cycle4,5 shows that stage II crack propagation takes place by a repetitive plastic rounding and closing of the crack tip. This mechanism has come to be called the "plastic blunting process".6 On the other hand, the extreme smallness of the phenomena occurring in stage I growth prevents elucidation of its mechanism by direct observation. Many suggestions have been proposed for this mechanism,&apos;"-" most of them involving some kind of random or systematic unslipping process. It is possible, however, that the mechanism of stage I growth is essentially the same as that of stage 11, but occurring on a much smaller scale.6 In this situation, it has been observed that cracks propagate along slip bands comprising dislocation cell structures and bounded on either side by regions of different dislocation density and distribution." Since there is evidence that such bands are softer than their matrix, * it may be expected that the plastic blunting process occurs within a single slip band during stage I growth, rather than on two gross bands as in stage 11, when the acting stress is greater.= The blunting of the crack tip allowed by coordinated

Jan 1, 1968

By J. A. Elias, A. J. Heckler

A simplified method is described for the prepauatior~ of comnposite samples for determination of one quadrant of an X-ray pole figure of a sheet material by the Schulz reflection method. Other special conzposite sampling methods were used to investigate specific areas of the pole figures. Plotting nets for these sampling plans were drawn by a computer. Step-scanning techniques were used to obtain detailed information about specific regions of a pole figure. Satisfactory pole figures were obtained for large-grained (0.5 to 1 in. diam) sheet materials. COMPOSITE sampling techniques for preparing one quadrant of an X-ray pole figure have been reported by a number of authors.&apos;-5 All of these methods are similar in that they use a composite sample prepared by stacking and holding together coupons cut from a sheet material. The schulze reflection technique is applicable to the composite sample so that only minor corrections are needed for defocusing and absorption effects. Also, a relatively large area of a sheet is sampled through the use of multiple cross sectional pieces. In general, sampling methods necessary to obtain a composite-type specimen have been formidable enough to discourage their use in routine pole figure determinations. This paper presents a versatile sampling technique which is simpler than those previously described along with an accurate method for obtaining stereographic plotting nets for any orientation of a composite specimen surface. A step-scanning technique is used to obtain more quantitative determinations of specific regions in pole figures, and the composite technique is satisfactorily applied to determine pole figures of very large grained sheet materials. REVIEW OF PREVIOUS TECHNIQUES Mueller and ~nott&apos; first proposed a technique for obtaining one complete quadrant of an X-ray pole figure by using a sheet sample and six different composite samples in which the individual coupons were aligned and held together by two types of jigs. Data from the seven different samples were overlapped to obtain a single quadrant. ~eieran~ and Lopata and ~ula~ obtained a complete quadrant of a pole figure from a single sample prepared by forming a cube from coupons cut from a sheet. These coupons were bonded together using an epoxy cement, and a cut was made at 54.7 deg to the normal, transverse, and rolling directions of the sheet, Fig. I. The normal direction of the cut face projects into one quadrant of a stereographic projection at equal angles to the three principal orthogonal reference directions mentioned above. X-ray reflection data obtained along a spiral path in physical space could then be plotted on a stereographic net offset 54.7 deg from the normal, as shown in Fig. 2. A rotation about an axis representing the surface normal is designated as a, and a rotation about an axis lying in the plane of the surface is 4. From Fig. 2 it is apparent that a tilting angle of no more than 55 deg @ is necessary to obtain the complete quadrant. ~eber~ devised an improved sampling technique by cutting coupons 45 deg to the rolling direction of a sheet, stacking and bonding these with an epoxy cement at an angle of 54.7 deg to a reference plane, and finally grinding and lapping a flat specimen surface parallel to this plane. An important consideration of this technique and that of Lopata and ~ula~ is variation in the diffracting power of the specimen surface due to the bonding medium. While any composite-type sample will have some such variation due to boundaries between coupons, the problem can be minimized by elimination of the bonding medium and by avoiding prolonged etching times which could lead to significant inter lame llar attack of the coupon faces. EXPERIMENTAL TECHNIQUE We were working concurrently though independently on a technique similar to ~eber&apos;s,&apos; but which does not

Jan 1, 1968

By R. E. Goodman, M. A. Mahtab

The, state of stress around a vertical wellbore in rock following nonlinear stress-strain laws is examined by means of finite element analysis. The wellbore, is considered an axisymmetrier body with axisymmetric loading. The initial vertcal and borizontal stresses are "locked" in the rock elements around the wellbore and a new state of stress is generated by the displacements which occur around the borehole. A point-uise variation of the elastic moduli is made on the basis of the new stress state and the triaxial data. The initial stresses tire now reintroduced along with tile changed moduli and original boundary constraints. This procedure is repeatrd until convergent stresses are reached. The effect of nonlinearity on streases is examined for a 6,000-ft wellbore in a scbistose gnezss and Berea sandstone using results of laborratory triaxial compression tests. The results sbou that the effect is restrited to one well radius from the bottom periphery of the bole. Beyond a distance of one-quarter radius, the elfect of nonlinearity on stresscs IS almost always less than percent for the cases conszdered, The consideration of a static pressure inside the well does not magnify the effect of monlinearity on borehole stresscs. INTRODUCTION The terms " wellbore" and "borehole" here designate cylindrical openings in the ground with vertical axis and a circular cross-section. A knowledge of the stress redistribution that occurs on excavating a wellbore is important in understanding (1) the behavior of the lined or unlined hole, (2) hydraulic fracture response, and (3) the effect of stress redistribution on drillability; also it is important in predicting initial stresses in the virgin ground, and in analyzing the response of measuring instruments placed in the borehole. Our knowledge of the state of stress around a wellbore has been restricted to homogeneous, isotropic, elastic material and derives chiefly from the analysis by Miles and Toppingl and the photoelastic work of Galle and Wilhoit2 and Word and Wilhoit.3 In this investigation the state of stress is examined for a nonlinear elastic material by means of finite element analysis. Many rocks possess stress-strain curves that depart notably from straight lines in their initial or final portions. While the literature contains abundant stress-strain data from triaxial tests (axisymmetric loading) on cylindrical rock specimens, there is little information on rock deformability under nonaxisymmetric loading conditions such as occur at each point around the bottom of a wellbore. Although there is some knowledge of the effect of intermediate principal stress on rock strength, there is virtually nothing known about its effect on rock deformability; therefore, we have assumed here that the effect of intermediate principal stress can be ignored. A schistose gneiss4 and Berea sandstone5 were selected as representative rocks for this analysis. The traditional graphs of deviator stress (ol -o3) vs axial strain were reworked to give the tangent modulus as a function of the deviator stress for varying values of the minor principal stress. The result is a nesting family of skewed, bell-shaped curves for tile gneiss (Fig. 1A) and the sandstone (Fig. 2A). A similar replotting of the lateral strain data defines the variation of Poisson&apos;s ratio (v) with the deviator stress and confining pressure. These curves, shown in Fig. 1B for the gneiss and in Fig. 2B for the sandstone, are not so well ordered as the tangent modulus curves. However, all of these display an increase of v with deviator stress application, but the rate of increase diminishes with confinement. The ET and v curves for the two rock types are tabulated in Tables 1A and 1A for use in a digital computer so that material properties corresponding to a given state of stress can be assigned by interpolation. STATEMENT OF THE PROBLEM Given a state of stress in the ground prior to drilling a hole (initial stresses), we assume that

Jan 1, 1969