Search Documents

By L. S. Darken

A quadratic formalism is developed lor the representation of the excess free energy, and of the activity coefficients of each component of a ternary system in the vicinity of a single component selected as solvent. This formalism, in contrast to the e formalism, is thermodynamically consistent at .finite concentrations. Numerous specific systems are used as illustrations, with cross checks where possible; these systems are Fe-C-X systems in the vicinity of 1600°C and include practically all those for which reliable data are available. It is shown that this formalism is usually quite adequate to rebresent available data up to a solute concentration of 20 to 30 at. pet. It is hoped that the formalism here presented may prove useful for the treatment and tabulation of data on ternary and multi-component solutions. SINCE, as discussed in a prior paper,' we do not as yet have an adequate basic understanding of the nature of binary metallic solutions, we cannot hope to be in better position for ternary and multicomponent systems. However, the accumulation of thermodynamic data on ternary and multicomponent metallic solutions has led to an urgent need for a rational formalism for the presentation and use of such data. Notable attempts to represent the isothermal isobaric thermodynamics of ternary solutions were made by Benedict et 01.' and by Wohl.3 These attempts were directed toward developing an analytic expression adequate to encompass the entire ternary triangle, primarily for organic systems; they generally involve the power-series approach, which has been shown' to be inappropriate, in general, even for binary metallic systems. On the other hand, interest in the thermodynamics of ternary and multicomponent metallic solutions has tended to focus either on general thermodynamic relations4-10 over wide compositional ranges, or on cases where one or more components have been at rather low concentration. This latter approach has been followed by Wagner,= by chipman,11 and by Alcock and Richardson,12 and more recently in many publications by the same and other authors. It has already been shown' that the thermodynamic behavior of binary metallic systems is, in general, relatively simple only in the terminal regions. Hence, it would appear that any approach, at this time, to a simple formalism for ternary and multicomponent systems must, a fortiori, be limited to the vicinity of the major component or solvent. The present treatment is aimed to be useful for a given solvent with moderate solute concentrations ranging up to 20 or 25 at. pet and in exceptional cases to considerably higher concentrations. In order to achieve even this modest objective, it is necessary that any appropriate formalism have built into it 1) a thermodynamic consistency (at least within the desired degree of approximation), and 2) an approach at infinite dilution to Raoult's law for the solvent and to Henry's law for each solute. Let us first formulate the condition of thermodynamic consistency for a ternary system. Under isothermal isobaric conditions, we may write for the molal value G of any extensive function dG - C1 dN1 + G2 dN2 + G3 dN3 [l] where the G's are the corresponding partial molal quantities and the N's are the atom fractions. Taking component 1 as the solvent, this relation may be rewritten in terms of the independent variables N2 and N3

Jan 1, 1968

By M. R. J. Wyllie

Eight laboratory-prepared aqueous base drilling muds representing common mud types, and 15 aqueous base drilling muds sampled in the field, have been used in an experimental investigation of the relationship between the streaming potential developed across a mud filter cake and the differential pressure causing filtrate flow. It has been established that the relationship is of the form, ES = kPy, where Es is the streaming potential, P the differential pressure and k and y are constants for any particular mud system at constant temperature. Beyond a certain minimum thickness this relationship is independent of filter cake thickness. Experiment has shown that both k and y are temperature dependent parameters and that in certain mud systems there is a critical value of P above which the simple relationship Es = kPY ceases to hold. The constant k tends to be a function of the mud resistivity at 77ºF, but the exponent y has not been found to be related to any commonly measured mud property. Certain field experiments are described which serve to demonstrate that the total SP (self potential) developed opposite a formation may be quantitatively accounted for as the algebraic sum of an electrochemical potential and an electrokinetic potential. A chart to permit the simple use of mean solution activities in the interpretation of electrochemical SP data is given and its use described. Various field methods for computing the electrokinetic SP component of the total SP are considered and rejected as being unduly complicated and inaccurate. It is shown that the employment of low resistivity drilling muds is essential if optimum accuracy is desired in the use of the S.P. curve for the purpose of computing connate water salinities by the electrochemical SP method. The use of such low resistivity drilling muds is accordingly recommended. INTRODUCTION In a previous paper' a quantitative analysis of the electrochemical compovent of the SP (self potential) curve was made and it was shown that the electrochemical* emf, E, could be expressed as E = ZtRT/F In a,/a,, where t is the transference number of the chloride ion, R is the gas constant, T the absolute temperature of the formation water and mud, a, the mean activity of the connate water (assumed to be sodium chloride) and a, the mean activity of the mud (as sodium chloride). The essential accuracy of this expression was confirmed by an investigation 2 of the SP kicks opposite sands of known water salinity. The logs examined had been routine runs in oil wells in a number of Mid-Continent states. Further confirmation of the original analysis of the electrochemical SP is implicit in the paper by Tixier3 wherein a relationship between the SP and the ratio of connate water and mud resistivities, based only on field data, is given. Allowing for the fact that mud and connate water resistivities were used instead of activities, and that electrokinetic effects were ignored, the constant reported by Tixier is essentially 2tRT/F. The error involved in the universal substitution of resistivities for activities will receive further reference below. Possible reasons for the fact that the electrochemical properties of shales do not appear to be significantly dependent upon either the nature or the quantity of cation exchange materials, such as clays, they contain has been considered at some length by Wyllie and Patnode.4 In the same paper, the probable effects of shale thickness and compaction on electrochemical efficiency, and hence the inherent difficulty of obtaining reliable laboratory-scale data, were also pointed out. However, it is

Jan 1, 1951

By A. T. Chatas

A description of the geometrical characteristics of spherical reservoir systems, a discussion of unsteady-state flow of such systems and examples of engineering applications are presented as backgmund material. The fundamental differential equation, a description of average spherical permeability and the introduction of the Laplace transformation serve as theoretical foundations. Engineering concepts are investigated to indicate particular solutions of interest, which are analytically obtained with the aid of the Laplace transform. These are numerically evaluated by computer, and presented in tabular form. INTRODUCTION A tractable mathematical analysis of unsteady fluid flow through porous media generally requires incorporation of a geometrical symmetry. The simplest forms include the linear, cylindrical (radial) and spherical. Most analytical endeavors have concentrated on cylindrical symmetry because it occurs more often in petroleum reservoirs. Nevertheless, some reservoir systems do exist that are better approximated by spherical geometry. Review of technical literature revealed but a single reference to unsteady spherical flow in petroleum reservoirs.l The motive and purpose of the present work was to remove this gap in technical information, and to provide the practicing engineer with some useful analytical tools. The mathematical details associated with the particular solutions of interest involved use of the Laplace transformation. Hurst and van Everdingen previously demonstrated the efficacy of this operational technique, and in many respects the present treatment was patterned after their earlier work.2 PRELIMINARY CONSIDERATIONS GEOMETRICAL CHARACTERISTICS Geometrically, a spherical reservoir system is defined at any instant of time by two concentric hemispheres whose physical properties of interest vary only with the radial distance. Every physical property is thus restricted to be a space function of only one variable: the distance along a radius vector emanating from the center. Such a system is composed of an outer region and an inner region, separated by a defined internal boundary. The inner region simply extends inward from this boundary, whereas the outer region extends outward from it to an external boundary. The position of the internal boundary is presumed fixed, so that the size of the inner region remains constant. On the other hand, the position of the external boundary at any given instant of time is determined by the distance into the system that a sensible pressure reaction has occurred, Thus, the external boundary may change position with time. It initially emerges from the inner region and advances outward to its ultimate position. When this ultimate position coincides with a geometric limit, the reservoir system is said to be limited. When it coincides with points subject to pressure gradients furthest removed from the internal boundary, yet short of a geometric limit, the system is said to be unlimited. In this investigation two different boundary conditions are imposed at the ultimate boundaries of limited systems. The first requires that no fluid flow occur across this boundary; the second that the pressure remain fixed at this boundary.3-5 UNSTEADY-STATE FLOW In a strict sense virtually all flow phenomena associated with a reservoir system are unsteady-state. The transient behavior of these phenomena requires accounting, however, only when time must be introduced as an explicit variable. Otherwise, steady-state mechanics may be used. Analytically, steady-state conditions prevail in a reservoir system only over that portion of its history when this relation is satisfied:

Jan 1, 1967

By Robert D. Pehlke, Arden L. Bement

A mass transfer coefficient for the removal of hydrogen from liquid aluminum by inert flush degassing has been determined experimentally at 700°C. A mathematical model has been derived, assuming transport control, which adequately describes the experimental data. The removal of hydrogen is shown to increase with decreasing bubble size, and with increasing flow of flushing gas. The presence of hydrogen in liquid metals in amounts which exceed the solid solubility at one atmosphere hydrogen pressure is highly undesirable, causing unsound ingots or castings, or resulting in difficulties during further fabrication. The removal of hydrogen from liquid metals can be accomplished by several methods, the principal two being vacuum treating and inert flush degassing. The design of a process involving either of these techniques is based on the mass transfer coefficient for hydrogen between the liquid metal and a hydrogen-dilute gas phase. In an effort to extend our knowledge in this area, an investigation of the rate of removal of hydrogen from liquid metals by inert gas flushing was undertaken. In view of the difficulties encountered with gas porosity in the casting of aluminum and aluminum alloys, and because of certain experimental advantages, liquid aluminum was selected for the study. PREVIOUS INVESTIGATIONS The removal of gases from liquid metals, particularly steel, has received considerable attention in recent years, The first definitive work of a theoretical nature was presented by Geller1 who considered the case where equilibrium exists between purging gases and metal bath during flushing, and also deviations from this ideal case. Hulme2 has surveyed a number of experimental studies and discussed some of the practical limitations of removing dissolved gases from molten metals. More recently a number of studies have been reported on degassing of steel melts by flushing with an inert gas.3'4 Considerable work has also been reported on aqueous and organic liquid-gas systems, and an excellent summary prepared by Hovis.5 The treatment of liquid aluminum by a flushing gas to remove dissolved hydrogen has been known for several decades.6 A number of gases have been used for this purpose, including nitrogen, argon, and chlorine, or mixtures of these gases. Tikkanen and Erkko7 have studied the relative density changes of aluminum castings poured from melts degassed with chlorine, nitrogen, and mixtures of the two gases. Their

Jan 1, 1962

By Eli T. Conner

(Wilkes-Barre Meeting, June, 1911.) My connection, under a commission from the Councils and Board of School Control of the city of Scranton, Pa., with a recent investigation of mine-caves and the resultant damages to surface-improvements, has led to the preparation, at the invitation of our Secretary, of the present paper. It is notorious that there are, in the anthracite-fields, frequent subsidences of the surface, due to the removal of the coal beneath. No particular attention is paid to such occurrences, unless they happen to injure surface-improvements. Many caves have happened in Scranton and its vicinity, which have excited but little remark, since they have generally caused no serious damage. In August, 1909, a cave occurred in the Hyde Park section of the city, generally known as the West Side, which nearly destroyed school-house No. 16 and considerable adjacent property. 'Fortunately, there were no pupils in the building at the time; but the thought of the possible result, bad the usual number of pupils and teachers been present, aroused the public to the gravity of the situation; and the School Board employed engineers to investigate the case. Two reports were made by separate sets of engineers, which differed in some particulars. These differences were pointed out in the public press and magnified, and the consequent agitation of the subject was taken tip and greatly exaggerated by some of the metropolitan newspapers, giving to uninformed people the very erroneous impression that the whole city of Scranton was in danger of sinking into the bowels of the earth. All this tended to affect the credit of the city and to depress real-estate values. The matter was considered at a joint meeting of the Board of Trade, the Councils, and the School Board. The Hon. J. Ben. Dimmick, a former Mayor, suggested that it

Sep 1, 1911

By U. N. Bhrany, N. Arbiter

F. C. Bond (Consulting Engineer, Processing Machin-Dept., AllisChalmers Manufacturing Co., Milwaukee) This is a very comprehensive paper. It deals with 1) size distribution functions, 2) energy-particle size relationships, and 3) criteria for grinding efficiency. The method of analysis used consists of tracing the increase and decrease of successive Tyler scale screen fractions as grinding progresses. The treatment is based upon the assumption that the work done in grinding is concentrated entirely upon the coarse size fractions in the mill. It is postulated that when a feed consisting entirely of particles of one size fraction is batch ground that finer fractions increase at a constant rate as grinding progresses, until finally each in turn begins to be ground and to decrease in amount. It is shown that under these conditions the empirical Gates-Gaudin-Schuhmann log-log size distribution plot can be logically explained. It is difficult to imagine any condition of ball mill operation, or any probability analysis, which would result in no grinding contacts on the finer particles present. However, the authors' postulate must be accepted if their conclusions are to be sustained. They conclude that the equation expressing the for mation rate of each size fraction can be resolved into the G.-G.-S. size distribution equation, and that the curvilinear portion of the log-logsize distribution line represents the only particle sizes which are being ground as well as produced. The author's statement that the G.-G.-S. Eq. 4 applies only to the finer portion of the size distribution while the Rosin-Rammler Eq. 6 applies over the whole size distribution range certainly supports the view that an exponential equation of the Rosin-Rammler type corresponds more closely to actual size distribution than does the power law G.-G.-S. equation. However, the G.-G.-S. equation is used to derive Eq. 8 and Eq. 9, which state that the energy required to grind a ton through a particular mesh varies inversely as the particle size to the exponent -a , whic h is the slope of the extended G.-G.-S straight line. This is stated to apply to both batch and continuous grinding. The G.-G.-S. slope exponent thus becomes very important in all crushing and grinding analyses by this method. The authors include some data which can be used to test their conclusions, as described below: We will consider grinding a feed composed of 14 mesh (1190u) to -200 mesh (74u). The average of 17 Allis-Chalmers grindability tests on quartz gives a work index (Wi) of 12.77. The -200 mesh product should have 80 pct passing 41.6u or P = 41.6. The equivalent 80 pct passing size of the 14 mesh feed particles is considered to be the particle diameter, as is true for spherical particles. Thus F is 1190u. In the work index equation 10 wi 10 wi 127.7 127.7 P f 41.6 1190 = 16.1 kw-hrperton. In Table I the ball mill exponent for quartz is 1.0 and from Eqs. 9 and 10: w=kw-hr/ton =Constan/(do/di) = 16.1 ton so that (1190) Constant/ Constant and the Constant equals 259. The tons per horsepower-hour value is 0.0464, in agreement with Table IV. The exponents for other minerals are listed in Tabl e 11. In Table IIA, the kilowatt hours per ton required to grind these minerals between the same sizes as the quartz, both by the authors' method and the work index equation, are listed. According to the authors kw-hrperton = 259/(16.1)".

Jan 1, 1961

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By J. E. Hilliard

AS is well-known, the length per unit volume, Lv, of any lineal feature (such as grain edges in a polycrys-talline specimen) can be estimated from a count of the number of point intersections with a random plane of polish. Specifically:' Lv = 2NA in which NA is the expected number of intersections per unit area of the plane. As with all other measurements in quantitative metallography, the estimation of Lv is subject to a sampling error, the magnitude of which will depend on the number of observations. The analysis is usually performed by counting the number of intersections falling within a closed test figure of known area applied to a succession of different areas on the plane of polish. The test figure can be drawn directly on the screen of a projection microscope or else superimposed on the image of the structure by means of a reticule in a focusing eyepiece of an ordinary microscope. The use of micrographs is very inefficient unless, as in the present study, these are also needed for some other purpose. Consider an analysis in which n applications are made of a test figure within which the expected number of intersections is N. The expected total number of intersections counted is then given by NT = nN If (N) is the standard deviation in N for one application of the figure and if the latter is applied to non-overlapping regions so that the N's are independent in the statistical sense, then the standard deviation in Nt is given by: a(NT) = a(N)vn Eliminating n from Eqs. [2] and [3] we obtain: Thus, k will be a constant for a given specimen and test figure. From Eqs. [ I] and [4] we find As is to be expected, therefore, the fractional error in the estimation of Lv is inversely proportional to the square root of the total number of intersections counted. It is very desirable to have at least a rough estimate of k before starting an analysis, since it will then be possible to decide in advance by use of Eq. [5J how many counts should be made in order to achieve the desired accuracy in Lv. Without such information the investigator, unless lucky, will either do needless (and tedious) work or else obtain a result with too large an error. It is very easy to calculate k for the limiting case of N — 0, i.e., when the test figure is vanishingly small. In this case there is a negligible probability of two or more intersections falling within the figure. There are therefore only two possible outcomes: a count of one with a probability of, say, p and a zero count with a probability of (I-p). Thus and Therefore: We might expect k to depend on N and therefore on the size of the test figure. Such a dependence was observed2 in the analogous case involving the estimation of grain-boundary area per unit volume by a count of the number of intersections of a test circle with traces of the boundaries on the plane of polish. In this case the theoretical value of k at the limit N — 0 was 1.41. It was also predicted that k would display an initial sharp decrease with increasing N. Experimentally, k was observed to decrease to a value of 0.62 at N - 6 and thereafter remain approximately constant. Furthermore, the value of k for N > 6 was found to be insensitive to the form of the structure provided no pronounced anisotropy or grain-size variation was present. In order to determinewhether k would exhibit a similar dependence on N in the estimation of grain-

Jan 1, 1967

By H. M. Manasevit, R. L. Nolder, L. A. Moudy

The observation that silicon-on-sapphire crystallo-graphic relationships exist which are unrelated to one another implies that there are regions of sapphire orientations containing these relationships which are separated by boundaries. The presence of these boundaries is demonstrated, initially on rounded cylindrical surfaces md more definitively on flat sapphire substrates with orientations along primary crystallograPhic zones. By means of reflection electron and X-ray diffraction techniques the crystallographic relationships have been studied. Two new relationships have been found, viz., (111) Si (1014) sapphire and (001)Si (1012) sapphire. Also found were multicrystalline relationships. The approximated boundaries between the crystallographic regions for each relationship are shown by means of a stereo-graphic projection. These studies also clarify discrepancies in the literature concerning the silicon-sapphire relationships on the (1120) sapphire substrate. SINCE the first report of the deposition of single-crystal silicon on sapphire in 1963' this epitaxial system has been the object of study by many semiconductor laboratories. In particular, there has been strong interest in the crystallographic relationships between the two materials. In 19652,3 we demonstrated that at least four distinct modes of orientation relationships can occur, namely, (100)Si (0112) sapphire, (111)Si 11 (1124) sapphire, (111)Si (1120) sapphire, and (111)Si (0001) sapphire. Each relationship could be obtained when the substrate plane was at the particular orientation indicated. Other workers as well as our laboratory concurred in general with the initial results. However, two discrepancies were observed. First it was observed by Joyce et a1.4 that on a (1120) substrate a (110) silicon deposit could be obtained. This was actually consistent with the (111)Si (1124) sapphire mode. Second, two other laboratories,5,6 as well as our own,7 have found a deposition on the basal plane which is different than our initial findings. Even though a (111)Si deposit was obtained, the crystallographic relationships were not identical to those previously found. Subsequent work in our laboratories has confirmed that the form found in these later investigations of the basal plane is by far the more predominant, and the original finding was an anomaly. The relationships which have been previously estab1.ished are summarized in Table I. It is seen that at present a rather small number of possible substrate orientations has been investigated. Therefore, it seemed appropriate that a systematic investigation of silicon epitaxy on various sapphire substrate orientations would be of value. It was expected that the answers to three questions could be obtained from such an investigation. These are as follows: 1) Since we know of at least four distinct modes of orientation relationships, over what range of substrate orientation do these modes occur and what is the form of the boundary expected between the different regions in which each mode occurs? 2) Are there other modes which occur on substrate orientations not previously investigated? 3) What is the explanation for the discrepancy between our laboratory's findings and those results reported by Joyce et al.4 concerning the relationships occurring on (1120) substrates? APPROACHES AND EXPERIMENTAL METHODS Two different approaches were used. Initially, depositions were carried out on rounded cylindrical surfaces which lay perpendicular to the (0112) plane.* By this means the presence of boundaries and their approximate positions could be determined. Finally and more definitively, depositions on flat sapphire substrates with a variety of orientations, mostly along primary crystallographic zones and about the (1120) plane, were studied by means of reflection electron diffraction and X-ray diffraction (using back-reflection Laue patterns and full circle goniometer techniques) in order to determine crystallographic relationships. The sapphire orientations investigated by these means

Jan 1, 1969

By F. F. Skiba, J. E. Warren

An idealized miscible displacement process in a three-dimensional, heterogeneous porous medium has been studied via experimental computation based on a "Monte Carlo"' model. The macroscopic dispersion which results solely from variations in the permeability of the porous system is related to the scale of the heterogeneity as well as the distribution function for the permeabilities. The interpretation of laboratory data is often subjective and may bt? misleading since ambiguities are not always recognized. Furthermore, an experiment performed on a conventional oil-field core does not yield a valid measure of Macroscopic dispersion for reservoir engineering purposes since the scale of heterogeneity that is significant in the laboratory is not significant in the reservoir. The analysis of field-performance data is discussed and it is shown that a unique characterization of the reservoir can not be obtained from this information alone. A possible method based on a minimum of data, for estimating the spatial distribution of permeabilities in a reservoir is postulated INTRODUCTION Many studies concerned with diffusion and dispersion in porous media have been undertaken because of the potential importance of miscible displacement processes for increasing the recovery of oil. The results from a large number of these studies have been published; in a recent paper Perkins and Johnstonl presented a comprehensive survey of the pertinent literature and excellent appraisal of the "state of the art". In this investigation, dispersion resulting from macroscopic variations in the properties of the porous medium is of primary interest. Dispersion, in this sense, is caused by fluctuations in the velocities of the individual fluid elements as they move through the porous system. This phenomenon is the result of inhomogeneities in the permeability and/or porosity, but it can be amplified by differences in the characteristics of the fluids involved. The objectives of this study are the following: (1) to determine the qualitative manner in which macroscopic rock properties affect the observed dispersion coefficient; (2) to evaluate the possible influence of macroscopic dispersion on laboratory experiments; and (3) to appraise the significance of macroscopic dispersion with regard to reservoir performance. Investigations of this type are essential since many miscible displacement processes, while seemingly understood from a mechanistic point of view, do not behave in a predictable manner under field conditions. Although a plausible explanation is that the process is dominated by the environment, little serious effort has been devoted either to describing the reservoir itself or to determining the probable effect of the reservoir properties on the process. This study is intended to be exploratory in nature; however, it is hoped that it will provide additional insight into the physical problem and that it will suggest areas suitable for further investigation. EXPERIMENTAL COMPUTATION MODEL To limit the observed dispersion to that which results solely from variations in the properties of the porous medium, an idealized process in which one incompressible fluid displaces another identical fluid is simulated. This simulation is achieved by means of a modified Monte Carlo method — by tracking a number of mathematical particles (tracer) through the system according to the steady-state, single-phase velocity distribution. Then the system is characterized by the distribution of the residence times of the particles. Since any number of particles can be injected simultaneously (a true delta-function distribution) and only pressure conditions are imposed on the boundaries, the model is completely determinate. Furthermore molecular diffusion, microscopic dispersion, adsorption and/or reaction are automatically precluded from the model. The heterogeneous medium is represented by a rectangular, three-dimensional array of uniform, homogeneous porous blocks which are randomly arranged; it is assumed that any two opposing faces of the system are held at fixed potentials

Jan 1, 1965

By G. L. F. Powell, L. M. Hogan

Large degrees of undercooling have been produced in bulk samples, 400 g, of copper and Cu-O alloys by melting in a slag of commercial soda-lime glass. The maximum degrees of undercooling obtained for copper, hypoeutectic and hypereutectic Cu-O alloy samples were 208°, 218°, and 97°C, respectively. No grain refinement as a function of undercooling was observed for the pure copper samples although in Cu-O alloys there was a marked decrease in grain size at degrees of undercooling greater than 150°C. The grain size change is the result of recrystallization during or immediately following the freezing process. THE first report of a large degree of undercooling in a bulk sample of metal was that by Bardenheuer and Bleckmann1 who produced 258°C undercooling in a 150-g sample of iron by melting in a glass slag. Garbeck2 extended this technique to nickel, cobalt, and copper and reported maximum undercoolings of 220°, 230°, and 60°C, respectively, while Fehling and schei13 undercooled a large number of metals as 2 to 20-g samples slagged in glass. Subsequently, one of the authors4 applied the glass slag technique to bulk samples, -500 g, of silver and obtained a maximum undercooling of 250°C. It was concluded that impurities which normally catalyze nucleation at small degrees of undercooling could be removed by oxidation and solution in a glass slag. Large undercooling of bulk samples of iron, nickel, cobalt, and their alloys has also been reported by other authors.8'10'11 Of the five metals mentioned above, copper was the only one which had not been undercooled by at least 200°C in the form of a large bulk sample. Since copper exhibits high liquid solubility for oxygen, as do iron, nickel, cobalt, and silver, it was considered that a much larger degree of undercooling than that observed by Garbeck for copper should be possible. This was obtained by modification of the technique applied previously to silver and the results of undercooling experiments with 400-g samples of copper form part of this paper. The undercooling behavior of Cu-O alloys was also studied, and the influence of a small oxygen content on the grain structure of undercooled copper was observed metallographically. 1) EXPERIMENTAL The copper was obtained as oxygen-free high-conductivity copper, the major impurities of which were Fe 0.01 pct, Zn 0.015 pct, and Si 0.02 pct. Melting was carried out in an open-ended vertical cylindrical fur- nace wound with Kanthal wire. A schematic diagram of the experimental setup is shown in Fig. 1. The undercooling technique previously used for silver4 involved preoxidation of the samples by melting the silver in contact with the atmosphere. The oxygen content was subsequently reduced by freezing the samples under a glass slag, whereby the oxygen was evolved as gas and the continuous glass cover over the surface prevented re-solution of oxygen when the sample was remelted. Since oxygen is not released as gas when a Cu-O alloy is frozen. the technique had to be modified for use with copper. Also, copper is oxidized during fire refining so that an in situ preoxidation treatment of the samples was not considered to be a prerequisite to large undercooling. Experimentation proved that this supposition was correct. Instead, it was found that care was necessary during melting of the copper under glass to minimize oxygen pickup. Samples of copper weighing 400 g were prepared by adding 50-g pieces of copper to a vitreous silica crucible partly filled with soda lime glass at a temperature of approximately 1000°C. Each piece was quickly immersed in the glass slag and held at this temperature until the oxide coating on the surface decomposed. This change was easily observed and coincided with the formation of gas bubbles in the glass

Jan 1, 1969

By R. V. Higgins, A. J. Leighton

Interest by petroleum engineers in the flow of three phases—oil, gas and water—in irregularly bounded porous media lies mostly in the performance calculation of water floods of reservoirs that have been partially depleted as the result of expansion of much of the originally dissolved gas. The authors present a method to forecast three-phase flow in complex geometry and explain the details by the use of a specific example of a five-spot water flood of a partially depleted stratified reservoir. In this example, the fluid and rock properties of a field given by Prats, et al,' were used. The computed results fit the field performance more accurately than Prats, et al, and Slider.' The time required for the high-speed digital computer to make the calculations, including the contributions from the different layered zones, is about one minute. INTRODUCTION The declining rate of discoveries of new oil fields in the United States makes the recovery of more oil from known reservoirs more attractive than previously. Water flooding has an excellent proved background for recovering additional oil economically. Accordingly, the main interest of this paper is in this type of recovery. In petroleum engineering studies, commercial interest in three-phase flow is mostly in the water flooding of reservoirs in which the oil has been partially produced by the expansion of dissolved gas. When water is pumped into these reservoirs, three-phase flow takes place. Although this paper is concerned with these conditions, the principles involved could be used for other conditions when and where they come to the fore. Many investigators have made contributions to non-empirical forecast methods using basic scientific engineering principles. Several of these used the oil in place at the start of the water flood and the oil remaining after a large quantity of water has passed through a core for key values in their calculations. Recently, Prats, et al,' and Slider have reduced assumptions by adding the third phase—gas—in their forecasting methods. Slider uses the mobilities in the immediate vicinity of inlet and outlet wells as an aid to simulate the resistance effect in the five-spat pattern of the flow of oil, gas and water. Prats, et al. minimize assumptions by using previously determined laboratory data for determining sweep efficiency in the five-spot pattern. In neither of these papers is the saturation profile continuously affected by permeability-saturation curves. Sheldon and Dougherty- recently described a method that employs continuously changing saturation profiles using permeability curves, has a minimum of assumptions and needs no prior sweep efficiency. The Higgins-Leighton method, described in this paper, has all of these features; however, many of the techniques are different from those of Sheldon and Dougherty. The Higgins-Leighton method, tested in May, 1961, is direct and easy to apply and requires very little computer time to calculate a forecast. The short computer time is especially helpful in the study of a reservoir containing many layers of different relative permeabilities. In the Higgins-Leighton method, the individual pressures do not have to be calculated, as the resistance to flow in each cell in the flow pattern is readily determined without the use of any iterative techniques. The saturation and permeability distributions are readily determined. These data and the shape factor, which is measured only once from the potentiometric model when mobility ratio is one, determine the resistance to flow in each cell. THEORY The authors showed in a previous paper4 that, as an aid to calculating performance, the reservoir can be divided into channels using the streamlines of a potentiometric model as a guide. See Fig. 1. This procedure also was used in this paper. The authors also showed that, by treating conduits as approximately one-dimensional and neglecting pressure gradients transverse to the main flow, the Buckley-Leverett equation may be expressed as The principles expressed by this equation are employed extensively in the three-phase flow, as they were in two-phase flow. In the three-phase flow, the channels (the size and shape of which are taken from a potentiometric model)

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', 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' 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', 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' 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 Z. F. Oszter

This paper discusses the materials handling problems encountered in cement plants in the areas of raw and finished product grinding and cement storage. It illustrates by examples taken from actual installations, how optimum feed rate is maintained in grinding circuits with properly applied bin dischargers (vibrating hoppers, Mexican feed boxes and aerated bin bottoms) for materials with high moisture content and/or difficult handling characteristics. Initial operating problems with the equipment will be discussed and maintenance requirements will be indicated for comparison purposes. In a modern cement plant, properly functioning bins and feeders are as essential as process control, for no process control system, however sophisticated, can perform efficiently if raw materials do not flow into the plant processing facilities at the required time and rate. Most materials used in the manufacture of portland cement - such as silica sand, iron ore fines, crushed limestone, gypsum, clinker or cement itself-have difficult handling characteristics under specific conditions. High moisture content or even dampness due to condensation may create bin hang-ups, arching and ratholing, all of which result in improper functioning of the feeders or proportioning equipment. The occurrence of these problems in existing cement plants spurred a search for the corrective devices discussed in this article. The equipment investigated included bin dischargers in the form of vibrating hoppers for moist raw materials, Mexican feed boxes for clinker and aerated bin bottoms for cement. Operating experience with different types of bin dischargers showed that the actual performance of some did not come up to the expected level under prevailing operating conditions. The aim of this paper is to record the results of this experience and there- by aid designers and/or operators in selecting the most suitable equipment for their requirements. BIN DISCHARGERS (GENERAL) A discharger is a device that causes stored solids to flow out of a hopper without metering the rate of flow as a feeder does. The investigation discussed here encompasses the following types of bin dischargers used in a cement plant: (1) vibrating hoppers (horizontal motion), (2) whirlpool hoppers (horizontal and vertical motion), (3) aerated bin bottoms (air fluidization), and (4) Mexican feed boxes. Most of the bins used in the cement industry are of the plug flow type. Raw materials with difficult handling characteristics such as high moisture content may cause flow problems when fed through small openings (12 x 18 in.), even if two walls of the bin are vertical.' The most appropriate solution to these flow problems would be the installation of separate mass flow bins for the different materials such as silica sand, iron ore fines, crushed limestone and gypsum, with each bin correctly designed from test data obtained by measuring the flow properties of representative samples. However, limited storage capacity and plant space usually preclude this approach. Quantitative design of bins for flow based on measured flow properties of representative samples of the bulk solid to be handled is no longer an "engineering art," thanks to the basic research done on the flow of bulk solids in bins during the last decade.2, Where space limitations exist, however, or where improvements are being sought for equipment already installed, bin dischargers offer a satisfactory solution to flow problems. They can be applied effectively under flat-bottom bins or mounted to existing plug-flow bins with very little modification to these facilities. In both the foregoing cases, plug flow will be maintained, but the flow channel will be considerably increased so that it is larger than the discharger's bin-side diameter, and material flow can be sustained if the discharger is correctly selected. The main purpose of a bin discharger is to draw the bulk solids across the whole area of the bin outlet, to make the opening fully effective. As long as the bin-side opening of the discharger is large enough

Jan 1, 1969

By Robert Maes, Robert de Strycker

The Fe-Ni-As phase diagram has been established by the study of about a hundred alloys, by microscopic observation, and by thermal analysis, with arsenic contents up to 50 pct. The iron and nickel arsenides present extensive solid-solution fields, owing to the substitution of nickel by iron and vice versa; the extent of the solubility field of each compound has been determined with an accuracy of ±1 pct. In the investigated range of compositions, the solidification reactions were established, and the temperatures of the invariant reaclions detevmined with a preciston of iZ°C in the most favorable and ±5°C in the least favorable cases. The isothermal lines of the liquidus surface have also been drazum, with an accuracy estimated at i5°C. Reactions in the solid state, which take place for the formation or the decomposition of certain phases , were investigated in detail. ThE constitution diagram of the Fe-Ni-As system is fundamental for the understanding of the properties of nickel speiss; these by-products of the extraction of certain nonferrous metals indeed often contain the three elements iron, nickel, and arsenic as main constituents. The Fe-Ni-As system has already been the object of earlier investigations. In 1932, Guertler and Savels-berg' have presented some elements of the ternary phase diagram, for arsenic contents up to 55 pct (all percentages in this paper are given in weight percent), including the vertical section FezAs-Ni&s2. This investigation is, however, incomplete and certain anomalies suggested the necessity to verify the results published by these authors: the section Fe2As-Ni5AsZ, for example, is presented as quasi-binary, with a field of complete miscibility in the solid state, even though these compounds do not have the same structure. Recently, ~useck' established an isothermal section at 800°C in the Fe-Ni-As system by X-ray diffraction and microscopic examination of water-quenched alloys. These techniques are sometimes inaccurate for the determination of the fields where alloys are liquid at the investigated temperature, and they may lead to erroneous conclusions when a compound exists at the investigated temperature but is not stable at room temperature and cannot be maintained by quenching. LIMITING BINARY PHASE DIAGRAMS The Fe-As constitution diagram has been the object of several investigations which have been reviewed by Hansen and Anderko.3 For the Ni-As phase diagram, a recent study has been effected by Yund.4 In this system, Heyding and calvert5 have determined the existence of a compound of unidentified structure at arsenic contents slightly higher than those corresponding to Ni5As2 and at temperatures lower than about 200°C; by analogy with iron and cobalt arsenides, this compound could correspond to the formula Ni2As, as suggested by Kulle-rud; although this is not definitely established. In the region of arsenic contents from 35 to 55 pct, Fried-rich7 detected anomalies in the solidification reactions; an interpretation of these anomalies was given by Hansen ,8 which assumed that the solidification reactions in practice are not in equilibrium, but are meta-stable. The main features of the Fe-Ni constitution diagram are the existence of a complete miscibility field, at least at high temperatures, for the fcc phase (which will be designated by My in the remainder of this work) and of a limited solubility field for the bcc phase (designated by Ma). EXPERIMENTAL METHODS The fundamental technique used in this investigation was microscopic observation, which allowed the determination of the reactions occurring during solidification of the alloys, and possible reactions in the solid state.

Jan 1, 1968

By A. McLean, D. W. Gomersall, R. G. Ward

An experimental study has been made of the solubility of nitrogen in liquid iron and liquid Fe-C alloys using levitation melting and a rapid quenching device. Iron alloy droplets were equilibrated with nitrogen gas at 1 atm pressure, quenched, and analyzed. Previous techniques for studying the Fe-C-N system have produced data which me in marked disagreement. This disagreement is due largely to errors caused by reaction between the molten alloy and the crucible material. With the levitation procedure, errors from this source have been eliminated and precise solubility data obtained for temperatures between 1450° and 1750°C. C-N interactions in molten iron have been expressed in terms of first- and second-order free energy, enthalpy, and entropy parameters. ALTHOUGH the solubility of nitrogen in iron base alloys is in general small, the effects of nitrogen on the properties of steel may be quite profound. For most purposes nitrogen in finished steels is undesirable, particularly in the low-carbon grades, since on cooling to room temperature the solubility limit of nitrogen in the steel may be exceeded and this can lead to embrittlement and loss of ductility on aging. On the other hand, nitrogen can improve the work-hardening properties and machinability of steels while in certain stainless grades nitrogen is important in order to stabilize the austenite phase. It is, therefore, desirable that one should be able to predict the solubility of nitrogen in liquid iron alloys. To do this, information is required concerning the interactions between nitrogen and the various alloying elements which may be present in liquid iron. There have been several investigations of these effects in recent years1"7 and the interactions between nitrogen and many elements dissolved in liquid iron are now known to a high degree of precision at steel-making temperatures. Unfortunately, a number of iron alloy systems which are of interest in steelmaking have been difficult to deal with by the experimental techniques generally used for this type of investigation. Among the most important of these are the Fe-C alloys. In the past, two methods have been widely used for determining nitrogen solubilities: the Sieverts' technique, in which the amount of nitrogen required to saturate a given mass of liquid metal at a particular temperature and pressure is measured volumetrically, and the sampled-bath technique in which liquid metal held in a crucible is equilibrated with a gas phase containing a known partial pressure of nitrogen, and samples drawn from the melt are quenched and analyzed. These two methods have been discussed in detail elsewhere.5,8 With the Sieverts' technique, errors may be introduced from the following sources: i) Gas adsorption on metal films which have condensed on the cooler parts of the reaction chamber. ii) Uncertainty in the determination of "hot volume" calibrations. iii) Crucible-melt interaction, particularly if a gaseous reaction product is formed or if the melt becomes contaminated with material from the crucible walls. The sampled-bath method may also suffer from errors due to reaction between the melt and the crucible material. In addition, there is the possibility that gas may be lost from the sample during solidification and cooling. In the present investigation, the solubility of nitrogen in liquid iron alloys has been studied by means of a new technique based on the use of levitation melting equipment and a rapid quenching device. In addition to the fact that problems of the type outlined above are avoided, this particular approach has the following advantages: i) The high-frequency current induces vigorous stirring within the levitated droplet so that gas-metal equilibration is rapidly attained. ii) The gas phase surrounding the melt can be changed very quickly and is easily controlled. For example, a droplet may be levitated in helium, deoxidized in hydrogen, equilibrated with nitrogen, and quenched, within a period of 15 min. ii) Melts can be readily under cooled or superheated, thus extending the effective temperature range of an investigation and allowing temperature-dependent data to be determined with a high degree of precision. Excellent reviews of levitation melting techniques and their application to physical-chemistry studies at high temperature have been published recently by Jenkins et al.,9 Peifer,10 and Rostron.11 In the present investigation a levitation melting technique has been used to obtain data for the solubility of nitrogen in pure liquid iron and liquid Fe-C alloys at temperatures between 1450° and 1750°C. The solution of nitrogen in liquid iron can be described by the reaction:

Jan 1, 1969

By David B. Brooks

Identified potential resources of lead and zinc are lower, relative to reserves, than the same ratio for other metals, Either there is little material below present grades or, if it exists, there are no data on it. The author examines both these factors in detail, and also introduces a third possibility. In their 1961 monograph, The Future Supply of the Major Metals, Bruce Netschert and Hans Landsberg noted that identified potential resources (that is, known but submarginal resources) of lead and zinc were much lower relative to reserves (that is, known and supermarginal resources) than was the same ratio for other metals.' Put another way, it seemed that for such metals as iron, manganese, copper, and aluminum, slightly more favorable economic conditions or modest advances in technology would greatly increase reserves by permitting the exploitation of vast quantities of leaner or deeper or otherwise less valuable mineral deposits. But for lead and zinc, Netschert and Landsberg concluded that either: "(1) there is relatively little material below present cutoff grades; [or] (2) such material exists, but for some reason there are no data.'12 This situation, which they called "the lead-zinc anomaly'' has been subsequently discussed by other authors. The purpose of this paper is to carry the discussion further by examining the evidence in favor of (1) and (2) and by advancing a third alternative. There are at least two reasons for devoting time to the lead-zinc anomaly. First, while it is unlikely that we shall have to turn to alternative lead or zinc resources in the near or middle-term future, the longer run resource picture for the two metals is not so certain.3 Both within the United States and for the world as a whole, new discoveries will have to come at a fair pace to keep up with growing rates of consumption. Thus, it is important to see whether there exist lower grade or nonconventional potential sources of two such useful metals that we could turn to at some time in the future. The second reason for studying the lead-zinc anomaly is also related to resource adequacy, but is of more general significance. Depending upon whether the anomaly is technologic, geochemical, spatial, or temporal in origin, we shall be able to draw some conclusions about the resource picture for other metals. Can we expect the apparent pattern of lead-zinc resources to be repeated? If so, will technologic advances in either exploration methods or mineral processing be able to alter the pattern? Or is the pattern unalterably set by geochemical distribution? Most important, are metals characterized by what seems to be an anomaly, different in any way from metals for which low-grade or nonconventional sources already have come to be exploited, as the porphyry deposits have for copper. In sum, investigation of the lead-zinc anomaly represents a case study combining both physical and economic aspects that may open up a neglected area in evaluating mineral resource adequacy. WHAT IS THE LEALLZINC ANOMALY? There are two dimensions to evaluations of the potential of any mineral resource? The first involves the degree of certainty that we ascribe to the estimate: how much do we know about the size of the resource as a result of geological mapping, exploration, and mine development? This dimension is expressed by the terms measured, indicated, and inferred reserves (and, if you wish, 'undiscovered reserves'), all of which terms carry with them the implication of exploitability under current cost-price conditions. The other dimension to evaluations of a mineral resource involves the economic limits of the estimate: how will the size of the resource vary with different (usually more favorable) cost-price conditions? (This economic dimension should be understood broadly to include all factors bearing on the cost-price relationship. In this sense it subsumes technology because any resource is exploitable at some cost; the role of technology is to bring cost down.) The latter dimension is expressed in the terms potential or sub-marginal 'ore,' and is logically quite independent of the former. For example, we have relatively good measurements of the tonnage of manganese in low-grade deposits in this country, though they are too lean to be considered reserves; on the other hand, we are much less certain of our measurements of new lead and zinc deposits in Missouri, though there is no question but that there is a sizable reserve. Thus, it may be uncommon but is not a contradiction to talk of measured submarginal ore or inferred potential ore.

Jan 1, 1968

By J. L. Gillson

Historically, rock deposits and sand deposits of titanium minerals came into production about the same time, although there may be some argument as to what is meant by production. Beach deposits of heavy minerals in India (Figs. 1-4) and Brazil (Figs. 5) were worked for monazite about the turn of century, but as there was then no market for titanium minerals, these were thrown away. The rock rutile deposits at Roseland, Va., Fig. 6, were worked to supply rutile for titanium chemicals and for coloring ceramics long before there was a titanium pigment business. The pigment industry started about the middle twenties, both in Europe and the U. S., and almost simultaneously the rock deposits at Ponte Vedra Beach near Jacksonville, Fla., were worked for titanium content. Since those days, production from both types of deposits has continued to grow at a rapid rate; many deposits of both types have been found, and reserves have grown to very large figures. In total tonnage of reserves, there is no doubt that the rock deposits are far ahead of the sand deposits; nevertheless there is a very large tonnage of commercial sands available. It is the quality of titanium mineral in the sand and the relatively lower costs of operating sand deposits that have kept them abreast, at least in annual tonnage produced, with the rock deposits. The principal titanium mineral used is ilmenite, but as soon as that mineral began to be sought as a titanium ore, it was obvious that there are ilmenites and ilmenites. Textbook ilmenite should have the composition FeOTiO2 and should analyze 52.6 pct TiO2 and 36.8 pct iron as Fe. The Indian ilmenite, for almost a generation the standard ore for manufacturing pigment in the U. S., was found to analyze about 60 pct TiO, and only 24 pct. Fe, and most of the iron is in the ferric condition. The whole process of pigment manufacture in the U. S. was built up on the use of a raw material of that grade, and the American chemical engineers who operate the pigment plants shuddered at the thought of using a rock ilmenite with 45 pct or so of TiO, and nearly 40 pct Fe. Intensive search was made around the world to find other deposits of rich black sand, like the Indian beaches, but although a few were found, there was some objectionable feature about each. A deposit in Senegal, south of Dakar (Fig. 7), was worked for a while, but an organic coating on the grains made attack by acid difficult. Modern practice would have included a scrubbing operation, in a caustic soda bath, to eliminate the organic coating. Brazilian deposits were numerous, but individually small, and shipping from them difficult. Deposits on the east coast of Ceylon had many attractive features, but the ilmenite analyzed only 54 pct TiO2 and could have been used only with a penalty. Sand deposits with 2 pct ilmenite, like those now worked in Florida, would not have been considered commercial ore, even if they had been known at that time. Most rock ilmenites are associated or mixed with hematite or magnetite, which accounts for the lower titanium and higher iron values than in the sand ilmenites. The Norwegians, English, and Germans, with cheap Norwegian rock ore at hand, learned to install in their pigment plants adequate capacity on the black side, as it is calltd, and counterbalanced the extra cost of plant, and larger amount of acid used, by the lower cost of ore. When World War II arrived, two of the largest pigment manufacturers in the U. S. had to learn how to use the Adirondack ilmenite, but one of them very gladly went back to sand ores when the Florida deposits came into large-scale production after the war. The other continues to use Adirondack ilmenite and finds it commercially attractive to do so. Rutile is not a raw material for titanium pigment manufacture by the sulfate process, since it is insoluble in sulfuric acid. In addition to its small consumption in chemicals and ceramics it began to be used in quantity in welding rod coatings. With the outbreak of World War 11, and the tremendous need for welding rods in shipbuilding and other structural steel construction, rutile suddenly became in heavy demand. The sand deposits on the eastern shore of Australia (Fig. 8A) which had been worked in a small way since 1934 were brought into production, and some stream placers in Brazil were worked and rutile concentrates shipped to American

Jan 1, 1960

By M. A. Goudie

The deposits occur in a large salt basin of Middle Devonian age. The potash, the final deposit in the salt basin, results from several interrupted cycles of evaporation and dessication. The deposits are extensive, and, at first glance, relatively undisturbed. With more and more wells being drilled, it has now become evident that salt solution has played a large part in changing the original deposits, resulting in some cases in partial to complete removal of the potash and the underlying halite. The most dominant factor in the removal of salt by solution appears to have been tectonic movement and consequent faulting, probably of relatively minor dimensions but of major importance. Evidence which indicates the tilting of the evaporite basin to the north and northwest is shown by the changing pattern of the basin during succeeding eras of potash deposition. The potash minerals of most importance economically are sylvite and carnallite. Reserve calculations indicate that 6.4 billion tons of recoverable high grade potash in K2O equivalent exist in the basin. The Devonian salt basin, which contains the Saskatchewan potash deposits, extends from just east of the foothills in Alberta, north as far as the Peace River area, across Saskatchewan and into Manitoba as far east as Range 10 west of the First Meridian and south into Montana and North Dakota (Fig. 1). The basin is closed everywhere except to the northwest. The known potash deposits are confined almost entirely to the Province of Saskatchewan, with the exception of a small area in western Manitoba bordering the Saskatchewan boundary. The following discussion will concern only the Saskatchewan part of the basin. The evaporite series in the basin is defined as the Prairie Evaporite Formation of the Elk Point Group, of Middle Devonian age. Recent work done by potassium-argon dating methods has indicated an Upper Middle Devonian (Givetian) age of from 285 to 347 million years for the potash. The Elk Point Group consists in ascending order of the Ashern, Winnipegosis, and Prairie Evaporite Formations. The Ashern formation, with an average thickness of 30 ft, sometimes called the Third Red Bed, consists of dolomitic shales and shaly dolomites. The Winnipegosis, is a reef-type dolomite, usually with good porosity, and in many cases oil-staining, although to date no production has been obtained. The thickness varies from 50 to 250 ft. The Prairie Evaporite formation, varying from 0 to 600 ft in thickness, consists of halite with interbedded anhydrite and shale, with considerable amounts of potassium salts in the upper part of the formation. The potassium salts are chiefly chlorides, although very minor occurrences of sulfates have been re- ported. The anhydrite beds do not appear to be continuous, although generally one or two bands of anhydrite underlie the lowest potash zone and are used as marker horizons. The shale occurs as seams interbedded with the salts, as large irregular inclusions in the salts and as very fine particles in intimate mixture with the salts. The Prairie Evaporite Formation is overlain by the Second Red Bed member, the Dawson Bay Formation and the First Red Bed Member of the Manitoba Group, listed in ascending order. The Red Beds are shales which vary in color from red to green, maroon, grey, grey-black, and reddish purples. They serve as marker horizons for coring the potash. The Second Red Bed averages 14 ft in thickness, the First Red Bed 35 ft. The Dawson Bay Formation, which everywhere overlies the First Red Bed and the Prairie Evaporite Formation in the area under discussion, is a reef type of carbonate, in some places limestone, in others limestone and dolomite, with vugular to pinpoint porosity averaging 130 ft in thickness. In some parts of the area, it has a salt section near the top of the formation, usually with interbedded shales and limestones. In other parts of the area, it is waterbearing and the salt is absent. Detailed mapping has indicated that the areas in which the Dawson Bay is water-bearing are areas which have been disturbed by faulting. Where the Dawson Bay is salt-bearing, the porosity has been plugged by salt. The total thickness of the salt varies from between 600 to 700 ft in the center of the basin to zero at the northern edge of the basin (Fig. 2).* The salt-free area in the center of the Province is believed to have resulted from removal of salt by solution. Evidence from several wells suggests that salt removal has been a continuing process from the time of deposition to the present day. One well drilled between the Quill Lakes for potash information encountered

Jan 1, 1961