Search Documents

By Ralph L. Erickson

The discovery of a Carlin-type gold deposit at Cortez, Nevada, in 1966 can be attributed directly to the use of geochemical exploration techniques. Most mineral deposits owe their discovery to geologic concepts, geologic analysis, and luck, but at Cortez, geochemistry as an exploration tool played the clearly dominant role in discovery. Of course, the economic significance of the discovery had to be determined by exploration and development drilling. The story of the discovery began in 1959 when field work was initiated on a new project, "Geochemical Halos Utah and Nevada," by R.L. Erickson and A. P. Marranzino for the US Geological Survey. The Cortez district was selected for geochemical work because it was an area with good geologic control where we could address the problem of how to prospect in barren outcrops for concealed ore deposits in potentially favorable structures (buried thrust zones) or favorable host rocks in the subsurface. Geologic mapping of the Cortez 15-min quadrangle, just being completed by Gilluly and Masursky (1965), showed that the quadrangle contained excellent exposures of both the upper and lower plates of the Roberts Mountains thrust fault, a major structural feature of north-central Nevada. Roberts (1960) had noted that a number of mining districts were associated with windows in the thrust. In 1959, Erickson and Marranzino did some reconnaissance rock sampling and spring-water sampling in the siliceous clastic rocks of the upper plate of the thrust. Results of this reconnaissance prompted a full-scale sampling program in 1960 in the upper plate rocks on the west flank of the Cortez window. Results of the investigation showed that anomalously high concentrations of metals occur in the upper plate rocks, and further, that the distribution of metals is fault controlled and shows a pronounced zoning pattern (central copper zone; intermediate zinc, copper, and lead zone; and outer arsenic zone). The anomalies were interpreted as primary leakage halos that originated from metal occurrences in the thrust zone or in carbonate rocks below the thrust and moved upward along normal faults that cut both upper and lower plate rocks. Erickson gave a talk about these anomalies in the summer of 1961 to the local AIME Section in Reno, Nevada; two short reports were published that year-" Geochemical Anomalies in the Upper Plate of the Roberts Thrust Near Cortez, Nevada" (Erickson et al., 1961) and "Hydrogeochemical Anomalies in Four Mile Canyon Near Cortez, Nevada" (Erickson and Marranzino, 1961). These releases prompted blanket staking in the area by several small companies. In 1963, geologic and geochemical mapping were started by the USGS in the lower plate carbonate rocks of the Cortez window west of the quartz monzonite stock at Mount Tenabo and north of the old townsite of Cortez. The results of the work showed anomalously high concentrations of arsenic, antimony, and tungsten in jasperoid, fracture filling, and shear zones in Silurian and Devonian carbonate rocks. The anomalous area was about 1.6 km (I mile) long and 300 m (1000 ft) wide. A brief report, "Geochemical Anomalies in the Lower Plate of the Roberts Thrust Near Cortez, Nevada" (Erickson et al., 1964a), was published in 1964. During this same time period, 1959-1964, several exploration or mining companies were active in the general area (chiefly mapping geology and acquiring property). In 1964, American Exploration and Mining Company (Amex) concluded that the entire Cortez district was worthy of an extensive exploration effort involving drilling as well as geologic, geophysical, and geochemical studies. To carry out the program, Amex formed a joint venture group with the Bunker Hill Company, Vernon Taylor, Jr., and Webb Resources. Their early efforts were directed to the upper plate rocks and to the lower levels of the old Cortez silver mine. The group also drilled some shallow rotary assessment holes adjacent to the area of the arsenic, antimony, tungsten anomaly in lower plate carbonate rocks described in Erickson et al. (1964a). Assay results from these holes offered little or no encouragement to the joint venture group. Splits of these drill samples were made available to the USGS. In order to enhance any metal content present in these rocks and to determine the mineral residence of any metals detected, heavy-mineral concentrates of each 3-m (10-ft) sam-

Jan 1, 1985

By F. Weinberg, E. H. McLaren

Dilute binary alloys have been solidified under controlled thermal conditions, and solute distributions, temperatures during freezing and melting, and the position and morphology of the solid-liquid interfaces have been examined. Binary Zn, Sn, and Pb based alloys were investigated with solute additions of AgU°, AU°, Sb'24, TlZo4, snH3, and Zns5. Approximately 150 cc of metal was used in each ingot. Appreciable solute segregation occurred in most of the ingots, depending on the solute concentrations, distribution coefficients. and freezing conditions of the alloys. Measurements were made of the general solute distribution by sectioning layers of material from the ingot circumference and counting the activity of each layer. Local segregation was detected by autoradiographic techniques. Freezing and melting curves for the same alloys were obtained by precision resistance thermometry. These curves are related to the corresponding measured solute distribution curves and those calculated from theoretical expressions. The position and morphology of solid-liquid interfaces during freezing were determined by both decanting and quenching techniques. The results indicate the manner in which freezing and remelting occurred under the conditions investigated. If an isomorphous dilute binary alloy is solidified or melted under equilibrium conditions, the concentration of solute in both solid and liquid would be uniform at all stages in the freezing process. The solute concentrations during freezing could be determined from the phase diagram of the alloy, and the freezing curves would have an alloy slope defined by the solidus and liquidus points. In practice, however, solidification under equilibrium conditions rarely, if ever, occurs, primarily because of low rates of diffusion in the solid and incomplete mixing in the liquid. As a result, the solute distribution and the shape of the freezing and melting curves obtained in a given system depends on a number of factors. These include, in addition to the extent of diffusion and mixing, the initial solute distribution and concentration, the distribution coefficient, the rate of freezing or melting, supercooling of the liquid, the morphology of the solid-liquid interface, and the effect of annealing. The purpose of the present investigation was to determine, experimentally, actual solute distributions and thermal curves for a series of dilute binary alloys, using tracer techniques to determine the solute distributions and precision resistance thermometry for the temperature measurements. In addition, observations were made of the positions and morphologies of the solid-liquid interfaces during freezing and remelting. The distribution of solute in a dilute alloy, assuming complete mixing in the liquid, and negligible diffusion in the solid, is given by1 where Co is the mean solute concentration of the alloy, C, is the solute concentration in the solidg is the fraction of the ingot solidified, and ko is the equilibrium solute distribution coefficient (assumed to be independent of concentration for dilute alloys). For the case of partial mixing, ko can be replaced by an effective distribution coefficient k. If solute transport in the liquid is assumed to be entirely due to diffusion, with negligible diffusion in the solid, then the solute distribution in the solid is given by Tiller et a1.' as where R is the rate of advance of the solid-liquid interface, D is the solute diffusion coefficient in the liquid, and x is the position of the interface from the start of freezing. The temperature at the interface TI is related to the solute concentration of the solid at the interface C, by the expression where Tm is the freezing temperature of the pure solvent and ml is the slope of the liquidus for the binary alloy being considered. Accordingly, by relating TL to the actual temperatures measured in the present experiments, solute distributions can be deduced from the temperature measurements. Using the above expressions, the measured solute distribution of a given alloy may be compared to the theoretical distributions, as well as to the distribution calculated from the observed freezing curve. EXPERIMENTAL The solvent materials used in making up the binary alloys were C. P. New Jersey zinc (99.999 pct), Vulcan Extra Pure Tin (99.999 pet), and Tadanac lead (99.998 pct). The solute elements were of 99.99 pct purity or better, and were activated in the

Jan 1, 1963

By J. T. Raleigh, D. L. Flock

This study investigated the nature of formation plugging with bacteria and attempted to relate its characteristics to physical rock properties. Fifteen core samples of four specific formation types were defined and plugged using a commonly occurring, uniform sized dead bacteria, Bacillus subtilis. Two of the formation types were fairly uniform-grained sandstones and two were heterogeneous carbonates. The injection rate and concentration of solids in the brine were held constant during the test runs. The pore geometry factor G was a significant petro-physical rock characteristic for correlation with depth of plugging, but lacked importance as a parameter for rate of plugging. Although the G factor may not completely describe porous rock geometry for plugging predictions, for reservoirs with multi-pore systems, there is evidence that the pore geometry factor might form a basis for establishing injection water specification in more llornogeneous reservoirs with only primary porosity development. INTRODUCTION GENERAL PURPOSE Injecting water into oil-bearing formations to increase recovery is common practice in the petroleum industry. One of the most serious operational problems involved with water injection wells is the tendency of the porous rock immediately adjacent to the injection wellbore to plug partially. This results in increased injection pressures and/or a decrease in the rate of injection. This plugging in injection wells can result from a number of causes, including bacteria contained in the injected water.1 The purpose of this work was to investigate the nature of injection plugging with bacteria and to attempt to relate its characteristics to physical rock properties. Many variables exist in research of this type, so, to reduce them, injection rate and concentration of bacteria in brine were held constant. The bacteria used for plugging was a uniform-sized, dead bacteria, Bacillus subtilis, a micro-organism commonly found in water, air and soils. Four specific porous rock types were plugged—two consolidated sandstones and two heterogeneous carbonates. The experimental work involved two distinct phases. The first concerned determining the physical rock properties and the second was the experimental work of actually plugging the selected core samples. PORE GEOMETRY FACTOR Capillary pressure is necessary to displace a wetting fluid from a capillary opening with a nonwetting fluid. The factors which govern capillary pressure are the surface or interfacial tension of the fluids involved, the system wettability and the equivalent pore radius. Thus: P,= 2ucososrd..........(1) Thomeer used the idea that the location and shape of the capillary pressure curve reflect characteristics of the pore structure of any porous media sample.' He presented the mathematical description of capillary pressure curves and of indicated differences in pore geometry of samples. where G = pore geometry factor PC = capillary pressure Pi = displacement pressure (V,)p, = fractional bulk volume occupied by the displacing fluid at any capillary pressure (V,), = fractional bulk volume occupied by the displacing fluid at infinite pressure or total interconnected pore volumes. The shape of the curve is defined by the pore geometry factor G. The curve shape is related directly to the sorting and interconnection of the pores of the sample. The location of an infinite number of curves can be defined by the same asymptotes. However, these curves differ in shape, each curve being described by a specific value of G. Thus, the description of a curve is unique when, in addition to the position of its asymptotes, its shape is defined. These curves can be used graphically to determine the characteristic parameters of any arbitrarily located curve. Those parameters that can be approximated from any capillary pressure curve are: (1) total interconnected pore volume, represented by (V,),- ; (2) extrapolated displacement pressure, the abscissa axis; and (3) pore geometry factor, by the curve shape. PROPERTIES OF BRINE AND BACTERIA The fluid used throughout this study was an aqueous solution of sodium chloride, 25,000 ppm by weight. The addition of salt to fresh water was intended to eliminate clay swelling in the sandstone cores and to retard leaching in the limestone cores. The viscosity of the brine solution at 100F, the temp-

Jan 1, 1966

By R. D. Burlingame, T. L. Joseph, Gust Bitsianes

DESPITE almost fifty years of commercial practice, the sintering of iron ore has received little fundamental study. Much of the theoretical work1-'has dealt with the constitution of sinter produced under widely varying conditions. While these studies have broadened our knowledge of the changes that occur in the sintering zone and in the freshly formed sinter during the early stages of cooling, they provide little insight into the changes that precede the formation of sinter. These preliminary changes merit study as a part of the overall process. Hessle. working with beds of Swedish magnetite concentrates, was one of the first investigators to study the sintering process in its entirety. On the basis of temperatures observed at various levels of the bed during sintering, he postulated a number of distinct reaction zones to account for the chemical changes leading to the formation of sinter. A more direct method of attack is that of arresting the sintering zone after it has progressed part way through the bed. A study of a vertical cross section through such a quenched bed provides direct information on the changes taking place at various levels. This method was used by McBriar et al.' to show that several well-defined zones of chemical change existed within beds that were typical of British sintering practice. The same general method of attack was developed independently in the present investigation to study partially sintered beds typical of American practice. Experimental Sintering Equipment The sintering operation was carried out on an experimental scale with the equipment shown in Fig. 1. The refractory-walled sintering chamber A was 11 in. deep and averaged 9 in. in diameter. Air was introduced through a tapered flow section B, which contained the orifice C for accurate metering of the incoming air. This section was located directly above the square ignition housing D, which in turn rested upon the sintering chamber A. The bed was ignited with burner E. The required suction for the operation was furnished by a fan F, which had an air capacity of 500 cfm (stp). Hot exhaust gases from the sintering chamber were cleaned in the dustcatcher G before entering the exhaust fan. In the study of partially sintered beds, it was essential to find some technique for removing the entire charge from the sintering pot without disarranging the unsintered bottom portion. This problem was finally solved by sintering the charge in a removable basket, which snugly fitted the sintering chamber. This basket was constructed of two thicknesses of window screen and was lined with a 3/16-in. layer of asbestos paper. The bottom of the basket consisted of two thicknesses of wire screen, which were fastened to the basket wall. For high fuel mixtures, additional insulation was provided by a somewhat thicker layer of asbestos cement. Preparation of Partially Sintered Mixtures The moist feed was carefully placed in the sintering basket, to prevent segregation of the particles, which varied widely in size and composition. A thermocouple was placed in the center of the basket with the hot junction halfway down, and the mixture was evenly distributed around it. During ignition and throughout the sintering of the upper half of the bed, the hot junction temperature increased very little. When the sintering zone reached the halfway point, as indicated by the sudden increase in the hot junction temperature, the charge was quenched. During quenching the suction was turned off and the orifice was tightly stoppered to prevent further influx of air. At the same time, nitrogen was admitted to the sintering chamber through the orifice tap. As soon as the nitrogen had displaced the air and products of combustion, the charge was removed from the sintering pot for immediate dissection. It is impossible to preserve the exact zone structure of the bed at the instant that combustion is arrested unless the downward transmission of heat is also immediately stopped. Fortunately, heat transfer is very slow in beds containing a stationary fluid, especially if the particle size is small. It follows that the minimum quantity of nitrogen should be used to displace the air and that static conditions be established as soon as possible. A very steep temperature gradient across the combustion zone for some time after the quench was evidence of in-

Jan 1, 1957

By K. G. Davis, P. Fryzuk

The diffusivity of silver in liquid tin has been determined, using the capillavy-reservoir technique, over the temperature range 250° to 500°C. The new value, D = 2.5 x 10'* exp(-2480/ RT) sq cm per sec, differs from that obtained by other workers in an earlier investigation. The analysis of data from the capillary-reservoir technique is discussed. In a recent investigation of the solidification of dilute alloys, values for the diffusion constant of silver in liquid tin were required in the analysis of the formation of impurity substructures. AS a result, measurements were made of the diffusion constants in the temperature range 250° to 500°C (melting point of tin 232°C), for the alloy concentrations used in the solidification experiments and at higher concentrations, to verify previous determinations.I)2 The capillary-reservoir method was adopted, using experimental procedures similar to those followed by Ma and Swalin,1 with the main exception that radioactive silver was used in the present investigation to facilitate solute-concentration measurements. EXPERIMENTAL PROCEDURE a) 100 ppm Samples. Glass capillary tubes of 2 mm inside diameter and approximately 5 cm- long were sealed at one end, evacuated, and filled with tin of 99.999 pct purity. The region of shrinkage near the mouth was cut off, and the tubes were then placed in a graphite holder and immersed, with the open end up, in an unstirred bath of alloy containing 100 ppm Ag 110, where they remained for periods of up to 30 hr. On removal from the bath they were cooled by an air blower. The bath was kept under a small positive pressure of argon, and the temperature controlled to within +1°C. A 10-hr diffusion period was used in the majority of the tests, scatter on runs of less than 5 hr being rather large. The procedure outlined above was chosen in preference to putting alloy in the capillary and pure tin in the bath, in order to avoid segregation when the tubes filled with alloy were first solidified. To minimize segregation when the diffusion period was complete and the capillaries again solidified, the earlier samples were held in thin-walled silica tubes which could be cooled very rapidly. Later tests were made in precision-bore Pyrex tubes, to eliminate effects caused by variations in the capillary diameter. No consistent differences in diffusivity as measured in the two types of tube were detected. After removal from the glass tubing, the samples were sectioned into 2.5 mm lengths and counted for y activity, using a scintillation counter with fixed geometry. Samples were also drawn directly from the bath and counted, so that values for C/C,, the ratio of the weight of Ag 110 in the sample to that in the bath, could be obtained. b) 5000 ppm Alloy. To check for possible effects of concentration, the silver content of the bath was increased to 5000 ppm. Complete mixing was found to have taken place in the capillary after a 10-hr period at 300°C. It appears that the greater density of the alloy was sufficient for buoyancy forces to cause instability in the alloy-tin interface, leading to rapid convective mixing. For the 5000 ppm alloy, therefore, the bath was of pure tin and the capillary tube was filled with alloy. With this arrangement, values of D consistent with those for the 100 ppm alloy were obtained, Fig. 1. CALCULATIONS OF DIFFUSIVITY The terminology used applies to a capillary of pure tin immersed in a bath of alloy. 1) Error-Function Method. under the present experimental conditions, the rod of liquid tin into which silver is penetrating may be considered semi-infinite. Assuming the concentration at the mouth of the tube to remain constant at Co, the concentration C at distance x from the mouth of the tube at time / is given by3 Plots of the inverse error function of (1 -C/Co) vs .v gave straight lines passing through the origin with slope 1/2-, x being corrected for shrinkage both on solidification and while cooling to the melting point (total correction about 6 pct at 500°C). Values for log D obtained in this manner are shown in Fig. 1. A least-squares fit to the relation D = Do exp(-Q/RT)

Jan 1, 1965

By J. D. Alexander, W. L. Martin, J. N. Dew

This paper presents data obtained using a flood-pot technique to determine the fuel available and the corresponding theoretical air requirements for in situ combustion of crude oils. Since the technique is relatively quick and easy, it is a practical and convenient tool for evaluating reservoirs as fireflood prospects. It is also a research tool which facilitates systematic study of the variables affecting fuel availability and corresponding air requirements. The understanding of these variables is of prime importance to those concerned with the technical and economic development of in situ combustion as an oil-recovery process. The experimental results show conclusively that the fuel available for in situ combustion is not a constant but, rather, varies with crude-oil characteristics, porous-medium type, oil saturation, air flux and time-temperature relationships. Thus, the fuel availability for specified field applications should be determined using actual reservoir crude and core material and the process conditions expected during in situ combustion in the reservoir. INTRODUCTION In situ combustion is a thermal process for recovering crude oil from reservoirs. The thermal energy released during the combustion of a small amount of the oil in place aids in the displacement of the remaining oil. Numerous articles have been published describing the in situ combustion process giving detailed results of laboratory and field experiments.10 In order to engineer an in situ combustion project, a number of important factors must be considered and determined. These factors include: (1) the amount of fuel consumed per unit of reservoir volume swept by the combustion zone, (2) the composition of the fuel consumed, (3) the amount of air required to consume this fuel, (4) the portion of the reservoir swept by the combustion zone, (5) the appropriate air-injection rates and pressures, (6) the amount of oil that will be recovered, (7) the rate of oil production and (8) the operating costs. Nelson and McNiell1 recently have described a procedure which utilizes laboratory combustion-tube data as a basis for the calculation of some of these design factors. Various authors have attempted to describe the in situ combustion process mathematically, and considerable progress has been made. Analytical solutions to the problem of heat transfer from a moving combustion front have been obtained for linear and radial systems."-' All of the published results involve the assumptions that: (1) fuel concentration is constant throughout the reservoir, or that fuel concentration is inversely proportional to the velocity of the front for a given rate of oxygen consumption; and (2) the fuel reacts instantaneously with injected oxygen, while liberating a constant amount of heat per unit weight of fuel at all temperatures. It seems both desirable and reasonable to test the validity of these assumptions experimentally. This paper presents laboratory data which were obtained by means of a "fire flood-pot" method for determining fuel availability and composition, and the corresponding theoretical air requirements for in situ combustion of crude oils under variable conditions. The mechanics of the method are similar to a conventional tube-run experiment.' The important differences involve the size of the reservoir model used and the method for providing the experimental environment. The new method subjects conventionally-sized core samples or unconsolidated sands to a programmed environmental sequence similar to that experienced by a similar volume of rock during the approach and passage of a combustion front in a long tube or in an oil reservoir undergoing in situ combustion. Restored-state samples can also be used. The small samples and relatively simple techniques involved allow an experiment to be set up, run and calculated in about three 8-hour days. This is a considerable improvement over long-combustion-tube techniques which can require several days to run and several more work days to set up and calculate. All the runs presented were run at 40-psig injection pressure. Pressure was not considered as a variable for these experiments, since we previously had found that it had only a small effect on fuel availability up to 600 psig.APPARATUS AND MATERIALS APPARATUS The fire flood-pot apparatus consists of a consolidated

By H. C. H. Darley

The influence of particle size and concentration on the development of chip hold-down pressure (CHDP) was studied in an apparatus designed to measure the change of filtration rate during the first second of the filtration process. CHDP is controlled by the "bridging" particles (i.e., particles in the 50 to 0.2 micron range), whereas filter loss is controlled by the colloid fraction. The results indicated that a fast drilling fluid with a low filter loss could be obtained if the concentration of bridging solids and the viscosity were kept very low. A novel group of drilling fluids, which we have named colloid emulsions, was developed to meet the above requirements. These colloid emulsions are made by dispersing art oleophilic colloid in oil and emulsifying about 5 per cent of this dispersion in water or a dense brine. Good results have been obtained in field trials with emulsions which use asphalt as the oleophilic colloid. INTRODUCTION In recent years it has been found that fast rates of penetration can be maintained in hard formations by using "clear" water as a drilling fluid. In some cases, however, water cannot be used because the very high down-hole filter loss creates various difficulties in drilling, in logging, or in producing the well. It then becomes necessary to "mud up", which results in a decrease in drilling rate. The reason that mud decreases drilling rates is a complex subject which has been well covered by such authors as Garnier and van Lingen, Cunningham and Eenink,' and van Lingen.3 One of the principal factors involved is the tendency of the drilling fluid to develop a filter cake on the bottom of the hole. The pressure differential which develops across this filter cake opposes the action of the bit in dislodging chips from the rock. This pressure differential is referred to by Garnier and van Lingen' as the chip hold-down pressure, which is denoted by the initials CHDP in this paper. The work described in this paper was undertaken to develop a drilling fluid which would minimize CHDP but which would still have a low filter loss—in other words, a fluid which would deposit a filter cake on the sides but would lay little or no filter cake on the bottom of the hole. EVALUATION OF CHDP EXPERIMENTAL Conditions governing the growth of filter cake on the sides of the hole are obviously different from those on the bottom of the hole. On the sides of the hole, cake will form continuously but at a decreasing rate, until eventually the rate of growth will equal the rate of cake erosion by the mud stream.' On the bottom of the hole, each time a bit tooth dislodges a chip the cake is removed, a fresh rock surface is exposed, and the filtration process must start over again. Thus the CHDP is determined by the amount of cake which can form in the interval between successive tooth strikes at a particular spot, normally a fraction of a second. Our first step therefore was to study filtration rates over the first second or so of the filtration process. These tests were conducted in the simple type of dynamic filter cell shown in Figs. 1 and 2. In this cell the mud was filtered through a train of rock cores while rotating blades kept the mud flowing over the surface of the core. The procedure was to evacuate the cores to 0.05 mm of mercury, saturate with brine, and assemble in the holder under brine. The permeability of the core train was then determined. A plastic seal was placed on the face of the core train, and a rip cord was attached from the seal to one of the blades. At the start of the test, 500 psi was applied to the mud, and the meniscus in the capillary measuring tube was adjusted to zero. When the blades were started, t,he seal was ripped from the face of the core, and filtration began. The filtration rate was at first measured by photographing the advance of the meniscus in the capillary tube. Using

Jan 1, 1966

By John M. Tinsley, Calvin D. Saunders, H. K. van Poollen

In the past few years much con-sideration has been given to the evaluation of the effect of hydraulic fracturing on the productivity of wells. Generally, these studies included the evaluation of fracturing materials, fracture extension and formation damage due to the use of various fracturing fluids. Only little consideration has been given to the characteristics, and in particular the flow capacity, of the fracture itself and its effect on well productivity. This paper presents the results of laboratory investigations pointed toward fie evaluation of the efficiency of various fracrures with special emphasis on the flow capacity of these fractures. Data presented in this paper are the results of both an electrical model study and physical testing. Under consideration are (I) effect of overflush, (2) premature production of well after treatment, (3) "tailing-in" with coarse sand near the end of the treatment, (4) effect of propping agent size and concentration, (5) reduction in effective frac-ture permeability caused by formation caused by formation fines, silt and clays, and (6) effect of various fluids on formation strength and competency. The results of this investigation indicate that the flow capacity of a fracture is affected by any or all of the various parameters mentioned above. The authors believe that a better understanding arid utilization of these factors should result in more efficient formation fracturing. INTRODUCTION Hydraulic fracturing has become almost a standard practice of many companies for stimulating production from old and new wells. Although most companies utilize this service, techniques of application vary widely between companies and areas. Probably too often when a well in an area responds favorably to a particular technique all future wells in the same area are treated in a similar manner. Possibly a modification of the technique would result in a further production increase. Variables, of which many are extremely difficult to evaluate from field results, hamper the selection of procedure changes. Attempts are being made by a number of organizations to analyze statistically treating techniques from production data. This is a very worthy and necessary approach but very possibly laboratory investigations may aid in evaluating some of the variables which tend to affect the results of a fracturing treatment. Some of the factors cannot be studied from practical field experience and only laboratory tests can show the possibilities which might exist. One of the factors which appears to be of major concern today is the flow capacity of the created fracture and how it can be changed. Papers on this general subject have been primarily concerned with the size of propping agent and the extent of fracture. Papers have also been written on the possible permeability damage to formations by fracturing fluids. In addition, it might be possible that another type of flow restriction is prevalent. This would be a restriction of flow through a sand-packed frac- ture caused by foreign materials integrating within the propping agent. This paper presents preliminary data obtained in an attempt to evaluate the effect of some factors affecting flow through sand-packed fractures. No attempt is being made to offer a fracturing technique adaptable to all areas and conditions, but to furnish data tending to show the possible effects which might be caused by variations in procedures and materials. PROPPING AGENT PERMEABLLITIES A hydraulically induced fracture containing sand as a propping agent may theoretically be classified as a packed-sand system. The flow of fluids through such packed systems has been the subject of much research. Although there have been numerous methods proposed for the evaluation of such systems, most writers agree on the general properties affecting their flow capacity or permeability. These properties include porosity, particle size, sphericity and the roughness of the particle. In some methods of evaluation . the particle size and sphericity terms are combined to produce an equation which is a function of the surface area of the particles. In this study the permeability of various fracturing sands was both measured and calculated. The apparatus for the permeability measurements consisted of a 52-in. Lucite tube with a 2.5-in. ID. A screen and drain plug were fitted in the bottom of this tube to retain and hold the sand in place while allowing fluid flow. Two pressure taps consisting of thin, highly perforated

By J. B. Cohen, O. F. Kimball

Aging of specimens of Au-20 at. pct Ni and Au-40 at. pct Ni has been followed with electron microscopy, resistance, magnetic measurements, and X-ray diffraction. The periodic structures observed by electron microscopy in Au-Ni alloys by Fisher and Embury are shown to be due to nickel-rich prepre-cipitates. Assuming this alloy system undergoes spinodal decomposition, our observations and available data on thermodynamic properties and elastic constants indicate that the gradient energy coefficient in this system is of the order of 10-5 ergs cm-1 . Shifts of X-ray peaks in this system have been found to be due primarily to residual stresses, not faulting. A large effect of applied stress during aging was detected by X-rays and electron microscopy for the Au-20 at. pct Ni alloy. The presence of a preprecipitation stage during aging of gold-rich Au-Ni alloys for short times below 300°C has been indicated by several previous investigations. Changes in electrical resistivity,1"3 elastic modulus,4 superparamagnetisrn,5 positions of X-ray diffraction peaks,6 and hardness7 have been studied. Electron diffraction and microscopic studies by Fukano8 and Fisher and Embury9 respectively showed that the precipitation process resulted in the formation of periodic composition fluctuations—perhaps spinodal decomposition. MOSS10 detected a satellite with X-rays in quenched specimens, but the wavelength was much smallerthan that found by Fukano. Recently Woodilla and Averbac11 detected a modulation similar to that found by Moss, but with electron diffraction. They suggested that the microstructures observed by Fisher and Embury may be simply due to a Moir'e effect from a back deposit of gold during electropolishing. Satellites in electron diffraction can be enhanced by small bends or tilts of a foil. Even if this were not the case, asymmetric satellites, as are found in Au-Ni are difficult to analyze. They are principally influenced by the large differences in atomic size in this system, and obtaining information on the associated composition fluctuations depends on assumptions as to whether or not there are isolated clusters or composition "waves", and in the latter case whether or not the entire specimen has decomposed. Assuming that decomposition was complete in alloys aged at low temperatures, Woodilla and Averbach suggested that the satellites indicated a composition which varied about the average value by only± 5 pct, much less than that predicted from the phase limits, Fig. 1; they attributed this small variation to a cessation of the reaction caused by the loss of mobile vacancies. In Ref. 12 it was shown that vacancies are trapped by the preprecipitate in this system. Furthermore, after reversion, aging proceeds at a much slower rate than immediately after a quench from the solution temperature,5'11,13 indicating the importance of vacancies in the kinetics. In some senses Au-Ni is a favorable system for studying preprecipitation and reversion mechanisms because of the simplicity of the equilibrium diagram14 and the availability of good thermodynamic15 and diffusion data.16,l7 The calculations necessary to relate the decomposition process to classical nucleation and growth or to spinodal decomposition, as well as those necessary for understanding reversion mechanisms, can, therefore, be made. As an example, the position of the chemical spinodal, calculated from the thermo-dynamic data of Sellars and Maak15 with the method described by Rundman and Hilliard,18 is included in Fig. 1. It is compared to the chemical spinodal calculated from the equilibrium miscibility gap by Cook and Hilliard's method.l9 Included is the coherent spinodal suggested by Golding and MOSS," using measured elastic constants for their calculations and the same ther-modynamic data. The spinodal estimated by Woodilla and Averbach,11 by determining the temperature limits for the observation of the satellites discussed above, is also given. Actually this last curve is open to question on the following grounds: As a phase boundary is approached the spacing of particles or the wavelength

Jan 1, 1970

By J. P. Denny

The constitution of the Ga-In system was determined by thermal methods. An experimentally determined metastable equilibrium line (an extension of the indium-rich liquidus) was obtained. The various alloys were studied metallographically using polished samples obtained by a casting method. These low melting alloys required a special dry-ice assembly to maintain a suitable temperature. RECENT interest in alloys that are liquid at room temperature has led to rather extensive investigation of gallium-base alloys. Widely distributed over the earth, gallium could be produced in substantially larger quantities than at present, if a significant demand existed.' One study' has established its presence in 12 out of 14 zinc blends, in all of 15 aluminum ores, in 4 out of 12 manganese ores, in 35 out of 91 iron ores, and in all of 7 magnetite ores. It occurs as a rule in minute amounts, however, leading to high extraction costs. Recent quotations run from $2.50 to $7.50 per g. During the course of the present investigation, portions of the system Ga-In have been redeter-mined, and the results of this study are presented herein. Thermal and metallographic methods have been employed. Lecoq de Boisbaudran, the discoverer of gallium, conducted the first investigation" on Ga-In alloys in 1885. The temperatures of incipient melting, and of completion of melting, were determined at four alloy compositions. In 1936, Hansen' constructed a eutec-tic-type phase diagram for the system Ga-In, based on . work. The existence of a Ga-In compound was regarded as improbable by Hansen, and subsequent investigations are in agreement. French, Saunders, and Ingle3 conducted a more complete study of the system in 1938, using thermal methods. Their phase diagram is a eutectic type, containing a unique concave-upward liquidus. The solid-solution range of gallium in indium was reported as 9.5 pct by weight, and that of indium as less than 1 pct, at the eutectic temperature. The eutectic composition, determined as being bracketed by the compositions showing a true horizontal at the eutectic temperature (16°C), was reported as 76 5-0.5 pct Ga and 24 i 0.5 pct In. Experimental Procedure The preparation of Ga-In alloys is simplified by the low melting points involved. Various compositions were prepared by melting in pyrex tubes, using a cover of distilled water or parafin to prevent the alloys from wetting the glass wall. In all cases, the melts were homogenized. by stirring. Where possible, both cooling and melting curves were determined. The extensive undercooling of gallium was found to prohibit a satisfactory cooling-curve analysis of gallium-rich alloys, however, and transition points on the gallium side of the eutectic could be determined only by melting curves. The inverse rate method of thermal analysis proved to be most satisfactory and was used to a great extent. Various heating and cooling rates were used, ranging from 0.2" to 5.0°C per min. Low temperature melting analyses were conducted within a constant temperature bath, maintained at about 70 °C. The alloys were solidified (under water or paraffin) within a pyrex tube, using dry ice; the tube was then sealed within a cold Dewar flask, the unit transferred into the constant temperature bath, and periodic temperature readings taken. The high temperature melting-curve determinations and all cooling-curve determinations were made in a vertical tube furnace. At near-eutectic compositions, the furnace was placed within a refrigerated room held at —20°C. Accordingly, the furnace on cooling approached —20°C asymptotically and permitted the determination of those phase transitions occurring below room temperatures. Temperatures were measured with a 30-gage iron-constantan thermocouple, immersed directly in the alloy. To prevent contamination of the melt, the leads and junction were coated with Lucite, applied by painting with a solution of Lucite in ethylene dichloride. Electromotive force measurements were made with a Leeds and Northrup precision potentiometer, type 8662. The couples were calibrated against the boiling point. of water and, at lower temperatures, against a calorimeter thermometer having a Bureau of Standards certificate. The melting points of gallium and indium used in the present investigation were determined as 29.-77° and 156.1°C, in good agreement with previously reported values of 29.78°C° and 156.4"C.' The spec-

Jan 1, 1953

By J. E. Owen

A model of a porous body is presented in which the pore space consists of a system of voids and interconnecting tubes. Relationships between porosity and resistivity formation factor are determined partly by calculation, partly by experiment. Con triction effects characteristic of the model are shown to be sufficient to account for high formation factors. It is shown that constriction may be combined with moderate amounts of tortuosity to give model pore systems exhibiting to a first approximation porosity and resistivity properties similiar to those of natural porous bodies. INTRODUCTION The relationship between the electric resistivity of a fluid-filled porous body and the geometry of its pore space is so complex that the calculation of the resistivity of a natural porous rock is a practical impossibility. Both the resistivity of a body and its porosity are measurable quantities, however, and previous successes at relating them have been reached by an empirical approach. Efforts at obtaining theoretically derived formulae relating them have generally been unsatisfactory. One of the reasons for this may lie in the pore geometry that has been assumed. THE TORTUOSITY CONCEPT A Parameter called the formation factor is useful in dis-cussing the resistivity of a fluid-filled porous body. This parameter is the ratio of the resistivity of a fluid saturated porous body to the resistivity of the saturating fluid. Formation factors are often available from measurements on cores or from electric logs, and many attempts have been made to correlate formation factors and porosities of geological formations. Whenever a successful correlation is found, the engineer working with electrical logs has a useful tool for the determination of porsities of pay section?. One of the more successful formulae applicable to these correlations is the familiar equation empirically obtained by Archie.' which F is the formation factor. $ is the porosity, and rn is an exponent called the cementation factor. When the for- mula applies, the cementation factor usually is found to be between 1.3 and 2.2. The values for formation factors experimentally obtained are often higher than simple pore geometry would lead one to expect. In an effort to account for such high values certain formulae have been derived based on a so-called "tortuosity concept." In deriving these formulae a synthetic porous body is usually assumed in which the solid material is an electrical non-conductor. and in which the pore system consists of three sets of fluid-filled tubes of uniform diameter connecting opposite faces of the body which, for convenience, is considered to be cubical in shape. The three sets of tubes account for the whole of the effective porosity of the body, and usually, it is specified that they do not interconnect. By considering that the pore tubes are not straight but tortuous, their resistance to the flow of electric currents can be made as high as needed to explain high formation factors. Such an explanation has some basis in fact, but it appears that the tortuosity concept is often incorrectly applied when other factors are largely responsible for observed high resistivities. Recently, Wyllie and Spangler have recognized that tortuosity as calculated by conventional formulae has little if any physical significance.' RESISTIVITY AND THE CONSTRICTION CONCEPT Any explanation of high formation factors which depends solely on tortuosity of uniform pore paths necessarily ignores the effect that variations in the cross-sectional area of the conducting paths have on the resistivity of a body. Although, as previously pointed out, the calculation of such paths for an actual body is impossible, it will he shown that a synthetic pore network can be devised which will yield to analysis, and lead to results in agreement with the experimental data represented by Equation (1). The porous body to be considered is assumed to be homogeneous and isotropic or, for present purposes, identical in its characteristics in the three directions parallel to its coordinate axes. It will he assumed to be built of identical unit cubes, each of which contains a single pore network connecting all faces of the unit cube. A unit of such a pore network is shown

Jan 1, 1952

By Leslie W. Austen

ROTARY kiln operators will agree that some of the most severe conditions a refractory must stand occur in the hot zone of a kiln burning Portland cement, dead burn dolomite, magnesite, peri-clase, and similar materials. Frequently the operator is faced with factors beyond his control which drastically shorten the life of refractories. Shutdown due to mechanical failure can be serious if the period is of sufficiently long duration to cause the dropping of coating or the loosening of the lining. A change in slurry can affect the coating and cause ring buildup. A change in type of fuel and its effect upon the flame can cause a shift in location of the hottest zone. Weekend shutdowns or any other interruption can cause the operator trouble and may damage the refractories, since stopping and starting a rotary kiln is certainly more difficult than stopping and starting a motor. Some operators have tried to set an estimate of damage for each shutdown in equivalent days of running time. Conditions affecting the refractory may be roughly grouped in four classes: chemical attack, mechanical stress, thermal shock, and abrasion. Chemical Attack: The drive to obtain maximum production through a kiln demands maximum operating temperatures, temperatures which are limited more by the ringing up or melting of the clinker. This can cause interface temperatures at the junction of coating and refractory which require the use of a basic kiln block to withstand the chemical attack. Chemical changes take place within the refractory itself, especially in chemically bonded or unburned kiln blocks. These changes cause the formation of the ceramic or burned bond. Migrating liquids or fluxes from the kiln charge have an effect within the refractory and lead to mineral or glass formation. The alkalies, sodium and potassium, migrate into the refractory as silicates, chlorides, sulphates or other salts. They may move under capillary action or may be caused to move by volatilization with condensation in the cooler portion. Mechanical Stress: Concentrated stress may be caused by several factors or combinations thereof. I—The rings of refractories must be kept tight and rigid within the kiln, and this alone demands considerable force to hold the blocks in place. So that the force will not be concentrated, the blocks should fit the circle as perfectly as possible, with the faces in contact overall. 2—As the kiln is heated, thermal expansion takes place at the hot end of the kiln block. Since this disturbs the plane face it too can cause a concentrated stress at the two ends of the block, and shearing stress can be set up within the brick itself because of the difference in expansion between the two ends. 3—If a lining becomes loose and moves in the shell very severe stress can be set up, and as the kiln rotates this load changes and gives the effect of repeated loading. Permanent expansion of the refractory can also cause severe loading. 4—Not least important, flexing of the kiln is frequently the cause of concentrated stresses. Thermal Shock: Thermal shock, the result of heating and cooling too rapidly, occurs on starting and stopping or when a large patch of coating drops, exposing the bricks. Again, its destructive effect is often the result of phase change, liquid to solid or the reverse; dense refractories loaded with glass-forming impurities are particularly susceptible. Thermal shock is a. problem with refractories set in the wall or roof of a stationary furnace, and becomes even more serious in a rotary kiln, the tendency to spa11 being magnified with movement and concentration of stress. Uniform rate of feed and loading insures both better coating and a more uniform stress. Abrasion: If the refractories do not take a coating, abrasion can become a most destructive factor. Movement of the lining in shell or movement of loose blocks causes abrasion, which is also most destructive if the refractories do not take a coating. An analysis of the problem of basic lining for the hot zone reveals, therefore, a number of desirable characteristics: high refractoriness, basic chemical reaction, resistance to spalling, good strength at all stages, ability to take coating, true sizing, volume stability, and abrasion resistance. Increased demand

Jan 1, 1953

By Leslie W. Austen

ROTARY kiln operators will agree that some of the most severe conditions a refractory must stand occur in the hot zone of a kiln burning Portland cement, dead burn dolomite, magnesite, peri-clase, and similar materials. Frequently the operator is faced with factors beyond his control which drastically shorten the life of refractories. Shutdown due to mechanical failure can be serious if the period is of sufficiently long duration to cause the dropping of coating or the loosening of the lining. A change in slurry can affect the coating and cause ring buildup. A change in type of fuel and its effect upon the flame can cause a shift in location of the hottest zone. Weekend shutdowns or any other interruption can cause the operator trouble and may damage the refractories, since stopping and starting a rotary kiln is certainly more difficult than stopping and starting a motor. Some operators have tried to set an estimate of damage for each shutdown in equivalent days of running time. Conditions affecting the refractory may be roughly grouped in four classes: chemical attack, mechanical stress, thermal shock, and abrasion. Chemical Attack: The drive to obtain maximum production through a kiln demands maximum operating temperatures, temperatures which are limited more by the ringing up or melting of the clinker. This can cause interface temperatures at the junction of coating and refractory which require the use of a basic kiln block to withstand the chemical attack. Chemical changes take place within the refractory itself, especially in chemically bonded or unburned kiln blocks. These changes cause the formation of the ceramic or burned bond. Migrating liquids or fluxes from the kiln charge have an effect within the refractory and lead to mineral or glass formation. The alkalies, sodium and potassium, migrate into the refractory as silicates, chlorides, sulphates or other salts. They may move under capillary action or may be caused to move by volatilization with condensation in the cooler portion. Mechanical Stress: Concentrated stress may be caused by several factors or combinations thereof. I—The rings of refractories must be kept tight and rigid within the kiln, and this alone demands considerable force to hold the blocks in place. So that the force will not be concentrated, the blocks should fit the circle as perfectly as possible, with the faces in contact overall. 2—As the kiln is heated, thermal expansion takes place at the hot end of the kiln block. Since this disturbs the plane face it too can cause a concentrated stress at the two ends of the block, and shearing stress can be set up within the brick itself because of the difference in expansion between the two ends. 3—If a lining becomes loose and moves in the shell very severe stress can be set up, and as the kiln rotates this load changes and gives the effect of repeated loading. Permanent expansion of the refractory can also cause severe loading. 4—Not least important, flexing of the kiln is frequently the cause of concentrated stresses. Thermal Shock: Thermal shock, the result of heating and cooling too rapidly, occurs on starting and stopping or when a large patch of coating drops, exposing the bricks. Again, its destructive effect is often the result of phase change, liquid to solid or the reverse; dense refractories loaded with glass-forming impurities are particularly susceptible. Thermal shock is a. problem with refractories set in the wall or roof of a stationary furnace, and becomes even more serious in a rotary kiln, the tendency to spa11 being magnified with movement and concentration of stress. Uniform rate of feed and loading insures both better coating and a more uniform stress. Abrasion: If the refractories do not take a coating, abrasion can become a most destructive factor. Movement of the lining in shell or movement of loose blocks causes abrasion, which is also most destructive if the refractories do not take a coating. An analysis of the problem of basic lining for the hot zone reveals, therefore, a number of desirable characteristics: high refractoriness, basic chemical reaction, resistance to spalling, good strength at all stages, ability to take coating, true sizing, volume stability, and abrasion resistance. Increased demand

Jan 1, 1953

By L. A. Panek

During the past three years, USBM has developed an apparatus and technique for direct measurement of existing pressure and change of pressure in mine rock. This relatively simple and inexpensive monitor is reliable for months after being installed. It is used to determine shift of ground pressure created by different sequences of mining, to ascertain the cause of rock failures, and to evaluate the need for artificial support. The technique has been employed to measure pressures in limestone, greywacke, concrete, diabase, and soft iron ore. When rock is subjected to a load it is deformed. Ordinarily this is observed in a mine as the displacement of one point with respect to another—the deflection of the roof, which may be observed as a convergence between roof and floor; or the extrusion of material from the rib, which may be observed as a decrease of the distance between the rib and the post of a timber set. The effect of excessive pressure may be a rockburst if the rock is strong, or it may be squeezing ground if the rock is soft. Some desirable effects of high stress (high in relation to strength) are the caving of roof in a longwall mining operation, the caving of ore in block caving, and the decrease in mechanical energy required to break down the mineral seam in a retreating pillar-robbing operation. In any case, whether the observable effect of rock load is desirable or undesirable, it is a displacement, and depends on the following four factors: 1) The structure—the size and shape of openings, pillars, and linings, the geologic bedding and jointing. 2) The mechanical properties of the rock—prin-cipally the strength, modulus of elasticity, and flow characteristics. 3) The load or applied stress—primary sources are the weight of superincumbent rock, which increases with depth, and unrelieved tectonic stresses; secondary sources are redistributed pressures caused by other nearby openings, especially large mined out zones (rock pressure depends partly on the rock structure created by mining). 4) Duration of load, related to the length of time the opening is exposed. CONTROL OF ROCK DISPLACEMENT Rock displacement can be controlled by control of these four factors. Consider now the means of exercising such control over these factors. Control of the structural features is obviously possible to a great extent, as such control is exercised largely by choosing the method of mining and the methods of natural and artificial support. Rock properties vary, even within a particular mine, but they are controllable only in the limited sense that control may be exercised by choosing the beds or zones to be mined so that rocks with undesirable properties will not occupy critical positions within the rock structure created by mining. Rock pressure is the most complex of the four factors through which ground control can be achieved because it is invisible, difficult to measure, and poorly understood. Rock pressure is controllable only to the extent that control is exercised on the rock structure created by mining. Considering openings within a particular mine, time of exposure varies, and is readily controllable because it is easily measured and easily understood — the longer an opening stands, the greater the likelihood of failure or excessive convergence. Control is exercised by choice of an appropriate sequence of driving openings of different classes, such as haul-ageways and rooms, which are required to remain well supported for different lengths of time under different conditions. Again, control is exercised through the method of mining. All controllable factors can be controlled by proper design of the mining method. The orientation and relative positions of the mine workings and the sequence of their excavation are likely to be much more important to ground control than is the design of artificial support. This implies that the major decisions in regard to ground control are made, knowingly or not, at the time the mining method is chosen. WHY MEASURE ROCK PRESSURE In addition to restrictions on the several factors, control implies the measurement of these factors in some sense, whether only qualitatively by visual observation, or by actual quantitative determination with a measuring instrument. Rock pressure is the most difficult of these factors to measure, largely because of the interaction between the measuring device and the rock. Nevertheless, the quantitative

Jan 1, 1961

By R. J. Lutton

Empirical data from 91 rock excavations have been used to construct a family of slope curves on a chart relating excavation height and inclination. The highest and steepest slopes from each excavation are plotted and connected by a line. The 31 1ines appear to establish one or two geometrical slope fields through which the curves are projected. The curves apparently represent a lumped eflect of many factors too complex to be considered individually. Important factors such as geological structure or ground-water conditions can dominate the picture, but where these factors do not overshadow others, the slope chart should be useful for estimating or checking optimum inclination. The U.S. Army Corps of Engineers (USAE) has been engaged since 1962 with the U.S. Atomic Energy Commission in a joint research program to develop nuclear excavation technology. The Corps is responsible for collecting the requisite data on engineering and construction problems associated with nuclear excavation. The ultimate objective is to develop a capability to use nuclear explosives as a construction tool on public works projects such as harbors and canals. The present study on empirical inclination-height curves for conventionally excavated slopes is an outgrowth of an investigation' sponsored by the USAE Nuclear Cratering Group on applicability of data from existing excavations in establishing guides for engineering judgment in crater excavation. Data were generously contributed by over 50 mining companies among other organizations. Inasmuch as several companies preferred that their contribution be anonymous, all data are treated so here. Some Factors Affecting Stability Many factors affect or are suspected of affecting the stability of excavated slopes, e.g., adverse structural orientation, degree of structural ordering, ground water, climatic conditions, mechanics of excavation, and plan-profile configuration. Such factors are lumped together in developing a slope chart. Obviously geological structure and ground water can overshadow all other factors. Where such is the case, the simplification of slope charts may give an erroneous picture (not conservative enough). A further complication results from the fact that the concept of stability may vary according to the use of the excavations. For example, slope adjustments that would be regarded as failure in a powerhouse excavation, might be tolerable in a mining operation. Mine slopes are continually being modified by further excavation, and instability that would develop in a similar permanent excavation slope over a period of years might not have time to develop in a mine. Inclination vs. Height Charts Charts of inclination vs. height for particular formations or conditions have been used in the past as an aid for designing excavation slopes. Some of these have been based upon empirical analyses of collections of data and experience. MacDonald" and Lane have presented slope tables and charts for excavated and natural slopes in relatively weak sandstone and shale. Coates' has shown an inclination-height chart for slopes in incompetent rock. Most organizations have developed slope design criteria which relate bench widths, heights, and inclinations to use in various rocks. Slope data collected in present studies can be similarly used to develop inclination-height charts; however, a plot of all slopes on such a chart shows very little because of the wide spread of the data. Reasons for the dispersion of points are: (1) many slopes are not carried at optimum inclination, particularly in mines where maintaining roadways and following irregularities of ore bodies are more important than achieving steep inclinations; (2) many engineers and geologists with varied backgrounds and viewpoints have been involved in designing these slopes; and (3) each slope is characterized by a unique combination of factors affecting stability and in turn optimum inclination. Significant curves can sometimes be constructed where the points can be grouped into related clusters from the same excavation. The essence of this approach appears when only two slopes from each excavation are considered—the highest slope and the steepest slope, connected by a line segment (Fig. 1). This reduces considerably the usable data because in many cases only one of these two slopes is available. Nevertheless, the concession seems justified. Construction of Slope Curves Ninety-one pairs of points and connecting line segments were assembled in this study. The line segments were assumed to approximate a geometrical slope field (Fig. 2), and a family of curves has been visually projected through to represent the slope field. It should be

Jan 1, 1971

By J. Downie, H. A. Kendall, G. H. F. Gardner

Some cases of the motion of two miscible fluids in uniform linear models are discussed. There is no bulk flow through the models, and the convection currents are caused solely by density gradients. Horizontal, vertical and inclined positions of the models are treated separately. Analytical formulas for the motion of a narrow transition zone between the fluids are obtained and compared with measurements, and the limitations are discussed. Observations of the ultimate stabilizing effect of diffusion on the motion me compared with the theory. It is shown that fluid motion in an inclined model with a constant rate of bulk flow can be directly deduced from the cases discussed, provided the viscosity ratio is small. INTRODUCTION In some reservoirs, convection currents caused by variations in the fluid densities may be significant. For practical purposes it may be desirable either to minimize the effect, particularly in miscible displacement, or to maximize the effect as in some schemes for the production of heavy crude oils by combustion. Experiments by Craig, et al,l have shown the effectiveness of gravity in reducing the efficiency of a miscible displacement when the mixed zone between the displacing and displaced fluids is narrow compared with the height of the system. The presence of a board transition zone from one fluid to the other has been shown by Perkins, et al,2 to nullify considerably the influence of a difference in density. Indeed, the theoretical discussion given by Perrine3 has shown that when the transition zone exceeds a certain critical length the efficiency approaches 100 per cent. Many additional results 495 have been published on miscible and immiscible displacements which show the importance of density difference and also mobility ratio. The aim of the present paper is to discuss some cases of gravity segregation of miscible fluids in linear models such as are often used in the laboratory to evaluate recovery processes, and to supply additional experimental evidence in support of the conclusions. An attempt is made to calculate simply the magnitude of the fluid motion caused by density pdients and mobility ratios, and to quantify the modifications in the fluid motion produced by the combined effect of molecular and convective dispersion. However, it should be pointed out that the uniform permeability of the models and the use of only a single pair of fluids (whereas several pairs of fluids may be employed in a recovery scheme) severely restricts direct conclusions about natural reservoirs. For conciseness, the motion generated by density gradients is emphasized for the case of no bulk flow through the models. This is not a drastic limitation because, as shown in Appendix A, a fluid motion with a constant bulk flow through any linear model is equivalent to a fluid motion with no bulk flow in the same model inclined at a different angle to the horizontal, provided only that the density is a linear function of the viscosity of the fluid mixture. HORIZONTAL MODEL In this case the model is assumed to have a rectangular cross-section with a horizontal breadth b, vertical height and a great length L. Throughout.the paper it will be assumed that x-axis points along the length, the y-axis along the width and the z-axis along the height. Thus, in the present case the x, y plane is horizontal. Initially, the transition zone between the fluids is assumed to be thin and parallel to the y, z plane. An example of this situation in a practical problem arises in a scheme for using air for cushion gas in the storage of fuel gas in aquifers. 6 With the passage of time the denser fluid sinks and spreads along the bottom of the model while the lighter fluid rises. The interface extends from the top plane of the model to the bottom. Experiments have shown that for practical purposes this interface may be treated as a plane surface at all times. The surface is significantly curved, however, for viscosity ratios greater than about 10. Other complications arise with large viscosity ratios because of the nonlinear dependence of the viscosity on the composition of the mixture, so that the simple results given here are restricted to fairly low viscosity ratios.

By D. L. Sikarskie, B. Paul

A theory is presented for the static penetration of a single rigid wedge into brittle material. The material considered is one which exhibits both crushing and chipping phases in the penetration process. If the wedge angle and three parameters are specified, the theory predicts forces and associated penetrations during both the crushing and chipping phases. For certain ranges of the parameters agreement with the limited experimental data available is promising, except for the initial phase of the penetration process where refinements on the proposed theory are required. The theory also predicts that for certain values of the wedge angle and other known parameters, the chipping process does not occur and penetration is due entirely to crushing. When viewed in detail, the drilling of rock by percussive means involves, among other things, system dynamics, actuating forces, dynamic stress levels and the actual penetration mechanism of the bit into rock. If rock drills are to be significantly improved, a thorough understanding of the entire system would be highly desirable. In this paper, one aspect of the system, namely the mechanism of penetration of the bit, will be studied in detail. It has been found1 that for the velocities encountered in percussive drilling (of the order of 20 fps) a static analysis adequately describes the penetration process, at least for some rocks. Hence, we will only be concerned with describing analytically the static penetration of a single wedge shaped tool. It will further be assumed that the wedge is long enough to permit a two-dimensional analysis. Numerous authors have studied the static penetration of a wedge into rock. The following papers give some of this work and provide additional references to other work. Cheatham2,3 has assumed the rock to behave plastically and has obtained force-penetration equations for both Coulomb-Mohr and parabolic yield criteria. Evans and Murrell4 have studied the penetration of two types of coal and have found equations relating the penetration characteristic (P/dSc) to wedge angle for the various strength coals tested. race' attempted to find a correlation between hardness as determined by indenting the material and other mechanical properties and concluded that the results of such a test are generally inconclusive. His paper however contains a very extensive bibliography on the indenting of many different materials. Gnirk6 also has a fairly complete literature review on the static penetration of rock, of which he makes use in his indexing studies. A wide range of behavior is found in the wedge penetration of different rocks under different external conditions. For example, a rock that is essentially elastic-brittle at standard pressure may become elastic-plastic at a high confining pressure.7 It has also been observed4. " that some rocks will merely be crushed and indented by a wedge, whereas others will crack and form chips, furthermore the existence or non-existence of chips depends in great measure on the geometry of the indentor, type of rock and the depth of penetration. Hence, to attempt a theory which embraces all possible behavior is not practical at present. We, therefore, confine our attention to predicting the force-penetration characteristic and volume removal behavior for a type of rock which exhibits both crushing and chipping phases in the penetration process. Such behavior is characteristic of the harder rocks such as granites and represents a more difficult problem than the behavior of softer rocks such as coal and certain sandstones. In this preliminary study an attempt will be made to describe only the essential features of this complex penetration process. A qualitative description of these essential features is obtained with the aid of Fig. 1.1, where a wedge of vertex angle 28 is shown at some intermediate stage of the penetration process. As the wedge advances, the rock is fragmented (i.e. crushed) in some local region surrounding the wedge, the shape of this region being unknown. Simultaneous with the fragmentation in this local region, essentially elastic stresses are assumed to be building up in the surrounding rock. When a certain penetration level, di+1, is reached, the stresses along some surface are sufficient to cause failure and a chip is thus formed. The process now repeats, i.e., a crushing phase followed by the formation of a chip.

Jan 1, 1965

By M. M. Rao, P. G. Winchell

The growth rates of bainitic plates were measured at 400°C in Fe-Ni-C alloys containing 0.10 atom-fract~on nickel and 0.0012 to 0.0075 atonz-fraction carbon. The growth rates are adequately represented by where xc is nearly the atom fraction carbon of the bulk austenite and PXCy is nearly the carbon atom fractlon in the ferrlte of radiids p In equilibrium with austetzite. The form of the equation is that predicted by a model in which carbon diffusion in austernite controls the gvowth, but the numerical constatnt is two orders of magnitude below that suggested by the model. THE growth of bainitic plates in steel is often assumed to be controlled by the diffusion of carbon away from the advancing plate tip. This hypothesis predicts that the growth rate will increase as the carbon content of the austenite, xCz, is reduced toward the carbon content of the saturated ferrite comprising the plate tip, PxCY The growth rate should vary approximately as (xCg- pxCy)-1. Experimental observation of the growth behavior at low carbon levels should provide a significant test of this model. An alloying element in addition to carbon is required so that low-carbon austenite can be experimentally observed while undergoing bainitic transformation. Nickel was selected. The presence of nickel complicates the interpretation of the data in two ways: First, diffusion of nickel during the transformation would make analysis very difficult. Nickel is assumed immobile during the transformation. Second, nickel affects the solubility of carbon in ferrite and austenite in equilibrium. This effect has been evaluated.' At the completion of our experimental work Goode-now et al.2 published data in essential agreement with the observations to be reported here. Since their discussion is abbreviated and their data are scanty in the region of interest, we believe the present work is of significance. I) THE MODEL OF BAINITIC PLATE GROWTH The rate of lengthening of a plate is assumed to be controlled by the diffusion of carbon from the advancing ferrite-austenite interface into the surrounding austenite. The precipitation of carbides is assumed to be a secondary process. For ease of analysis the carbon-atom ratio,* pxCy, of austenite in equilibrium with ferrite which is convex with minimum curvature radius p, and the carbon-atom ratio, PxCY, of that ferrite in equilibrium with austenite are assumed independent of location on the ferrite-austenite interface. Since these carbon contents vary with the radius of curvature of the ferrite, p, their assumed positional independence must be held as an approximation. The consequences of these assumptions have been developed approximately by zener3 and Hillert,4 and the resulting equation for a platelet has been applied to bainite by Speich and cohen5 and Kaufman, Radcliffe, and Cohen.8 The Zener-Hillert equation* for plates is: The analysis of Hillert is supported by that of Hor-vay and cahn7 which involves no mathematical approximations but does include the assumption that the a/y interface coincides with an isoconcentration line. The solutions of Horvay and Cahn for elliptic paraboloids are replotted in Fig. 1. The shape of the paraboloid is expressed in terms of the ratio of the principal radii of curvature at its tip, A =p1/p2, which is also the ratio of the minor to the major axis of the elliptic cross section. The Zener-Hillert equation for plates is also plotted. The agreement is within a factor of two for (pxyaCr - xyC )/(xyC - PxCaY) between 0.5 and 100. This is the range of interest here and in most other work on bainite. The original form of the Zener-Hillert equation was the form given above with the right-hand side replaced by (pxCya -xCy)/(PxCya). This replacement is not appropriate here. 11) THE EXPERIMENTAL PROCEDURE Alloys were prepared and three kinds of experiments carried out. Continuous-cooling-transformation experiments were carried out on wires by measuring temperature and resistance during continuous cooling. Isothermal-transformation experiments were carried out on wires by measuring electrical resistance as a

Jan 1, 1968

By P. H. Wendland

Calculations are presented for the design of a silicon photodiode in which the incident light beam makes multiple passes between the detector surfaces. Total internal reflection is used for this "light-trapping" effect. By this means, the optical path length can be extended to several millimeters, while the electrode separation remains less than 102 cm, as required for nanosecond response time. Data are presented for a Schottky barrier photodiode constructed on a multiply reflecting silicon base wafer. It is shown that the long-wavelength response is considerably extended in such structures without a corresponding sacrifice in high-speed response. The development of efficient and powerful lasers at 1.06 p has stimulated interest in detectors which operate at this wavelength. In typical silicon photodiodes, for detecting 1.06 p radiation, the requirements for high speed and high sensitivity are mutually exclusive. Since the absorption coefficient is only 25 cm-', a lo-'-cm path length is required to absorb 92 pct of the incident 1.06 p radiation. If the electrode separation is greater than 10 cm, however, the carrier transit time will be greater than 1 nsec. This problem can be solved by allowing the incident light beam to make multiple passes between the electrodes. The optical path length can then be extended to several millimeters, as required for complete absorption, while the electrode separation remains less than 10' cm, as required for nanosecond response time. In a typical photodiode geometry, one ohmic contact and one rectifying contact are formed on the two opposite surfaces of a base wafer, and the wafer thickness determines the electrode separation. The objective of the multiple reflection design is to allow all 1.06 p radiation to enter the detector front surface and to form the back detector surface so that no 1.06 radiation can exit. Total internal reflection at the back detector surface is well-suited for light trapping of 1.06 p radiation because the relatively large dielectric constant of silicon leads to a critical angle of 16.5 deg for total internal reflection. LIGHT TRAPPING It is well-known that, as light passes from one medium such as air into another medium such as glass or silicon, the angle of refraction is always less than the angle of incidence. In the limiting case, where the incident rays approach an angle of 90 deg with the normal, the refracted rays approach a fixed angle +, beyond which no refraction is possible: this is called the critical angle. It follows from Snell's law that where = critical angle, n - index of refraction of air, n' - index of refraction of the medium. Applying the principle of reversibility of light rays, all internal angles of incidence greater than +, will produce total internal reflection and "light trapping". The index of refraction of silicon at 1.06 p is 3.5,' and the critical angle is thus 16.5 deg. Fig. 1 shows these relationships for silicon. This very small critical angle in silicon is significant because all incident angles between 16.5 and 90 deg will produce total internal reflection and "light trapping". This effect can be implemented with a "prismlike" geometry, so that incident light can be introduced into the sample without loss and "trapped". PHOTOSIGNALS A precise knowledge of the absorption coefficient at 1.06 in silicon is of critical importance to the design of fast and efficient silicon photodiodes for 1.06 radiation. Dash and newman2 show a value of 25 cm-l, and our measurements have corroborated this value. Assuming that the collection of photoinduced minority carriers is perfect, the quantum efficiency of a photodiode is dependent only on the absorption coefficient. It then follows from Lambert's law that where QE is the quantum efficiency in pct, a is the absorption coefficient, d is the optical path length, and the reflectivity at the surface is assumed to be completely suppressed by an optical interference layer. Fig. 2 gives the maximum quantum efficiency for 1.06p radiation of a silicon photodiode with optical path length d, using Eq. [2]. The ultimate response time of a fully depleted photodiode to an incident light pulse can be considered to be the arrival times of all photoinduced carriers at the contacts, i.e., the minority carriers at the junction interface and the majority carriers at the oppo-

Jan 1, 1968

By E. H. Snow, B. E. Deal

Instability effects in semicanductor devices have long been attributed to the motion of charges on or within oxide layers on the surface. These effects are of critical importance in metal-insulator-semiconductor MIS) field-effect devices. For this reason, the capacitance-voltage or the conductance-voltage characteristic of these devices can be used as a sensitive detector to study charge transport and polarization effects in the oxide or insulator. A review is given of the results of such studies of thermally grown SiO2 as well as of other insulating films of importance in silicon technology. Three types of effects are distinguished. The first of these is the drift of mobile cations within the dielectric, examples being thermally grown SiO, which is often contaminated with sodium ions, and a variety of other glasses in which the mobile ions are a part of the glass composition. The second effect is a dipole-type polarization which occurs in phosphosilicate glass films obtained by reacting P2O5 with thermally grown SiO2. The third effect involves the transfer of charges between the dielectric and the silicon electrode. This occurs in silicon nitride and other deposited dielectrics. It is concluded that MIS studies have provided a powerful technique joy the study of charge transport and polarization effects in insulating films. The knowledge gained from these studies has led to an understanding of surface effects on conventional transistors and diodes as well as making possible stable MIS transistors. THE metal-insulator-semiconductor field-effect transistor is conceptually the oldest type of active semiconductor device.' The earliest attempts at making this device were frustrated because of high surface state densities at the interface between the semiconductor and the gate insulation.' However, by using a silicon substrate with thermally produced silicon dioxide as the gate insulation, this problem was solved and metal-silicon dioxide-silicon devices with good characteristics were made as early as 1960. 3 Yet it was still over 5 years before these devices became a commercial reality. This delay was largly due, not to surface states, but to stability problems associated with polarization effects within the insulating layer which caused the threshold voltage of the device to drift under temperature and bias treatments. The solution to these problems has not only made possible stable MIS devices, but it has added immensely to our understanding of failure mechanisms in conventional bipolar transistors and has added to the reliability of ali types of planar devices. In this review, we shall first describe the effects of various types of polarization phenomena on MIS device characteristics. Then, since thermally grown SiO, is by far the most important insulator used in these devices, we shall review historically the type of instabilities which have been observed in thermal oxides, the attempts at understanding and eliminating them, and the present status of the problem. We shall then turn our attention to the various deposited insulators which have been used, including lead glasses, phosphosilicate glass, vapor-deposited silicon oxide, and silicon nitride. Interestingly enough, many of these materials show polarization effects which are quite different from those generally observed in thermally grown SiO2. THE EFFECTS OF POLARIZATION PHENOMENA ON MIS CHARACTERISTICS The simplest MIS device and the one which has been most frequently used in the study of polarization effects is the MIS capacitor. Two modifications of this structure with single- and double-layer dielectrics are illustrated in Figs. l(a) and (b), respectively. The capacitance of this structure as a function of voltage applied to the metal gate electrode is plotted in Fig. 2 for the case of an n-type silicon substrate. When the silicon surface is accumulated (positive bias) the measured capacitance is just that of the insulating layer C. When the surface is inverted (negative bias), the capacitance is that of the insulator and a silicon depletion layer in series CoCs/(Co + Cs). Indicated on the horizontal axis of Fig. 2 is the voltage VT at which the silicon surface becomes strongly inverted. This voltage corresponds to the threshold or turn-on

Jan 1, 1969