Search Documents

By R. Q. Shotts

Nitric acid oxidation rate analyses of three coals, previously studied microscopically by the Bureau of Mines, revealed three components. Relative quantities agree with those found for the four components given by the Bureau and results are consistent with current ideas of coal constitution. Possible multi-component composition for bright coal and a reactivity-rank relation are suggested. THE physically dissimilar components of bituminous coals often are easily recognized mega-scopically. Under the microscope, reflected light or light transmitted through thin sections reveals the presence of the different components, even when these are intimately mixed. Optical methods for the quantitative estimation of the relative abundance of the various components, both by means of thin sections and by particle count, have been fully described.'. ' It has long been recognized that there are chemical and physical differences between the various petro-graphic components of bituminous coals, although analytical differences usually are small.:'. ' Only in the case of fusain have chemical differences been used for quantitative determination of a component. C. C. Hsiao and associates, at the Mineral Industries Experiment Station of the Pennsylvania State College, have described a method of analysis which is based upon the differences in the rate of nitric acid (8N) oxidation, fusain, and the other components of coal."," The reproducibility of their method and its applicability in checking microscopic determinations of fusain content have been supported by several independent investigations.'. " The writer has proposed the use of differences in oxidizability for the estimation of other components." "' The results of the oxidation of whole coals and of float-and-sink fractions of coals were reported. In most cases the plots of the logarithms of the percent dry, non-fusain, organic residue from oxidation, against time, revealed the presence of at least two distinct components. Both components appeared to oxidize according to a first order law, but the reaction constants for the components were distinctly different. One or more of the dull density fractions were found to contain but one component, and some of the lower rank coals oxidized in such a way as to suggest the presence of three components. A suitable way to check the identity and significance of the components delineated by oxidation would be to analyze a sample of coal both by the nitric acid oxidation procedure and by a microscopic method. The writer was wholly unfamiliar with either of the microscopic techniques commonly used, and to make such a comparison it was necessary to rely upon microscopic analyses made by someone else. It is hoped that some laboratory which is equipped to make both types of analyses will some day make them upon identical samples. During the past 20 years, four Alabama coals have been analyzed petrographically and the results published by the United States Bureau of Mines. They are: 1—Flat Top mine, Mary Lee bed; 2—Empire mine, Black Creek bed; 3—Wylam No. 8 mine, Pratt bed, all in the Warrior field; and 4—Soot Creek mine, Fairview bed, in the Coosa field."-" Of these, only the Flat Top mine is still operating. Because of the closing of these mines, it first appeared necessary to rely upon the indirect and unsatisfactory procedure of sampling the beds in other mines located as near to the closed mines as possible. Upon investigation, however, it was found that the Bureau of Mines still had, in storage, the very same samples which had been used in the published petrographic studies. The Bureau very generously furnished about 2000 g each of the Pratt, Mary Lee, and Fairview bed coals, largely lumps but with some fines. The blocks of coal, when received, still were covered by the paraffin coating which had been placed on the polished surface, in the case of the Mary Lee coal almost twenty years ago. Procedure The procedure for oxidizing the coal sample and removing the alkali-soluble humic acid has been described. In the present study, oxidation periods of 1/6, 1/3, 1/2, 3/4, 1, 2, 3, and 4 hr were used. All oxidations were made in triplicate. After the paraffin had been removed in boiling water and the coal washed carefully with cold benzene, the entire sample of approximately 2000 g, obtained from the Bureau of Mines, was crushed to pass a No. 4 sieve. About 200 g of this material was pulverized to pass

Jan 1, 1954

By Louis B. Slichter

Prospecting is certainly the world's biggest and best gambling business. It is a game where the chips cost many thousands and where many millions, even billions, can be won. An attractive feature of this gamble is the fact that the players are free to rig the odds as favorably as possible. Only the limited bounty of nature, the restrictions of the laws of the land, the competition from other players, the limited sagacity of the player himself- only these, and other factors, restrict the possibilities for large winnings. Today probably about half the land surface of the earth remains to be searched for signs of ore visible at the surface. It is fortunate that so much land still remains for exploration by relatively economical methods. Looking to the future, it is a good guess that the costs of finding ore will increase in response to increasing technical difficulties. To continue to achieve profits from prospecting, it will become necessary to analyze the business more skillfully, in a more and more sophisticated manner. Whenever feasible, qualitative concepts must be replaced by their quantitative equivalents, improvements, and extensions.

Jan 6, 1960

By C. B. E. Douglas

In an interesting article on evaluation of the water problem at Eureka, Nev., W. T. Stuart describes how a difficult water problem, or one phase of it, may be evaluated by means of a small scale test. Test data are plotted by a method rendering, under certain conditions, a straight-line graph that can be projected to show how much the water table will be lowered by pumping at any specified rate for a given time. A formula is then used to determine the size of opening, or extent of workings, necessary to provide sufficient inflow to enable pumping to be maintained at that rate. At first glance this might seem the answer to a miner's prayer, but a word of caution is in order. It may not be the whole answer. Moreover, results obtained by the method described are reliable only for conditions approximating those assumed. Even where conditions do not meet this requirement, however, it may be possible to draw helpful inferences from the results, perhaps enough to facilitate another approach to evaluation of a problem.

Nov 1, 1955

By H. W. Schreiber, H. R. Babitzke, R. E. Rodriquez, Patterson J. A., C. D. Masters, J. J. Schanz, R. B. McCammon, E. A. Noble, P. M. Krishna

A five-year effort by the ASTM subcommittee E10.14 on uranium resource evaluation has culminated in the approval of E901-82, Classification System for Uranium Resources (ASTM Standard E901-82, April 1983). The standard classification was prepared with participation of representatives from the US Department of Energy, the US Department of Interior's Geological Survey and Bureau Mines, utilities, and consultants. The classification a resource is based on level of knowledge and the economic producibility of the resource at a given time. Alphanumeric identifiers are used rather than words to avoid confusion with existing terms, definitions, and classifications. An example demonstrates the dynamic nature of the process of classification and attests to the necessity of addressing both the aggregate needs of broad-scale resource planning and the specific needs of individual property evaluation. Problems that remain infixing the classification a given uranium resource include uncertainty in estimating the quantity of undiscovered resources and resolving differences that may exist in deciding when the drill-hole spacing is adequate to determine the tonnage and grade of discovered resources.

Jan 1, 1985

Earth Sciences: The standard of living of the world's people depends upon human initiative, investigation, and know-how in exploiting the earth's resources. The Earth Sciences-geology, mineralogy, geophysics and geochemistry, meteorology, geography-are basic sciences dealing with those problems concerning the origin, constitution, and utilization of the earth's materials. The development of these fields in one department provides the greatest opportunity for direct contact and stimulus for scientists and teachers working in related fields. In the Earth Sciences, geology and mineralogy are obviously closely related. As rocks are the materials of which the earth is made, so in turn minerals are the constituents of the rocks. Conversely, a study of earth structures and the interplay of physical and chemical forces in the earth is necessary in order to understand the origin and source of minerals. Building on the basic core of knowledge of earth processes, earth structures, and earth history, the geologist will continue to push out the frontiers of stratigraphy, structural geology, pale-biology and petrology with growing insight into principles that control surface and subsurface rock conditions. As the field of geology has advanced, several lines of specialization have developed. The economic geologist

Jan 1, 1950

By A. Dorsey Lyon

INTRODUCTION A GREAT DEAL of attention has recently been given to the metallurgy of zinc, and much of that which has been said and written on this subject has been in the nature of a criticism. For instance, we are quite accustomed to hearing it said that the metallurgy of zinc has stood still, that it is antiquated, out of date, etc. We are, however, inclined to agree With George C. Stone1 When he says that, although there have been no such spectacular and radical changes as Were worked in the metallurgy of copper and iron by the introduction of converters, there has nevertheless been a decided and marked progress in the metallurgy of zinc, and especially during recent years. However, it must be admitted that, due to the fact that the chemistry of zinc is so different from that of the other metals, the retort process is the only process by which zinc is produced commercially; and, moreover, that losses which occur in the zinc industry-that is, in the mining, milling, and smelting of zinc ore-are very great, for it is stated' that probably less than 50 per cent. of the zinc mined, even

Jan 7, 1914

Henry Carlisle: This is June, 1961, and my friend Arthur Bunker is on the other side of the tape recorder. I, Henry Carlisle, am pursuing the hobby of chronicling the interesting parts of the mining careers and the lives of well-known mining engineers. Now, Arthur, supposing we skip, for the time being, where you were born and brought up, and start right in with where you went to school. Arthur Bunker: I graduated from Taft School and then went to Yale, Sheffield Scientific School, and took a course in electrical engineering. I was quite fascinated by higher mathematics and physics, and this offered more in mathematics and physics than any other one of the engineering courses at that time. The school was only a three-year course, and so I graduated from there in 1916.

Jan 8, 1963

By M. R. J. Wyllie

The relationship between the electromotive force (E.M.F.) across a shale barrier and the concentrations of sodium chloride solutions on either side has been investigated. It is shown that the action of a shale barrier is analogous to a glass membrane separating two acid solutions of different hydrogen ion concentrations. The shale behaves as a sodium electrode and is responsive to the activities of the sodium ions in the two solutions in such a way that the potential can be calculated by means of the Nernst equation. This conclusion is confirmed by laboratory experiments. In a borehole the total E.M.F. of a shale cell is the algebraic sum of the ~otential across the shale and a boundary potential. The relationship between total E.M.F. and the resistivity ratio of two sodium chloride solutions is indicated for a number of formation temperatures. The E.M.F. thus predicted is then compared with the .elf potential read from an electric log and good agreement is demonstrated. Based on both the self potential and resistivity curves of the electrical log. a method is given for calculating connate water content in a bed having in-tergranular porosity and containing both connate water and hydrocarbons. INTRODUCTION The first paper on electrical well logging by C. and M. Schlumberger and E. G. Leonardon in 1934' attributed the self potential curve principally to streaming potentials, i.e. to electroki-netic effects. Almost immediately great difficulties were encountered in reconciling many of the curves they obtained with this interpretation. and a ~econd paper' by the same authors soon appeared. In this second paper self potentials were attributed to the combined effects of streaming potentials and electrochemical potentials, the electrochemical potential being considered the result mainly of the interaction of fluids of differing salt concentrations, i.e. a boundary potential, and partly of potentials set up at the faces of impermeable materials. Some experiments involving a gray clay for the impermeable material were quated. The Schlumbergers and Leonardon deduced from the equation for a simple boundary ~otential that the electrochemical potential, as opposed to the electrokinetic potential, could be expressed in the form E=Klog- .......1 pe where K is a constant, pm the mud resistivity. p, the resistivity of the connate water in a porous bed. However, no general expression for the constant K was obtained. Although the literature between 1934 and 1943 contains a number of quotations of their results, the valuable work of the Schlumbergers and Leonardon was not extended so that the electrochemical potential has been generally attributed wholly to boundary potentials between the mud in the borehole and the connate waters in porous formations. Unfortunately, however, the fundamental premise of all these papers, that a boundary potential can give rise to current flow in a borehole, is thermodynamically untenable. As will be shown. the fact that the electrochemical potential can be fairly accurately express as E = K log pm/pc, a form in which a boundary potential may also be written, is partly fortuitous. The boundary potential is indeed an integral part of the expression for the electrochemical potential in a horehole, but in magnitude it represents only about 20% of the total potential. In 1943 an important step in the elucidation of electrochemical potentials was made by Mounce and Rust3 who showed that if a wall of shale separated two compartments which contained saline solutions of different concentrations, and if the two solutions were themselves brought into contact in the pores of a porous inert membrane (such as unglazed porcelain) a current flowed through the shale and saline solutions. The direction of positive current was from the shale into the more dilute solution. The paper of Mounce and Rust, while repeating some of the observations of the Schlumbergers and Leonardon, seems to be the first to show that the shale was the seat of a genuine electrochemical effect capable of causing current flow. In the same paper Mounce and Rust pointed out the similarity between the fundamental conditions of their experiment and the conditions which existed when a bed of shale in the ground was simultaneously in contact with a porous sand containing saline connate water and mud fluid of salinity different from that of the water in the sand. Since it is now generally recognized that the S.P. curve measures ohmic potential changes in the mud fluid in the well bore resulting from changes in current flow, it is apparent that currents having their origin in the electrochemical interaction of mud filtrate and connate waters with shale beds are a very important portion of the total S.P. The work of Mounce and Rusta and others appears to indicate that, in general, the electrochemical portion of a particular kick on a S.P. curve far exceeds any electrokinetic potentials resulting either from streaming potentials or Dorn effects. The Dorn effect, or sedimentation potential. arises when small particles are allowed to fall through certain fluids under the influence of gravity. a difference of potential being observe? between two electrodes placed at different levels in the stream of falling particles. The Dorn effect is unlikely to affect seriously the S.P. curve as now measured. A successful analysis of the electrochemical aspects of the S.P. log should

Jan 1, 1949

By A. C. Emery

The physical and chemical properties for the metalliferrous feed materials to the Imperial Smelting Process blast furnace are defined and the means of achieving these properties by the sintering or briquetting routes discussed. Modern sinter plant performance data is included with a plant flow sheet showing draughting patterns. The necessity for ensuring an adequate SO2 content in sinter machine gases when roasting sulphide materials to enable conversion to sulphuric acid in a contact acid plant is described. The ability of the process to consume limited quantities of oxide materials, residues from other processes and other secondary sources is discussed. The argument for a briquetting route is developed out of this limitation and the merits of cold and hot briquetting examined. Development work on the hot briquetting route for this purpose is briefly described and several flow sheets outlining possible applications for hot briquetting to different raw material situations including the preparation of low lead charges for the furnace are presented.

Jan 1, 1977

By T. B. Counselman

PRobably the earliest concentration of iron ore in this country was carried on in the northeastern magnetite areas. Magnetic concentration was relatively simple and gave a concentrate that, after agglomeration, and in spite of beillg high in phosphorus hecause of the dry concentrating methods then in use, was satisfactory blast-furnace feed. Production from the northeastern magnetitc area languished after the discovery and development of the Lake Superior ores, where there was seemingly an unlimited tonnage of high-grade ore, first the Michigan ores, and subsequently the still larger deposits of Minnesota. For years only the higher grade, direct shipping ores of the Lake superior district were mined. Those of Michigan, and of the Vermilion Range in Minnesota, were hard ores and came principally from underground mines. The softer ores of the Mesabi came partly from underground mincs, but mostly from the enormous deposits that could be worked by open-pit methods. The costs of ore from these open-pit deposits, even with the high initial cost of stripping, were so low that the underground mines of the Mesabi could hardly compete. The so-called Old Range ores, of Michigan and of the Vermilion Range of Minnesota, were hard and lumpy, and opened UP the burden in the blast furnace, therefore making possible harder driving with greater output. Also, some of these Ores were of such grade and physical characteristics that they were suitable for charge ore in the open hearth. These ores, therefore, commanded a price premium. and could compete with the softer, finer, but cheaper, Mesabi ores. The big tonnage reserves, however, are on the Mesabi. and it is these big open-pit deposits that contitue- our stock pile of iron ore, from which we are drawing so heavily during this war period. On the western end of the Mesabi the ore was found to be high in silica, and not suitable as furnace feed. As early as 1902, attempts were made to concentrate this so-called wash ore and by 1910 the first washing plant was put into operation. The Process was very simple, consisting essentially of a sizing and washing operation, with scalping screens, logwashers and tables, which later were replaced by classificrs. Ore particles finer than I00 mesh were deliberately wasted, because these fines were undesirable in the furnace. Fortunately, these tailings were impounded, and already are being reworked. Before long some of the low-grade ores began to require more elaborate treatment than the simple washing operation and, tables having been discarded because of their low capacity, in favor of classifiers, jigs were next introduced. These were able to discard coarser pieces of tailing than classifiers. Jigs are in active use today, various improvements having been made in their design and application. Various other machines and methods of concentration have been tried from time to time, and have in turn been discarded for more improved processes., The latest and most important process that has been 38

Jan 1, 1944

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By T. B. Counselman

PRobably the earliest concentration of iron ore in this country was carried on in the northeastern magnetite areas. Magnetic concentration was relatively simple and gave a concentrate that, after agglomeration, and in spite of beillg high in phosphorus hecause of the dry concentrating methods then in use, was satisfactory blast-furnace feed. Production from the northeastern magnetitc area languished after the discovery and development of the Lake Superior ores, where there was seemingly an unlimited tonnage of high-grade ore, first the Michigan ores, and subsequently the still larger deposits of Minnesota. For years only the higher grade, direct shipping ores of the Lake superior district were mined. Those of Michigan, and of the Vermilion Range in Minnesota, were hard ores and came principally from underground mines. The softer ores of the Mesabi came partly from underground mincs, but mostly from the enormous deposits that could be worked by open-pit methods. The costs of ore from these open-pit deposits, even with the high initial cost of stripping, were so low that the underground mines of the Mesabi could hardly compete. The so-called Old Range ores, of Michigan and of the Vermilion Range of Minnesota, were hard and lumpy, and opened UP the burden in the blast furnace, therefore making possible harder driving with greater output. Also, some of these Ores were of such grade and physical characteristics that they were suitable for charge ore in the open hearth. These ores, therefore, commanded a price premium. and could compete with the softer, finer, but cheaper, Mesabi ores. The big tonnage reserves, however, are on the Mesabi. and it is these big open-pit deposits that contitue- our stock pile of iron ore, from which we are drawing so heavily during this war period. On the western end of the Mesabi the ore was found to be high in silica, and not suitable as furnace feed. As early as 1902, attempts were made to concentrate this so-called wash ore and by 1910 the first washing plant was put into operation. The Process was very simple, consisting essentially of a sizing and washing operation, with scalping screens, logwashers and tables, which later were replaced by classificrs. Ore particles finer than I00 mesh were deliberately wasted, because these fines were undesirable in the furnace. Fortunately, these tailings were impounded, and already are being reworked. Before long some of the low-grade ores began to require more elaborate treatment than the simple washing operation and, tables having been discarded because of their low capacity, in favor of classifiers, jigs were next introduced. These were able to discard coarser pieces of tailing than classifiers. Jigs are in active use today, various improvements having been made in their design and application. Various other machines and methods of concentration have been tried from time to time, and have in turn been discarded for more improved processes., The latest and most important process that has been 38

Jan 1, 1944

By R. W. Heckel

A method is described whereby the relationship of both the "at-pressure" powder compact density and the "zero-pressure" compact density to the applied pressure may be obtained from continuous measurements of punch movement during a single compaction operation. The rate of density increase with applied pressure for the metal powders is found to be proportional to the volume fraction of pores for pressures exceeding a lower limit which varies from 15,000 to 30,000 psi, depending on the powder. 1 HE compaction of powdered materials is carried out primarily to increase the density of the material. If the ultimate goal of the over-all process is the attainment of minimum porosity, compaction is responsible for most of the densification. For example, loose powdered metals have porosities of about 65 to 75 pct. Axial loading in a die or hydrostatic pressure may effectively reduce the porosity of most powders to 20 pct or lower with pressures of 50,000 psi and greater. The relationship between the density of a powder compact and the pressure needed to achieve that density is most commonly obtained by pressing several compacts, each at a different pressure. Measurement of the densities of the individual compacts when they are removed from the die then provides the density-pressure relationship. However, this method suffers from several disadvantages. Because a large number of compacts and pressing operations are required, it is inherently slow. Furthermore, since density measurements are made after the specimens are removed from the die, data obtained refer only to pressures greater than those necessary to form a coherent compact. The third disadvantage is the inability of this method to provide information about the compact density at the applied pressure in the die. Because the determination of density-pressure relationships has been a slow, tedious process, investigators have sought to express the compaction behavior of powders by other means. The general terms "comparability" and L'compressibility" have been used to rank powders qualitatively according to their compaction characteristics. To restore a quantitative nature to the subject and to prevent confusion, Schwarzkopf1 has proposed that compac-tibility be defined as the minimum pressure needed to produce a given green strength, while compressibility should be used to indicate the extent to which the density of a powder is increased by a given pressure. Probably the most widely used of the compaction parameters is the "compression ratio,'' which is generally defined as the ratio of the compact density obtained by pressing at a given pressure to the apparent density of the loose powder.2 However, the description of the compaction behavior of a powder by these parameters provides only limited information about the process because of their applicability to just one specific condition such as a given green strength or a given pressure. The method of obtaining density-pressure relationships which will be described in the present paper was designed to eliminate the shortcomings of the previous techniques while maintaining equivalent accuracy and precision. EXPERIMENTAL TECHNIQUES The general principle employed makes use of the fact that the linear movement of the punch during a

Jan 1, 1962

By M. R. Tek, R. L. Nielsen

The scaling laws as formulated by Rapoport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions. Two hypothetical systems — a gas reservoir subject to water drive and the laboratory scaled model of that reservoir — were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these par ticular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system. This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. INTRODUCTION Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models. To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoportl and others.2-5 Rapoport's scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems. It will be recalled that these scaling laws presuppose three conditions: (1) the relative permeability curves must be identical for the model and prototype; (2) the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; and (3) boundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because: (1) The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media. (2) The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent. (3) Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible. The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible. The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. The numerical treatment employed was similar to that of Douglas, Peaceman and Rachford6 and it included both capillary and gravitational forces as

By Walter Crafts, D. C. Hilty

The liquidus diagram for the iron field of the Fe-S-O system has been derived experimentally. The solubility of oxygen in molten Fe-S alloys has been measured at several temperatures and found first to decrease slightly and then increase rapidly with increasing sulphur content. IT is well-known that the types, distribution, and amount of inclusions in any given steel are closely allied with the oxygen and sulphur contents of the metal, so that the term "deoxidation" as conventionally applied in steelmaking refers not only to the control of oxygen and oxide inclusions but to collateral modification of sulphide inclusions as well. Inclusions commonly occur in steel as mixtures of oxides and sulphides, sometimes closely associated and sometimes apparently completely divorced. Their effects on the properties of the metal may vary widely depending on their types and distribution as well as on their quantity. The mechanism by which inclusions form, however, and the means for controlling it are exceedingly complex and not well understood. Nevertheless, as described by Benedicts and Lofquist' and by Wentrup' all oxide and sulphide inclusions, excluding mechanically entrained matter, must result from specific modifications of the basic equilibria of the Fe-S-O system. In an investigation of the mechanism of deoxidation and inclusion formation being carried out at the Union Carbide and Carbon Research Laboratories, it became evident that quantitative data for the fundamental Fe-S-O system were essential to a better understanding of the problem. Consequently, the liquidus surface of the portion of the Fe-S-O diagram of interest in steelmaking has been determined, and the solubility of oxygen in liquid iron containing sulphur has been measured. Fe-S-O Diagram The portion of the Fe-S-O diagram of specific interest in steelmaking is the Fe-FeS-FeO corner. Of the contiguous binary diagrams for this partial system, only those for Fe-FeO and Fe-FeS are known with reasonable certainty.3 Several attempts to derive the FeO-FeS diagram have been made, chief among them being the work of Giani as reported by Oberhoffer,4 and, more recently, that of Vogel and Fulling." Although both Giani and Vogel and Fulling published diagrams for the FeO-FeS system, their observations indicated that the oxide component approximated Fe3O4 rather than FeO, so that there is doubt whether the system is truly binary. Benedicts and Lofquist and also Wentrup published qualitative sketches of the ternary diagram based on the assumption that it was of the type described by Marsh for a ternary system in which two of the contiguous binary systems are simple eutectic systems and the third exhibits a solubility gap in the liquid. Apparently, however, the only previous effort to establish the diagram experimentally was that of Vogel and Fülling whose results are reproduced in Fig. 1. Although defining approximately the limits of the miscibility gap and the location of the ternary eutectic, point E in Fig. I, Vogel and Fülling's diagram must be considered only semiquantitative. It was derived from a limited number of observations mainly along two quasi-binary sections, and the directions of the tie-lines and the location of the critical point, K in Fig. 1, were estimated by computation. Experimental Procedure The experimental technique employed for the present investigation involved melting mixtures of ferrous sulphide and ferric oxide in contact with iron. Since it was considered that in this part of the system the vapor pressures of sulphur and SO, are so low that they would not affect the results significantly, the melting was carried out in a current of argon to sweep out products of any initial re-

Jan 1, 1953

By R. H. Jacoby, V. J. Berry, R. C. Koeller

The experimental phase behavior of several field gas-condensate systems, one field volatile oil system, and a series of synthetic systems having gas-oil ratios from 2,000 to 20,000 scf/bbl stock tank oil was measured. The data were used to calculate the depletion performance of the field systems at their respective reservoir temperatures and of the synthetic systems at various temperatures. The results of the performance calculations were used to prepare correlations of total stock tank oil and separator gas in-place, and the amount of each recovered by primary depletion. Data needed to use the correlations are initial gas-oil ratio, initial tank oil gravity, reservoir temperature and reservoir pressure. The correlations presented should be useful for estiinating recovery from reservoirs producing volatile oils or rich gas condensates having gas-oil ratios from about 2,000 to 30,000 scf/bbl of stock tank oil. INTRODUCTION Prediction of depletion performance for gas-condensate reservoirs usually requires extensive laboratory work or lengthy phase-behavior calculations. For many small gas-condensate bearing zones the expense of such work is often not justified. Furthermore, where large holdings are involved, it is useful to have estimates of future oil and gas production in advance of the results of a laboratory analysis, even if such work is ultimately planned. For these reasons, it would be of value to have a method for rapidly estimating the depletion performance of gas-condensate reservoirs from the type of field data usually available. Few data on the reservoir depletion performance are available for the purpose of developing correlations. The increasing number of gas-condensate discoveries and the advent of high-speed computers have stimulated the development of methods for predicting performance. Such calculated performance data have been shown to be reliable and are becoming accepted as a basis for reserve estimates and for engineering studies. It therefore seems satisfactory to use results of performance predictions themselves for the purpose of developing the desired correlations. Such an approach is described in this paper. Several rich natural gas-condensate systems (gas-oil ratios of from about 3,600 to 60,000 scf/bbl) and one natural volatile oil system were studied (2,363 scf/bbl). In addition, to provide a more systematic basis for correlating purposes, a series of related, synthetic reservoir fluid mixtures was also studied (GOR's from about 2,000 to 25,000 scf/bbl). The purpose of this paper is to show how the recovery of oil and gas by primary depletion may be estimated from surface gas-oil ratio, reservoir temperature and pressure, and tank oil gravity. procedure Experiments Companion samples of gas and oil were collected from the high pressure separators in nine different fields at the time of their discovery. These samples were each recombined in accordance with the respective producing GOR to reproduce the reservoir fluids. Eight of these fluids, designated F-2 through F-9, were rich gas-condensate systems at the respective reser-

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By M. P. Tixier, M. Gondouin, G. L. Simard

In the quantitative interpretation of the SP logs, the electrochemical component is generally taken equal to — K log Rm /Rw where K has the theoretical value corresponding to solutions of pure sodium chloride. This method may be misleading when relatively large quantities of salts other than NaCl are present in form-rrtion waters—as is generally the case for low salinities —or in gyp-base muds. In such cases recourse is made in the field to changing the K value or to adding a correction term to the equation on the basis of local experience. An investigation has been made in the laboratory of the influence of HCO3, SO4=, Ca++, and Mg++ on the amplitude of the SP deflection, which included in par-ricular determinations of the activity coefficients of Ca++ and Mg++. The theory, the experimental techniques and a tentative method for applying the results are described. So far, applications to actual field cases where chemical analyses of waters were available have provided an excellent confirmation of the proposed method. The case of the SP in very salty brines is also considered and it is shown that the use of activity data for sodium gives satisfactory results. HISTORICAL In early times of electrical logging, the SP curve was used exclusively as a tool for the location of permeable beds and the definition of their boundaries. Later, with the introduction of methods of quantitative analysis, attention was called to the possibility of deriving from the SP log some information on the formation water resistivity, which constitutes one essential element for the computation of water saturation from log data. The first laboratory experiments' in 1932 had shown that, as a result of electrochemical phenomena arising at the contact of the niud and the formation, the amplitude of the deflection of the SP curve was, at least partly, a function of the salt concentration of interstitial water and accordingly of its resistivity. It was indicated that the amplitude of the electrochemical component of the SP should be equal to K log R1/R2, R1 and R2 being the resistivities of the mud and of the formation water respectively. However, the exact value of the K coefficient was not determined at that time because the surface clay used for a laboratory model was far from an ideal shale membrane. Consequently, the way was not yet open towards quantitative interpretation. In the early forties, Tixier reported2, on the basis of field log observations, the following equation for the case of NaCl solutions, SP = -71 log f1c1/f2c2, c1 and c2 being the concentrations and f1, f2 the mean activity coefficients. The local variations observed in the empirical values of K, chiefly in the Rocky Mountains, led him to advance the opinion that, besides the relative amounts, the nature of the salts present in the solutions had a bearing on the magnitude of the SP. At about the same time Mounce and Rust- ad published laboratory results which emphasized the part played by the shale in the generation of the SP. On the basis of laboratory tests and theoretical considerations, Wyllie in 1948' published the following equation of the electrochemical SP for the case of NaCl solutions, SP= - 71 log (a1)m/(a2)m

Jan 1, 1958

By Paul Billingsley

The Georgetown mining district is located in Deerlodge county, Mont., about 20 miles west of Anaconda. It lies along the divide between the headwaters of Warm Springs creek, draining eastward, and Flint creek, flowing west and north through the Philipsburg valley. The actual divide is formed by Cable mountain, a narrow ridge running north and south, terminating to the south in low hills overlooking Georgetown lake, and to the north in a high, rugged plateau that falls off to the west to the wide valley of Flint creek. On the eastern, or Warm Springs, slope of Cable mountain is the Cable mine, with the Hidden Lake mine a few miles north. On the western slope is the Southern Cross group, with smaller mines, the Twilight, Cincinnati, Red Lion, and Big Bill, at intervals to the northward. The Georgetown mines lie in the hills south of Cable mountain; the Philipsburg mines, of which the Granite Mountain, Bimetallic, Hope, and Trout are the chief, form a separate group to the north, overlooking Flint Creek valley and Philipsburg. (See Fig. 1.) The region ranges in elevation from 5,000 to 9,000 ft. It is, with the exception of the lower valleys, well timbered with small fir and pine, although Cable mountain and the higher summits to the northward rise slightly above timber line. The main valleys have the smooth contours of glaciation, and their mouths are frequently choked with morainal deposits. General Geology. In a very general way the valley of Flint creek coincides in this region with the boundary between the pre-Cambrian slates of western Montana and the folded Palaeozoic rocks of the Rocky Mountain system. The great limestone formations of the Carboniferous and Devonian, lying at flat dips, extend from Anaconda well toward the head of Warm Springs creek, but on the divide the Cambrian appears, closely folded and intruded

Jan 1, 1914

By C. B. Manula, R. Venkataramani, T. V. Falkie

A linear programming model has been developed to solve a production planning problem for a multi-plant stone operation serving several customers. The factors considered are production capacities, inventory restrictions, demand requirements, and available stone reserves. The objective function to be optimized is the total annual net returns from the quarry operations. For the purposes of model testing, annual stone sales from each plant were divided monthwise, productwise, and assigned to a market. The concept of assigning demands to market sectors and not to any particular plant creates the self-correcting interplay between the plants and the markets. Along with the inventory restrictions, and ending inventory requirements, this interaction enables the model to minimize the external diseconomics. The model was tested against an operating prototype and a number of alternative problems were also considered. Finally, the basic linear programming model has been enlarged by the addition of another dimension and the problem has been derived in a decomposition format.

Jan 1, 1972