Search Documents

By D. E. Baldwin

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

Jan 1, 1967

By Klaus Schwerdtfeger

The rate of formation of silicon monoxide from silica spheres has been measured at 1500°C in H2, H2-H2O, H2-He, H2-Ar, and CO-CO2 gas mixtures. Indications are that in the hydrogen-containing gases the reaction rate is mainly determined by convective diffz~sion. In the CO-CO2 mixtures the rate approches a constant value which is independent of sphere diometer (and gas velocity). This suggests that a surface reaction becomes rate controlling in CO-CO2 gas. WHENEVER silica is exposed at elevated temperature to a strongly reducing atmosphere there occurs, inevitably, a volatilization of silica by the formation of gaseous silicon monoxide. Because of the importance of this reaction in metallurgical processes and laboratory experiments, its kinetics were studied in the present work. EXPERIMENTAL The spheres made of vitreous silica were suspended by a molybdenum wire in the hot zone of a gas-tight high-purity recrystallized alumina furnace tube of 3.1 cm internal diam and 100 cm length. Molybdenum wire resistance heating was used. A purified H2, H2-H2O, H2-He, H2-Ar, or CO-CO2 gas mixture was passed through the tube, from top to bottom. The amount of silicon monoxide formed was determined by weighing the silica balls before and after the experiments. The silica balls were made from commercially available translucent fused quartz glass. In experiments with hydrogen-containing gases the silica stayed vitreous during the experiment, whereas in experiments with CO-CO2 mixtures, which were performed over longer times, the balls were always crystallized to cristobalite (as determined by X-ray examination) at the surface. Constant-pressure head capillary flowmeters were used to control the composition of H2-He, H2-Ar, and CO-CO2 gas mixtures. The oxygen and water vapor were removed from the gases containing hydrogen by passage through columns containing platinized asbestos (450°C) and an-hydrone. The CO was purified from O2 and CO2 by passing it through activated copper (200°C) and through ascarite. The H2-H2O gas mixtures were prepared by passing purified Hz through a mixture of oxalic acid dihydrate and anhydrous oxalic acid held at a constant temperature in a water bath. The equilibrium water vapor pressures of oxalic acid dihydrate-anhydrous oxalic acid mixtures were taken from the available data.1 The furnace system was checked for sufficient gas tightness at regular intervals. As an additional check, the H2O and CO2 contents were measured in the exit gas and were found to agree closely with the metered contents when H2-H2O or CO-CO2 mixtures were used. When purified H2 was used, H2O pressures found in the exit gas were always below 0.1 mm Hg. These contents originate either from an incomplete reoxidation of silicon monoxide in the cold part of the furnace or from small leaks which cannot be avoided. If the latter is the case, such small uncontrolled H2O contents decrease the partial pressure of SiO in equilibrium with SiO2 by less than about 7 pct, which is considered small in comparison to other experimental uncertainties. The temperature of the furnace was controlled automatically in the usual manner. The reported temperature was measured with a Pl-Pt/l0pct Rh thermocouple and was accurate to ± 5°C. TRANSPORT-CONTROL LED RATE OF REDUCTION The gaseous convective diffusion from the surface of a sphere can be predicted from the semiempirical relationship obtained by Ranz and Marshall2 in a study of the evaporation of water and benzene drops in air. The following equation applies:

Jan 1, 1967

By S. H. McFarlane, K. H. Zaininger, G. W. Cullen, C. C. Wang, G. E. Gottlieb

Heteroepitaxy of silicon on stoichiometric spinel has been studied. Both boron-doped (p-type) and arsenic-doped (n-type) single-crystal silicon films have been grown by the pyrolysis of silane on sioichiometric spinel (MgAl2O4) substrates prepared by a flux technique. The spinel crystals are free from strain and subgrains, and exhibit no thermal instability at temperatures employed for silicon deposition. The growth parameters of epitaxial silicon have been studied. Effects of substrate orientation, surface preparation, and growth conditions on the silicon Properties have been determined. The silicon-spinel composite was characterized by X-ray diffraction, electron diffraction, electron microscopy, and optical techniques. Information on epitaxial orientation relationship, defect structures and film-substrate interface has been obtained. A parallel orientation relationship was found between silicon and spinel. The silicon films exhibit no subgrain structures. However, growth pyramids, microtwins, stacking faults, and dislocations are the prevalent defects. The deformation of epitaxial silicon was measured by an X-ray technique. It was found that the epitaxial silicon is strained under compres-sive stress due to the thermal expansion difference between the film and substrate. The stress was estimated using the classical plate theory and beam equation. Electrical properties of the epitaxial silicon have been measured, including resistivity, carrier concentration, and Hall mobility. MIS and vertical through diffused junction structures have been fabricated wing silicon-spinel composite. Information on MIS characteristics, minority carrier life time, and junction properties has been obtained. THE achievement of single-crystal growth of large area silicon semiconductor films on oxide insulating substrates is of technical importance to many solid state electronic devices. The semiconductor-insulator composite structure is also of scientific interest because the heteroepitaxy is determined by the spatial relationship between the atomic arrangement in the substrate and that of the atoms in the appropriate crystallographic plane of silicon. The degree of crystalline perfection of the silicon deposits depends also on the physical condition of substrate surface. Epitaxial growth of silicon on commercial flame fusion Al-rich spinel (MgO:3.3A12O3) has been reported by several investigators.'-4 The lattice of spinel of this composition is distorted. The material cracks and exsolves the excess alumina in the lattice at temperatures (-1200°C) used in fabricating silicon devices. Stoichiometric MgA12O4 spinel exhibits neither lattice distortion nor thermal instability. Results are given on the epitaxial growth and characterization of silicon on stoichiometric spinel. EXPERIMENTAL Fabrication of Silicon-Spinel Composite Structures. The MgA12O4 spinel substrate wafers were prepared from bulk single crystals grown at about 1200°C by a flux technique5 using a PbF2 flux. The crystals are free from strain and subgrains, and exhibit a typical dislocation density of less than 50 lines per sq cm. The crystals are thermally stable without cracking and exsolution of alumina to at least 1500°C. The spinel crystals were oriented for cutting by the Laue X-ray back-reflection method.6 Substrate wafers (about 25 mil thick) of (Ill), (loo), and (110) orientations have been obtained bu cutting the X-ray oriented crystals using a standard type diamond wheel. An accuracy of better than 1/2 deg is maintained throughout the operation. The spinel substrate wafers were mechanically lapped and polished to produce the flat and smooth surface which is required for silicon epitaxy. Lapping was carried out with about 30-µ boron carbide abrasives to obtain a flat coplanar surface. The lapped surface was further polished using successively finer grades of alumina, ending with the 0.3-µ grade. After polishing the wafers generally have a flatness of ±0.4 µ per cm, as revealed by interferometry. Fig. 1 shows (111) oriented MgA12O4 spinel wafers of various shapes.

Jan 1, 1970

By J. A. Stachura

Continuous mining machinery provides the coal industry with one way to compete for a larger share of the total energy market. Various types of machines are discussed and some of the problems with continuous miners, encountered by operators, are reviewed. Equipment manufacturers are working with mine personnel to provide solutions for problems that arise. While coal production over the past 25 or 30 years has been on a horizontal plane, coal's share of the total energy market has declined. To participate more effectively in this total energy market, it is necessary to produce coal more efficiently. It is the obligation of all management, employes, and mining departments to gear the deep mining industry to the rapid progress and changing of today's modern industry. This can be accomplished in the near future with the selection of the proper type of continuous miner best suited to each operator's individual situation. In most mining operations there is tremendous incentive to undertake the continuous mining program. It can reduce the size of the mine greatly by permitting a minimum of working places; it makes pillar recovery work more efficient from the standpoint of overall cost, amount of coal recovered, and safety. The work force can be reduced materially permitting closer and more efficient supervision. It simplifies maintenance because equipment can be more readily standardized. The trend of the coal market favors the use of continuous mining machines. Although there appears to be a general feeling that continuous mining is still a relatively new program and will be slow in replacing conventional mechanical equipment, the fact is that tremendous strides have been made since the first machines were installed in 1948. This program is advancing at approximately the same rate that mobile loading machines replaced hand loading. From 1948 to 1955 there were approximately 450 continuous mining machines in service. In October 1959, a survey revealed that there were more than 700 continuous mining machines in service. Many operators have expressed a desire to undertake this program, but they feel that they could not do so at this time because of one or more of the following reasons: 1) the thickness of their coal seams, 2) seam characteristics, 3) soft bottoms, 4) bad roof conditions, 5) size consist, 6) insufficient flexibility in machines, 7) difficult ventilation problems, and 8) high maintenance costs. With the realization on coal about the same today as it was in 1948, or slightly less and since coal is still failing to participate to a greater degree in the total energy market, it is not surprising that the coal industry is desperately exploring more economical methods for deep mining. The manufacturers are aware that the coal industry is willing to invest in continuous miners if the equipment is built for maximum flexibility, will produce higher tons per man, and assure long life between overhaul programs. CONTINUOUS MINING MACHINES Before discussing details regarding the selection of a continuous miner, let us have a preview of some of the continuous mining machines which are available to the coal industry today. Jeffrey Manufacturing Co.: The machine shown is the Jeffrey 76 A.M. Colmol. This is their most widely used miner, and has been particularly successful in central Pennsylvania and in high-wall mining in western Kentucky. One of the outstanding features of this auger-type miner is its portability. The entire mining range can be changed from its lowest point to the maximum height without stopping the mining operation. Jeffrey 76 B.M. Colmol: This machine is similar to the 76 A.M. model; however, it is built bigger and stronger for a mining range of 50 1/2 to 72 in. This is the model that is now available (Fig. 2). Jeffrey has added two arms to the top row, omitted the odd arm in the bottom row, thus permitting a 50 pct larger throat opening. This eliminates one

Jan 1, 1961

By Robert W. Lawson

INTRODUCTION The two basic types of mine-air compressors are the positive-displacement compressors and the dynamic compressors. Positive-displacement compressors con¬fine successive volumes of air within a closed space and increase the pressure by reducing the volume of that space. Reciprocating and rotary compressors operate on this principle. Dynamic compressors utilize an impeller to acceler¬ate the air, with the increased air velocity converted to pressure through a stationary diffuser. Centrifugal compressors and turbocompressors operate on this principle. Reciprocating Compressors Reciprocating compressors include both single-stage and multistage units. For mine service the most com¬mon is the two-stage intercooled compressor having an atmospheric pressure intake and a discharge pressure of 689 to 861 kPa (100 to 125 psi). The cylinders of a reciprocating compressor may be either single acting or double acting. In the single-acting cylinder, compression occurs on only one side of the piston, while, in double¬acting cylinders, compression occurs on both sides of the piston. Fig. 1 shows a cross-sectional view of a typical two-stage double-acting water-cooled reciprocat¬ing compressor. Fig. 2 illustrates the sequence of opera¬tions occurring in the cylinder of a reciprocating com¬pressor. Most modern mine compressors are double acting and are driven by electric motors. Induction motors are used on low-powered compressors, and synchronous motors are used on larger compressors. Steam-driven compressors once were common and still are applied where a reliable source of low-cost steam is available. For this type of compressor, steam is a desirable power source since a governor that senses the discharge pres¬sure can be used to control the steam inlet. This can be used to regulate the compressor speed to compress only air sufficient to meet the mine's requirements (up to the capacity of the unit). Unlike a constant-speed electric¬driven compressor, no other regulation is required for the cylinders of a steam-driven unit. Rotary Compressors Rotary compressors include sliding-vane, liquid¬piston, two-impeller straight-lobe or cycloidal, and helical-lobe units. Each of these has specific advantages and applications. Of these, the helical-lobe or "screw" compressor has the greatest application for supplying "100-lb air" [689 kPa (100 psi) ] to small- and medium¬size underground mines. Fig. 3 illustrates the operation of a sliding-vane compressor. Fig. 4 illustrates the operation of a helical lobe or screw compressor. Fig. 5 illustrates the operation of a two-impeller straight-lobe compressor. Fig. 6 illustrates the operation of a liquid¬piston compressor. Centrifugal Compressors Centrifugal compressors deliver 0.94 m3/s (2000 cfm) or more of air, and they are best suited to supply¬ing large volumes of compressed air for base loads. The efficiency of a centrifugal compressor declines rapidly when its operation deviates from its design point. If the power cost is not a critical consideration, a centrif¬ugal compressor is attractive as a result of the low initial and installation costs and the high degree of reliability. Fig. 7 illustrates a typical four-stage centrif¬ugal compressor, and Fig. 8 is a cross section of this type of unit. COMPRESSOR SELECTION The process of selecting a compressor requires sev¬eral considerations, including the capacity, the control system, the location, the power supply, and the cooling provisions. Compressor Capacity To determine the required capacity, the total air quantity and pressure requirements must be calculated. The air consumption of various tools, as listed in Table 1, must be included in this calculation. Some of the other considerations include: 1) The correction factor for altitude, as listed in Table 2, must be included in the calculation. While pneumatic tools require additional compressed air as the altitude increases, the compressor produces less. 2) A correction factor must be incorporated for the use of multiple tools. Table 3 lists appropriate correc¬tion factors for rock drills.

Jan 1, 1982

By J. H. O’Dette

BLISTER formation on wrought aluminum products has been a matter of concern to aluminum fabricators for many years. The fact that blisters make their appearance during the final stages of fabrication means that much time and effort have been expended before it is evident that the material must be scrapped. This and the loss of production are ample reasons for embarking on a research program to find out more about blistering and possible methods of control. Many excellent papers on this subject have appeared in the technical literature in the last few years. In one of the more recent papers, Kostronl gives a fairly complete review of the problem of gas in aluminum and blistering in general. From the work of the various investigators, a general mechanism for the formation of blisters has evolved and may be summarized briefly as follows: Hydrogen gas is usually present in the cast ingot, in the atomic or ionic state, in excess of the solid solubility. At room temperature the gas atoms are relatively immobile, but at elevated temperatures they become mobile and diffuse to free surfaces where they associate to form molecules. If these free surfaces are discontinuities within the ingot then the associated gas molecules collect and build up considerable pressure until equilibrium is reached with the dissolved gas around the discontinuity. Since the solubility of the gas is proportional to the square root of the partial pressure it is evident that large pressures may be built up at the discontinuities. When the metal is rolled to sheet the pressure of the gas in these, discontinuities is sufficient to deform the surface of the sheet and form a blister, particularly when the strength of the material is reduced during the annealing. Although the mechanism of the formation of blisters was fairly well known, the actual impor- tance of the individual roles played by the gas and the discontinuities was not so well established that positive methods for control of blistering were at hand. It was the object of the investigation to learn more about the nature of the discontinuities giving rise to blistering and to develop methods of control. The investigation was carried out in three stages as follows: a) A study of the more commonly used methods of removing hydrogen from aluminum, such as fluxing and holding, and their effects on blistering. b) Location of the discontinuities giving rise to blistering. c) A study, by indirect methods, of the movement of gas in the metal and its effect on blistering, and the development of methods of blister control during fabrication. Removing Gas and Reducing Blistering by Fluxing or Holding-—The work of removing gas and reducing blistering by fluxing or holding was carried out on metal of Alcan 2s (1100) composition with an artificially high gas content. The high gas content was obtained by introducing water into the melt prior to casting; the water reacted with the aluminum forming A12O3 and hydrogen. The procedure was to raise the gas content by introducing water into an aluminum melt and using half the melt to cast one ingot immediately. The remaining metal was treated again with water to make up for the reduction in gas content during holding, while the first ingot was being cast, and then given a fluxing or holding treatment and cast into an ingot. Samples for gas determination were cast from the molten metal before and after the gassing operation; before, after, and during fluxing, and while casting. Samples were also cut from sections of the ingots. Gas determinations were made on these samples using the well known vacuum fusion technique. All ingots were hot rolled at 500°C to 0.250-in. slab which, in turn, was cold rolled to 0.040 in. sheet. The sheet was coiled and annealed for 2 hr at 500 °C and then cut into 5-in. squares which were inspected for

Jan 1, 1958

By G. R. Speich, S. A. Kulin

The isothermal formation of martensite at subzero temperatures has been studied in an austenitic stainless steel. The amount of martensite formed isothermally in a given time was found to follow a C-curve behavior with decreasing temperature. Since isothermal rnartensitic transformation occurs at temperatures above the true M. temperature of the olloy, the observed M. temperature is dependent on the cooling rate. It is believed that a general martensite theory must include the basic concepts of the strain embryo theory, and also the important role of thermal fluctuations. UNTIL recently the isothermal component of the martensitic transformation had been considered as a trivial effect in theoretical considerations of this type of transformation. One of the postulates used in defining a martensitic transformation was that the reaction was temperature dependent but not time dependent. Cohen and his associates1'2 were the first to report the isothermal formation of mar-tensite. They obsel-ved the time dependency of the martensite reaction at room temperature in plain carbon steels and certain tool steels. They found that a small amount of additional transformation occurred after the specimen reached the quenching temperature and that the rate of this isothermal transformation decreased rapidly with time. These authors considered this phenomenon to be a "tailing-off" of the main athermal reaction. Kurdjumov and MaksimovaL 5 studied the isothermal formation of martensite at subatmospheric temperatures in steels containing 0.6 pct C-6 pct Mn, 1.6 pct C, and in an iron-base alloy containing 23 pct Ni and 3.4 pct Mn. In each case it was found that both the initial rate and the total amount of isothermal transformation exhibited sharp maxima with decreasing temperature. They stated that it was possible to completely suppress the martensitic transformation in the 0.6 pct C-6 pct Mn steel by drastic quenching in liquid nitrogen. They interpreted their results in terms of a thermally activated process. The maximum in the rate of nucleation vs. temperature curve was assumed to be a result of the temperature dependence of the rate of growth of the martensite plates. It has been demonstrated," however, that martensite plates can form with great rapidity even at temperatures approaching absolute zero. In addition," attempts to confirm the reported suppression of the martensitic transformation by rapid cooling of a 0.6 pct C-6 pct Mn steel have proved unsuccessful. On this basis, and on other theoretical grounds, the validity of some of Kurd-jumov's results and conclusions is open to serious question. Cohen, Machlin, and Paranjpe' have proposed a theory of martensite formation based on the nucleation of strain embryos and a growth mechanism involving cooperative atomic displacements. According to this theory thermal fluctuations can make a strain embryo supercritical and thereby produce martensite isothermally. Das Gupta and Lement' have investigated the isothermal formation of martensite in a 17 pct Cr-0.7 pct C steel at temperatures below the ambient. They found that partial suppression of the martensitic reaction occurred on quenching into liquid nitrogen. That part of the reaction which can be suppressed is taken to correspond to the isothermal rather than the athermal part of the rnartensitic transformation. The initial rate of isothermal transformation for this steel was found to go through a maximum with decreasing temperature. This is in conformity with the findings of Kurdjumov and Maksimova for the steels used in their investigation. Das Gupta and Lement conclude that isothermal martensite formation occurs by the nucleation of new plates and is always preceded by some athermal transformation. An iron-base alloy containing 14 pct Cr-9 pct Ni was used in this investigation to study further the kinetics of the austenite to martensite transforma-

Jan 1, 1953

By G. R. Wilms

A study has been made of the structural changes in polycrystalline high purity aluminum during the tertiary stage of creep under conditions of constant tensile load. It appears that there is no basic modification in the particular mechanism of deformation that characterized the preceding secondary stage of creep, but that the changes in crystalline structure produced during the tertiary stage are only a consequence and not a cause of the accelerated strain. Intercrystalline fissures were observed on entry into the tertiary stage, and these might supplement the accelerated rate of strain owing to local stress concentrations. There is evidence that the formation of these fissures may be associated with the presence of surface imperfections, and it is suggested that a more generalized treatment of the results may explain certain creep characteristics of commercial alloys. IN previous work'-4 a study was made of the changes produced when specimens of annealed aluminum were deformed at various temperatures and rates of strain, especially under conditions leading to creep under constant tensile load. It was shown that, as the temperature was increased or as the rate of strain was decreased, elements of the substructure produced within the grains became coarser, and with the strain rate constant at a given temperature, their size tended to an equilibrium value. The preliminary process of the breaking down of the grains into a substructure was identified with primary creep.' The subsequent equilibrium, condition was associated with secondary creep. The question then arose whether the tertiary stage corresponded to any further significant change in the substructure. The question is also of interest in view of the opinion, expressed for instance by Andrade,5 that the tertiary stage is merely an accelerated rate of strain due to the reduced area of section of a specimen under continuous tensile load, and the other opinion, as suggested by the observation by Sully, Cale, and Willoughby6 of a tertiary stage in com-pressive loading, that this third stage is a real effect in its own right. The object of this paper is to record some observations on the tertiary stage of high purity aluminum. Experimental Procedure The material was aluminum of 99.98 pct purity as previously used. The specimens had a rectangular section of l/8 x l/2 in., and extensions were measured over a 1 in. gage length. After machining to shape, the specimens were annealed to give a standard grain size of 0.2 mm and then electropolished. The general procedure involved metallographic and X-ray diffraction examination at progressive stages of elongation (marked by arrows on the curves shown in Figs. 1 and 2), the specimen being extended to fracture under constant load at a given temperature. The X-ray tests were made by the back-reflection method from a stationary specimen, using CoK radiation. The longer edge of the micrographs in the paper corresponds to the direction of applied stress. Results The previous work showed that, in the secondary stage of creep in aluminum, three basic types of structure could be developed according to the temperature and rate of strain. These were: 1—slip, which predominated at the lower temperatures and higher rates of strain and produced the finest and least perfect substructure; 2—cell structure, which superseded slip as the temperature was raised and the rate of strain decreased and which was characterized by a coarser and more perfect structure and an absence of slip lines; and 3—a stage termed "boundary microflow," which predominated at still higher temperatures and lower rates of strain when the cell size had increased and become comparable

Jan 1, 1955

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 W. C. Knopf, G. Samsel

For several years International Minerals & Chem-ical Corp. has used a radiation method to assay potassium content of products from potash and feldspar beneficiation. The procedure is rapid, accurate, and inexpensive and is easily handled, if necessary, by nontechnical personnel. Data are presented to show the accuracies that may be expected if the technique is used for control of pilot plants or large-scale plants. Natural potassium is composed of three isotopes —K, K and K. Five other isotopes, artificially produced, have been detected. In natural potassium only the K isotope, with a half life of 1.4x10 years, is radioactive. Its modes of decay are beta and gamma rays, positrons, and x-rays. Only the more abundant beta rays, however, are important in radiation measurements using geiger tubes. Barnes and Salley' measured radiation of potassium salts in solution in 1943. Gaudin and Parnell2 first reported the use of solid salts arranged in a vessel surrounding the ß counter in 1948. Later, Scheel used the method to assay potassium content of various crude and finished potash salts. Experimental: Dip-type Geiger-Muller tubes (TGC-5) supplied by Tracerlab Co. were used in all radiation counting measurements. The tube was encompassed by a glass container about M in. larger in diameter than the tube. This container, which served as a sample holder, was closed at the bottom with a rubber stopper to facilitate removing the samples. Samples were not weighed, but about 40 g of sample were needed for each test. Electrical impulses from the tube were counted on either a Tracerlab scaler or a rate meter. To obviate the need for constant attendance, preset features were incorporated into the scalers, so that at a predetermined time or count the units were automatically deactivated. When rate meters were used, counting rates were read directly from instrument or from a recorder used in conjunction with the rate meter. Each tube was calibrated, according to the manufacturer's instructions, to determine the correct operating voltage to insure reproducible results. A linear calibration plot of percent K2O vs counting rate was also obtained for each tube by recording the counting rate plus background count of several standards of different K2O content chemically analyzed by the chloroplatinate method.' These standards were of approximately the same particle size distribution as samples to be assayed. Specific background count caused by cosmic rays and other surrounding radiation was checked periodically by counting a sample of pure NaCl. Fig. 1 is a representative calibration plot. Assays were determined directly from the standard plot or calculated by subtracting the background count from the counting rate, then multiplying by the slope of the line. Chemical assays of the test samples were made by the Sherril method," a flame photometric technique, or the chloroplatinate method. Data: The radiation assay method has been used for three years by International Minerals & Chemical Corp. in the laboratory and in pilot plants processing ores containing potassium. These dry concentration plants process sylvinite and langbeinite ores at Carlsbad, N. M., and feldspar ores at Top-sham, Maine. Sylvinite ore is predominantly a mixture of syl-vite (impure KCl) and halite (impure NaCl). Contaminants include sulfates, such as langbeinite (2MgSO - K2SO4), kieserite (MgSO H2O), poly-halite (MgSO4.K2SO, 2CaSO, . 2H2O), and slimes. Langbeinite ore is predominantly langbeinite and halite with small amounts of sylvite, sulfates other than langbeinite, and slimes. Several variables, in addition to correct calibration and standardization of tubes, were found to be important in obtaining reproducible results. Particle size distribution of standard samples should be

Jan 1, 1959

By Henry A. Wriedt, John Chipman

Equilibrium in the reaction of hydrogen gas with oxygen in liquid nickel, iron, and their alloys has been studied at temperatures of 1500° to 1700°C. The equilibrium con^stant, 0/p, [% O], is greater in nickel then in iron by a factor of 45 at 1594 C. In the alloy steel range, the activity coefficient of oxygen is slightly increased by the presence of nickel. NICKEL is more resistant to oxidation than is iron, in both the solid and the liquid state. The solubility of oxygen,' however, is greater in nickel than in iron. The two metals and their liquid alloys constitute an interesting series of solvents in which the chemical behavior of oxygen can be studied by methods similar to those used for liquid iron.' Experimental Method The furnace assembly, including the tube for preheating the entrant gas, was the same as that described by Gokcen and Chipman. " controlled mixture of water vapor, hydrogen, and argon preheated to approximately the melt temperature impinged upon the surface of the liquid metal. The charge, weighing about 35 g, was held in an alumina crucible at constant temperature in the gas stream for a prolonged period, then quenched in a stream of cold helium, sampled, and analyzed for oxygen by the vacuum fusion method. The temperature of the liquid metal was measured by a Leeds and Northrup optical pyrometer which had been checked against a similar instrument certified by the National Bureau of Standards. The freezing point of electrolytic iron in hydrogen, taken as 153j°C, was used for frequent checks. This point and the emissivity determined by Dastur and Gokcen" established the transmissivity of the optical system. The emissivities of Fe-Ni alloys at 1535" are known from the experiments of Smith and Chipman,h nd it was assumed that the variation with temperature is parallel with that of iron. Recorded temperatures are probably accurate within k5"C. Preparation of Gas Mixtures—Since the ratio H,O:H, had to be much higher for nickel than for iron, it was found convenient to construct two separate systems for gas preparation. In both cases, the ratio of argon to the sum of the other two gases was at least 4, the proportion of hydrogen and water vapor being varied to produce H,O:H, ratios from 0.2 to 70. The primary purpose of the argon was to diminish the errors caused by thermal separation' in the gas phase. In addition, its presence lowered the saturation temperature required for a given H,O:H, ratio and minimized the difficulty of preventing condensation. Both systems were modifications of the one used by Dastur and Chipman.' For H,O:HI ratios greater than 2, welding-grade argon purified by passage over magnesium chips at 620°C was mixed with hydrogen produced in the electrolytic cell shown in Fig. 1. The design was based on a cell used by Bodenstein and Pohl," the important features of which area nickel wire cathode and sodium hydroxide solution as the electrolyte. The rate of hydrogen evolution was controlled by a calibrated ammeter and a rheostat in the circuit, the current efficiency being taken as 100 pct. Precise control of the argon flow rate at about 300 ml per min was achieved by means of a capillary flowmeter having a controlled pressure drop. The molar flow rate was measured volumetrically before each run, and a constant flowmeter setting was maintained thereafter. The A-H, mixture passed through a water saturator of the type previously described," the saturation temperature being maintained automatically with a fluctuation of k0.04"C. The tube between saturator and furnace was heated by a resistance winding. The system was of all Pyrex construction with the exception of the stainless steel tube containing magnesium chips. This was connected with the glass tubing at the upstream end by

Jan 1, 1957

By C. Wachmann

The theory presented ill this paper makes possible examination of the phenomenological aspect of secondary oil recovery by alcohol flooding. Limitutions of the theory are contingent on three primary assumptions: (1) chemical equilibrium, (2) no fingering, and (5) no dispersion. The second assumption implies one-dimensional flow. In a geometric sense the theory extends the Buckley-Leverett construction to systems of two slightly miscible fluids and to systems of three components. This is made possible by considering single-phase fluids to be made up two phases, namely of an excess of one phase and a negative amount of the other phase. The procedure is discussed in the paper. The theory also permits direct determination of the relationship between the aqueous phase saturation and the concentration of water in the aqueous phase in a zone of variable composition. The velocity of planes of constant water concentration and of constant saturation can be computted likewise and the existence of saturation discontinuities inferred. Consequently the behavior of such a Zone, which may exist at some time between a displaced and a displacing fluid, May be determined. mathematical details are given in Appendices. Problems of experimental verification are pointed out on the basis of a supporting experiment. Implications of the theory and effects of violations of primary assumptions are discussed briefly. The major conclusion is that, if the primary assumptions are met, with pure alcohol and alcohol-rich .slugs piston-like displacement of oil and water is possible. INTRODUCTION The recovery of oil from porous media by alcohol displacement has been studied for several years. Extensive experimental investigations into alcohol flooding have been carried out in the laboratory although little has been done in the field. A large amount of practical information is available, the extent being indicated by the first three papers cited 1-3 and their respecrive bibliographies. This has permitted the Formulation of inferential theories. A more analytical approach in which theories are developed primarily by reasoning from known physical principles and relationships has, however, lagged behind. In fact, an analytically derived theory of alcohol displacement has not been available. The theory presented in this paper on oil and water displacement by alcohol from a porous medium is a result of the analytical approach. The theory is idealized, as certain assumptions had to be made to make the mathematics tractable and reduce complexity of argument. In particular, major assuptions concerning conditions of equilibrium, interfacial stability, longitudinal mixing, etc., led to the exclusion of fingering, for example, from all considerations. In many laboratory experiments reported in the literature fingering could not be avoided. There is, therefore, disagreement between theories based on experimental evidence available and the analytically derived theory presented here. This theory presented is based on known physical data and makes use of the Buckley-Leverett construction and its suggested extension. The mathematical derivation of the theory is given in Appendices. Three examples relevant to alcohol displacement are discussed to explain the use of this theory: (1) system of two immiscible fluids (the Buckley-Leverett case); (2) system of two slightly miscible fluids involving an extension of the Buckley-Leverett construction; and (3) system of three components, including injection of pure alcohol or alcohol-rich fluid as a special case, which involves a more general extension of the Buckley-Leverett construction. In these discussions strict assumptions were made to impose equilibrium conditions and exclude mixing and fingering. Effects of relaxing these assumptions and consequent extensions of the theory are indicated. The complexity of the alcohol displacement process makes experimental verification of the theory essential. Some relevant evidence is offered in this paper.

Jan 1, 1965

By P. A. Beck, Hsun Hu

It has been known for many years that in cold drawn polycrystalline aluminum the recrystallization texture is practically identical with the deformation texture.l,2,3 V. Goeler and Sachs4 stated that in cold rolled polycrystalline aluminum, too, the deformation texture is retained upon recrystallization. This behavior was interpreted by Dehlinger9 as an indication that each recrystallized grain inherits the orientation of the matrix in which it grew. However, Burgers and Basart showed6 that in deformed aluminum single crystals the orientation of the recrystallized grains strongly deviates from that of the matrix. Barrett demonstrated8 that even in deformed polycrystalline aggregates the local orientation within each deformed grain changes on recrystallization. No explanation has as yet been offered of how the aggregate texture can be retained, while the local texture is everywhere changing. According to Burgers and Louwerse6 the orientation change on recrystallization in compressed and rolled aluminum single crystals corresponds to a rotation around [112] axes lying in the active slip planes and perpendicular to active slip directions. This orientation change was rationalized by them on the assumption that the orientation of the recrystallized grains is determined by the orientation of the available nuclei, and that these nuclei are small fragments of the deformed crystal, rotated in the deformation process. Barrett pointed out,? however, that these assumptions were by no means proven, and that the total volume of material necessary for the nuclei is so small, that it must be available in any orien- tation whatever. It is, therefore, unlikely that an "oriented nucleation" theory could satisfactorily explain the orientation changes in recrystallization of a deformed crystal. Furthermore, in contradiction to Burgers and LOU-werse's conclusions, Barrett's results for the orientation change could be best expressed in terms of rotations around. [Ill] axes. Thus agreement is lacking not only in the interpretations but apparently even in the experimental results contributed by the various investigators. The orientation relationships on coarsening* in a highly oriented polycrystalline matrix were first investigated by Burgers and Basart.5 They produced a fine grained highly oriented matrix by the recrystallization of a compressed or rolled single crystal of aluminum, and studied the orientation relationships in three successive generations of crystals: the "deformation texture" in the single crystal, the recrystallization texture, and the coarsening texture. Burgers24 showed that the preferred orientation of each successive generation was almost always widely different from that of the preceding one, but that the coarse grains often approximately reverted to the orientation of their grandparent, the deformed single crystal. This observation, and the idea that the orientation of the coarse grains is determined by that of the available "nuclei," led to the concept that the coarse grains are nucleated by the remnants of the deformation texture.12,15 On the other hand, Bowles and Boas16 found, in a carefully conducted experiment, that in recrystallized fine grained copper of cube texture coarse grains grow in orientations which can be derived from the orientation of the recrystallized matrix by rotation around a [lll].] direction: the coarsening texture here does not appear to be related in any simple manner to the original rolling texture from which the cube texture had been obtained by recrystallization. It may be assumed that in a recrystallized material, even with a strong texture, there always are present some grains of practically any orientation. These may serve as "nuclei" for the coarse grains, provided that the conditions are favorable for their growth. The existence of a coarsening texture may be considered as an indication that grains with certain relative orientations with respect to the preferred orientation of their neighbors can grow much faster than others. This "oriented growth" theory was briefly suggested by van Arkel, as early as 1936,13 and it was mentioned by Burgers." Recently, C. G. Dunn produced direct

Jan 1, 1950

By J. A. Marshall

"INTRODUCTION The purpose of this paper is to investigate the potential for use of UAVs in underground mines and present a prototype design for a novel autorotating UAV platform for underground 3D data collection. The mining industry has recently shown increased interest in the use of unmanned aerial vehicles (UAVs) to assist in everyday operations [20, 11, 4]. Above ground, small UAVs are in some cases a more efficient, inexpensive, and safer alternative to manned aircrafts currently used for photography, inspection and security [6, 10]. For example, attaching camera, infrared and LiDAR payloads, UAVs can provide a low-cost method of obtaining highly accurate 3D photogrammetric data and aerial photography. UAVs are now commonly used in open pit mining operations for applications that include stockpile surveying, 3D pit modelling, facilities management, accident reporting, progress monitoring, and environmental assessment [11]. There are numerous companies offering UAV services for mining applications; e.g., UAV Geomatics, Leica Geosystems, Sky Futures, SenseFly and FlyTerra, to name just a few. Common services include stockpile and open pit volume computations, environmental assessment, aerial mapping and photogrammetry. Although the aforementioned services are thosemost commonly provided by UAV surveying companies, they may also add value to other daily mining operations. For instance, SenseFly offers pre and post-blast monitoring in order to identify the presence of misfires and wall damage. This data can also be used to reconcile the blast results with expected results. GEM Systems offers a UAV equipped with a magnetometer for mineral exploration surveys. Companies such as Barrick have even used UAVs for solids modelling at tailings dams and stability monitoring [11]. UAVs in mining have so far been mostly limited to surface applications. Harsh underground environments pose many obstacles for flying UAVs. The confined space, dampness, reduced visibility, air movement, and lack of control signal propagation hinders most operators from being able to fly a drone underground. It may be that truly practical uses for UAVs underground will require either autonomous or semi-autonomous flight capabilities. Although there are many difficulties with flying underground, the potential benefits from a working system could greatly improve mining operations. The potential benefits of deploying UAV platforms underground include access to unreachable and dangerous locations and removing personnel from unsafe situations. These benefits have the potential to greatly improve mine monitoring and mine safety. Research has shown that current UAV technologies exist that allow for autonomous indoor flight. Extensive research has been done to develop UAV systems that are capable of performing on-board simultaneous localization and mapping (SLAM), which can allow them to navigate and map a foreign environment autonomously [5, 7, 1, 13]. Grzonka [7] successfully used an open hardware quadrotor to autonomously navigate and map an office building. The research outlines the localization, mapping, path planning, height estimation and control of the autonomous quadrotor. Other research has been done that exploits autonomous UAVs for search and rescue. Kassecker [9] proposed a software and hardware framework for a quadrotor capable of indoor and outdoor urban search and rescue and Rudol [19] developed a system for human body detection and geolocalization using an autonomous UAV.The use of autonomous UAVs in search and rescue has the potential to improve situational awareness and surveillance for a rescue team. Similar research has been done that uses autonomous UAVs for indoor exploration. Pravitra [15] outlined a strategy for autonomous exploration of miniature air vehicles (MAVs) within indoor environments and Rudol [19] developed a method for using an unmanned ground vehicle (UGV) and a UAV for cooperative indoor exploration. These studies show that aut"

Jan 1, 2017

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas'-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller's results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By J. T. Norton, N. J. Grant, A. M. Gervais

An investigation of the creep deformation of coarse grained specimens of high purity aluminum in the temperature range 800' to 1150°F, permitted the formulation of a theory explaining the formation of kink bands. The X-ray Laue back-reflection technique was used in conjunction with metallographic studies to determine the crystallographic elements involved in kinking and to measure the rotations of the bands. THAT metals deform by the formation of kink bands was first shown by Orowan1 in 1942. He manually compressed cadmium single-crystal wires, the basal planes of which were parallel to the wire axis and to the direction of compression. He found that kink bands formed quite suddenly. Similar experiments were performed by Hess and Barrett2 on zinc with a hydraulic testing machine to permit slower more careful compression. They observed that kink bands formed progressively with increasing compression. Recently Washburn and Parker3 observed the formation of kink bands in zinc single-crystal specimens during tension creep experiments. Each of these authors observed the formation of kink bands but did not make a complete analyses of the crystallographic relationships. In addition to these studies on hexagonal metals, extensive investigations of kink bands formed during deformation of aluminum in tension at room temperature have been reported. Two types of deformation bands have been observed during the deformation of aluminum at room temperature. These are: 1—kink bands: these exhibit a sharp curvature of the lattice along certain surfaces, and 2—bands formed by secondary slip along certain regions. As has already been shown by Honeycombe,4 it is important to point out the difference characterizing the two types of deformation bands. In the first case. for moderate deformation, only one system of slip planes normally is active; in the second case, however, two slip systems are active in the deformation bands. The latter have a rumpled surface because of the action of the two slip systems. In the present work only the first type of band was studied and to avoid any confusion the term kink band will be used exclusively in the discussion which follows. Honeycombe' observed that very narrow kink bands (width, 0.01 mm) formed along the (110) planes during the deformation of high purity aluminum at room temperature. He observed narrower ones in commercially pure aluminum and found also that their width decreased with decreasing temperatures. Chen and Mathewson5 reported the formation of kink bands during tensile testing of single crystals of high purity aluminum at room temperature and pointed out that some rotation occurred about a [211] direction. The influence of the bending moment on the formation of kink bands was illustrated by the experiments of Rohm and Kochendorferl0 who used a Polanyi apparatus to obtain pure shear during tensile testing of high purity aluminum single crystals. Under such conditions only one slip system was active and no kink bands were observed; consequently, no noticeable asterism was detected after deformation. The influence of the bending moment will be discussed later. Calnan8 gave a detailed account of the work performed on kink bands obtained during deformation of aluminum at room temperature. Because no explanation of the process of kink band formation in creep was given in the previous investigations and because some of the above observations are vital to the following discussion the brief review of the literature has been given. It is proposed in the present work to provide an explanation of kink band formation during creep deformation at elevated temperatures. Results and Discussion Kink bands were observed to form during tensile creep testing of high purity (99.995 pct) aluminum specimens, having a very large grain size, in the temperature range 800" to 1150°F. These observations were made in grains which occupied the whole cross section of the specimen. The specimens were originally round (diameter, 0.187 in.; gage length, 1 in.) and had parallel flats milled to permit easier metallographic and X-ray observations. The final

Jan 1, 1954

By E. Topanelian

In recent years considerable progress has been made in the development and application of mathematical techniques for the solution of certain problems involving economic "strategies". Such a problem might involve, for example, the scheduling of shipments of a commodity from a number of sources to a number of destinations. The object would be to schedule the shipments in a manner so as to satisfy the destination requirements and at the same time minimize the transportation costs. The solution to such a problem is not necessarily intuitively obvious. The "obvious" solution is frequently far from optimum. If the shipments are to be made from, for example, only two sources to four destinations, the optimum schedule is readily found. However, if shipments are to be made from, say, 10 sources to several hundred destinations, even a competent and experienced scheduler may spend considerable time in finding a reasonable answer. Even then he is not sure that he has the optimum solution. Furthermore, he has no way of knowing how far from optimum the answer is. Consequently, he does not know whether he should accept this solution or seek a better one. Prior to the advent of large high-speed digital computers, little more could be done with such problems because of their great size and multiplicity of possible solutions. A problem involving 20 sources and 50 destinations would require choosing, from a very large number of possible combinations, the optimum combination of 1,000 variables. The best one could do was to utilize intuition, extrapolation from past experience, and other non-exact approaches. With a high-speed computer, however, such problems can be solved providing a reasonable computational procedure (or algorithm) can be utilized. The purpose of this paper is to describe such a procedure (linear programming) and to apply this procedure to a production scheduling problem. THE RESERVOIR PROBLEM A simple reservoir model is used throughout this analysis to illustrate the use of linear programming. Even the simplest reservoir behavior problem is nonlinear in both space and time, but if the geometry is fixed for a particular study and the time variable is quantized, the resultant system may be described by linear constraints on the variable production rates. Crude oil is available from five separate sources and is delivered to a pipeline. Sources 1, 2, 3 and 4 are ideal reservoirs and therefore are subjected to reservoir flow restrictions. Source 5 is labeled "Outside Source" and crude oil from this source is assumed to be available in unlimited quantities without considering reservoir conditions. The general problem is to determine the schedule of crude oil production from five sources which, over an eight-year period and subject to certain restrictions, will result in maximum profit. "Profit" is defined to encompass all economic factors and expenses involved in producing and selling crude oil to a pipeline facility. Such economic factors are lumped to the extent that the term represents the dollars profit from one barrel of crude oil from any of the five sources entering the pipeline system. Table 1 presents the assumed profit per barrel for oil from any of the five sources over each of the four time periods. In this chart the eight years is broken into four equal intervals of two years. The selection of four time periods of two years each and the values representing the potential profit are completely arbitrary selections for the purpose of illustrating this idealized reservoir problem. In any practical application of this work a careful economic study would have to be made in order to estimate potential unit profits for time periods as far in the future as eight years. In this particular study it is assumed that all of the potential unit profits (C) are known and fixed before the production schedule is attempted. For each reservoir arbitrary values are assigned for

By W. B. Nelligam, J. Tittman

Refinements in radiation logging techniques during recent years have involved increasing usage of scintillation ditectors. These detectors produce voltage pulses whose heights are related to the energies of the gamma rays which initiate them. Analysis of the gamma-ray spectrum, as indicated by the pulse heights, yields information about the chemical elements composing the formations surveyed. Refined scintillation counter techniques can furnish chemical information concerning earth formations in situ, from a study of the gamma-ray spectra emitted by the formation either naturally or as a result of neutron bombardment."' Accompanying the rising interest in gamma-ray scintillation spectroscopy, there has been increased activity in the development of accelerator-type neutron sources (in contrast to encapsulated chemical-mixture sources). Such neutron generators are attractive for several reasons: (1) they greatly reduce radiation hazards to personnel; (2) there is a great reduction in contamination danger if they are lost in the hole; (3) they can produce larger neutron intensities than can conveniently available encapsulated sources; and (4) they are capable of being pulsed, thus permitting new techniques in logging. In the past, both accelerator and encapsulated neutron sources have been used by others in conjunction with scintillation-detector pulse-height analysis. The results have not been too encouraging, due to the interference among different gamma-ray spectral lines and to the fact that the gamma-ray peaks were not too clearly distinguishable above the large and ill-defined background "noise".'.' This paper is a status report on laboratory studies of a technique using a borehole accelerator as a neutron source, which gives an improved scintillation spectrum, thus permitting more accurate chemical analyses of the formations penetrated. The Schlumberger-accelerator neutron source' is presented; the origins of inelastic and of thermal-neutron, capture gamma rays are dis- cussed, and results are given for some laboratory measurements performed in borehole geometry. THE NEUTRON GENERATOR-TUBE ACCELERATOR In the attempt to develop an accelerator neutron source which would have the desirable properties outlined in the "Introduction" and which would operate satisfactorily under borehole logging conditions, a small-diameter, cylindrical, neutron generator tube has been developed. This tube utilizes the principle of accelerating deuterons (heavy hydrogen nuclei) by a high voltage so that they bombard a tritium (heavy, heavy hydrogen) target. The resulting reactions produce large numbers of neutrons of 14 Mev energy. Furthermore, the tube is permanently sealed and, thus, the use of pumping techniques in the sonde is avoided. The tube consists of a pressure control, an ion source in which the deuterons are stripped of their electrons, an accelerating gap down which the deuterons are sped by the high voltage, a secondary electron suppressor and a tritium-loaded target. Obviously, the tube requires auxiliary circuitry in a sonde for the control of the ion source, tube pressure, high voltage and other operating variables. Neutron yields, both continuous and pulsed, have been produced under simulated field conditions in the range between 1 and 10 times those conventionally used in neutron logging applications. For the experiments discussed in this paper, neutron-pulse repetition rates in the range between 500 and 5,000 pulses/sec are adequate. Furthermore, by making suitable adjustments in operating conditions, one can vary the pulse width from a minimum of several microseconds up to dc. As will be seen later, the present application to well logging does not require that one have available pulses of neutrons which are appreciably shorter than the time it takes for fast neutrons to slow down to thermal energy in water. As a consequence, much of our interest has been directed towards the operating characteristics of the tube with pulses longer than about 10 microseconds in duration. Our experience to date shows that tubes can be operated in this manner with reasonably constant, average neutron outputs.

By Ken-ichi Hirano, B. L. Averbach, Morris Cohen, N. Ujiiye

The influence of plastic deformation in compression on the self-diffisivity of a iron has been measured in the temperature range of 742º to 885°C. The diffusivity is enhanced in proportion to the strain rate, but this dependence decreases with increasing temperature. The magnitude of the strain is relatively unimportant in this connection. Strain rates from 0 to -2 x 10-3 sec-1 were investigated, with a corresponding increase in diffusivity up to 3 x 103 times. The results are analyzed in terms of vacancy diffusion and the excess vacancies introduced during dejornzntion. It is concluded that grain boundaries are the main vacancy sinks in polycrystalline iron and that the vacancy lifetime is therefore dependent on the grain size. A typical vacancy lifetime at 800°C with a nzean grain-boundary intercept of 0.3 cm is of the order of 1 sec. LATTICE imperfections are considered to be the dominant mechanism of diffusion in close-packed crystals. Vacant lattice sites are usually the most important defect, but the enhancement of diffusion along grain boundaries and other short-circuiting paths has been well documented. Inasmuch as plastic straining produces a variety of defects, it would be expected that the application of stresses sufficient to cause plastic flow during a diffusion run might result in an increase in the self-diffusivity. Enhanced self-diffusion in a iron under compressive creep has been observed by Buffington and Cohen' and subsequent coworkers.2 Lee and Maddin3,4 have also found an enhancement in the self-diffusivity of silver single crystals undergoing torsional deformation. On the other hand, Darby, Tomizuka, and Balluffi5 have not observed any significant diffusion increase in single crystals of silver deformed either in tension or in compression, although the results in this case may have been affected by the occurrence of recrystallization. With poly crystalline silver, however, Forestieri and Girifalco6 have found en- hanced diffusion during compressive creep. A small decrease in the self-diffusivity of zinc single crystals under compressive loading has been reported,7 but the applied stresses in these experiments were too low to cause appreciable plastic deformation. For the cases in which enhanced diffusion has been found, a number of common effects have been noted. Usually, the ratio of the diffusivity during straining, Ds, to the diffusivity in the unstrained condition, D,, varies linearly with the strain rate, and does not seem to be a function of the plastic strain, at least up to ? = 0.2. Values of Ds/Du of the order of 100 have been reported'-4 but the enhancement decreases with increasing diffusion temperature. The activation energy for diffusion decreases markedly with increasing strain rate and approaches a level of roughly one-half the unstrained value. Most of the current explanations of these observations assume that excess vacancies are generated by the deformation, but rather large concentrations of excess vacancies and/or long lifetimes are required by these theories.',',',' Since the original work on iron, there has been an important improvement in the compressive-loading technique.2 The early data were also analyzed without taking into account the motion of the coordinate system in the specimens undergoing plastic deformation. This problem has now been solved by a number of investigators (see Refs. 11, 12, and the appendix of the present paper). Consequently, the strain-diffusion runs on iron for the Ds measurements were repeated and extended over wider ranges of temperature and strain rate. In addition, the original diffusivity data9 on the unstrained iron have been corrected for the divergence of the beam in the counting system and new data points have been added, resulting in somewhat different values for the self-diffusi-vities in a iron.10 The new values for Du are used in this paper. A comparison of Ds and Du confirms the large enhancement in the self-diffusivity of iron during compressive creep, as previously reported.',2

Jan 1, 1963

By David E. Nicholas

INTRODUCTION In this paper, a numerical process for selecting a mining method, with the emphasis on underground mass mining techniques, such as caving, induced caving, and stoping, is proposed. In the past, selection of a mining method for a new property was based primarily on operating experience at similar type deposits and on methods already in use in the district of the deposit. Then, the chosen method was modified during the early years of mining as ground conditions and ore character were bet- ter understood. Today, however, the large capital investment required to open a new mine or change an existing mining system make it imperative that the mining methods examined during the feasibility studies and the method actually selected have a high probability of attaining the projected production rates. Although experience and engineering judgment still provide major input into the selection of a mining method, subtle differences in the characteristics of each deposit, which may affect the method chosen or the mine design, can usually be perceived only through analysis of measured characteristics. The parameters that must be examined when choosing a mining method include: 1) geometry and grade distribution of the deposit; 2) rock mass strength for the ore zone, the hanging wall, and the footwall; 3) mining costs and capitalization requirements; 4) mining rate; 5) type and availability of labor; 6) environmental concerns; and 7) other site-specific considerations. This paper encompasses a detailed look at the first two parameters since they, plus mining costs, have the greatest impact on the selection of a mining method. The proposed method selection process is for a project where drilling has defined sufficient geologic reserves, but little or no underground development has been done. Since each deposit has its own characteristic geanetry/grade distribution, and rock mechanics properties, mining method selection should be at least a two-stage process. In Stage 1, the deposit is described in terms of geometry, grade distribution, and rock mechanics properties. Using these parameters, the mining methods can be ranked to determine which are most applicable; they can then be considered in general terms of mining and capitalization cost, mining rate, type and availability of personnel, environmental concerns, and other site-specific considerations. In Stage 2, the most likely mining methods are costed out, based on a general mine plan. Mining and capitalization costs are used to determine a cut-off grade fran which a minable reserve can be calculated; economic comparisons can then be made to determine the optimum mining method and economic feasibility. During the mine planning phase of Stage 2, rock mechanics information would be used to provide realistic estimates of underground opening size, amount of support, orientation of openings, and caving characteristics, and open pit slope angles. If ground control or operational problems should be encountered with the methods being considered, modifications could be made. Although planning on paper extends start-up time, it is cheaper to err on paper than to find the error after mining has begun. METHOD SELECTION - STAGE 1 The main purpose of Stage 1 is to select those mining methods which should be considered in greater detail. The simplest way to do this is by defining those characteristics required for each mining method and then determining whether the characteristics of the deposit are suitable. However, no one mining method is so restrictive that it can be used for only one set of characteristics, as indicated by the classification system proposed by Boshkov and Wright (1973). In the mining method selection proposed, geometry, grade distribution, and rock mechanics characteristics are ranked according to their acceptability for ten general mining methods.

Jan 1, 1981