Search Documents

By J. S. Aronofsky

The results of polentiometric model studies and numerical computations are described. The purpose of these studies was to determine the influence of the mobility ratio on flooding efficiencies during water encroachment in petroleum reservoirs, in systems of regular well geometry. The direct line drive well pattern having relative spacing distances of one and one-half to one was investigated for water to oil mobility ratios of 10. 1. and 0.1. The results of these studies indicate that the pattern Sweep efficiency i.; very dependent upon the mobility ratio. INTRODUCTION The almost universal occurrence of water in the immediate neighborhood of oil-hearing sands lends considerable emphasis to the desirability of having a rational analysis of water encroachment into petroleum reservoirs. This problem has 1)become increasingly important with the accelerated use of water injection for secondary recovery operations. The water encroachment problem has been studied from both the theoretical and experimental viewpoints by Muskat.' He has formulated, in a precise manner, differential equations to express flow in a water encroachment system, and rigorous analytical solutions have been developed for the linear, radial, and spherical cases.' However, it becomes exceedingly difficult to solve the encroachment problem rigorously for any general two-dimensional system in which the shape of the two-fluid interface is not immediately evident from the geometry of the system. The difficulty lies in the fart that a rigorous analytical solution requires that the shape of the two-fluid interface be known simultaneously with its instantaneous posititon and with the pressure distribution on both sides of the interface. Although this paper concerns the study of the two-fluid system of oil and water. the methods described in the paper are generally applicable to other systems of incompressible and immiscible fluids in which the displacing fluid sweep-the displaced fluid down to a residual value. In each region it is assumed there is only one mobile fluid and its mobility is constant throughout the region. The Same procedure. with minor alteration, is applicable for gas injection systems provided the assumption can be made that both the displaced and displacing gases can be described by steady state formulas. For the particular two-fluid problem considered in this paper, there is a relative permeability. K, and a viscosity. p. associated with the water zone, and a different permeability and viscosity associated with the oil zone. The K/p (mobility value associated with a particular volume element of the porous medium. at any instant of time. will depend upon whether the element is ahead of. or behind the advancing interface. It is assumed that all volume elements ahead of the interface are characterized by the K/p value for oil, while all the elements behind the interface are characterized by the K/p value for water. That is. the K/p value in the oil zone is constant throughout the oil zone. and the K/p value in the water zone is also constant but different from the value in the oil zone. The symbols are defined at the end of the paper. The well pattern chosen for initial study is the direct line drive; the swept-out area at breakthrough was determined by two different methods which are described in the section on theory. These are: 1. Potentiometric analyzer 2. Numerical computation The use of the potentiometric analyzer for studying water flooding problems is adequately described in the literature, particularly for cases where the mobility in the water zone is assumed equal to the mobility in the oil zone." ' 5

Jan 1, 1952

By P. S. Rudman

Crystal-structure occurrences in binary phase diagrams of fifth- and sixth-period elements, zirconium through palladium and hafniun through platinum, are surveyed with the conclusion that average .group number GN is the dominant correlative. This is interpreted in terms of curves of enthalpy differences between competitiue structures that are continuous func -tions of GN. It is further reasoned that the absolute formation energies may also be continuous functions of GN and if so can be simply obtained from the regatives of the heats of sublimation of the elements. It is shown that from such a heat of formation-GN plot the heats of mixing of alloy systems can be quantitatively predicted. For example, for sixth-period alloys very large exothermic heats of mixing are predicted for systems that span GN = 6 (i.e., Hf-Pt, Ta-Os, ...). Thermodynamic data with which to test these predictions do not appear to be available. If in the latter transition (A subgroup) elements we eliminate from consideration the fourth period where magnetic terms contribute to crystal-structure stability, as shown in Table I, there is a very regular sequence of structures. Using the group number (GN) assignments of Table I, which differ from the usual assignments of 8 to the elements listed as 8 through 10, it is perhaps trivial to say that the phase occurrence is a unique function of GN. However, the interesting aspect of alloys from among these elements is that the alloy structures also appear to correlate with GN. The GN for an alloy is given as GN = SixiGNi, where xi and GNi = atomic fraction and GN of component i, respectively. It is difficult to credit anyone in particular with having discovered this correlation, so many people have implicitly or explicitly employed it in the last 10 years. However, Bloom and Grant' in their correlation of s-phase occurrence with total number of d + s electrons were certainly one of the first to appreciate the correlation and Hume -Rothery2 appears to be responsible for the use of average GN as defined above. Only in the last few years have the phase diagrams of enough alloy systems been determined for this correlation to be surveyed with any degree of completeness. We shall here, in an abbreviated format, summarize the available phase-diagram literature and discuss some consequences of the GN correlation. BINARY DIAGRAM SURVEY In Figs. 1 to 3 the phase occurrences are summarized. In Fig. 4 the occurrence limits for the more common structures are summarized. That the GN of an alloy falls within the occurrence limits of a particular structure is a necessary but not a sufficient condition for the appearance of that structure. The preceding figures certainly illustrate that GN is a dominant correlative for structures of these elements. Many of the diagrams have complications, with ordering and other phase changes occurring in intermediate phases, that cannot be clearly accommodated in such a schematic survey and have hence been left out. Such omissions cannot, however, negate the GN correlations discussed here. Relative atomic size is also an important crystal-structure determinant and since atomic size varies quite regularly with GN for these elements it is possible that the apparent GN correlation is basically a size correlation. In Fig. 5 relative size occurrence limits are given for intermediate phases and extended terminal solutions. The relative volume, ?v/(v) - 2(vA - vB)/(vA + uB), has been taken as the measure

Jan 1, 1968

By Richard H. Jahns

Gem materials, comprising those minerals and closely allied natural substances used for personal adornment, for the fashioning of ornamental objects, or for other decorative purposes, have been valued in human commerce for thousands of years. Amber, jet, natural glasses, and at least 15 mineral species were sought and used by primitive man, and most of the materials in common use today were known well before the time of Christ. Their variety is remarkably great, and the traceable relationships among their respective sources, modes of occurrence, properties, and histories of production, preparation, and valuation are extremely complex. Their uses have been affected by numerous shifts in fashion, and from time to time by superstition, by decrees of church and state, by personal interpretations of religious significance, and by various presumptions concerning medicinal value and influence on human behavior (Ball, 1935, 1950; Kunz, 1913; Pogue, 1915, Whitlock, 1946). Thanks to their durability, portability, and relative scarcity, some gem stones have been employed as media of exchange, and the investment value of the more precious stones has long been recognized and exploited, especially in Asiatic and European countries. Demands for gem materials have stimulated trade among various peoples for centuries, and they also were responsible for much early geographic exploration. During more recent times they have prompted numerous geological studies and a considerable amount of research on the growth and properties of crystals. Apart from decorative uses, some naturally occurring gem materials also are employed directly in numerous fields of industry. Substances that are manufactured for gem or industrial uses are termed synthetic gem materials; most of these correspond in composition and crystal structure to certain minerals found in nature. Imitation materials are those whose use is based upon their resemblance to substances of different composition and greater intrinsic value. Gem stone is a designation commonly restricted to those materials suitable for personal adornment, and the complementary term ornamental stone is applied to the remaining materials not used industrially. This distinction is arbitrary and somewhat inconsistent, owing to numerous geographic differences and historical shifts in end uses. Further, there is a complete gradation between the category of ornamental stones and that of dimension stone. Substances that lose their identity in use are excluded from the gem category; thus a transparent crystal of beryl would not be classified as a gem stone if it were processed as a source of beryllium metal. The adjectives precious and semiprecious have been employed as a means for distinguishing different gem materials on the basis of recognized value. Diamond, emerald, ruby, and sapphire consistently have been regarded as precious stones, and in cut form they commonly are referred to as the noble gems. Opal, pearl, turquoise, and chrysoberyl also have been classed as precious at various times. In recent years, however, it has been increasingly recognized that no satisfactory distinction can be based upon assessed value or any group of specific properties, and now all gem stones normally are referred to as precious. Indeed, the tendency is to extend this term to include many of the so-called ornamental stones, and the designation semi- precious has been essentially discontinued. Unlike most mineral commodities of economic value, gem materials are little altered in processing for end uses. Some stones are used in essentially their natural form, and most are reshaped and polished without significant changes in their fundamental properties. Thus the preparation of gem and ornamental stones is intended chiefly to emphasize or enhance

Jan 1, 1975

By Herbert E. Hawkes

“For his Pioneering leadership in the science and technology of mineral exploration, especially in the development and world-wide application of geochemical methods as major exploration tools, and for his lecture ‘Geothermal Hydrogen’.”

Jan 1, 1981

By E. A. Cisney

The formula for calculating the crystallographic angles of a tetragonal lattice is: C°S =where $ is the angle between planes (HKL) and (hkl). The angles in Table I have been calculated for indium

Jan 1, 1960

By C. C. Templeton, S. S. Rushing

Methods previously developed for the study of air-liquid displacements in microscopic capillaries (inner diameters of 3 to 40 microns) have been used to investigate oil-water displacements in capillaries initially filled with water. Displacement calculations assuming perfect displacement and no capillary pressure hysteresis yielded oil effective viscosities smaller than the macroscopic viscosities. For a given liquid pair, the oil effective viscosity decreased both with decreasing capillary size and with increasing oil-water viscosity ratio. This behavior can be explained by the existence of an annular water film (20 A to 260 A thick) on the capillary wall. When the capillary was first filled with oil, the ratio of the oil effective viscosity to the normal oil vircosity was highest for the first water displacement and decreased with subsequent displacements. Sometimes the oil effective viscosity ratio during the initial water displacement was greater than unity. INTRODUCTION In a previous paper1 a technique was described for studying air-liquid and liquid-liquid displacements in very small capillaries of uniform diameter, in the hope that such microscopic data would further the understanding of the nature of multiphase fluid flow through porous media. That paper contained comprehensive data for air-liquid displacements in Pyrex capillaries, with a few data for oil-water displacements in capillaries initially filled with water. The purpose of this paper is to present more complete results for oil-water displacements in capillaries initially filled with water, and to describe for the first time such observations in capillaries initially filled with oil. In this way the effect of the wetting history of the system upon the displacement process may be studied. METHODS The basic techniques employed were described in the previous paper.1 The measuring procedure and the working equations will be briefly summarized here and a few modilications will be pointed out. For the present study, temperature control within 0.l °C was obtained by placing the microscope and its immediate accessories in a thermostated air bath made of a steel frame covered with plexiglas. The water and oils used in this work were the same as in the previous paper,1 with the addition of "Medium Mineral Oil, U.S.P. supplied by the Harshaw Scientific Co. When the first liquid was introduced into the horizontal capillary, the air-liquid static capillary pressure Poc was measured along the observable length of the capillary. From the relation Poc = 4y cos 0/d, the capillary diameter d could be calculated from Poc if the "microscopic" air-liquid boundary tension, To = y cos On, was taken as equal to its known "macroscopic" value. This calculation involves assuming that cos 6 = 1 and that capillary size or the resulting high interfa-cial curvature has no effect on surface tension, or y cos 6, as was verified by our work, and also by that of Cohan and Meyers on air-liquid-solid systems. For several capillaries these calculated diameter values were compared with values measured visually with a filar micrometer eyepiece. For a given capillary the average values for the two methods agreed within 1 or 2 per cent, but the visual values were less reproducible than the calculated values. Further the calculated values were much less subject to personal bias and took into account any slight ellipticity that might be present. Hence diameters calculated from air-liquid Poc data were used throughout this paper. The second liquid was introduced in such a way that a single spherical interface separated the two liquids. Displacements were always made with a constant pressure, Pt, between the ends of the capillary. The times (t) at which the interface passed certain

Jan 1, 1957

By David Cornell

Back pressure test data on natural gas wells are short time test data of unsteady state nature. Performance curves from which unsteady state effects have been eliminated, called stabilized performance curves, are necessary for understanding the behavior of gas wells and for scheduling gas fields for periods up to 20 years. A calculation procedure has been developed whereby three day, seven day, 30 day, or other duration stabilized performance curves may be calculated from back pressure test data. Such curves are essential for correct gas deliverability calculations and gas field scheduling. I NTRODUCTION The manner in which a back pressure test is made will affect the slope of the back pressure curve because the radius of drainage changes with the length of time that the well is allowed to flow. The reciprocal of the slope of the back pressure curve is defined as n in Equation 1 and varies from 0.5 to 1.0. If the flow is viscous and the radius of drainage constant, n will have a value of 1.0. Because of turbulence in the formation near the wellbore, the pressure drop is usually greater than the viscous flow pressure drop and the exponent n is frequently less than 1.0. In addition, since the radius of drainage moves away from the well-bore during the test, the gas flows a greater and greater distance as the test proceeds, with a resulting increased pressure drop at the high flow rates, and, again, n for the test tends to be less than 1.O. The effect of the changes in the radius of drainage may be calculated and the results used in preparing stabilized performance curves for fixed values of the radius of drainage. Such curves can be made for a defi- nite time period, such as three day, seven day, or 30 day stabilized performance curves. These cukes will show the effect of turbulence near the wellbore and can be used to predict the performance of the gas well at any flow rate or for any production period. In order to discuss the importance of the radius of drainage quantitatively, it should be defined. Muskat" defines the radius of drainage as follows: "... the radius of drainage gives the distance from the well to which the approximately steady state condition has been established ... it being assumed that no gas is removed from any point until the radius of drainage has passed that uoint." An eauivalent definition of the radius of drainage is that radius at which the extrapolated portion of the steady state portion of the pressure gradient gives a pressure equal to the shut-in reservoir pressure. This definition is also implied in the back pressure equation. The radius of drainage is, then, the distance into the reservoir from which the gas appears to originate at any time. It is, however, only a useful concept and does not actually exist. Fig. 1 illustrates the definition of the radius of drainage. At the end of each flow period the radius of drainage may depend on three factors: (1) the duration of each flow period; (2) the rate of gas production; and (3) the rate and duration of the previous flow periods. In general, back pressure tests are conducted with various lengths of production time and at various flow rates. Each individual back pressure test, therefore, requires an individual analysis. The calculations of the change in the radius of drainage have been made, however, for the case of equal flow periods for all four production rates and several sets of increasing flow rates. Fig. 2 illustrates the results of such calculations made by the method of Cornell and Katz2 for the case of 8, = 50,000 and the flow rate increasing in the ratio 1.0:2.0:3.0:4.0 for each flow period. The increase in the radius of drainage during the test is apparent.

Jan 1, 1956

By J. E. Miller, W. E. DeVaney, L. Stroud

During a recent phase study of a natural gas, two stable equilibrium liquid phases were observed at temperatures below —200°F and pressures above 200 psi. This paper reviews the published literature on the occurrence of multiple equilibrium liquid phases and presents analytical data for the vapor and two equilibrium liquid phases of the liquefied natural gas at five experimental conditions. In addition, data for 30 conditions of two-phase equilibria are included. INTRODUCTION The low-temperature phase behavior of gases, most of which contained helium, has been investigated in the laboratories of the Helium Activity for many years. Since 1952, experimental studies of these systems have been continuous as part of the research program at Amarillo, Tex. Because of their value to private industries interested in participating in the Helium Conservation Program, several "Open File" reports containing phase equilibria data for heliurn-bearing natural gases have already been released by the Helium Activity. A paper containing information on the general phase behavior, operating criteria and extensive vapor-liquid data for two helium-containing systems was recently published.2 Additional publications presenting experimental data on the phase relationships of various gas systems are now in process and will be available in the near future. PREVIOUS EXPERIMENTAL WORK Although the formation of multiple liquids has been reported for various systems, to our knowledge this paper is the only substantiated evidence of a vapor-liquid-liquid equilibria in a naturally occurring gas. In 1940, Vink, Ames, and others13 reported the presence of two liquid phases in a hydrocarbon system consisting of mixtures of crude oils, solvents and natural gas. Eilerts and co-workers; published data on the recombined fluids from a gas-condensate well. This condensed gas, containing approximately 76 per cent methane and 24 per cent ethane-plus, exhibited two distinct liquid phases. Weinaug and Bradley 14 observed "unusual" phase behavior in a reservoir mixture. These workers postulated that the anomalous phase behavior was due to the "imminent formation of a second liquid phase". Botkin, Reamer, Sage and Lacey 1 studied two California crude oils that exhibited multiple phases. Kohn and Kuratal0 recently reported two equilibrium liquid phases in the methane-hydrogen sulfide system. Roof and Crawfordl1 and Eakin, et a16 also have reported experiments with binary systems that formed two stable equilibrium liquid phases. APPARATUS AND PROCEDURE A U. S. Bureau of Mines Phase Equilibrium Apparatus was used in conducting this study. The apparatus and procedures employed in its operation have been previously described3 and will not be repeated in detail ill this report. Briefly, the apparatus consists of a windowed cell which can be maintained within ±0.5°F for temperatures between room temperature arid -320°F. Pressure within the cell can be maintained within 0.1 per cent of gauge reading up to 800 psig. Equilibrium vapor and liquid samples are obtained in special containers4 for analysis by a mass spectrometer. Although the accuracy of the analyzer is about ±0.1 mol per cent, the reproducibility of the phase-equilibrium apparatus is considered relatively poor. Values reported from methane and nitrogen are considered accurate to ±l.0 and ±0.6 mol per cent. Data for ethane-plus in the vapor are accurate within 0.2 mol per cent; liquid-phase data for this aggregate compolnent are accurate within 1.5 mol per cent. For helium in the vapor phase, the analytical data are accurate within 0.2 mol per cent; liquid-phase analyses for this component were obtained by the charcoal adsorption method described by Frost8 and are accurate within 0.006 mol per cent. All of these references to the accuracy of reported values are conservative estimates based upon a statistical treatment of reproducibility data obtained with the apparatus.

By Robert G. Nisle

In conducting a pressure build-up test on an oil well, it is often necessary to shut-in the well for a short time prior to initiation of the test. The effect of such a short tern shut-in on the resulting pressure build-up curve is derived. The derivation is based on the Kelvin point-source solution and is a straight forward extension of the fundarnental theory. The magnitude of the effect is evaluatecl and illustrated by an example. It is shown that if the production subsequent to the short term shut-in is at least 10 times the duration of the shrct-in, the effect will he less than 10 per cent. INTRODUCTION It often happens that before a pressure build-up test can be made an additional valve must be mounted on the wellhead, or on arrival at the well, it might be learned that the well had been shut-in several hours the previous day because of a shortage of tankage. In such cases the question is How will such an interruption of production affect the results of a pressure buildup test This work explains how to take account of such a short term shut-in in order that the test can proceed without further delay. The solution presented here makes use of the Kelvin point-source solution used by Horne1 in his treatment of the pressure build-up in wells. Use is also made of the super-position theorem. MATHEMATICAL DERIVATION The well is assumed to have been in production at a constant rate for a relatively long time (several days or more). It is then shut-in for a short time, after which it is again placed in production at the previous constant rate. It continues to produce at this rate until shut-in for a pressure build-up test. The Kelvin point-source solution is expressed by Eq. 1. The sequence of production rates is listed in Table I. TABLE 1 — SEQUENCE OF PRODUCTION RATES Time interval Production Rate At time AT the pressure will be the resultant of the sequence of production rates shown in Table 1. An expression for the super-posed effects may be derived from this table. (Effect of + q for T + ? T) + (Effect of - q for T + ?T-t,) + (Effect of + q for T + ? T - t,) + (Effect of - q for ? T) = Resultant........(2) The desired result is obtained by substituting the effect defined by Eq. 1 for the terms in Eq. 2. The result is Eq. 3. When the logarithmic approximation is made for the exponential integrals and the terms collected the result is Eq. 4. Let T — t2 = T = flowing time following short term shut-in t2 — t1 = 6 = duration of short term shut-in Then T - t1 = r + 6 Substitution of these quantities in Eq. 4 gives p(rw ?T) = P,, -When d<<t - AT then d—&apos;-, T— approaches 1. r + AT In this case the effect of the short term shut-in may be ignored. In order to get some idea of the numerical magnitudes involved take the worst possible case, when AT = 0, that is at shut-in. If the error is made less than some maximum acceptable value at this time, it will be less than the maximum for all subsequent times. d + r Let the difference between the value of y-— and 1 be less than 0.1. Then —6+r — 1 < 0.1 S < O.l T Therefore, if the flowing time following a short term shut-in is at least 10 times the duration of the short term shut-in, the error introduced by ignoring the short term shut-in will not exceed 10 per cent of T+AT A T

Jan 1, 1957

By R. Schuhmann, E. J. Michal

Experimental measurements are reported for the oxygen pressures of iron-silicate slags in equilibrium with solid silica. CO2-CO mixtures were bubbled through the slags in silica crucibles to find equilibrium values of p co2/p co over the entire composition range of silica-saturated slags at 1250°, 1300°, and 1350°C. From the data, the activity vs. composition relationships were calculated for O, FeO, and Fe. SLAGS composed primarily of silica and iron oxides have particular practical importance in copper smelting,&apos; in acid steelmaking, and in accounting for the service behavior of silica refractories. These slags all contain both ferrous and ferric iron, so that, strictly speaking, they are ternary or three-component liquid solutions. For convenience, the components may be taken as the three oxides, FeO, Fe2O3, and SiO2. The constitution of the FeO-Fe,O,-SiO, system at slagmaking temperatures has been studied in a recent investigation described in a separate paper.&apos; This work showed the compositional ranges of iron-silicate melts at temperatures up to 1450°C and also the temperature-composition relationships for melts in equilibrium with various solid phases, including iron, wüstite, magnetite, and silica. Broadly speaking, the compositions of iron-silicate melts at a given temperature should be thought of as varying over a two-dimensional field. In the first place, the silica/total iron oxide relationship varies from maximum silica for silica-saturated melts to minimum silica for wüstite-saturated melts or, at higher temperatures, to zero silica for straight iron oxide melts. In copper smelting, for example, the converter slags approach the low silica limits while the reverberatory-furnace slags have intermediate silica contents, perhaps approaching silica saturation in some circumstances. Liquids formed in the pores of silica refractories presumably are close to silica saturation. The second type of composition variation is in the ferric/ferrous relationship, which depends on how oxidizing or how reducing the process conditions are. Relatively reducing slags, such as acid-steelmaking slags, are lower in ferric oxide and ap- proach saturation with metallic iron. Oxidized slags, such as copper-smelting slags, contain larger percentages of ferric oxides and approach or reach saturation with magnetite. Since these two kinds of compositional variations are closely related to two general properties of interest for all metallurgical slags, basicity, and oxidizing power, it is hoped that the thermodynamic study of iron-silicate slags eventually will throw further light on the significance and nature of these long-recognized but as yet poorly understood slag properties. Measurements of oxygen activity as a function of slag composition, using CO-CO2 gas mixture to fix oxygen activities, afford an experimentally convenient approach to the study of the thermodynamic properties of iron-silicate slags. In a preceding investigation,³ his approach was used in a systematic study of iron-silicate slags in equilibrium with solid iron. The present paper gives the results of a similar study of the iron-silicate slags in equilibrium with solid silica. The experimental work consisted in the determination of CO2: CO ratios in gas mixtures at equilibrium in the system: gas: liquid slag: silica Determinations were made at temperatures from 1250° to 1350°C for slags ranging in composition from iron saturation to magnetite saturation. Also, equilibrium CO,: CO ratios and slag compositions were determined for the univariant systems: gas: liquid slag: silica: iron gas: liquid slag: silica: magnetite Experimental Methods The experimental procedure was essentially the same as that used in the previous study of iron-silicate slags saturated with solid iron," but in order to work with the silica-saturated slags an all-silica crucible assembly was used. A liquid slag was prepared from pure prefused stock materials and brought to equilibrium with the silica crucible at the desired temperature. Then a mixture of CO2 and CO was continuously bubbled through the melt

Jan 1, 1953

By R. V. Smith

The theory for unsteady radial flow of gas, as simplified by Aronofsky and Jenkins, has been reviewed and compared with actual well performance. This study indicated that the theory should be modified by the introduction of an empirical "rate of flow" function "Y(q)". The expansion of the theory on the flow of gas to include a Y(q) term bridges the gap between the theory of unsteady-state flow of gas and actual gas-well behavior. Apparently, the Y(q) term is a function only of the rate of flow for a given well. The completion factor or skin effect is associated with the Y(q) function in such a manner that at least two or more sets of drawdown or pressure build-up tests are needed to separate the Y(q) function and the completion factor. Since the Y(q) function used in this report represents an energy loss dependent on rate of flow near the well-bore of a gas well which is in addition to the loss required by Darcy&apos;s law, the Y(q) function is related to the exponent of the back-pressure curve for a gas well. Techniques and equations are presented which permit the estimation of stabilized and "short-time" deliver-abilities. INTRODUCTION The published theory of radial unsteady-state flow of gas through porous media may be divided into papers1-" which use Darcy&apos;s law as a premise and other papers4-7,16-10 which recognize that Darcy&apos;s law may not hold for flow into gas wells. Miller8 compared theory based on Darcy&apos;s law with actual gas-well behavior and concluded that flow into gas wells deviates from Darcy&apos;s law. Houpeurt4,5 in a thorough investigation of flow into gas wells concluded that the flow rate into a gas well is proportional to the difference between the squares of upstream and downstream pressures elevated to a power between 1.0 and 0.5. Houpeurt also concluded that the deviation from Darcy&apos;s law is caused by irreversible, kinetic energy exchange between the flowing gas and the porous media. Tek9 presented a complete evaluation of the flow problem in gas wells. Several investigators state that the inclusion of gas compressibility into flow equations based on Darcy&apos;s law does not begin to explain the flow problems suggested by gas-well behavior. For many years, natural-gas engineers have realized that gas wells by actual test exhibit behavior which cannot be explained entirely by Darcy&apos;s law. Many tests on gas wells have indicated that a value of 0.85 for the exponent in the back-pressure equation is a fair average for many wells. However, Darcy&apos;s law either for steady or unsteady-state conditions indicates that the exponent should be 1.0. The purpose of this paper has been to examine this deviation from Darcy&apos;s law and to present a reasonable mathematical means of evaluating the deviation. The method presented in this paper combines previous theory with an empirical function to provide a description of gas-well behavior. All equations in this report are for constant terminal-rate-of-flow conditions. While present methods of testing wells usually are not based on this concept, they are believed to approach the constant terminal-rate conditions closely in many tests. However, several gas wells were tested by holding the rates of flow constant in order to approach more closely the constant terminal rate of flow; the results confirm the concepts presented in this paper. REVIEW OF THEORY FOR UNSTEADY-STATE FLOW OF GAS The basic partial differential equation for the flow of natural gas has been given in technical literature; the equation for a cylindrical reservoir with the well in the center is

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

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

By I. S. Levy, J. L. Brimhall, B. Mastel

Specimens of&apos; Hastelloy X-280, a low-cobalt version of the solid-solution- hardened nickel-base alloy Hastelloy X, were given a series of thermomechanical treatttzents. They were then Lensile-tested at room temperature and 1350°F. The Microstructures of the specimens tested at 1350°F were studied by electron Microscopy. Microscopy showed that dislocation tangling, caused by the cold working, mas still in evidence after aging-for 24 hr to temperatures as high as 1550oF. (Cr, Mo)23 C6, (Fe, Mo)6 C&apos;. as well us an unidentifiable phase precipitale due to the combination of the thermomechanical treatments and subsequent tensile strain. The room -temperature and 1350°F tensile properties show a general increase ill strength and decreccse in ductility with increusing cold-work level. This is due to the dislocation taangling-. PerLurbatiot~s to this general trend occur- due to solid-solution weakening at the higher cold-work levels as a direct result of the preciplitation reactions. THE requirements for corrosion resistance and high-temperature strength place nickel-base alloys among the more promising materials for advanced nuclear reactor applications. However, when these alloys are irradiated at service temperatures, they generally suffer severe embrittlement&apos;" and reduction in rupture life.3 Certain experimental treatments resulted in improvements in the tensile and stress-rupture properties of Hastelloy X-280, a solution-hardened nickel-base alloy.&apos; To provide a basis for understanding these improvements, unirradiated tensile specimens were examined by electron microscopy. The results are reported in this paper. EXPERIMENTAL DETAILS Material. Hastelloy X-280, a low-cobalt version of Hastelloy X, is a member of the solution-hardened class of nickel-base alloys. The Hastelloy X-280 was obtained from 4- and 1-in: thick plates in the mill -annealed condition. Table I gives the chemical composition of the material as reported by the National Spectrographic Laboratories. ~x~erimental Treatments. Plates of suitable size, 12 by 3 by 1 in., or by + in., were cut under a water jet from the stock material and examined ultrasonic -ally for homogeneity of grain size and absence of voids. The acceptable plates were solution-treated at 2150°Ffor 1 hr and air-cooled. These plates were then reduced in thickness by cold rolling 0, 10, 15: or 20 pct along the original rolling direction. The plates were machined into buttonhead tensile specimens with a gage diameter of 0.125 in., with a gage length of 1.125 in., and with their tensile axes parallel to the rolling direction. The specimens were examined for quality of finish with a penetrant dye, ultra son ically checked for voids, and finally given a dimensional tolerance check with an optical comparator. The acceptable specimens were aged for 24 hr at either 1000" or 1500 F and air-cooled. These specimens were then held at 540°F in a water loop for approximately 18 days to qualify as controls for similar specimens irradiated under the same conditions of time, temperature, water chemistry, and flow rate. Table U summarizes the treatments. These treatments were investigated to establish whether prior precipitation and stabilization of carbide phases at strain sites within the matrix could reduce the amount of carbides available for precipitation at grain boundaries during service. Tensile-Test Conditions. Specimens were tested to fracture at room temperature or 1350 F in air on an Instron tensile-testing machine at a crosshead speed of 0.05 in. per min. The time for the furnace to reach equilibrium at 1350°F was approximately 30 min; soak time prior to testing was 15 min; and test duration was approximately 5 min.

Jan 1, 1968

By William C. Hagel

Oxygen-18 exchange between gaseous oxygen, held at a pressure of 125 mm Hg in a PL-IORh chamber, and splzeres of a Fe2O3 containing three or less grains was determined from 9000 to 1250°C. Isotope equilibrum on crystal surfaces appears to be quickly established, and the rate-determining factor is self-difiusion confbvrning to the relation D = 2.04 esp(- 77,900/RT) sq cin sec-&apos;. Changing sphere diameters caused no detectable variation in diffusion coefficients. Oxygen ions diffuse faster than cations at low temperatures, but a crossover may occur at about 1150°C. INVESTIGATIONS by Davies et al.&apos; and Himmel et al2 have suggested that the mechanism of iron oxidation above 600°C involves a predominance of cation migration in the wüstite (FeO) and magnetite (Fe301) layers and oxygen-ion migration in the hematite (a Fe2O3) layer. Evidence for little cation migration during hematite growth consists of: 1) the disappearance of radioactive-silver markers, presumed to vaporize from the oxide-gas interface, and 2) a large discrepancy between parabolic rate constants for magnetite oxidation and those calculated from Lindner&apos;s3 data for cation self-diffusion, D:, using Wagner&apos;s theory.4 Reversing the calculations gave anticipated self-diffusion coefficients, D:, for oxygen ions of 3.6 x 10-9 sq cm sec-1 at 1100°C and 5.9 x 10-10 sq cm sec-1 at 1000°C. Since then Kingery et a1.5 reported three values for oxygen-ion diffusion in polycrystalline, sintered hematite. These were obtained near 1200°C, at an oxygen pressure of approximately 150 mm Hg, and could be represented by the Arrhenius expression, Dg = 1 x 10 11 exp(-146,000/RT) sq cm sec-1. Extrapolation to temperatures below 1100°C infers that oxygen ions are less mobile than cations. This contradicts the aforementioned oxidation evidence and also coble7s study of the initial sintering of hematite. From the latter, apparent diffusion coefficients for the slowest rate-controlling species could be calculated, and fair agreement was found with an extension of Lindner&apos;s data for cation diffusion. Fig. 1 depicts such comparisons graphically. Three points have been included for the determinations of Himmel et al. on cation diffusion in large, natural hematite crystals surrounded by pure oxygen at atmospheric pressure. No dependence on crystal log raphic direction was noted. Since the value of 10-l4 sq cm sec-&apos; at 1000°C actually represents an upper limit, one might interpret that there is a small decrease in D: when single crystals are employed. On the whole, however, cation diffusion appears relatively well-defined,* while the few results tension to lower temperatures. After having gained experience measuring oxygen-ion diffusion in certain glasses7 and in the corundum oxide, a Cr2O3, the purpose of this work was to determine D; values for hematite and thus possibly assist a mechanistic understanding of the kinetics of the oxidation and reduction of solid iron oxides. I) EXPERIMENTAL PROCEDURE A) Material. Hematite powder of moderately high purity was received from the Columbian Carbon Co. An accompanying lot analysis listed 0.03 H2O, 0.008 weight percents, with the remainder Fe2O3. Subsequent treatments were ball milling, drying, prepressing, and hot-pressing for 60 min at 1350°C under 5000 psi. Since graphite dies coated with a slurry of alumina were used for hot-pressing, the surrounding atmosphere at temperature was sufficiently reducing to form some magnetite on the outer surface of each cylinder. After removing any attached alumina, a stabilizing anneal of 50 hr at 1300°C in oxygen was added. The final product was nonmagnetic and exceeded 97 pct of the density (5.24 g cm-3) usually given for FezO3. On polishing cylindrical ends, counts were made of the number of grain boundary intercepts per

Jan 1, 1967

By W. van der Knapp

A theoretical and experimental ana1ysis is given of the change in volume of a porous medium due to changes in external and internal pressures. The result enable one to deduce directly the effect of large increlents in stresses. It is shown that the three-dimensional representation of the volume of a nonlinear porous system is a cylindrical rules surface, the generating lines of Which have a slope according to the solid rock compressibility. Analytical expression for such a ,surface is very simple. It follows that the bulk and pore compressihilities depend only on the difference between pore fluid and external hydrostatic pressure, i.e., the effective rock frame pressure. From measurements on a number of sandstones as well as limestones it can be conduded that for practical values of this effective frame pressure the pore compressibility falls in the range between the compressibility of water and that of undersaturated crude oil. In general, pore compressibility is found to be higher the lower the porosity. For a particular limestone reservoir pore compressibility and porosity could he related by means of a simple empirical formula. Results from the described study find a direct applicxtion in material balance calculations and in prob1ems of liquid flow through porous media. Since the elastic constants of the rock bulk material enter into equatiorls for velocity of acoustic waves, the results are also of importance for calculating rates of wave propagation in these media. It is shown for sandstones with 1.5 to -10 prr cent porosity that laboratory measured porosities may differ from those lender reservoir conditions by about I per cent of their value. The difference for low porosity limestotnes can be of the order of 10 per cent. The subsidence caused by elastic deformation of reservoir rock is shown to he very small. General formulas relating pore and bulk volume to external and internal pressure can be derived, as has been shown by Biot1, Gassmann2 and Geertsma.2 Their theories describe the differential change in volume caused by a variation in pore fluid tension and external hydrostatic tension. It has been known for a long timei that the stress-strain relations of a packing of spheres are essentially noninear. Our experiments show that for almost all reservoir rocks which arc elastically deformable the relation between volume and loading is also a nonlinear one, resembling that of a packing of spheres very closely. The behavior of these media is often described by using the linear stress-strain relations for small incremental stresses and strains in the vicinity of a pre-stressed condition." The present study differs from these considerations in that it enables one to deduce the effect of large increments in the stresses directly. In analogy to the stress tensor for a homogeneous body the stress tensor for a porous body can be written in the conventional way: where the tensor components, o, and t, represent the total forces acting on the faces of a unit cube of porous material. This tensor can be split into two parts. a hydrostatic and a deviatoric part, as follows2,3,1 where the component of the hydrostatic stress system, a. equals The hydrostatic stress causes a change in volume whereas the deviatoric stress system gives rise to a change in shape only. One of the difficulties in describing changes of volume occurring in reservoir rocks is that the relations connecting these volume changes with external and internal stresses are essentially non-linear. A complete description requires the bulk volume, V,,, as well as the pore volume,V,,. to be known as a function of both the changes in pore fluid tension (v) and the hydrostatic tension on the outer boundary.

By N. Jaffe, J. E. Dorn

This investigation is concerned with the possibility that the creep resistance of metals might be reduced as a result of recovery at the creep temperature when the applied stress is reduced. For this purpose the effect of stress on the recovery of the creep behavior of pure polycrys talline aluminum was investigated in the cross-slip temperature range where the activation energy for creep is 28.000 cal per mole. Before recovery all specimens were given the same substructure by creeping them under a stress of 7800psi to a strain of 0.218. Recovery was then accomplished by holding the specimens at the creep temperatures of 327" and 337°K under a series of reduced stresses for various periods of time. When the original stress was reapplied the specimens crept more rapidly than when the original stress was first removed, illustrating that recovery did take place. The extent of recovery could be indexed in terms of the UNDER certain conditions of test, the creep rate decreases over the primary stage of creep often by as much as a factor of 10~B. Since the external variables of stress and temperature are held constant, such pronounced and significant changes in the creep rate must be attributed to various substructural modifications that are introduced as a result of creep itself. At the secondary stage of creep, over which the creep rate remains constant and the substructure no longer changes, a balance must be achieved between the rate of increase in creep resistance due to creep straining and the rate of recovery of the substructural modifications. The general validity of these concepts on substructural changes has already been shown,&apos; including the corollary that a more creep-resistant substructure is introduced under a more rapid creep rate arising from a higher applied stress.&apos; In spite of these well-documented facts, very little progress has been additional strain undertaken by the recovered specimen to achieve the same strain rate as the unrecovered specimen. It is significant that the higher recovery stresses gave the most rapid rates of recovery. Over the high range of reduced stresses, where some creq took place during recovery, the activation energy for recovery equaled that for creep, namely 28,000 cal per mole. But for stresses during recovery below 0.80 of the original creep stress, no creep was detected during recovery and an activation energy for recovery of about 72,000 cal per mole was obtained. Whereas the activation energy of 28,000 cal per mole suggests that cross-slip of dislocations is responsible for this recovery, the mechanism for the 72,000 cal per mole process is not known. This activation energy, however, is about twice that for self-difhsion. made in determining, in satisfactory detail, precisely what substructural changes so effectively modify the creep rate.&apos; Up to the present, electron microscopy has failed to uncover, in a satisfactory quantitative manner, the mechanistic processes involved. Since, however, creep, at least at moderate and high temperatures, arises from one or more types of dynamic recovery processes, it should be possible to shed additional light on the total problem by investigating the recovery of creep resistance. The importance of such an investigation becomes apparent when the fact is recognized that there has yet been no direct experimental proof that the rate of softening due to substructural changes as a result of recovery equals the rate of hardening due to creep-induced substructural strengthening over the secondary stages of creep. Furthermore the problem of what happens to the creep-strain induced strengthening when the material recovers either at zero or some reduced operating stress has significant engineering consequences. It was for these reasons that the present investigation on the recovery of creep behavior was initiated. Inasmuch as a substantial background has now developed on the creep of high-purity aluminum, this material was again selected for the current investigation. Furthermore extensive evidence is now

Jan 1, 1963

LEAD and tin mixed together make an alloy through an attachment of natural affinity which they have with each other, so that when they are mixed it is difficult to recognize by the sight which one it is, if the essences of their qualities are not surpassed by more than half. If, indeed, one can be recognized, it is the tin which gives off a certain sharp odor to the sense of smell, besides having a greater hardness and whiteness. The best of these is the one that is the purest in its kind. In the alloy of tin it is the one that is the whitest and hardest, or, when it is bent or some thin end of it is held tightly by the teeth, is heard to crackle as frozen water does. Because this happens more with this than with other metals, I do not wish now to inquire into it. It is enough that it be a sign for recognizing when tin is pure or mixed. Pewterers say that it should be joined in alloy with lead and that it is better for working when it contains four or six pounds to the hundred because it makes it softer to the hammer and more fluid in casting. But I do not like to hear this, even if all they say is true, for I see that the tin that comes from England, when worked as well as in cakes that show it to be pure, is much more beautiful and better in all works than is that made in Venice. And I believe that it is not better for any work when mixed than when pure, except for the advantage of the master who, if he does nothing else, sells lead that is worth little for tin that is worth much. I find that when mixed it serves in only two things. One is when it is wished to make a solder for copper, for use when it is not worth while to solder the bottoms and large pieces of vessels with low-grade silver as should be done. The other is when master potters calcine it in order to make glazes and to apply white to their vases. I consider lead to be useless in all other works of any kind of metal, since it unites with none of the other metals of good company except tin.

Jan 1, 1942

By J. T. Norton, A. L. Mowry

S. J. SINDEBAND*—(1) Discussing the properties of the powders used, Mr. Rostoker mentioned a silicon powder as being between 150 and 325 mesh. We always had much difficulty in measuring particle size of silicon powder by screen analysis, because of its tendency to agglomerate. I would be interested in the method used to measure these particle sizes. (2) Mr. Rostoker maintains that the measurement of magnetic properties is more accurate in determining complete homogenization of an alloy, than is the X ray diffraction technique. I wonder whether this is really true. Even when diffusion is complete, magnetic properties will change during continuation of the heat treatment, due to change in density and pore shapes, as has been demonstrated by Mr. Rostoker himself. Which effect is faster, the diffusion of the alloy components or the attainment of final density, will, in my opinion, depend very much on initial material and sintering conditions used. For magnetic applications, of course, the magnetic measurements are most suited for determining the point at which equilibrium conditions are reached. (3) I was glad to see that Mr. Rostoker had used sample sizes which are similar to those I had used previously, and I believe that ring samples of about 2 in. od will become more or less standard samples for the determination of magnetic properties of powder metallurgy products. I should like to ask how the density was changed for the samples on which the effect of porosity was investi- gated, and, further, why the sintering temperature used for this series of experiments was so low. It is somewhat surprising that a 24-hr treatment at 850°C acts so much more to spheroidize the pores than a treatment for 1 hr at 1100°C. But the effect seems to be indisputable. I have made an investigation of the influence of sintering temperature for longer times on the maximum permeabilities of iron powder compacts which was presented at the International Powder Metallurgy Conference in Graz, in 1948, and I hope that these results will be published here within a short time. The results of this investigation are confirmed by Mr. Rostoker&apos;s findings that the change of time from 1 to 24 hr changes the permeability behavior considerably in the direction as indicated by Polder and Van Santen for spheroidiza-tion of pores. I do not quite understand Mr. Rostoker&apos;s explanation of the influ-

Jan 1, 1950

By John Henry Harden

THIS shaft, located in the northern anthracite coal-field about 2300 feet southwest from the court-house at Wilkes-Barre, in the County of Luzerne, Pa., is the property of the Lehigh & Wilkes¬Barre Coal Company. The sinking and cribbing to the solid rock (see isometrical view on Plate IX) was conducted by Mr. W. W. Kenrick. Dimensions of the earthwork to the foundation of the cribbing, 51&apos; 1/2" long, 17&apos; 4" wide, 31&apos; deep. Total amount of excavation in cubic yards, 998.5. Dimensions within the permanent cribbing, 45&apos; 4" long, 11&apos; 16" wide, 31&apos; deep, divided into six parts by buttons 8" x 12", each division, 7&apos; x 11’ 16". Permanent cribbing 12" x 12"; temporary cribbing 8" x 10"; backing, 3" thick ; clay puddle to ex¬clude surface water, 12" thick; conductors, 8" x 5"; lining, 1" thick from top to bottom of the shaft. In the latter part of the year 1871, Mr. H. Murray & Son contracted with the company for the continuation of the work at the rate of four hundred dollars per yard in depth, to a point 100 feet below the Hillman seam of coal ; for the succeeding 100 feet the price to be increased twenty-five dollars per yard, and so on until the bottom rock of the Baltimore seam of coal was found. The size of the shaft within the cribbing timbers is 45 feet 4 inches by 11 feet 6 inches, having six divisions of equal dimensions, 7 feet by 11 feet 6 inches, giving room for what is known as the" standard diamond" car, whose capacity when level full is 100 cubic feet, but when loaded by the miner 4 inches above the top contains 112 cubic feet.* Two divisions of this shaft are intended for hoisting coal from the Hillman seam, two from the Baltimore seam, one for ventilation, and one for pumps. The works were planned for two pairs of direct-acting, high-pres¬sure engines, twenty-four boilers, and for an estimated output of 2500 tons of coal per day from the two seams. * Size of the "standard diamond" car, 8&apos; 10” x 4&apos; 3" x 2’ 8” = 100.1.4; loaded 4 inches above the top, 8&apos; 10” x 4’ 3” x 3’ “0"= 112.7.6.

Jan 1, 1877

By N. A. Gokcen, S. Fujishiro

THE pressures of Mn(g) in equilibrium with Mn7C3 and graphite have been measured by McCabe and Hudson&apos; and Butler, McCabe, and paxton2 by means of graphite, zirconia, and Ta-Mo Knudsen cells. The details of corrections in their measurements were somewhat different from those of the authors;3&apos;4 hence, further independent results at higher temperatures seemed to be warranted. The stoichiom-etry of Mn7C3 in equilibrium with graphite and its stability up to the vicinity of 1373 °K have been satisfactorily established by Kuo and persson5 (See also Hansen and Anderko).&apos; The method of investigation has been described elsewhere in detail.3&apos;4 The carbide was prepared in the manner described by McCabe and Hudson. In order to avoid possible reactions with refractory oxides, only graphite cells were used for the effusion experiments. The authors&apos; results are the runs M-1, M-2 and M-3 listed in Table I, together with those of McCabe and co-workers. The heat capacity, c;, of Mn7C3 is not known; therefore McCable et al. plotted their data as log P vs 1/T and the results were only expressed first&apos; as, and later2 as logP=- -14220+6.06 without further thermodynamic calculations. However, it is possible to estimate closely the necessary heat capacity as follows: When H; - &98 of CrC, or MnO, is plotted vs x, corresponding to various compounds, with temperature as the parameter, it is seen that the relationship is represented by a line of very small curvature. This relationship is also followed by other compounds which do not exhibit solid state phase transformations. In view of the fact that chromium and manganese have the atomic numbers of 24 and 25 respectively, it is reasonable to assume that H; - H;ae vs x for MnC, follows the same type of lines as those for CrC, for which adequate data are available.7 Similar correlations appear to be valid for oxides of two neighboring metals for which reliable data are available. The procedure is thus to move the lines for CrC, vertically (parallel to the axis of HT - H298) until they coincide with the points for MnC1/3, the only carbide of manganese whose H; - HO298 has been

Jan 1, 1963