Search Documents

By R. L. Vaughn, H. R. Cooper

Effective procedures for automating X-ray assaying by use of a computer are described. These methods have advanced the reliability and accuracy of X-ray data, resulting in improved assay measurements for operations and process control. Assaying by X-ray fluorescence has been increasingly used in the mining industry, both for analyzing solid H.R. COOPER and R.L. VAUGHN, Members AIME, are Systems Manager-Process Control and Senior Systems Analyst, respectively, Scientific & Engineering Computer Center, Ken-necott Copper Corp., Salt Lake City, Utah. TP 68B236. Manuscript, Sept. 16, 1968. Discussion of this paper, submitted in duplicate prior to Dec. 15, 1969, will appear in SME Transactions, March 1970, and AIME Transactions, 1970, Vol. 247. samples for smelter, mill, and mine quality control and for continuous stream analysis of mill pulps. The important advantages of X-ray methods in laboratory assaying are reduced analysis cost, and for on-stream analysis, immediate assay information on concentrator streams. Of the numerous reports concerning X-ray analysis in mining operations, one is referenced which can be noted particularly in conjunction with this report.' Continuous X-ray analyzers have been installed at Kennecott and other mining companies for operations guidance and trend information for a number of years. However, complete operational acceptance of X-ray data has been slow, due* to the difficulty in maintaining

Jan 1, 1970

Mr. John J. Porter, Chicago Ill. (communication to the Secretary†):—Mr. Baker's account of his experiences with stock-distribution has been particularly interesting to me, as it

Jan 1, 1905

By John C. Calhoun, Raymond G. Loper

Theoretical calculations for the decline of pressure and the variation of instantaneous producing gas-oil ratio with increased cumulative production have been made for reservoir systems under various producing drawdown values and for various spacings. These variables have been included by a modification of previously Published techniques. This consists essentially of utilizing a radial system flowing gas and oil in a succession of steady state calculations. The results of the calculations do not show any variation in the predicted behavior of the internal gas drive reservoir when either the drawdown or the spacing are changed. The prediction of performance on an internal gas drive reservoir has been presented by various authors. ',"' These predictions are in the form of a pressure decline and an instantaneous producing gas-oil ratio as a function of cumulative oil withdrawal. In making these predictions a basic assumption has been that the saturation of the reservoir is uniform throughout at any time, or that the pressure differential imposed on the reservoir system is zero. Because of these assumptions the validity of such predictions has been questioned. Their application to actual systems, where saturations are not uniform and where producing differentials exist, has been a hazardous procedure. The authors cited have used different methods of approach to the derivation of equations for making their predictions and for arriving at a solution of both equations, but basically the same relationships are utilized in each case. In all instances the instantaneous producing gas-oil ratio is calculated from the equation: In utilizing this equation an average reservoir pressure and an average reservoir saturation must be known, the former to define ß, v, s, and uo/ug and the latter to define kg/ko. Predictions can be made, therefore, only for the conditions of a uniform reservoir saturation and a uniform reservoir pressure. This restriction of requiring a uniform oil saturation can be resolved by applying two-phase steady state flow for the radial system to the reservoir's history. Radial steady state flow is defined when the ratio of phases flowing,

Jan 1, 1949

By John C. Calhoun, Raymond G. Loper

Theoretical calculations for the decline of pressure and the variation of instantaneous producing gas-oil ratio with increased cumulative production have been made for reservoir systems under various producing drawdown values and for various spacings. These variables have been included by a modification of previously Published techniques. This consists essentially of utilizing a radial system flowing gas and oil in a succession of steady state calculations. The results of the calculations do not show any variation in the predicted behavior of the internal gas drive reservoir when either the drawdown or the spacing are changed. The prediction of performance on an internal gas drive reservoir has been presented by various authors. ',"' These predictions are in the form of a pressure decline and an instantaneous producing gas-oil ratio as a function of cumulative oil withdrawal. In making these predictions a basic assumption has been that the saturation of the reservoir is uniform throughout at any time, or that the pressure differential imposed on the reservoir system is zero. Because of these assumptions the validity of such predictions has been questioned. Their application to actual systems, where saturations are not uniform and where producing differentials exist, has been a hazardous procedure. The authors cited have used different methods of approach to the derivation of equations for making their predictions and for arriving at a solution of both equations, but basically the same relationships are utilized in each case. In all instances the instantaneous producing gas-oil ratio is calculated from the equation: In utilizing this equation an average reservoir pressure and an average reservoir saturation must be known, the former to define ß, v, s, and uo/ug and the latter to define kg/ko. Predictions can be made, therefore, only for the conditions of a uniform reservoir saturation and a uniform reservoir pressure. This restriction of requiring a uniform oil saturation can be resolved by applying two-phase steady state flow for the radial system to the reservoir's history. Radial steady state flow is defined when the ratio of phases flowing,

Jan 1, 1949

By R. G. Davies

Above 500°C, where dislocation climb is rate controlling, it is observed that the activation energy for creep is independent of the apblied stress, although it varies from 62 kcal per mol at 15 pct A1 to 73 kcals per mol at 20 pctA1 due to the change in short-range ovdw with composition. Below 500°C, when the nonconservative motion of jogs upon dislocations is thought to be rate controlling, the activation energy is stress dependent. A sudden increase in the creep rate, for which no explanation is availab2e, occurs at the highest stress at all temperatures. CREEP has been the subject of numerous research publications during the past decade with most of the attention focused upon high temperature steady state creep, i.e., above 0.5 TM, where TM is the absolute melting temperature, as it is the most amenable to theory. Many of the proposed theories involve the climb of dislocations out of their slip plane, which allows them to circumvent obstacles or to annihilate dislocations of opposite sign on other slip planes.' For pure metals it is observed that the activation energy for this creep is equal to the activation energy for self diffusion.' The position is much more complicated in alloys as there are many possible microcreep mehanisms, as well as dislocation climb, which could be rate controlling. There has been little systematic investi- gation of alloys, but Weertman (lead and indium base alloys)3 and Lawley, Coll, and Cahn (Fe-A1 alloys)4 suggest a microcreep process as the controlling factor. A Cottrell-Jaswon atmosphere mechanism is proposed for the lead and indium alloys, and a mechanism involving crystallographic order for the Fe-A1 alloys. Lawley et al. showed that the formation of a super-lattice, either of the FeAl or Fe3A1 type. decreased the rate and increased the activation energy for creep. From the measurements of lattice parameters,= resistivity,' magnetic properties,7'8 internal friction: and dilatometry10 the ordered region is thought to commence at -19 at. pct Al. However, a recent investigation by the author indicates a value

Jan 1, 1963

The experience of the Committee in connection with securing papers for the Butte meeting in August has shown that authors should begin early on papers contemplated for the New York meeting next February. So many papers come in at the last moment that it becomes exceedingly difficult to have them all read, indorsed, and delivered to the printer in time. Authors should remember that those serving on this Committee are doing so without pecuniary compensation, and it is hardly fair to expect them to devote so much of their time at the last moment from their private affairs. So many papers have to be sent back for corrections of one kind or another, that it is desirable that they be sent to the Committee at least three months previous to the date of the meeting. The Committee desires to be kept informed of contemplated papers on metallurgical subjects particularly. It is hoped that the Institute Transactions can be made an up-to-date record of metallurgical advance of all kinds. For the present, at least, the activities of the Committee include the refining and preparation of the non-metallic products such as phosphate rock, precious stones, cement, lime, gypsum, building stones, etc. Papers covering any of these subjects are desired. General papers covering geological, mining and metallurgical phases of an ore deposit can be sent either to this Committee or to one of the other Committees on Mining Geology or Mining Methods. C. W. GOODALE,Chairman. DARSIE C. BARD, Secretary.

Jan 8, 1913

By Riki Kobayashi, D. L. Katz

The conditions at which methane and water form solid hydrates have been extended from 4,000 to 11,200 pounds per sq. in. The curve at high pressure had been in doubt because of Villard's report of a critical decomposition temperature. The paper explains the nature of the hydrate curve at high pressure and its dependence on fluid phase behavior. Extropolation of the methane hydrate curve indicates a pressure of 40,000 pounds per sq. in. will cause hydrates to form at a temperature of 100" F. INTRODUCTION The conditions under which natural gases form hydrates when liquid water is present have been studied by several investigators up to pressures of about 4,000 pounds per sq. in. The early work by Villard on hydrates of methane and ethylene13, ethane'', and propane14 included references to critical decomposition temperatures of the hydrates. Villard's work raised the question as to whether the pressure-temperature curve for methane-water hydrate continued above a temperature of 70.7°F. or whether this was the maximum temperature at which hydrates would form. This paper extends the experimental data on methane hydrate to 11,200 pounds per sq. in. A discussion of phase relations in the hydrate region for the hydrocarbons and of the methane-water system at high pressures is included. PREVIOUS WORK Villard suspected that methane hydrates could not be formed at temperatures above 21.5°C. (70.7°F.)la irrespective of the pressure of the liquid-vapor system. However, he was limited in his studies by his apparatus in 1888 to pressures not greater than 300 atmospheres (4,410 pounds per sq. in.) so that there has remained some doubt as to the existence of such a maximum temperature in this region. For the case of ethylene hydrates, however, Villard does specifically state the occurrence of what he terms a "critical temperature for the decomposition of the hydrates"13. Fig. lB shows a pressure-temperature projection of the three-phase locus, hydrate-ethylene rich liquid-vapor (H-L2-V) terminating at the critical locus in which L² and V assume similar properties. Throughout this paper the following phase designations will be used: H for hydrate, Ll for the water rich liquid. L2 for the

Jan 1, 1949

By Riki Kobayashi, D. L. Katz

The conditions at which methane and water form solid hydrates have been extended from 4,000 to 11,200 pounds per sq. in. The curve at high pressure had been in doubt because of Villard's report of a critical decomposition temperature. The paper explains the nature of the hydrate curve at high pressure and its dependence on fluid phase behavior. Extropolation of the methane hydrate curve indicates a pressure of 40,000 pounds per sq. in. will cause hydrates to form at a temperature of 100" F. INTRODUCTION The conditions under which natural gases form hydrates when liquid water is present have been studied by several investigators up to pressures of about 4,000 pounds per sq. in. The early work by Villard on hydrates of methane and ethylene13, ethane'', and propane14 included references to critical decomposition temperatures of the hydrates. Villard's work raised the question as to whether the pressure-temperature curve for methane-water hydrate continued above a temperature of 70.7°F. or whether this was the maximum temperature at which hydrates would form. This paper extends the experimental data on methane hydrate to 11,200 pounds per sq. in. A discussion of phase relations in the hydrate region for the hydrocarbons and of the methane-water system at high pressures is included. PREVIOUS WORK Villard suspected that methane hydrates could not be formed at temperatures above 21.5°C. (70.7°F.)la irrespective of the pressure of the liquid-vapor system. However, he was limited in his studies by his apparatus in 1888 to pressures not greater than 300 atmospheres (4,410 pounds per sq. in.) so that there has remained some doubt as to the existence of such a maximum temperature in this region. For the case of ethylene hydrates, however, Villard does specifically state the occurrence of what he terms a "critical temperature for the decomposition of the hydrates"13. Fig. lB shows a pressure-temperature projection of the three-phase locus, hydrate-ethylene rich liquid-vapor (H-L2-V) terminating at the critical locus in which L² and V assume similar properties. Throughout this paper the following phase designations will be used: H for hydrate, Ll for the water rich liquid. L2 for the

Jan 1, 1949

By R. J. Borg, C. E. Birchenall

Vapor pressures of eight compositions of Mg-Cd have been determined by the Knudsen effusion technique. Measurements are made at several temperatures for each alloy and the results interpolated to a common temperature for the calculation of the activities. It is not possible to obtain values for the partial molar enthalpies because of the narrow temperature range spanned by the measurements. The results show less precision than those obtained for pure Cd which are also reborted. These alloys are extremely susceptible to corrosion which provides a possible explanation for the lack of precision. Nevertheless, the values of the activity are in good agreement with the results of electromotive force measurements previously reported for the Mg-Cd system. MAGNESIUM and cadmium are mutually soluble in all proportions in the solid state and form a single phase which has the hexagonal close-packed structure. X-ray investigations of this system1,2 have revealed a rapid variation in the relative dimensions of the unit cell occurring at approximately 50 at. pct. The rate of change of the c/a ratio with composition increases markedly for alloys containing more than 50 at. pct Cd. The purpose of this investigation was to determine the activity of Cd in the neighborhood of 50 at. pct and to discover whether or not the thermo-dynamic properties undergo a similar rapid change. The activities are calculated from vapor pressure data obtained by the Knudsen effusion method. Unfortunately the precision of the data from which the thermodynamic values are calculated is insufficient to establish unequivocally a departure from the normal trend. The results are reported primarily because they substantiate the only previously reported activity measurements3 and because they include values for pure Cd which in contrast to those of the alloys show excellent precision. The Cd-Mg system is known to contain several ordered phases. However, this investigation was made at temperatures considerably above the transition temperatures, and the alloys may be regarded as nearly random solid solutions. EXPERIMENTAL Apparatus—The effusion cells are machined from solid tantalum rod. The lids consist of threaded, flanged rings of tantalum onto which are welded shallow cups of 4 mil molybdenum foil. The orifice, which is located at its center, is made by grinding off a small dimple produced with a sharpened punch. Such an orifice has been mounted in plastic, sectioned, and the thickness of the orifice edge measured in a microscope. The maximum thickness is 0.00016 in., and the Clausing factor is consequently equated to unity. The area of the orifice is determined by tracing its projection on the screen of the metallograph at a known magnification and measuring the area of the tracing with a planimeter. Excellent temperature control is achieved with a 4 by 18-in. alundum core, gradient-wound, furnace

Jan 1, 1960

By J. F. Martin, J. E. Friedline, L. M. Melnick, G. E. Pellissier

A method has been developed for determining gases in steel in which the gases are extracted by vacuum fusion and analyzed by mass spectrometry. This method is especially applicable for determining small amounts of oxygen and nitrogen. The lower limit of detection for both oxygen and nitrogen is 0.0001 pct with an average deviation of about 0.0001 pct for oxygen and 0.00008 pct for nitrogen. Hydrogen analyses by this method are in good agreement with analyses by vacuum tin-fusion. A complete analysis of all gases extracted from a steel sample can be accomplished in 40 min. It has long been known that certain of the gaseous elements exert particular effects on the physical and mechanical properties of steel. With the recent advent of new experimental steelmaking practices it is now possible to obtain steels with extremely low contents of such elements as oxygen and nitrogen. Present methods for determining these two gases are of insufficient sensitivity as applied to steels made by these newer practices, and therefore a new technique was considered necessary. A review of the literature showed that the conventional vacuum-fusion method1 was the procedure most suitable for modification to obtain a sensitive method for determining the gaseous elements present in steel simultaneously. However, for trace amounts of gases it is impractical to increase the sample size because of the physical limitations of the apparatus and the increased degassing time which would be necessary. Also, when only one of the desired gases is present in trace amounts, quantitative separation by conventional vacuum fusion is almost impossible. Since small analytical mass spectrometers are now commercially available, it was decided to use this instrument in conjunction with conventional vacuum fusion for both qualitative examination and quantitative determinations of the gaseous elements in steel. This technique has already been used with success for the analysis of copper and lanthanum.2"4 APPARATUS A vacuum-fusion unit was first designed and built, Fig. 1. This unit was designed especially for the determination of small amounts of nitrogen and/or oxygen in steel since an acceptable method (vacuum tin-fusion)5 was already in use for the determination of hydrogen in steel. The vacuum system was con- structed of Pyrex glass throughout. Mercury check valves were used in place of stopcocks in the high-vacuum section of the system. Ferromagnetic samples are normally introduced into the vacuum system through a mercury lift. In case of a break in the lift, the mercury is prevented from entering the furnace by a self-locking, mercury safety check valve, Fig. 2. This check valve, developed at our laboratory, is a modified T/S ball-and-socket joint placed in the sample inlet system between the mercury and the furnace. A sudden entry of mercury or air pushes this valve into the closed position, completely blocking the top of the mercury lift. The check valve is held in the closed position by the vacuum of the apparatus. A separate entrance tube also was placed in the system, so that samples which would amalgamate with mercury, or which are nonmagnetic, could be introduced into the furnace. A slightly modified Guldner-Beach furnace is used in this system.8 The optical flat was connected to the system by a ground-glass joint to facilitate cleaning, a quartz delivery tube was substituted for Pyrex, and the crucible was used uncovered. Fusion takes place in a high-purity graphite crucible insulated by — 200-mesh high-purity graphite powder. The crucible set (crucible and funnel) is suspended in a quartz crucible. The furnace envelope and the sample inlet system are exhausted by 2 single-stage, mercury-diffusion pumps in series. The first of these is a high-speed pump, which will work only with a very low backing pressure. This pump is used to exhaust the gases from the furnace as rapidly as possible. The second mercury-diffusion pump is used as a backing pump for the first; the capacity of this pump is much lower. However, it will work efficiently against a back pressure of 3 mm of mercury. To assure positive pumping, a modified Naughton and Uhlig condenser7 is used. With the use of proper valving the second pump is also used as a circulation pump in the catalyst system.

Jan 1, 1959

By L. Kaufman

An equation is derived relating AF a", the difference in free energy between austenite and martensite, to temperature and composition in the iron-chrmnium and iron-chromium -nickel systems. This equation is used to calculate at M,for low-carbon stainless steels. A method for utilizing the energetic equations to calculate M, in stainless steels is suggested. In addition, an explanation is given for the anomalous effect of chromium in forming a y loop concurrent with the lowering of M, in iron-base alloys. It has been demonstrated that the effect of alloying elements on the Ms temperature is strongly dependent upon the influence of the alloying element on the thermodynamic properties of the y and a phases. Notwithstanding the fact that alloying elements affect the Ms temperature by changing the lattice parameters, elastic constants, and inter-facial energies of the a and y phases, their primary role is a thermodynamic one. Thermodynamic analyses have been presented for the iron-carbon'-4 and iron-nickele systems in which equations are derived for Fa and FY, the free energies of the b.c.c. (or b.c.t.) and f.c.c. phases. These equations, which can be obtained with the aid of the equilibrium diagrams, can be used to help predict the Ms vs composition curves. This paper deals with a calculation of the thermodynamic properties necessary to compute AFffLY for iron-chromium alloys. In view of the lack of direct information concerning the activity of chromium in the a and y phases of iron, the regular solution approximationa,g was utilized to bridge the gap. After the necessary equations were derived for the iron-chromium system they were used to calculate AFff— Y for stainless steels by combining them with the previous results obtained for iron-nickel alloys. The results of these calculations offer a means of rationalizing the anomalous behavior of chromium in iron-base alloys. The latter depends on the regular solution approximation in contrast to the explanation offered by zener , which is based upon the "magnetic specific heat effect." THERMODYNAMICS OF THE IRON-CHROMIUM SYSTEM If y represents the atomic fraction of chromium in an iron-chromium alloy, then the free energy of the b.c.c. phase Fa, the f.c.c. phase FY, and the difference in free energy between the f.c.c. and b.c.c. phases of the same composition are given through [3] and Fmy represent the free energies of mixing of the a and y phases; while @era -' and uFea 4Y represent the difference in free energy between f.c.c. and b.c.c. chromium and iron respectively. The values of AEFe0ff4Y in the temperature range between 0" and 900K ( .

Jan 1, 1960

By W. J. Garmoe

Distinct areas of mineralized and non-mineralized brecciated rock are found in the Ontario Unit of the United Park City Mines. These breccias contain an appreciable fraction of the present ore reserves and probably will become even more important as deeper levels are developed. Three ages of breccia have been recognized in the Ontario Mine to date: ( 1) Pre-mineral breccia that may or may not be mineralized, (2) early post-mineral breccia with mineralized fragments, and ( 3) late post-mineral breccia that is younger than and cuts through ( 1) and (2) above. These age ,relations are known from detailed underground mapping. Locally, both pre-mineral and post-mineral breccias have sufficient lead-zinc-silver content, either as matrix or fragmental material, to constitute ore. Ore-grade breccias occur throughout a vertical range of over 1000 feet and are found in both the Humbug and Deseret formations from above the Ontario 1300 level to below the Ontario 2200 level. The mineralized breccias, which are difficult to distinguish on mine maps, are not well delineated by underground workings, and their horizontal outline and vertical extent are not known beyond the immediate mining areas. The greatest amount of truly mineralized breccia is exposed on the 2200 level.

Jan 1, 1968

By S. Pearson, L. Rotherham

Measurements have been made of the variation of internal friction with temperature for spectroscopically pure silver, and for o series of solid solutions of silver with cadmium, indium, and tin, using a Ke-type torsion pendulum apparatus. Some experiments have also been made to investigate the effect of nonmetallic impurity on grain boundary relaxation in silver. The effect of the alloying elements is to increase the grain boundary viscosity, and to raise the activation energy for grain boundary relaxation from 22,000 cal per mol for pure silver to 43,000 cal per mol for the solid solutions; the same value being obtained, within the limits of experimental error, for all the alloying elements and solute concentrations investigated. Results of the experiments show exactly the same trend as those obtained previously for a similar series of copper solid solutions. They are in agreement with the general theory of grain boundary relaxation developed by Zener and Ke, but do not seem to be in agreement with either of the mechanisms so far put forward to explain grain boundary slip. EXPERIMENTS described in this paper were made as part of an investigation into the effect of alloying elements in substitutional solid solution on grain boundary viscosity, measured by the internal friction method developed by Ke.&apos; The apparatus and experimental procedure have already been described.&apos; Preliminary work was done with spectro-scopically pure silver. Experimental Work Materials Used—All the materials were supplied in the form of % mm diam wire by Messrs. Johnson Matthey and Co. Ltd. The silver and all the alloying elements used were spectroscopically pure. Composition of the alloys is given in Table I. They were all melted and chill cast under the same conditions as the copper alloys. The ingots were rolled to 1/4 in. sq rod with intermediate annealing at 600°C in an atmosphere of cracked ammonia when necessary. These rods were then drawn to wire by standard procedure with intermediate bright anneals at suitable steps in the drawing operations. Specimens from the annealed wires were sectioned longitudinally. mechanically polished, and examined under the microscope. They were then etched and reexamined. Grain sizes were reasonably uniform for all wires, and no traces of a second constituent were observed. The wires were weighed before and after testing so as to measure any loss of the alloying element due to evaporation. The loss was less than the accuracy of measurement—< 0.1 pct of the total weight— for the Ag-In and Ag-Sn alloys, but appreciable loss of weight occurred in the case of the Ag-Cd alloys. The percentage loss is given in Table 11. X-ray diffraction photographs were taken of all thc wires after test to determine whether there was any marked preferred orientation. Results of the examination are given in Table 111. Experiments on Spectroscopically Pure Silver-— Experiments were made on specimens in which the grain diameter was small relative to that of the wire, and on one specimen in which the grains were grown by the usual strain-anneal method to extend across the full diameter of the wire. These large grained specimens will be referred to as single crystal specimens. The wires were annealed for 11/2 hr at 700°C, and experiments were made at two frequencies of vibration, 1.5 and 0.4 vibrations per sec. The variation of internal friction with temperature is shown in Fig. 1. Experiments on a second small grain specimen gave approximately the same results as for the first wire. It will be seen from Fig. 1 that the peak in the curve for the small grained wire, which is presum-

Jan 1, 1957

By B. K. Larkin, H. R. Bailey

A model Of heat flow in an underground combustion process is studied. This model includes convection effects and thus is more general than previous studies which considered conduction as the only mechanistn for heat transfer. Both linear (tube run) and radial (field application) geometries are considered. The effects of ignition heaters, vertical losses and finite source width are considered for the linear case. The results are in the form of equations and are presented in graphical form for a number of cases. Convection effects increase frontal temperatures about 25 per cent over those computed for conductive transfer for typical field operating conditions. This increase in temperature is a result of heat being transported toward the front by the injection gas. Even greater temperature increases are realized as the per cent oxygen in the injection gas decreases. It is well known that compression costs are of considerable importance in estimating the economic feasibility of underground combustion. By assuming an ignition temperature for the combustion fuel, predictions of limiting conditions on fuel density and injection rate necessary to sustain the combustion zone are made. For typical field conditions, at least 0.75 lb/cu ft of fuel are needed with air as the injection gas. If the injection gas is 10 per cent oxygen and 90 per cent nitrogen, this figure is 0.69 lb/cu ft. INTRODUCTION The possibility of increasing oil recovery by underground combustion has been considered for many years. Recent field tests indicate that the underground combustion process is technically feasible. The economic feasibility depends to a large extent on the amount of air which must be injected to sustain combustion. Prediction of the success of employing underground combustion in a particular reservoir must be based on existing field tests, laboratory tube runs"&apos; and solutions of mathematical or analog models. Vogel and Krueger5 devised an electrical analog of the heat transfer problem in an underground combus-ion process. They considered the problem of heat conduction from a cylindrical source with increasing radius assuming no vertical losses. Bailey and Larkin and Ramey solved the corresponding problem including vertical losses by considering a mathematical model of the conduction process. Ref. 9, 10 and 11 also present information related to the subject study. The present paper generalizes these results to include both conduction and convection mechanisms for heat flow. A model of the conduction-convection process is described. The partial differential equations governing this model are written and solved for a number of cases of interest. The formulas are evaluated and the results are presented graphically. A CONDUCTION-CONVECTION MODEL Description of the Model A number of models are employed in this paper, each of which approximate a particular physical situation depending on the geometry, i.e., linear or radial flow, the type of vertical heat losses, the nature of the source, etc. The geometry for both the linear and radial case is shown in Fig. 1. A particular linear model is described in this section, and variations from this model are described in later sections of the paper. Consider a homogeneous, semi-infinite porous prism with lateral sides insulated as shown in Fig. 1. The left end, x = 0, is heated for a time t, with a plane heater supplying a constant heat flux and during this time; gas (nitrogen) is passed through the pores of the prism at a constant mass flux. At time t, the heater is turned off and air is injected in place of the nitrogen. This causes combustion to take place, and a uniform source of width 1 is formed which moves at a constant velocity v, to the right. The frontal velocity v, is related to the gas velocity v,, through stoichiometric considerations

By B. W. Christ, P. M. Giles

Mossbauer effect measurewents have been made on I-mil-thick foils of commercial 1 wt pct C steel and Fe-2 wt pct C alloy. The experimental method required about 3 to 5 vol pct of a phase in the nzultiphase steel sample for detection. Room-temperature Md&apos;ssbauer patterns obtained on austenitized atid quenched samples exhibit fifteen, and possibly twenty-one, lines. A sharp parama&tetic singlet and a quadrupole doublet, poorly resolued from the singlet, are attributed to austenite. Remaining lines are due to tnartensite. Accurate evaluation of austenite line paranzeters is not feasible if significant amounts of other phases such as carbides or martensite occur simultaneously with austenite. This is demonstrated by comparison of hyperfine interactions determined for austenite in multiphase high-carbon samples with those reported for Fe-C austenite in a nearly 100 pct austenitic sanple. Lines from carbides are incompletely resolced from austenite lines, as demonstrated by comparison of austenite line positions with carbide line positions calculated frow published values of hyperfine interactions. One martensite line overlaps an austenite line in the pattern for commercial 1 wt pct C steel. Results of this study suggest that the usefulness of e M6ssbauer pectroscopy for quantitatizle analysis of austenite in bulk samples of quenched and tempered high-carbon steels is restricted by poor resolution. Use of Mossbauer spectroscopy for phase identification and for evaluation of atomic and electronic structures appears quite feasible, however, The Mijssbauer effect has been widely discssed,&apos;-and e Mossbauer effect measurements have been reported on materials of metallurgical interest.7"20 In particular, it has been proposed that e Mossbauer patterns of commercial steels and laboratory-made Fe-C alloys, in the quenched condition, are composed of lines originating in two phases, Fe-C austenite and Fe-C martenite.-&apos; Evidence accumulated in this study demonstrates that three absorption lines found in the central region of the e Mossbauer pattern obtained on quenched steels are attributable to retained austenite. This interpretation is supported by parallel decreases in the intensity of these three lines caused by subambient cooling of commercial 1 wt pct C steel samples after water quenching to room temperature. A second result of this study is to clarify effects of line resolution and sensitivity in the Mossbauer patterns of multiphase steels on the accu- rate determination of austenite line parameters. Experimental line widths (full width at half height) are generally 1.5 to 3 times larger than the natural line width of 0.19 mm per sec. At least two lines, and sometimes more, from a single phase such as cementite (Fe3C), other carbides, martensite, and austenite fall in the energy band, i0.85 mm per sec. hhis band width is employed simply for convenient reference. It represents approximately the energy interval between the + 112 to 112 transitions in ferrite and is expressed as the velocity needed to Doppler shift a 14.4kev 7 ray to the aforementioned ferrite energy levels. This energy band is referred to as "the central region of the Mossbauer atttern:: in this paper. Hence, due to the large number of lines from different phases in a multiphase steel falling in a relatively narrow energy band, absorption lines from the different phases may overlap. Analysis of available data, presented below, indicates that this occurs to a significant extent for phases which commonly occur in quenched and quenched and tempered high-carbon steels. One consequence of limited resolution in the Mdssbauer patterns from multiphase steels is difficulty in accurate determination of such line parameters as position, width, and intensity. In fact, it appears that quantitative analysis for retained austenite in quenched and tempered high-carbon steels is not practical with the present experimental method. Line resolution is influenced to some extent by sensitivity. We point out below that atom or volume fractions of less than about 3 to 5 pct are not detected by the present experimental method. Thus, the presence of a multiplicity of phases does not always lead to impaired resolution. Finally, we report in this paper room-temperature MGssbauer parameters determined for austenite in a freshly quenched, commercial 1 wt pct C steel and in a freshly quenched laboratory heat of an Fe-2 wt pct C alloy. These parameters are compared with others reported in the literature. Three types of hyperfine interactions are detectable in a Mossbauer effect measurement: isomer shift, quadrupole interaction, and magnetic dipole interaction. These interactions are evidenced by one, two, and six line patterns, respectively.&apos;-4 More than one type of interaction has been reported in certain metallurgical phases thus far studied by this method. Isomer shift is the experimentally measured displacement of line position from some arbitrarily defined reference position. In the case of a multiline pattern, isomer shift is given by the displacement of the centroid (center of gravity) of that pattern from the reference position. All isomer shifts measured at finite temperatures contain a second-order Doppler effect characteristic of that temperature. The isomer shift is related to the total s electron density at the nucleus, becoming more negative with increasing s electron density. The first-order quadrupole effect arises from the interaction between the nuclear quadrupole moment and any axially symmetric electric field gradient in

Jan 1, 1969

By A. S. Odeh

This paper analyzes the effect of limited entry to flow at the wellbore on the steady-state productivity of a well. Wells that have been opened to flow along a fraction of their productive interval are termed wells with limited entry. Previous work treated the cases of a partially penetrating well, a well producing from the central portion of the productive interval and a well in which several intervals equally spaced were open to flow. In this paper the open interval can be located anywhere within the productive interval. Thus, in a sense, it generalizes previous work. The finite cosine transform was used to arrive at a solution for steady-state flow of a slightly compressible fluid. The solution was programmed for a CDC 1604 computer. Numerical vaIues for rd = 660 ft, r, = 1/4 ft, and range of sand thickness of 20 to 200 ft are presented in graphical form. The effect of rd and r, values on the result is shown in a table. The correct calculation of skin and damage ratio in the presence of limited entry to flow is explained and illustrated by examples. Moreover, the paper shows how to calculate the net decrease in productivity due to the combined effect of limited entry and perforations. INTRODUCTION In some wells only a fraction of the productive interval is open to flow. Location of this fraction is usually dictated by formation characteristics and reservoir behavior. For instance, if a gas cap exists, the open interval is located away from the gas-oil contact to prevent any possible gas coning. Wells that intentionally have been opened to flow along a fraction of their productive formation are tened wells with limited entry. Obviously, unintentional completions of this type also exist. Limited entry to flow decreases well productivity. Magnitude of the loss depends on the fraction of the formation open to flow, on the thickness of the sand, on the location of the open interval and on the ratio of rd /r, where r , is well radius and rd is the drainage radius of the well. The use of pressure buildup data on producing wells to calculate the condition of the formation around the wellbore is an accepted practice. van Everdingenl and Hurst2 introduced the concept of she skin factor s considered to be due to a thin layer of different permeability immediately around the wellbore. These authors dealt with the case of a well of complete radial geometry, i.e., a well with open-hole completion that completely penetrates the formation. The presence of a low-permeability skin results in a loss of productivity, as does limited entry. Therefore, if pressure buildup data obtained on a well with limited entry are used to establish the presence or absence of skin (i.e., formation damage), and a correction is not made for this loss of productivity, the calculations would result in an erroneous skin value. They might indicate the presence of formation damage when in reality there is none, or they might indicate a value larger than the true value. This could lead to an incorrect basis for planning remedial measures. Muskat3 studied the problem of partially penetrating wells for the case of incompressible flow. He presented equations and figures which allow estimation of loss in productivity. Brons and Marting,4 using equations based on Nisle&apos;s work,5 studied the loss of productivity for three cases. The first was for a partially penetrating well; the second was for a well producing from only the central portion of a productive interval; and the third was for a well in which several intervals equally spaced were open to flow. Their work was for steady-state depletion-type reservoirs wherein the well radius of drainage is established and the fluid is considered to be slightly compressible. Considered in this paper is the problem of wells with limited entry in which the open intervals are located anywhere within the productive sand. The finite cosine transform is used to arrive at a

Jan 1, 1969

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system&apos;s behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell&apos; have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine&apos; have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.&apos; A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By A. W. Schillinger

Drilling was not entirely satisfactory in the search for native copper in the Osceola amygdaloid footwall, for oreshoots mined were more continuous than drilling indicated. The Geophysics Dept of Michigan Technological University developed a portable induced polarization drill hole probe with a 3-electrode configuration and a sampling radius of 5 ft ±. A 20 ma current of alternately reversed polarity is pulsed at 3.5 sec intervals and readings are taken every 2 ft. Potential is read during the "on" cycle for resistivity determinations. During the "off" cycle, after 10 ms delay, the 1P potential is sampled for 10 ms and recorded. Resistivity (ohm-ft) and "S" values (mv/v) are plotted on semi-log paper. A lithologic log is plotted based on resistivity and anomalous S values noted. The low resistivity amygdaloid contains sills of high resistivity basalt. Verified by core logs, contacts between amygdaloid and basalt could be picked to ± 0.5 ft from resistivity alone. Lithologic logs of long-steel holes, based on resistivity, then replaced sludge logging. An S value of 30 mv/v was established as the lower limit of diagnostic 1P effects in mineralized amygdaloid and 50 mv/v as a threshold value for ore. Of 914 footwall zones drilled, 315 were probed yielding 203 anomalies. Success increased from 34% for drilling to 72% for drilling and probing. As checked by mining, 46 probe results were 78% correct. Failures were due to the instrument not detecting copper or to the misinterpretation of an anomaly. No false anomalies were found. During the past ten years, the Calumet Div. of Calumet &Hecla, Inc has conducted extensive exploration and mining in the footwall of the Osceola amygdaloid lode in the native copper district in Michigan&apos;s Upper Peninsula. In an attempt to delineate more clearly the extent and nature of high-grade footwall copper mineralization, 876 diamond drill and long steel holes aggregating slightly more than 47,000 ft have thus far been drilled from the accessible hanging-wall workings in the Osceola No. 6, 13, 14 and 15 shafts. Early in this program it became obvious, through inspection of footwall workings in progress, that due to the erratic and bunchy occurrence of the native copper in the footwall zone, the exploration drilling was failing to outline completely the extent of the oreshoots. Early drilling re-

Jan 1, 1964

By J. F. Wolff

EXPLORATION of Lake Superior iron ores is done principally by drilling. The soft iron ores are churn drilled and the harder ores are diamond drilled. In exploratory work in the harder formations, where the diamond drill is used, the sample comprises both core and cuttings. The care given in handling and preserving the core depends on the im-portance that may be attached to certain characteris-tic strata in different formations. The following instructions for the care of core samples are issued by the Oglebay Norton Co. for diamond drill work done for them on the Gogebic Range of Michigan: The driller shall be careful upon removing the core from the core-barrel to keep the core in exact order. Not only should the pieces of core be kept in exact relation to each other, but care should be used not to turn them end for end. Core should be placed in the core-box from left to right in each row. Each drill should be furnished with small, hinged cover core-boxes, and each drill should be well supplied with empty boxes. Core remaining at the drill when the drillers are not working should be locked up. The practice of picking sample pieces of core and carrying them away from the drill to be examined in better light, or shown to someone else, is unnecessary and prohibited. The reason for this is that these pieces are seldom returned to the core-boxes and valuable information or "markers" may be lost. If the inspector suspects that a critical point is reached in the drilling and he cannot determine at the drill, he should remove the entire core-box and see that the drill is immediately furnished with an empty core-box. It is customary to split the core longitudinally so that one-half may be ground to a pulp for analysis and the other half may be kept for future inspection. The core should be filed in permanent core-boxes, properly labelled, if the core is to be kept in exact order as taken from the core-barrel. Otherwise, it is often filed in tin boxes inn the way that sludge samples are filed. To obtain the average analysis, the core and cutting samples obtained in the same run are analyzed separately and the analyses are combined in the proper ratio.

Jan 9, 1922

By D. Y. F. Lai, R. J. Borg

Tracer diffusion of Fe59 has been measured in the a-stabilized Fe-1.8 at. pet V alloy from 700° to 1500°C. The activation energies are obtained in both the presence and absence of magnetic order. Furthermore, it is established that diffusion in the alloy is identical to that in pure iron and consequently the values of Do and Q accurately represent the temperature dependence for self-diffusion. The purpose of this investigation is to obtain an accurate estimate of the temperature dependence for self-diffusion in bee iron in both the presence and absence of magnetic order, and, in so doing, to establish the temperature range of the magnetic effect.&apos;" As the temperature interval suitable for diffusion measurements is severely limited in both bee phases of pure iron because of the intervening fee ? phase, the experiments were performed on an a-stabilized alloy containing 1.8 at. pet V. This alloy is bee over the entire range from room temperature to the melting point. Although there have been several independent investigations of self-diffusion of iron in a, iron,1, 3-6 there still exists considerable disagreement regarding the values of Do and Q for the paramagnetic region. The two systematic studies of diffusion in 6 iron6, 7 previously reported are also only in fair agreement; but in view of the extremely small temperature range available for diffusion studies, i.e., 1390o to 1535oC, this is not surprising. It is comparatively easy to obtain accurate values of Do and Q for the a-stabilized alloy inasmuch as measurements can be made over the entire temperature range -700o to 1500°C. However, in order to assume that these same values apply to pure iron requires careful comparison of the data in the a, and 6 regions in both the alloy and pure iron. We have made several measurements in the appropriate temperature ranges and are unable to establish any systematic difference between the diffusion coefficients of iron in pure iron and in the alloy. We therefore conclude that the values obtained for the alloy are truly applicable to pure iron; the complete evidence favoring this conclusion will be discussed later in this paper. EXPERIMENTAL The experimental methods will be given here only in barest detail since they have been thoroughly de- scribed elsewhere.l, 7 The alloy was prepared by induction melting and chill casting under argon. Diffusion samples were machined from the ingot and annealed in hydrogen for several days at 900°C to give an average grain diameter of 1 to 2 mm. The penetration profiles of the tracer were established by a sectioning technique, the residual activity being counted after the removal of each section. The tracer used is Fe59 which emits two high-energy ? rays of 1.098 and 1.289 Mev, respectively; these were detected by a ? scintillation counter equipped with a pulse-height analyzer. For the measurements in the temperature range -700o to 1130°C the samples were vapor-plated with Fe59, encapsulated in quartz under vacuo, and annealed in resistance-heated furnaces which are controlled to ±1°C. The specimens diffused at higher temperatures are prepared as edge-welded couples, the two halves being separated by a thin washer of the alloy to prevent sintering. The diffusion anneal is then carried out by inductive heating under a dynamic vacuum. The temperature is monitored pyrometrically. RESULTS The diffusion coefficients are obtained from the penetration profiles in the usual way using the error-function complement relationship. The results over the entire temperature range are shown in Fig. 1. In the linear region, 900o < t > 1500°C the least mean squares (lms) values of the diffusion coefficients are given by D = 1.39 exp[-(56.5 ± 1) x 103/Rt] cm2/sec [l] The average departure of the measured diffusion coefficients from the values given by Eq. [1]In order to determine whether or not the slope is truly constant over the entire range from 900" to 1500°C, the data are arbitrarily divided into two groups, the first containing values between -900" and 1133°C and the second between -1133o and 1500°C. The lms values for the two groups are given by Eqs. [2] and [3]: D = 0.519 exp[-55.7 x 103 /RT](900° to 1133°C) cm2/sec [2] D = 1.45 exp[-56.7x103/RT] (1133° to 1500=C) cm2/sec [3] Thus, there is no significant difference between the high- and low-temperature segments of the linear region. This not only assures us of the consistency of the values obtained by induction heating as compared to those obtained from the resistance-heated

Jan 1, 1965