Search Documents

By W. M. McKewan, R. G. Olsson

In order to improve the understanding of the iron oxide reduction mechanism, the rate of equal molar diffusion of hydrogen and water vapor was measured through disks of porous iron and alumina. The measuretneizts were made at atmospheric pressure at 800°, 900- , and 1000°C. It is shown tint the reduction of dense hematite pellets with hydrogen at these temperatures and at atmospheric pressure is to a great extent limited by gaseous diffusion between the bulk-gas phase and the iron-wustile interface. THE investigation of the reduction of iron oxides with hydrogen has been the subject of numerous papers. Although the effects of temperature, pressure, and gas composition on the rate of iron oxide reduction have been studied, the actual mechanism of the reduction process has not been definitely determined. As a means of improving the understanding of the mechanism, the rate of equal molar counterdiffusion of hydrogen and water vapor was measured through porous iron and porous alumina disks at reduction temperatures. The alumina was used as a basis of compari-son because its pore structure is homogeneous and independent of temperature. The porous iron was prepared by the reduction of dense hematite disks with hydrogen at the temperature of the subsequent diffusion experiment. The hematite disks were made from pellets that were identical to those used by McKewan1,2 in his study of iron oxide reduction with hydrogen. In these investigations the weight loss and rate of movement of the various interfaces were measured as a function of the reduction temperature and gas composition. With these reduction data and the presently reported diffusion data it is possible to determine the importance of H2-H2O transport between the bulk-gas phase and the Fe/FeO interface to the reduction of dense hematite. APPARATUS AND PROCEDURE A porous iron or alumina sample was held in the base plate of a sealed canister containing porous granular wustite, as illustrated schematically in Fig. 1. The canister was suspended from an automatic recording balance in a vertical tube furnace through which flowed hydrogen and water-vapor mixtures of known composition at atmospheric pressure (0.96 atm). Hydrogen from the gas stream diffused through the porous sample, reducing the wustite in the interior of the canister. The water vapor that was formed diffused back to the gas stream. The rate of diffusion was determined from the recorded weightloss data. The furnace, automatic balance, and gas-metering equipment have been described previously.2 The main body of the canister and the tubing extension were made of stainless steel while the base plate was made of high-purity iron. The main body of the canister was % in. in diam. It held approximately 6 g of porous granular wustite that was prepared by reducing 0.1- to 0.2-cm-diam specular hematite particles at 900°C with an H2-H2O mixture of a composition near the iron-wustite equilibrium value. The upper end of the 1/8-in.-diam tubing extension was sealed with a segment of plastic tubing and a tubing clamp. The temperature at the top of the furnace was sufficiently high so that no water vapor condensed inside the extension. In the experiments with porous iron, iron oxide samples were prepared from reagent-grade hematite that was pressed and sintered in a manner previously described.' The resulting pellets had a density of 5.15 ± 0.05 g per cu cm. Each pellet was centerless-ground

Jan 1, 1967

By Will Mitchell, Ford F. Miskell, C. L. Sollenberger

THE taconites in general are hard, tough ores, difficult to grind. Liberation of iron mineral constituents usually is accomplished by grinding the ore through at least 100 mesh, and often it has been found necessary to grind substantially through 325 mesh to achieve satisfactory recovery and grade in the concentration process. Because of the fineness of grind required and the enormous tonnages of material contemplated for treatment in the future, costs resulting from grinding media wear and power consumption, together with capital investment required for comminution, approach astronomical figures. Economy in any one of these elements per ton of material ground could very well reflect a considerable saving to the ferrous industry in yearly costs. With this in mind, the Research Laboratories of Allis-Chalmers Manufacturing Co. have launched a program to discover means of effecting this economy. The initial phase of the work as described here deals with a heat treatment of crushed raw ore, followed by thermal shock in cooling, for the preparation of rod mill feed. Several investigators 1-5 have noted an improvement in the grindabilities of ores treated in this manner. One investigator2 subjected low grade iron ores to heat treatment in an electric furnace. Basing his conclusions on screen analyses of crushed products, he observed that treated ore was more easily crushed than the untreated. He reported that observation of the treated and untreated ore through a microscope revealed cracks following the grain boundaries in treated ore, whereas no cracks were present in untreated ore. However, few if any have quantified the improvement in terms of total hp hours saved, reduction of wear of grinding media, or reduction of capital costs of grinding equipment involved. By means of Bond rod mill and ball mill grindabilities6 and by comparative wear tests, the conditions maintained during this investigation have been definitely evaluated. Heat treatment experiments have been made on a batch scale in a muffle furnace, followed by continuous scale experiments in a rotary kiln under various conditions of temperature, atmosphere, retention time, and quenching to determine the combination that would give the greatest improvement in the grindability of a taconite without affecting adversely the magnetic susceptibility of the magnetite. From these data, relative costs have been calculated for grinding both untreated and heat treated ore on a basis of plant capacity of 120 tons per hour. Heat requirements for a rotary kiln, as well as the kiln capacity required for the treatment, have been estimated. For these tests 50 tons of ore were obtained from the Aurora mine of the Erie Mining Co., Hibbing, Minn. The ore was reported' to have been selected from the magnetic portion of the lower chert horizons. An average chemical analysis of this portion is shown in Table I and an approximate mineralogical analysis is shown in Table II. No attempt was made to differentiate between the various silicates because the compositions vary widely on the different horizons and depend on the degree of oxidation. Minnesotite, however, is the predominant silicate. The siderite grains were dispersed throughout the ore but were essentially associated with the rock-forming minerals. In a series of preliminary tests in which batches of -6 mesh ore were heat-treated in a muffle fur-

Jan 1, 1952

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By A. M. Gaudin, C. S. Chang

Adsorption was measured for barium ion and laurate radical, using radioactively marked agents, over wide range of concentrations. Laurate adsorbed in absence of barium fails to float. With barium, flotation-effective laurate coating needs not exceed 5 pct of monolayer. Barium adsorption then is several fold larger. No barium laurate coating is obtained. IN general, fatty acids or their alkali salts do not cause clean quartz to float. However, the presence of multivalent ions, such as the alkaline-earth metal ions and heavy-metal ions, can activate quartz and thus make it respond to soap flotation. The mechanism of activation and soap flotation of quartz has been studied.'-' The problem was attacked by making flotation tests rather than by careful adsorption measurements on suspensions of ground quartz in solutions containing the activator and the collector. Attempts have been made to determine the adsorption of barium ion and calcium ion on quartz. The results have only qualitative value as the analytical methods were not adequately precise. Artificial radioactivity and progress in radioactivity detection have introduced a new analytical tool, which can have a sensitivity several orders of magnitude greater than that of other methods."-" In the experimental work described in this paper the radioactive tracers, barium 140 and carbon 14, were employed. Lauric acid was synthesized, having a radioactive carbon atom and was used with Ba" to study the adsorption of barium and laurate ions from aqueous solution on ground quartz. Flotation tests, using the vacuum flotation technique,' were run in conjunction with the adsorption experiments. The object of this investigation was to understand better the reaction mechanism of barium activation and soap flotation of quartz. Methods and Materials Quartz Sample Preparation: Lumps of hand-picked quartz, washed free of stains, were reduced in size to —20 mesh in laboratory crushers. The product was boiled in concentrated hydrochloric acid to remove the abraded iron. The cleaned product then was ground dry in a laboratory porcelain mill, using siliceous pebbles, to pass through a 200 mesh screen. The ground product was deslimed by sedimentation at 37 microns. This —200 mesh +400 mesh material was subjected again to boiling with concentrated hydrochloric acid, followed by washing with distilled water and further desliming. The cycle was repeated until no trace of iron could be detected in the wash liquor. During the final washing, the pH of the wash water was checked to insure the complete removal of the acid from the mineral. The quartz then was rinsed with conductivity water. The resulting product was dried in a drying oven and stored in mason jars. The impurities (spectrographic) were as follows: Al, 0.01 to 0.001 pct; Li, Mg, Mn, Na, each less than 0.001 pct; Fe and Cu not detected. The specific surface of the quartz sample was determined by the krypton gas adsorption method.' The average of two measurements was 881 sq cm per g. The calculated specific surface for 37/74-micron quartz particles, using a surface factor of 1.96; is 900 sq cm per g. This is in good agreement with the measured value of 881 sq cm per g. Beta-Counting Method: The radioactive isotope of barium used in the experimental work was barium 140. It has a sufficiently long half-life, 12.8 days, emits high-energy beta particles (up to 1.0

Jan 1, 1953

By J. R. Hornaday, A. E. Morris, N. A. Parlee, D. C. Carmichael

Rates of absorption and effusion of hydrogen in solid iron were measurede by a Sieverts type of apparatus. With clean a iron these rates are diffusion controlled down to 420°C and are represented by the equation D = 6.70 x e-'100' crn2 sec-I Certain surface impurities can affect these rates at all temper-atuves. Thin oxide films lower, the apparent diffusivity of annealed iron markedly. Cold-worked iron (75 pct) annealed at 430°C has practically the same solubility and diffusivity as anealed iron. THE mechanism of absorption and evolution of hydrogen in solid iron alloys is not fully understood despite the large number of investigations on the subject. While there appear to be many instances of simple solid-state diffusion control of rate, there are many cases where investigators have either concluded that surface reactions are rate controlling or that some "nondiffusible" hydrogen resides in discontinuities. The purposes of this research were 1) to develop a sound method for studying this absorption and evolution behavior, 2) to examine it in pure annealed iron, and 3) to determine the effects of important variables, such as cold work, surface condition, and composition. This paper relates the progress toward these objectives. Much of the present knowledge of the movement of hydrogen through iron alloys comes from permeability studies which involve the measurement of the rate of passage (under steady state conditions) through membranes of the order of 1 mm in thickness. This method is open to question on various grounds. If surface reactions are involved, two different surfaces operating under two different conditions must be considered. The presence of defects in the membrane even of microscopic dimensions may lead to erroneous results. A few people have studied evolution of hydrogen into vacuum, under carefully controlled (unsteady state) conditions, and some of the best diffusivity data come from these sources. Most of these methods measure rates of evolution only and thus do not make a rigorous test for reversibility of behavior. It appeared that an "unsteady state" method that would permit the measurement of both rates of evolution and absorption under controlled temperature, pressure, and other conditions would offer several advantages. Accordingly the modified Sieverts apparatus described below was designed and built. APPARATUS This apparatus is shown in Fig. 1. The sample bulb is detailed in Fig. 2. Hydrogen of 99.5 pct purity and argon of 99.99 pct purity were used. The hydrogen was deoxidized and dried by means of a Deoxo unit followed by a "Dri-erite" tower. The argon was only dried. The gas burette was graduated in tenths of a cubic centimeter, and, with the aid of a travelling magnifier, readings could be estimated to the nearest 0.01 cc. The level of mercury in the burette was adjusted by means of a leveling bulb operated by a fine screw adjustment. An electrical system, consisting of two slightly separated platinum electrodes inserted in one leg of the manometer and a small light bulb was used to obtain rapid accurate adjustment of the pressure. The capillary tubing of the Vycor sample bulb was 3 mm bore as was all tubing back to the gas burette. The samples were discs made from vacuum-melted Ferrovac-E iron of greater than

Jan 1, 1961

By A. M. Gaudin, C. S. Chang

Adsorption was measured for barium ion and laurate radical, using radioactively marked agents, over wide range of concentrations. Laurate adsorbed in absence of barium fails to float. With barium, flotation-effective laurate coating needs not exceed 5 pct of monolayer. Barium adsorption then is several fold larger. No barium laurate coating is obtained. IN general, fatty acids or their alkali salts do not cause clean quartz to float. However, the presence of multivalent ions, such as the alkaline-earth metal ions and heavy-metal ions, can activate quartz and thus make it respond to soap flotation. The mechanism of activation and soap flotation of quartz has been studied.'-' The problem was attacked by making flotation tests rather than by careful adsorption measurements on suspensions of ground quartz in solutions containing the activator and the collector. Attempts have been made to determine the adsorption of barium ion and calcium ion on quartz. The results have only qualitative value as the analytical methods were not adequately precise. Artificial radioactivity and progress in radioactivity detection have introduced a new analytical tool, which can have a sensitivity several orders of magnitude greater than that of other methods."-" In the experimental work described in this paper the radioactive tracers, barium 140 and carbon 14, were employed. Lauric acid was synthesized, having a radioactive carbon atom and was used with Ba" to study the adsorption of barium and laurate ions from aqueous solution on ground quartz. Flotation tests, using the vacuum flotation technique,' were run in conjunction with the adsorption experiments. The object of this investigation was to understand better the reaction mechanism of barium activation and soap flotation of quartz. Methods and Materials Quartz Sample Preparation: Lumps of hand-picked quartz, washed free of stains, were reduced in size to —20 mesh in laboratory crushers. The product was boiled in concentrated hydrochloric acid to remove the abraded iron. The cleaned product then was ground dry in a laboratory porcelain mill, using siliceous pebbles, to pass through a 200 mesh screen. The ground product was deslimed by sedimentation at 37 microns. This —200 mesh +400 mesh material was subjected again to boiling with concentrated hydrochloric acid, followed by washing with distilled water and further desliming. The cycle was repeated until no trace of iron could be detected in the wash liquor. During the final washing, the pH of the wash water was checked to insure the complete removal of the acid from the mineral. The quartz then was rinsed with conductivity water. The resulting product was dried in a drying oven and stored in mason jars. The impurities (spectrographic) were as follows: Al, 0.01 to 0.001 pct; Li, Mg, Mn, Na, each less than 0.001 pct; Fe and Cu not detected. The specific surface of the quartz sample was determined by the krypton gas adsorption method.' The average of two measurements was 881 sq cm per g. The calculated specific surface for 37/74-micron quartz particles, using a surface factor of 1.96; is 900 sq cm per g. This is in good agreement with the measured value of 881 sq cm per g. Beta-Counting Method: The radioactive isotope of barium used in the experimental work was barium 140. It has a sufficiently long half-life, 12.8 days, emits high-energy beta particles (up to 1.0

Jan 1, 1953

By E. R. Buchanan, R. E. Reed-Hill, F. W. Caldwell

Poly crystalline zirconium tensile specimens containing a sizable fraction of grains unfavorably oriented for slip were deformed at 77°K to strains as high as 9 pct. The contributions of the various twinning modes to the observed strain were deter~nlned quantitatively using measured twin volume fractions and Schrnid orientation factors. The results indicate that at very small strains twinning can account for almost all of the observed strain. At higher strains slip becomes the dominant means of deformation. The most significant form of twinning was {1121}. PLASTIC deformation in polycrystalline aggregates is poorly understood even today. This is primarily because slip cannot be easily observed or studied in the interior of a metal. On the other hand, twinning can be discerned inside a specimen after deformation. Twinning studies can therefore often yield quantitative information about polycrystalline deformation not otherwise obtainable. This investigation uses the above aspect of twinning in determining the contributions of twinning to the strain during tensile deformation at a cryogenic temperature (77°K). Polycrystalline zirconium was used as the object of study. There are several reasons why this information should be of value. In contrast to many metals, both the strength and ductility of zirconium increase with a corresponding decrease in temperature. This may be related to the fact that twinning becomes more profuse and complex with decreasing temperature. Secondly, it has been discovered that certain room-temperature properties of zirconium can be greatly improved by about 1 pct prestrain at 77°K. These include an increase in tensile elongation in critical directions by over 50 pctl and a large mechanical hysteresis effect.2 Both effects are believed to be related to twins formed by prestraining: plastic twin-boundary movements accounting for the increased ductility, and anelastic twin-boundary movements for the internal friction. Thirdly, hexagonal metals which twin primarily on only the (1012) mode (such as beryllium, zinc, and magnesium), unlike zirconium, are very limited in ductility at cryogenic temperatures. Therefore, the significance of the various types of twinning in zirconium at cryogenic temperatures is of more than academic interest. Three twinning mechanism-s have beenobserved in quantity at WOK, {1012), {1121), and {1122). Each must be considered as a unique deformation mechanism because of their markedly different twinning shears and basically different angles between twinning and basal planes, as shown in Table I. It is therefore necessary to separately consider the role of each of these twins during plastic deformation at 77°K. Fortunately, each twin possesses unique morphological features which makes its identification possible on a properly prepared metallographic section. The use of the polarized microscope3 can also greatly aid in identification of the various forms of twins. EXPERIMENTAL PROCEDURE Arc-melted sponge zirconium plate, previously describedY4 was machined into 1/8 by 1 1/4 in. threaded-end, round tensile test specimens. These specimens were deformed to the desired strains in a liquid-nitrogen cryostat mounted in an Instron testing machine. The tensile axis was the transverse plate direction. This specimen direction resulted in large numbers of grains in orientations unfavorable for deformation by prism or basal slip. Deformed specimens were sectioned, mounted in epoxy cold-mount resin, ground on metallographic papers, electropolished, and anodized to provide surfaces sensitive to polarized light. Standard quantitative metallographic procedures5

Jan 1, 1965

By W. A. Backofen, B. B. Hundy

THE pole figure interpretation in a recent paper' on the torsion texture of copper was questioned in discussion,2 nd a simpler alternative interpretation was proposed. In the hope of reaching a positive decision in favor of one or the other, further studies were undertaken with 70-30 brass and Armco iron. This additional work has not answered all the questions about torsion textures, but the new pole figures do offer support for the original interpretation and show that the texture is more complicated than the alternative suggests. From the previous paper,' the appropriate method of describing the texture in a twisted bar is known to be a pole figure plotted on a plane tangent to the surface of the bar. The axes of the pole figure are chosen to coincide with directions in the tangent plane which are parallel and at right angles to the bar axis (the longitudinal direction, L.D., and the transverse direction, T.D.), for these directions define symmetry axes of the texture. Therefore any ideal orientation that might be selected to represent the texture gives the crystallographic plane parallel to the tangent plane and the crystallographic direction parallel to the axis of the torsion specimen. A unique texture was assumed in the original interpretation, and the four ideal orientations, (112) [111], (111) [112], (110) [001], and (112) [131], were used in its description. Of the four orientations, the first three provide possible slip directions along the direction labeled T.D. in the pole figure, which is one of maximum shear stress in a specimen being twisted. The alternative, however, simply considered the texture to be made up of two components, one a [Ill] alignment along the specimen axis with a random orientation about the axis as in the wire texture, and the other a superimposed orientation (110) [0011. Experimental evidence for this simpler interpretation was derived from two observations: 1—twisting did not alter the orientation of a single-crystal torsion specimen of a brass in which an octahedral plane was perpendicular (and a [Ill] direction parallel) to the axis, and 2—deformation bands with approximately the (110) [001] orientation were seen to develop in the crystal during the test. The new pole figures were obtained with a brass and an iron torsion specimen of 3/4 and 1 in. gage length diameter, respectively. After annealing, polishing, and twisting, seven rods, about 0.030 in. in diameter and with different orientations, were prepared from material near the surface of each specimen; these were required for the method of constructing pole figures by using a Geiger counter as developed by Norton.5 Since the shear strain gradient in plastic torsion is linear, it was readily computed that the material supplying the brass rods had been subjected to a shear strain of about 3.2, and that supplying the iron rods had received a shear strain of about 2. All details concerning specimen design, testing procedure, rod preparation, and use of the Norton technique can be found in ref. 1. Only this procedure was followed in working out the (111) and (100) pole figures for brass, Figs. 1 and 2. However, X-ray photograms were also involved in the construction of the (110) pole figure for Armco iron, Fig. 3. The photograms were made according to the conventional transmission diffraction technique using a small sheet taken parallel to the specimen surface and etched to a thickness of about 0.0025 in. They were especially important in locating the boundaries

Jan 1, 1954

By F. E. Jaumot, A. Sawatzky

IT as long been recognized that there is a need for accurate systematic data on diffusion in solids. Originally, such data were needed to check existing theories but, more recently, it would appear that it could be used most advantageously as a guide to, or a check for, much needed extensions of the theory of diffusion. The awareness of the conditions which are important to experimental accuracy, the availability of the materials for accurate experiments, and the techniques have improved to the extent that it would now seem to an experimentalist that, unless extensions of the theory are made which take into account the fundamental parameters of the atoms, the electrons, and the lattice in some detail, the field of diffusion will rapidly become an empirical or semi-empirical field. The present series of experiments was designed to provide accurate data on the diffusion of the elements of the IB and IIB subgroups in silver; that is, to provide data on the diffusion in silver of three face-centered-cubic elements with roughly similar chemical properties, but with different atomic or ionic sizes and, similarly, for three noncubic elements. The elements of these two subgroups have different valences, so that the data will not give any significant information on differences due to the crystal structure of the solute. The present paper reports data on the diffusion of copper and mercury in single crystals of silver; the data for the other elements of the IB and IIB subgroups were reported previously by the authors1'2 and by Tomizuka and Slifkin.8 Experiment The general technique used was identical to that previously described1,2 and will not be repeated here. In all cases, the sectioning technique and high specific activity isotopes" were used. The short half-life of Cu64 made decay corrections more important than for the other isotopes, and necessitated around the clock shifts while data were being taken. Each section was counted at least twice (as was true for all longer half-life isotopes) and a monitor was counted along with each sample. Data for the diffusion of mercury in silver are presented at a greater number of temperatures than for the other solutes. This is a result of two circumstances. First, six (of 30) samples had to be discarded because of pitting or recrystallization which occurred during the diffusion runs, even though the samples had sustained high temperature prediffusion anneals. Either of these was a very rare phenomenon in the other experiments, and their occurrence in this case is considered in more detail below. Second, the curves of the logarithm of the activity of the sections vs the square of the penetration deptht were, in general, quite good, but frequently seemed to show one point to have high activity and the next point to have slightly low activity in terms of the best straight line. When this occurred, there was remarkable regularity of the high-low sequence. In any event, the results, as far as diffusion coefficients are concerned, were not affected by this amalgam effect. In fact, the mercury data had the highest internal consistency of any of the data. Results Diffusion of Copper in Silver—-Diffusion runs were made at eight temperatures In the temperature range from 700° to 950°C, for times ranging from 5 to 45 hr. A linear curve of the logarithm or the activity vs the square of the penetration depth was obtained for all samples. Table I summarizes the results for individual samples and Fig. 1 gives a plot 02 the logarithm of the diffusion coefficients as

Jan 1, 1958

By Ottar Dragge, C. Danielsson, Bo Kalling

THE desulphurizing of pig iron has been accomplished with a number of different additions. The oldest and still most commonly used agent is soda, the extensive use of which commenced about 1925, when it was used principally for cupola furnace iron. More recent experience' seems to show that better results can be obtained with sodium hydroxide. The well-known desulphurizing properties of lime have also been exploited in different technical processes. Another material with even more powerful effect is calcium carbide.' The desulphurizing ability of manganese, when added to the ladle in sufficient quantity, should also be mentioned in this connection. During recent years increasing attention has been paid to the desulphurizing properties of metallic magnesium." An addition of a suitable alloy of magnesium is now in use purely for the purpose of sulphur elimination. Of the desulphurizing agents mentioned, lime is by far the cheapest, provided that the reaction can be brought about rapidly and completely. Therefore, a method that makes full use of the desulphurizing ability of lime may be able to compete with other processes. A method developed at the Dom-narfvet Iron and Steel Works (Sweden) will be described, which enables pig iron to be rapidly desulphurized to very low sulphur contents by using a burnt lime powder. as the desulphurizing agent. Lime in Older Processes In cases where lime has been used for the desul-phurization of pig iron, it has generally not been used alone, but mixed with other substances such as fluorspar, to obtain the formation of a molten slag during the process. This method has been tried by Tigerschiold,' who treated the iron with a lime-fluorspar mixture, the stirring of the iron being brought about inductively with low frequency alternating current. Very good results were obtained. A process of this type has also been suggested by R. P. Heuer, U. S. A. The principles of this method, which has been tested in Great Britain by Newell. Lanener. and Parsons." re that a mixture of lime and fluoispar is added to the hot metal in the ladle, while a powerful stream of nitrogen gas is blown into the bath to produce the required intermixing. The results of the tests were unsatisfactory, however. A similar process has been developed at The Steel Co. of Canada, according to a statement by H. M. Griffith.' Here the tests were carried out in a carbon-lined ladle provided with carbon tuyeres in the side wall for blowing nitrogen into the bath. The addition consisted of about 20 lb of a mixture of burnt lime and fluorspar per ton of pig iron. Good results appear to have been achieved. The sulphur content of the pig iron is stated to have been reduced from 0.025 to 0.050 pct down to 0.006 pct. Various methods of desulphurizing pig iron have been tried using lime powder without fluxing material for fusing. Eichholz and Behrendt7 have experimented with blowing a powdered limeicoke mixture with air into the ladle. Their results were, however, not conclusive and the experiments do not appear to have been continued. Similar experiments have been carried out at Domnarfvet, using nitrogen instead of air in order to avoid oxidation. But these attempts were not particularly successful. It appears to be difficult to achieve the required agitation by this means. The strong cooling effect of the gas on the iron is also a serious drawback. A method in many respects similar to that tried at Domnarfvet was tried by Eulenberg and Krus at the end of the 1930's. Here again desulphurization was carried out with lime alone, brought into contact with the molten iron in a rotary furnace. The temperature was kept at the required level, 1400" to 1500°C, by the introduction of a pulverized coal burner in one end of the furnace. The speed of rotating was not given. A paper by Bading and Krus states that, in one of the first experiments, the sulphur content in 56 tons of pig iron was brought down from 0.186 to 0.035 pct in 117 min, but that a considerable shortening of the time would be possible. According to later reports by Eichholz and Behrendt,' it should be possible by this process to achieve a desulphurization speed of 0.35 pct S per hr for a consumption of 6 to 10 pct limestone and 2 to 3 pct coke, as fuel exclusively. The final sulphur content is, however, not stated. Domnarfvet Method After a number of different procedures had been investigated, the tests at Domnarfvet were directed to desulphurization with lime in a rotary furnace. Before going into the practical details of the method, the theoretical aspects will be discussed briefly. If the pig iron does not contain alloying elements other than carbon, the reaction can be expressed most simply by the usual equation: FeS + CaO + C = Fe + CaS + CO [I] 4H,. ~ 34,000 cal That this reaction can be carried through to a very complete desulphurization of pig iron has been shown by OelsenD in a discussion in connection with the Eulenberg and Krus' method. He mentions two laboratory tests, in one of which the sulphur content in the pig iron at 1400°C was reduced from 0.540 to 0.006 pct after treating with 3.35 pct lime. The pig iron had a low manganese content, but other analysis is th. given. Mention also should be made of the recently published investigations by Fischer and Cohnen'" dealing with the influence of the carbon content of the iron on desulphurization with lime, although in this case fluorspar was added also. The tests show that efficient desulphurization is possible with lime in the steel bath, provided that the carbon content is sufficiently high. The temperature employed in these tests was considerably higher (1620") than that normal for treatment of pig iron. The author concludes that the product S% X C% - 0.011 at the temperature in question.

Jan 1, 1952

By Albert G. Anderson

I'm sure many of you have wondered what the term "Cross Flow Filtration" means. I shall explain it so that you more clearly understand it, and I hope you will be able to make it a useful tool for you in the future. In the simplest form of cross flow filtration contaminated fluid. or feed stream, is fed at a high velocity down the center of a porous tube, clear fluid weeps through the porous wall of the tube while the suspended solids concentrate in and are swept away by the feed stream. (See figure A). The high velocity flow of the feed stream tends to keep the inside of the tube free from a build up of contamination. By utilizing this cross flow approach to filtration, we have been able to concentrate such things as paint pigments, titanium dioxide, iron oxide, precipitated metals such as uranium yellowcake or silver, cadmium, some metal hydroxides, as well as gelatinous particles such as yeast cells. You will recognize that cross flow filtration lends itself to a recirculation loop system due to the necessity of maintaining a high velocity on the inside of the tube. This means that we generally have a storage tank from which we are recirculating, and the solids will concentrate in the storage tank. (See figure B). The storage tank design becomes an integral part of the total system and is very important to the total function of the system. The cross flow filter system requires only low pressure, and will work with centrifugal pumps, diaphram pumps or in some cases no pump at all. As an example, if we wanted to filter a 10 GPM (37.8 liters/min.) stream that contained 1% of submicronic particulate, we would flow the 10 GPM (37.8 liters/min.) into the storage tank and recirculate through the filter at a flow of anywhere from 30 GPM to 100 GPM (113.6 to 378.5 liters/min.) depending upon the filterability of the particles from the liquid. As the contaminated fluid recirculates, the solids will build up in the storage tank and loop until they reach a level, which is dependent upon the design of the tank. Then the settling rate in the tank becomes equal to the amount of the solids being fed to the tank, and as a result you reach an equilibrium in the amount of solids being presented to the inside of the porous tubing. We have in effect created a dynamic settling tank or thickener. A flow of 10 GPM (37.8 liters/min.) of clear filtrate comes off for either recycle or disposal with periodic removal of thick mud from the bottom of the tank. The thick mud can be either processed through a centrifuge, a rotary drum filter. a flat bed pressure type filter, or other means. Often the small volume of filtrate from the mud de-watering stream can be recycled back into the recirculation tank so that any fines in the feedback are again concentrated, thus closing the loop on the continuous process 100% removal system. Cross flow filtration will not do efficient molecular separation. We work primarily with discrete solids but we do separate them totally from the liquid carrier fluid. I know some of you will raise your eyebrows when I say that we remove 100% of the particulate, but in measurement tests, we have measured the filtrate cleanliness with instruments that record as low as .05 ppm and the instruments read zero, so that it is less than .05 ppm.

Jan 1, 1980

By A. G. Weber, H. J. Welge

A published calculation method for predicting incompressible, multidimensional fluid displacement has been adapted to the problems of water and gas coning in oil wells. Since depth and radial distance from the wellbore are the two key dimensions affecting the shape of a gas or water cone, coning calculations are well suited to the use of two-dimensional methods. The alternating direction implicit procedure (ADIP) was used for relaxation calculations of two-phase potentials in a two-dimensional grid. From the potentials and the capillary pressure relationship, saturation and pressure distributions were calculated which trace cone growth with time. Predictions have been made for well producing histories both before and after core breakthrough. To check validity of the method, two-dimensional calculations have matched the coning behavior and produced water-oil ratio history of a laboratory sand-packed model. They have also matched the coning behavior of several producing wells, for which the calculations were compared with produced water or gas cuts and logs showing water or gas cone movements. INTRODUCTION The application of two-phase, two-dimensional calculations using ADIP to various reservoir flow problems has been described in the literature.1-3 When adapted for computer solution, this method has proved to be a powerful tool for simulating well and reservoir behavior. This paper discusses the method as applied to well coning calculations. Single-well and coning calculations comprise an especially difficult class of two-dimensional problems which require special techniques for computer calculation and determination of reservoir characteristics. Refs. 4 through 8 describe previous approaches to the coning problem. Several examples of water and gas coning calculations, including studies on both laboratory models and producing wells, are presented here. The two-dimensional method accounts realistically for the most critical parameters affecting coning behavior, including production rate, formation stratification, horizontal and vertical permeabilities, depth of well penetration, gravity and capillary forces. The method considers the different densities and viscosities of the two phases and the relative permeability and capillary pressure characteristics of the rock and fluids. In addition to tracing cone growth in the vicinity of the wellbore, the method calculates the overall movement of the fluid interface throughout the well's drainage volume. Incompressible fluid flow is assumed to occur between the producing interval and the well's limit of drainage. Calculations can be made for the producing history both before and after cone breakthrough. A typical two-dimensional grid or array of blocks used to solve a coning problem contains about 400 blocks. The well's cylindrical drainage volume can be represented by about 20 radial subdivisions and the formation thickness by 20 vertical subdivisions. The grid spacing is normally smaller near the withdrawal interval to define the cone shape accurately. For this work we used an IBM 7074 digital computer having a core memory of 10,000 ten-digit words. A typical study, covering 5 to 10 years of well producing history, required from three to six hours of computing time. MATHEMATICAL SIMULATION OF CONING BEHAVIOR BASIC METHOD In coning calculations, the reservoir volume drained by the producing well is represented by a two-dimensional system of blocks as shown in Fig. 1 for water coning studies. The horizontal dimensions of the blocks increase with radial distance from the well axis in geometric progression, i.e., the block size is small near the wellbore and large near the well's drainage radius (re). Vertically, the blocks are bounded by horizontal planes located at different depths through the

Jan 1, 1965

By E. P. Miesch, P. B. Crawford

A study was made of the effect of permeable and impermeable lenses in a reservoir on the production capacity of a well. Both steady-state and unsteady-state data were obtained. An electrical resistivity model was used to obtain the steady-state data and thermal models. were constructed to obtain the unsteady-state data. The productivity of a well is affected very greatly only when the lenses are close to the well. The effect of circular lenses on the productivity ratio can be correlated with the distance from the center of the lens to the center of the well divided by the radius of the lens. When this dimensionless distance is equal to six or greater, the effect of the lenses on production capacity will be negligible. The pseudo steady-state productivity of a heterogeneous reservoir can be predicted using steady-state data. INTRODUCTION Many analytical solutions of reservoir behavior assume that reservoir rock is uniform and homogeneous. Although this assumption is used, all of the data from core analyses and well logging indicate that the reservoirs are heterogeneous. Very little work has been done on the performance of heterogeneous reservoirs. The work of Landrum,l et al. showed that transient phenomena in oil reservoirs could be studied with thermal models. Pickering and Cotman2 used thermal models to study flow in stratified reservoirs and investigated the effect of inhomo-geneities in oil reservoirs on transient flow performance. Loucks3 made a mathematical study of the pressure build-up in a system composed of two concentric regions of different permeability. Root, Silberberg and Pirson4 studied the effect of the growth of the flooded region on water influx predictions using a thermal model consisting of three concentric cylindrical regions of different thermal properties which simulated the aquifer, the flooded region and the unflooded portion of the original hydrocarbon region. Tomme, et a1.5 made a mathematical study of vertical fractures. The object of this investigation was to study the effect of highly permeable and impermeable lenses in the vicinity of the wellbore on the pressure depletion history of the well. Steady-state data were obtained for both conductive and nonconductive lenses that completely penetrated the formation. The lenses were symmetrically located at various distances from the wellbore. The unsteady-state data were obtained on seven thermal models. EXPERIMENTAL EQUIPMENT AND PROCEDURE STEADY-STATE DATA The steady-state data were obtained from an electrical resistivity model 30 in. in diameter and approximately 11/2 in. deep. The outside of the model was lined with a 30-in. diameter copper strip, which served as the outer boundary of the reservoir. The bottom was covered with a sheet of plexiglass so that it would be nonconductive. The model was filled with a slightly saline solution. The well size was varied from an 0.064-in. diameter copper wire to a 10-in. diameter copper cylinder. Readings were taken with an impedance bridge using AC current to prevent polarization at the contacts. Copper and wax lenses were used to represent infinitely conductive and nonconductive lenses, respectively. The resistance was first measured for each well diameter with no lenses in the reservoir. Then the conductive and nonconductive lenses were spaced symmetrically at various distances from the well and the resistance read from each lens location. The diameters of the conductive lenses were 3, 1.022 and 0.624 in., and those of the nonconductive lenses were 3, 2.25 and 1.563 in. The ratio of each resistance reading to the resistance reading with no lenses in the reservoir

Jan 1, 1965

By C. S. Matthews, M. J. F. Rosenbaum

The purpose of thir work wax to lcarn it~lzut infori~lation could he obtained from various typs of pilot water floods and to attempt to find the optunum pilot patter11, for a revervoir which had previously been depleted by a solution gas drive. The study was made in the laboratory with mathemetical methods a dynamic analog and a potentiotnetric analog. Results werp tested against the field llistorics of a nrrnlber of pilot water floods. At a reasonable valrre of currzulative injection, the total production rate for the one-injector five-spot should reach about 6.5 per cent of injection rate, and for a four-injector five-spot, about 9 per cent. Accurate estimates of ultimate recovery cannot be made on the basis of such snzall prorluction rates. However, with a pilot composed of nine ir1jector.s and 16 producers the production rate is approximately 50 pcr cent of injection rate at a reasonable value of camulative injection. Sonle inforn~ation for extended performance predictions might he obtained from such a large pilot. These conclusions were drawn on the basis of results obtained for unit mobility ratio, and a sturly using tlre potentiometric analog was made of the effect of other mobility ratios to determine the range of applicability of these predictions. For the four-injector, five-spot pilot with the ratio of production to injection rate (before water breakthrough) is about twice that for with it is about two-thirds; and with M0= 10, it is about one-third For high mobility ratios, it was found that the production rate increased considerably as water-cut increased. These result can be used to modify, qualitatively, the inter.pretntions based on curves for the unit rnobilit\. ratio CaSeS. It was found that the maximum ratio of production rate to injection rate obseriled in field pilot floods was of rhe scime order as that prerdicted by these methods. The time required to reach thisr maximum did not generally agree with the time predicted for a homogeti~orir reservoir. The differcrlce between predicted and observed time of response gives an indication of the permeability profile and of the condition of the producin,g wells. Pilot water floods of the pattern type are generally carried out in reservoirs which have been depleted by solution gas drive and are at low pressure. Under these conditions, oil and water can be considered incompressible. It is assumed that, as the water is injected, an oil bank forms ahead of it and that there is a distinct interface between the water zone (or bank) and the oil zone (or bank) and between the oil zone and the region ahead of the oil zone. It is further assumed that only gas is mobile in the unflooded (gas) region, only oil is mobile in the oil bank and only water is mobile in the water bank. The saturations and the mobilities associated with each zone are assumed uniform. We idealize our reservoir to be homogeneous, horizontal and of constant thickness. Effects of gravity within the producing layer are assumed negligible. If the actual time-dcpendent flow problem is approximated by a acries of steady-state problems. the potential and stream function in the oil bank and water hank satisfy Laplace&apos;s equation in two dimensions. We can therefore use a poteiitiometric analog of this system. Potentiometric models have yielded uscful results in this laboratory&apos; and clsewhere in the study of a variety of secondary recovery problems. For the case where M = I, we generally prefer to use theoretical mcthods as well as a simpler dynamic analog. Except where otherwise noted, the ratio, side of five-spot/wellbore radius. is taken to be 3,600. a figure which corrcsponds to a normal-size wellbore in a 10-acre well spacing. THEORETICAL EVALUATION OF VARIOUS PILOT PATTERNS, Mw0 = 1 <&apos;he theoretical models which we used to examine the performance of various pilots are shown in Fig. 1. Image theory was used to determine the ratio of production rate to injection rate as a function of the volumc of the flood. The ratio of production rate to injection rate was chosen because this is an easily measurable quantity which is characteristic of a pilot

By D. A. Flanagan, P. B. Crawford

A study has been made of the feasibility of storing liquid meihane at low pressures in undergrohd caverns. Methane liquefies at — 258°F at atmospheric pressure. It is shown that the methane evaporation rates will rapidly decrease and cool the surrounding rock so that at the end of one month they would be between 50 Mcf/hr for a cavern of 25-ft radius and 700 Mcf/hr for a cavern of 100-ft radius. At the end of 10 years, the evaporation rates would be 18 and 100 Mcf/hr, respectively, for caverns of the same radii. The evaporation rates may be reduced by a factor of 2 to 10 by the application of insulation. The cost for mining the caverns is estimated to be $.75 to $1.25/Mcf of storage. This is substantially less than surface storage; it is believed to be safer and to result in lower maintenance, savings in space and savings in strategic materials. INTRODUCTION During the past few years, there has been an increasing interest in the economic feasibility of liquefying methane. Methane liquefaction is being considered for tanker transport; in addition, liquefaction is being reconsidered for shaving peak gas demands. Several articles have described natural gas liquefaction and the progress of the tanker in making trial runs from the United States to Great Britain to determine the feasibility of tanker transport. At the present time, pipelines are not designed to supply peak gas loads during extremely cold periods such as are oiten encountered in the North and Northeast. Gas is being stored in underground reservoirs en-route to its destination, but in many instances satisfactory storage in porous reservoirs has not been practical, especially along the Eastern seaboard where few petroleum reservoirs have been found. In England and other foreign countries, it is unlikely that satisfactory porous structures could be found, and it may be desirable to mine or excavate the rock to obtain storage. By storing the methane near the consumption point, product availability can be increased during periods of need. Methane liquefies at - 258OF at atmospheric pressure, and 1 cu ft of liquid methane will make about 600 cu ft of gaseous methane. Methane liquefaction for shaving peak gas demands was conducted in Cleveland, Ohio in the early 1940&apos;s. The Cleveland plant was designed to liquefy 4 MMcf/D of natural gas.&apos; Development of several low-temperature processes has resulted in some improvements in the design of low-temperature storage vessels. Liquefied-methane storage vessels can be designed for various evaporation rates. The desired rate will depend on the reason for liquefaction. For shaving peak gas demands, the desired evaporation rates would be low, possibly between 0.2 and 1.5 per cent/day. For methane tanker transport and supply, the evaporation rates must be high and must equal the total daily rate of supply. It is apparent that the cost of storage facilities will depend on the objective and will vary appreciably. Surface storage of methane may cost between $3,000 and $20,000/MMscf, depending on the size and design. The failure of the surface storage vessels at Cleveland, with the resulting injury and loss of many lives, has caused considerable emphasis to be placed upon developing new, safer and improved methods of storing liquid methane. Storage of liquefied propane and butane in underground salt domes has been very satisfactory. The per-barrel cost of underground storage for any sizeable capacity is very small compared to the cost of surface storage for LPG mixtures. The larger the capacity of the underground storage, the less is the per-barrel cost. Numerous advantages appear to exist if liquefied methane could be satisfactorily stored underground, the principal advantages being increased safety, lower initial cost, lower maintenance cost, savings in space and savings in strategic materials. The purpose of this paper is to determine the feasibility of storing liquid methane underground in mined caverns. If a large spherical cavern is dug underground and is filled instantaneously with liquid methane, the surface of the sphere may approach the liquid methane temperature almost instantaneously. The temperature distribution for a sphere of radius r. in an infinite medium, initially at one temperature with spherical surface kept at T, from time t = 0, is given by Eq. 1.&apos;

By W. S. Row

This mine is operated at 4000 to 4500 tons daily through a single shaft, with one rock hoist and one senice hoist. Latest shaft construction is concrete with wooden dividers. Economics of drifters and air-leg drills are compared, as well as costs of shrinkage and sub-level stoping. Economy of large haulage equipment and underground crushing is shown. THE property of Kerr-Addison Gold Mines is located in the Larder Lake district of Ontario about 3 miles from the Quebec boundary on the highway between Kirkland Lake and Noranda. It is roughly 320 miles due north of Toronto and 365 miles north of Buffalo, N. Y. The property was staked in 1906 and after several development efforts finally came into continuous production at 500 tons per day in 1937. In 1939 the daily tonnage -was increased to 1200 and in 1941 to 2100 tons. At the end of the war, because of the shortage of labor, tonnage was down to 1100 tons per day, rising to 2200 in February 1947. A further plant expansion inaugurated in 1946 was completed by December 1948, bringing the capacity of the mine and mill to a production rate of 4000 tons per day. Current production, with 880 men employed, is 4400 tons per day. Changes in mining methods and equipment, required by the recently completed expansion program, have yielded interesting comparative results, some of which are described in this paper. Hoisting of ore and waste from the mine, as well as all service facilities for handling men and supplies underground, are provided for in a single vertical five-compartment shaft. Ore and waste are hoisted in 12½-ton skips in balance by a 14-ft diam hoist with drums 115 in. wide, grooved for 2?-in. wire ropes. This hoist, with a total rope pull of 85,000 lb, was designed to operate at a speed of 2600 ft per min to a depth of 4000 ft, hoisting 15 tons of ore in each skip. It is amplidyne-controlled and powered with dual drive dc motors of 2250-hp each, through flexible couplings to two pinion shafts, one on each side of the main gear. Each motor has a peak load rating of 4500-hp at point of greatest duty. Direct current power to the hoist motors is supplied from a motor generator set consisting of a 5000-hp synchronous motor driving two 1750-kw dc generators. The loading pockets used so far are at depths of 1500 and 2725 ft respectively. Consequently the present hoisting speed is 2000 ft per min and the skips capacity only 12½ tons. Now that the shaft has been deepened the new pockets at 3934 ft will be increasingly used so that the speed and the skip size will be increased up to original specifications as required. Some difficulty in meeting the required hoisting cycle is expected because of the shape of the skips. As now in use the cross-sectional dimensions of the inside of the skip are 4 ft 6 in. x 3 ft 8 in. To hold 121/2 tons, therefore, they are 14 ft 10 in. deep. Since the skips are the conventional overturning or Kim-berley type, considerable hoisting time is lost in retardation in entering the headframe dump and in acceleration when leaving. It is expected that to obtain the specified tonnage from this hoist when the skip size is increased to 15 tons it may be necessary to resort to bottom dump skips. This type of skip could reduce considerably the present excessive acceleration and retardation times. Analyses of the skip-hoisting operation are made monthly, and every effort is made to reduce controllable delays. For a nine-month period in 1951 such an analysis shows that the hoist was in operation 92.5 pct of the available time, which was allotted as shown in Table I. Overall hoisting costs, cage and skip, for the year 1951, averaged $0.168 per ton of ore milled. The mine is serviced entirely by means of a single cage, capacity 35 men or 15,700 lb, operating in con-

Jan 1, 1954

By T. A. Pollard, C. R. Sandberg, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler&apos;, and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl&apos; have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By C. R. Sandberg, T. A. Pollard, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler&apos;, and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl&apos; have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By S. F. Reiter

The paper presents isothermal recrystallization curves for 0.08 and 0.15 pct C steel at subcritical temperatures following small amounts of plastic deformation. The effects of deformation, temperature, and aging on nucleation and growth rates ore described. The free energy of activation for grain boundary migration in steel is given. SEVERAL excellent reviews of the literature have appeared concerning the recrystallization of metals.&apos;-&apos; The present investigation follows the approach advanced by Mehl, Stanley, and Anderson,6-7 in which the rate of recrystallization was analyzed in terms of N, the rate of nucleation, and G, the rate of growth of recrystallization nuclei. Two lots of low carbon, capped steel of the analysis given in Table I were studied. Each lot consisted of a 150 lb coil which had been hot rolled to 0.083 in. and then cold rolled to 0.042 in. at the mill. Strips 0.930 in. wide were sheared perpendicular to the rolling direction. Both steels were normalized before studying their recrystallization characteristics. The strips were cleaned, painted with a magnesia-acetone paste, and made into packs of equal weight, wrapped in 0.002 in. copper foil. The packs were placed in a salt bath at 900°C for 30 min and air cooled. A relief anneal followed in a second salt bath for 15 min at 650°C. The relief anneal was found necessary from early tests in which a longer incubation period and slower rate of recrystallization were observed in relief-annealed lot A steel than in similar material which was strained and recrystallized directly after being normalized. This effect, which indicates the presence of transformation and/or cooling stresses in steel air cooled from above the A, temperature, has also been observed by Samuels8 and Masing.9 Figs. 1 and 2 show the microstructure of lot A and B materials and illustrate the rather uniform No. 8 ASTM grain size produced by this heat treatment. Winlock and Leiter10 observed that strip specimens which had their sharp edges removed elongated more uniformly than those which were not polished. Similarly, when the sheared edges were removed on a belt grinder, it was found in the present investigation that such samples recrystallized more uniformly than did unpolished strips. Therefore, all strips were carefully rounded prior to their extension. The approximate strain limits for the production of large recrystallized grains are from 6 to 12 pct extension." It was found that for the purpose of this investigation, 8 and 9 pct elongation were suitable deformations. The strain rate employed was 0.01 in. per in. per min and produced a yield point elongation of 4 pct. Winlock and Leiter found that mild steel of No. 8 ASTM grain size gave the same yield point elongation when extended at 0.012 in. per in. per min. All of the lot A and B strips extended in tension developed a straight, stretcher strain line at each grip when the upper yield point was reached. The lines were parallel and made an angle of 55" with the edge of the strip. They approached each other with increasing strain and met near the center of the sample at the end of the yield point elongation. Immediately thereafter, a small drop in load was observed and then the load increased in a regular manner with increasing extension. The grips were initially 8 in. apart. After extension, the 6 in. gage length was carefully cut into 1 in. samples. The remainder of the strip was discarded. After a flash pickle in hot 50-50 hydrochloric acid, six samples, each of which had been taken from a different strip, were placed in a basket and submerged in a lead pot for isothermal recrystallization. Although no recovery effect was observed, strain aging did occur after extension. Therefore, samples were always recrystallized within 24 hr after their cold deformation. After recrystallization, the samples were etched with a solution comprised of one part by volume of nitric acid with three parts of water. Bromide printing paper was exposed directly at low magnifications and later used with a mask to measure the desired quantities. First, the average diameter of the largest grain visible in each sample was determined using dividers. Next, the number of recrystallized grains per unit area was counted and recorded as n. Then, for each sample, the combined area of the recrystallized grains was measured by transcribing the grain outlines to standard graph paper. Many determinations of the area of the recrystallized grains were repeated five times and indicated a standard error that was not greater than 25 pct. The average area for six samples was divided by the area of the mask to yield the percentage recrystallized. Recrystallization of 0.08 Pct C Steel The progress of recrystallization at 670°C after 8 pct elongation of lot A steel is shown in Fig. 3, a through f. The shapes of the growing crystals are approximately equiaxed, as is assumed in the

Jan 1, 1953

By A. D. McMaster, M. L. Gimpl, N. Fuschillo

Low-power, high-speed, radiation-resistant, monolithic thin-film integrated circuits require thin-film resistors of high sheet resistance which are compatible with the processing requirements for monolithic silicon devices. For this purpose the structure and electrical resistance as a function of temperature were investigated for Nichrome-silicon monoxide cermet films on passivated single-crystal silicon substrates. A suitable material was found to be a mixture of 50 vol pct Nichrome-50 vol pct silicon monoxide. This mixture, when evaporated onto 300°C substrates. showed a low temperature coefficient up to 400°C, at both 270 and 1100 ohms per square. Resistors deposited on room-temperature, 300°, and 500°C substrates were heat-treated at 600°C in air for short periods of H- seed, low-power, radiation-re sis tant, monolithic thin-film integrated circuits&apos;-3 require thin-film resistors of high sheet resistance and low temperature coefficients of resistance which are compatible with the processing requirements for monolithic silicon active devices. The present investigation is directed towards the replacement of the Nichrome in the sandwich-type SO2-Nichrome-SiO resistors, currently used for DTL monolithic thin-film integrated circuits, with codeposited Nichrome-SiO cermet resistors. These Nichrome-SiO resistors are provided with an overcoat of SiO and are deposited onto thermally oxidized silicon wafers used in monolithic mi-crocircuit technology. The cermet resistors are shown to provide a lower temperature coefficient of time totaling up to 45 min. Electron nzicrogvaphs and diffraction patterns were taken of the as - deposited and heat-treated films. These results were then compared with electrical-resistance measurements for each condition. The room-temperature and 300°C substrate films showed a large decrease in resistance. The 500°C substrate film remained stable throughout the entire cycle. Electron-microscopy results show that the increased stability after heating to 600°C is largely due to crystallization of the films which begins in the vicinity of 400°C. The slight oxidation of the films when heated in air had a negligible effect on the resistance. The resistance of the film was a function of the substructure. resistance and a more versatile compatibility with silicon monolithic processing. Previous work3, on chromium-SiO coevaporated cermet resistors on glass and ceramic substrates without overcoats of SiO has been reported. For the present application, the thin films should possess the following ideal requirements: a) a temperature coefficient of resistance as near zero as possible up to 300°C; b) have a sheet resistance higher than 200 ohms per sq to conserve space on the chip; c) be easy to evaporate; d) be oxidiation- and radiation-resistant; e) show long-time stability at 100°C; } be stable in N2 atmospheres up to 600°C for short periods of time (5 to 45 min). Requirements a to e are well-known. The last requirement was optional and was due to a mask elimi-

Jan 1, 1967