Search Documents

By R. P. Abendroth

The hydrogen-water vapor ratio at which Fe-Cr alloy, chromium oxide, andiron chromite coexist in equilibrium was determined between 900" and 1200°C. A thermogravimetric method was used to determine equilibrium conditions. The results fit a straight-line relationship in the temperature region studied, and are given by Reduction experiments were also performed to confirm the results of the equilibrium investigation. ThE oxygen pressures at which Fe-Cr alloy, chromium oxide, and iron chromite coexist in equilibrium have been previously determined by Boericke and angert,' Morozov and Novokharski,' and Katsura and uan. Only one determination (at 1300°C) was made by Katsura and Muan, but it agrees with the results of orozov and Novokharski. The results of Boericke and Bangert, however, differ appreciably from the results of these investigators. Previous studies have assumed that the equilibrium metallic phase is pure iron, but Dahl and Van vlack have shown that the iron contains from about 1 wt pct Cr at 1000°C to over 2 wt pct above 1300°C. The chromium oxide also contains a small amount of iron in solid solution. In the present study, hydrogen-water vapor mixtures were equilibrated with the condensed phases, using a therrnogravimetric method to determine equilibrium conditions. The reaction can be written EXPERIMENTAL General Procedure. The starting material was a sintered pellet of Fe2O3-Cr2O3 solid solution with a hole in the center, and was placed on a fused silica hook. This assembly was raised into the preheated hot zone of the furnace in a helium atmosphere, hooked onto a fused silica hangdown suspended from one arm of an Ainsworth Model RV-AU-1 recording balance, and the starting weight determined. A flowing hydrogen-water vapor atmosphere was then exchanged for the helium by evacuation, and the sample reduced until the weight loss indicated the sample composition to be in the alloy-Cr2O3-chromite field. The tem- perature was adjusted incrementally until constant sample weight was achieved for several hours, to within 0.02 mg. A hydrogen-water vapor atmosphere of different composition was then admitted, and the same procedure carried out. At the end of a series of determinations, the sample was examined by X-ray diffraction to verify the presence of the desired phases. Microscopic examination of the silica hook showed no interaction with the sample, nor did it lose any weight. Several criteria were used to insure equilibrium besides constancy of weight. For a given hydrogen-water vapor composition, equilibrium was approached from both oxidizing and reducing sides by varying the furnace temperature slightly. The resulting slow weight loss or gain was observed for several hours. Constant weight could be re-established by returning to the original furnace temperature. The last criterion used was varying the relative amounts of the phases by further reduction or oxidation, and observing any changes in temperature required for constant weight for a given hydrogen-water vapor atmosphere. None were observed. This procedure was essentially the same as approaching the equilibrium from oxidizing and reducing sides, but larger weight excursions were carried out. Sample Preparation. Reagent-grade Fe2O3 and Cr83 powders were mixed in the desired proportions and heated in air at 1250°C for 2 hr. The mixture was re-ground and heated in air overnight at 1250°C. X-ray diffraction showed complete solid-solution formation as a result of this procedure. The solid solution was then pressed into l/2-in.-diam pellets using Carbowax 4000 as a binder. The hole was drilled in the center, and the pellets were sintered 24 hr at 1250°C in air on a bed of Fe2O3-Cr2O3 of the same composition, contained in an alundum boat. After cooling, the pellet surfaces were abraded with 310 paper to remove any surface compositional differences, such as loss of Cr2O3. Chemical analysis of the sintered pellets was 67.16 wt pct CrP3 and 33.02 wt pct Fe203. Atmosphere Generation and Control. The hydrogen-water vapor atmospheres were generated by passing Matheson ultrahigh-purity hydrogen, with no further purification, through two water bubblers contained in a constant-temperature water bath. Since the water-vapor dew points required in this study were below room temperature, the bath was insulated, and was cooled by thermoelectric-immersion devices. The bath temperature was controlled to 0.0l0C. Since rather high flow rates of about 900 ml per min were used through the furnace tube, an independent check of the dew point was made to insure saturation of the hydrogen by the water vapor. Although the dew point could only be determined to within 1/2"C, the determined dew points agreed with the water-bath temper-

Jan 1, 1967

By G. Rowan, M. W. Clegg

The basic equations for the flow of gases, compressible liquids and incompressible liquids are derived and the full implications of linearising then discussed. Approximate solutions of these equations are obtained by introducing the concept of a disturbed zone around the well, which expands outwards into the reservoir as fluid is produced. Many important and well-established results are deduced in terms of simple functions rather than the infinite series, or numerical solutions normally associated with these problems. The wide range of application of this approach to transient radial flow problems is illustrated with many examples including: gravity drainage of depletion-type reservoirs; multiple well systems; well interference. INTRODUCTION A large number of problems concerning the flow of fluids in oil reservoirs have been solved by both analytical and numerical methods but in almost all cases these solutions have some disadvantages — the analytical ones usually involve rather complex functions (infinite series or infinite integrals) which are difficult to handle, and the numerical ones tend to mask the physical principles underlying the problem. It would seem appropriate, therefore, to try to find approximate analytical solutions to these problems without introducing any further appreciable errors, so that the physical nature of the problem is retained and solutions of comparable accuracy are obtained. One class of problems will be considered in this paper, namely, transient radial flow problems, and it will be shown that approximate analytical solutions of the equations governing radial flow can be obtained, and that these solutions yield comparable results to those calculated numerically and those obtained from "exact" solutions. It will also be shown that the restrictions imposed upon the dependent variable (pressure) are just those which have to be assumed in deriving the usual diffusion-type equations. The method was originally suggested by Guseinov,1 who (a) postulated a disturbed zone in the reservoir, the radius of which increases with time, and (b) replaced the time derivatives in the basic differential equation by its mean value in the disturbed zone. In this paper it is proposed to review the basic theory leading to the equations governing the flow of homogeneous fluids in porous media and to consider the full implications of the approximation introduced in linearising them. The Guseinov-type approximation will then be applied to these equations and the solutions for the flow of compressible and incompressible fluids, and gases in bounded and infinite reservoirs obtained. As an example of the application of this type of approximation, solutions to such problems as production from stratified reservoirs, radial permeability discontinuities; multiple-well systems, and well interference will be given. These solutions agree with many other published results, and in some cases they may be extended to more complex problems without the computational difficulties experienced by other authors. THEORY In order to review the basic theory from a fairly general standpoint it is proposed to limit the idealising assumptions to the minimum necessary for analytical convenience. The assumptions to be made are the following: 1. That the flow is irrotational. 2. That the formation is of constant thickness. 3. Darcy's Law is valid. 4. The formation is saturated with a single homogeneous fluid.

By E. J. Harris

Mathematical solutions to complex mine ventilation problems are possible, but often the airway network is so complex that the mathematical solution becomes tedious and impractical. A fluid network analyzer, designed and built for analyzing mine ventilation problems has been in service at the Bureau of Mines' Pittsburgh station for approximately nine years. Using this analogue as a model, the mine ventilation network is simulated electrically by a combination of series and parallel circruits laid out to conform with actual mine airways. A tungsten filament lamp, referred to as a Fluistor, is used to simulate mine airway resistance. With the analogue ventilation model established, voltages of the proper amplitude to represent ventilating pressures are impressed across the circuit at points where mine fans and airshafts are located. Comparison of mathematical solutions of complex systems against analogue results showed a maximum variance of 3% for pressures and 2% for quantities. An electrical fluid network analyzer, designed and built especially for analyzing mine ventilation distribution problems, has been in service at the Bureau of Mines, Pittsburgh, for approximately nine years. It is a nonlinear, low voltage, fluid network analyzer of the type developed by the late Malcolm S. McIlroy, Professor of Electrical Engineering, Cornell University, who cooperated with G. E. McElroy, of the Bureau of Mines, in adapting the instrument to mine ventilation systems. Several modifications have been made since the original installation, but a considerable part of this paper is drawn from G. E. McElroy's original description of the analogue.' The analysis of water or gas distribution systems2,3 led to the development of this type of network analyzer. Several similar units are now employed by utility companies for this purpose; however, the Bureau of Mines unit is the only one designed specifically for mine airflow problems. Other airflow analogues employing a similar principle have been used at the Central Research Station of the Netherlands State Mines,4 in England,5 and in South Africa.6 Information on these devices indicated they were somewhat inflexible because commercial lamps with a sufficient range of resistance are difficult to obtain. Computers for mine ventilation analysis have been developed in Germany which instead of lamps use a variable resistance to adjust for turbulent flow laws. These units are expensive, but excellent results have been reported in their application. THE ANALYZER Theory of Application: As airflow generally follows the law of turbulent fluid flow, resistance to flow is nonlinear; consequently, the problem has been to find a nonlinear resistance element of suitable range that can be used for electrically simulated airflow. Tungsten filament lamps operated on alternating or direct current approximate the square-law resistance characteristics of mine airflow over a large range below maximum or rated voltage; that is, the voltage drop varies approximately as the square of the current. Consequently, the heart of the network analyzer is a nonlinear resistor known as a Fluistor, which is simply a custom-made low voltage, tungsten filament lamp that is available in a progressive series of relative resistance values ranging from 0.05 to 500 in nominal 5% steps. However, variations in manufacturing large groups result in differences of 1 to 3%, but series arrangements required for high-loss branches can be matched within about 1%. Utilizing a combination of Fluistors and load circuits, the mine ventilation system is duplicated electrically. Intake load circuits are connected from power intake to primary points of the circuit network; segments of unregulated flow along intakes and returns are represented by Fluistors, regulated splits and leakage paths are represented by load circuits and Fluistors of proper capacity; and mine exhausts are connected to ground from the last point of the network to complete the circuit. For the special purpose of representing booster fans or natural draft conditions, boosters or separate source circuits are provided that can be inserted between any two points of a network to increase voltage to the required value. Physical Layout: The analyzer consists of three 42

Jan 1, 1963

By Warren W. Jamar, Philip D. Pearson

IN recent years, there have been many changes in mining operations in the Lake Superior district. To follow these trends on the Menominee Range of Michigan, information has been assembled from all of the iron informatior]mining operations in the area. There are 15 operating underground mines in the Iron River-Crystal Falls area of the Menominee Range. Within the past two years, two idle properties were reopened pasttwoyears,and are now producing, and a third and fourth are being reopened. Also, there are two siliceous openpits operated by independent companies outside this immediate area. Six companies operate the underground mines, employing some 1850 employees. Table I shows pertinent facts about these properties. During 1949, the largest mine in Iron County shipped 571,287 tons, and one of the newer mines shipped 39,378 tons, with a range total to 3,535,-373 tons. Since the Menominee Range was opened in the 1870's, the mines in Iron County have shipped 85,890,922 tons. From an operator's viewpoint rather than a geologist's, the ore is' classified as semi-hard, composed of hematite and limonite. It is not as soft as the ores of the Marquette Range nor is it as hard as the hard ores of the Marquette and Vermillion Ranges. The ore bodies have slate hanging walls and slate footwalls. In most cases the hanging walls and footwalls are soft and high in sulphur. The sulphur comes from pyrite, and these slates will ignite when piled more than 6 to 8 ft high. Ore is mined on this range by: 1—Sub-level stop-ing; 2—Shrinkage stoping; 3—Sub-level caving; 4— Block caving; 5—Top slicing. The predominance of these methods is in the order named. Underground drilling is important in the mining cycle. Some changes made and trends toward future changes fall into four categories: Drill bits, drill steel, drill machines, and compressed air pressures. Several types of bits have been tried and are in use. They include the detachable tungsten carbide, insert bit; the intraset steel bit, which is a conventional steel rod with tungsten carbide insert; the one-use bit; and the multiple-use bit. For many years, detachable multiple-use bits have been standard. In tests conducted recently to im- prove drilling efficiency, this bit was used as the basis for comparison. Under existing conditions, a multiple-use bit can be resharpened about three times before it is discarded. A thorough test of 2-in. tungsten carbide threaded bits was conducted under various ground conditions. 1—The drilled footage ranged from 48 to 600 ft per bit; averaging 357 ft per bit. Under these same conditions, a multiple-use bit ranged from 8 to 80 ft per bit. 2—The bit cost was greater for the insert bit in each case. 3—The average drilling speed for the insert bit was 12 in. in 62 sec and for the multiple-use bit was 12 in. in 64 sec. In a second test, 2V4 -in. tungsten carbide bits with the large 1 3/16-in. thread were used in moderately soft ground on a 152-lb drifting drill on a long feed jumbo. 1—The drilled footage ranged from 450 to 5000 ft per bit, averaging 1810 ft per bit. Under these same conditions, a multiple-use bit averaged 64 ft per bit. 2—The bit cost was reduced by the use of the insert bit. 3—No increase in drilling speed was recorded. 4—Minimum footage obtained was caused by thread failure. To improve the thread life, thread size on the rod was increased and the bit was attached to the rod with a pipe wrench. After this, bits failed in equal proportions because of cracked skirts, broken inserts, and gage loss. Tungsten carbide bits are now used in the operation where this second test was conducted because labor costs were lowered as a result of reducing the number of bits being changed by the miners. At the same operation and under the same conditions as the second test, 1Y4-in. insert bits with standard 1-in. threads were used. These bits did not drill much more than 300 ft before thread failure and were then welded to the rods and used until total failure. Sometimes this footage was considerable, sometimes it was not. Chisel-type 2%-in. insert bits were

Jan 1, 1952

By W. G. Lidman, H. J. Hamjian

' I HE fabrication of many materials from powders involves a sintering process. A mass of powder will sinter because of the excess free energy over the same mass in the densified state caused by the higher total surface area of the powder. An understanding of the kinetics and mechanism of sintering should assist in improving the properties of such materials. The present investigation conducted at the NACA Lewis laboratory deals with the sintering of chromium carbide. Dry sintering (sintering at a temperature below the melting point) was divided into two stages by Shaler:' the first stage, during which the particles preserve much of their original shape and the voids are interconnected, and the second stage, during which densification occurs and the pores are isolated. The mechanism of forming interfaces between particles, or welding together of particles, has been investigated by Kuczynski2 and may be described by any one or a combination of the following mechanisms: viscous flow, evaporation and condensation, volume diffusion, or surface diffusion. The mechanism by which pores are closed or eliminated (densification) during sintering, is of interest. Grain growth observed during sintering may be attributed to the variation in the surface energies of individual grains, causing some grains to grow at the expense of others. Grain boundary migration occurs presumably by a diffusion process, therefore the rate of grain growth would be expected to increase exponentially with increasing time and temperature. Thus, for practical sintering times of less than 1 hr, a certain minimum temperature may exist at which major structural and property changes will occur. Densification and kinetics of grain growth during sintering under pressure of chromium carbide were investigated to provide additional information which will aid in describing more accurately the sintering process and the mechanisms involved. This material was selected for this study because of the current interest in high strength, oxidation resistant refractory materials, such as carbides, which are sintered to produce solid, dense materials from powders. Sintering under pressure is a process where the heat and pressure are applied to the compact simultaneously, specimens for this work were prepared by sintering under pressure at different temperatures and for various time periods. Experimental Procedure Preparation of Specimens: Chemical analysis of the commercial chromium carbide used in this investigation was as follows: Cr 86.19 pct, C 12.14 pct, and Fe 0.2 pct. X-ray diffraction powder patterns gave characteristic diffraction lines of Cr3C2 crystal structure. Powder particle size was determined microscopically and the average initial particle size was 6 microns with 85 pct between 2 and 10 microns. Specimens sintered under pressure were formed in graphite dies3 heated by induction. Sintering temperatures were measured with an optical pyrometer by sighting into a 3/8-in. hole drilled 1 in. deep at the midsection of the graphite die. A load of approximately 1 ton per sq in. was applied to the powder. The die assembly was heated in 20 min to the highest temperature (2500°F') at which no increase in grain size could be observed, and less than 2.5 min were required to heat from this temperature to the maximum temperature (3000°F). Sintering temperatures and times for the specimens of this investigation are indicated in Table I. Analysis of Specimens: Specimens polished with diamond abrasives were etched to reveal the grain boundaries with a 1:1 mixture of 20 pct potassium hydroxide and 20 pct potassium ferricyanide heated to 160°F. Representative areas of each sample were photographed at 1000 diameters. The largest diameters of all well-defined grains were measured, but only the measurements of 15 of the largest grains were averaged in order to determine an index of grain size on the assumption that they were among the first to begin growth. Densities were determined from differential weighing of the samples in air and water. The reported density values are considered correct within ±0.01 g per milliliter. Results and Discussion Metal compacts have exhibited grain growth when sintered at temperatures about two-thirds of the absolute temperature of their melting point.' Grain growth also occurs during the sintering of chromium carbide and is illustrated by the micrographs shown in Fig. 1. These micrographs were prepared from specimens sintered for 90 min at temperatures ranging from 1371°C (2500°F) to 1648°C (3000°F). Average grain size and density measurements of specimens investigated are presented in Table I. The relationship between grain size and sintering tem-

Jan 1, 1954

By J. D. Price, G. V. Woody

MANY producers of by-product coke have spent considerable time and given considerable thought to the use of a substitute for low-volatile coal as an admixture with high-volatile coking coal for charging high-temperature by-product coke ovens. Generally speaking, there are two principal reasons for this, one being occasioned by the lack, in certain localities, of low-volatile coal delivered at a cost sufficiently low to permit its use, and another being the desirability of finding a substitute product with a delivered cost lower than the cost of the low-volatile coal. It was for the first of these reasons that the Colorado Fuel and Iron Corporation was particularly interested in this subject. This paper deals with the experimental work and the results obtained by that company in producing and using low-temperature char. OVERCOMING DEFECTS IN COKE FROM COLORADO COALS Colorado coking coals are all of the high-volatile type and were all laid down during he Cretaceous period, which makes them something like 80 million or so years younger, geologically, than eastern coals. They behave, when coked alone, something like the Pittsburgh-seam coals. They make a very brittle, highly cross-fractured, fingery coke of fair shatter value but quite low in resistance to abrasion. Mixed with a low-volatile coking coal such as Pocahontas or certain Oklahoma coals, they make a coke that is equal to virtually any of the eastern cokes. The inferior quality of Colorado coke has long been recognized, and a considerable amount of time, effort and money has been spent upon investigations of means for improving its physical properties. While it has been found that some improvement in quality of coke can be secured by the adjustment and control of such variables as pulverization, moisture content, oven temperature and bulk density of charge, and through the addition of certain inert materials such as pulverized coke breeze, pitch, high-volatile noncoking or semicoking coal, the benefits so gained have been comparatively small. There are no low or medium-volatile coking coals available within a reasonable freight-cost distance. Therefore a substitute for low-volatile coal was developed through low-temperature carbonization-or, more correctly speaking, through the partial devolatilization-of high-volatile coals. The blending of this lower volatile char with the high-volatile coking coal distinctly improves the quality of coke made from Colorado coal. There is nothing especially new about this idea. In 1908 a Japanese patent was taken out by Kotaro Shimomura, whose claim was: The method of making a nonfingery coke out of bituminous coals, without using natural coals of low-volatile matter; according to which a certain coal is heated at a temperature between about 300°C. and about 600°C. so as to leave about 15 to 25 per cent volatile matter in the coal; this coal is mixed with the

Jan 1, 1944

By Aureal T. Cross

THIS paper is a review of the published literature on research in coal geology, principally exclusive of resource studies, which appeared or became available during 1950 and the latter part of 1949. This report is not to be construed as being complete. The papers referred to in the bibliography are those among many more, which were read either in full or in abstract. Undoubtedly other papers were published which either escaped the author's notice or were not available to him. Those which were seen in abstract only (about one fourth of those listed) were not available in time for the inclusion of more than a notice. An outline of all papers listed in the bibliography has been arranged by subjects and reasonable subdivisions with some papers cited under more than one subject. Most papers are indexed according to the principal subject of discussion or research only as to an unusual or noteworthy section of the entire report. There will likely be some disagreement as to the quality or merit of some of the papers selected and the specialist may be supercritical of the outline or organization of papers in his field. It may be that attention has occasionally been drawn to papers reporting old information or conclusions of questionable value. Conferences and Meetings One of the best indications of the growing interest in coal geology problems in the United States is the increasing number of times this field has been the focus of attention at conferences and meetings. Notable among these are the joint meeting of the Society of Economic Geologists and the Geological Society of America at El Paso, November 1949, at which the principal thesis was concerned with low rank carbonaceous fuel deposits, especially of western United States. Among the papers given which are already available were those presented by Barghoorn,'" Parry? Roe? and Parks."' At the annual meeting of the Botanical Society of America in New York, December 1949, a joint meeting of the Paleobotanical and Microbiological Sections was held for which a symposium on Microbiology in Relation to the Geologic Accumulation of Organic Complexes was organized. Publication of the six papers presented by Ralph G. H. Siu, Elso S. Barghoorn, Irving Breger, Claude E. ZoBell, James M. Schopf, and A. C. Thayson is anticipated. At the regular meetings of the Paleobotanical Section at the same time, several other papers of interest reported on coal ball studies, partial coalification of petrified wood, and floras. In Chicago, April 1950, a symposium on Applied Paleobotany was held by the Society of Economic Paleontologists and Mineralogists in conjunction with the American Association of Petroleum Geologists. The five papers presented at this meeting dealt with the use of Paleozoic plant microfossils for stratigraphic work, J. M. Schopf, Devonian-Missis-sippian fossils of the black shales, Aureal T. Cross, Mesozoic plants of stratigraphic value, Th. Just, plant microfossils of the Tertiary, L. R. Wilson, and studies of the Brandon lignite, Elso S. Barghoorn. Early publication of these in the Journal of Paleontology is expected. The Nova Scotia Research Foundation and the Nova Scotia Dept. of Mines sponsored an excellent 3-day conference in June 1950, which dealt with several aspects of coal geology. Papers on coal classification, P. A. Hacquenbard, structure and sedimentation problems in Nova Scotia, T. B. Haites, new techniques of thermal analysis, W. L. White-head, geochemical investigations of Nova Scotia coals, Irving Breger, the role of fossil plant spores in coal correlation and the stratigraphy of the coal-bearing strata of the Appalachian Region, Aureal T. Cross, were given. Some discussions of these papers by those in attendance were recorded, and the entire proceedings is being prepared for publication. In September 1950, an unusual 3-day field conference was held by the Ohio and West Virginia Geological Surveys under the sponsorship of the Coal Geology Committee. This study of the stratigraphy, sedimentation, and nomenclature of the Upper Pennsylvanian and Permian coal-bearing strata of southeastern Ohio, southwestern Pennsylvania, and northern West Virginia was augmented by two discussions on associated rocks (clays and shales) and stratigraphic nomenclature at Wheeling and Morgantown, West Va. An extensive guidebook was prepared, and transcriptions of the Morgantown meeting were made. As a follow-up of the September field conference, a round-table discussion was held on this general topic at a special open meeting of the Coal Research Committee in conjunction with the November meeting of the Geological Society in Washington. Short prepared statements to invite discussion were given on each of several topics by L. M. Cline, Carl 0.

Jan 1, 1953

By W. J. Leas, L. A. Rapoport

As a preliminary, consicleration is given to the conventional definition of relative permeability and to the conditions governing the simultaneous flow of oil and gas through porous media. For the conditions of flow prevailing throughout most of a gas drive reservoir, the oil and gas can reasonably be supposed to be in capillary equilibrium with each other. Under these conditions, and these conditions only, the relative permeability to liquid can be expressed as a function of saturation. The relative permeability to liquid in that case is dependent upon the distribution of fluids which itself is shown to be related to the capillary pressure, and, in turn, to the saturation. As a consequence, relative permeability to liquid can be expressed in terms of the volume and surface area of a network of liquid channels bounded by the rock and the gas phase. While the volume of this network can be evaluated accurately, the surface area cannot. However, for any such volume, maximum and minimum values of the corresponding surface area can be calculated from capillary pressure data. It is then possible to establish for any saturation the limits within which the value of the relative permeability to liquid must lie. As a consequence of the theoretical development, the validity of an experimental method for measuring relative permeability to liquid which utilizes a stationary gas phase is demonstrated. In this method capillary barriers are cemented to the ends of the core sample to permit the maintenance of capillary equilibrium between the two phases. At the same time, this procedure eliminates undesirable secondary phenomena such as end effects, fissure effects, etc., the presence of which adversely affect the results of other laboratory methods. The results obtained by theoretical calculations, and experimentally, are discussed. In view of the overall precision that can presently be obtained in reservoir calculations, the agreement between the calculated and measured relative permeability to liquid data can be considered satisfactory. In conclusion, for reasons of economy and simplicity, the procedure of calculating limiting relative permeability to liquid curves from capillary pressure data is indicated for general engineering purposes. It is shown that the above procedure can easily be extended to the cases where connate water is present. Its use for reservoir studies is particularly recommended in conjunction with the method for measuring relative permeability to gas' which simultaneously yields the capillary pressure data necessary for the calculations. THEORETICAL Definition of Relative Permeabilities — Basic Equations for Heterogeneous Flow The equations by which the relative permeability concept is defined and upon which the formulation of all of the gas-oil flow problems rests at the present time are expressed as: V, = — Grad PL = — Grad PL ....(la) PL µL, k k KG VG = —K - Grad PG k- Grad Pc .... (lb) Mo where ,. and G refer to liquid and gas; V is the volumetric rate of flow per unit gross area. µ the viscosity, Grad P the potential gradient. and k the specific permeability of the porous medium.* (For horizontal flow, Grad P becomes the pressure gradient; otherwise, gravity must be included.) According to these expressions, each of the constituent phases is considered similar to a homogeneous system where the volumetric rate of flow is proportional to the pressure gradient, and for each of which the constants of proportionality, k, and kG, are termed effective permeabilities. by analogy to the specific permeability as defined by Darcy's law in its original form. In order to obtain a convenient basis of comparison, the effective permeabilities are referred to the specific permeability, k, of the considered porous medium, with the help of the relations: k1. = KI. k kQ = Kn k..........(2) K,. and Kr. are defined as the relative permeabilities to the liquid and to the gas phase respectively, and frequently expressed in per cent of specific permeability. It may be seen that Equations (1); which appear to be a direct generalization of Darcy's law, correspond to the assignment at any given time of a set of "local" permeabilities to each point of the porous medium, and represent in a differential form the two fluid flow system as a simple superposition of the individual single flow systems. The above interpretation implies that the effective or relative permeabilities are independent of pressure or rate of flow,

Jan 1, 1951

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By J. B. Newkirk, W. G. Martin

Oxidized cobalt powder is known to have a magnetic hysteresis loop which is asymmetric with respect to the magnetization axis. The experiment described herein shows that the orientation relationship between the basal plane of hexagonal cobalt and the oxide which forms upon it at 400°C is {111}ux//{100.1 }co <110>ux//<11.O> This orientation relationship allows a magnetic interaction between the antiferromagnetic oxide and the ferromagnetic substrate which could account for the offset hysteresis loop. Oxidized nickel powder has a symmetrical hysteresis loop and so is apparently not influenced by any magnetic interaction between metal and oxide. The orientation of the oxide (cubic) was found to be identical with that of the nickel substrate when the oxide forms on a polished surface parallel with {111} Ni. IN 1956 W. H. Meiklejohn and C. P. Bean discovered that fine particles of cobalt which had been prepared in a certain way have a magnetic hysteresis loop which in a strong field is asymmetric relative to the magnetization axis.&apos; The cobalt powder exhibited this unusual magnetic property only after it had been oxidized in air or oxygen; it lost the shifted hysteresis loop when the oxide was reduced in hydrogen. X-ray and neutron diffraction&apos; experiments showed the presence of cobaltous oxide (COO) and hexagonal cobalt in samples exhibiting the biased hysteresis loop and specifically showed no indication of any other compound with the exception of mercuric oxide.* The magnetic hysteresis loop becomes symmetrical above the Nee1 temperature (paramagnetic state) of COO. Therefore, it was concluded that the anomalous magnetic behavior is associated with the influence of cobaltous oxide upon the metallic cobalt. It has been proposed that crystallographic coherency may exist between the cobalt and a plane of some antiferromagnetic material which has unbalanced spin distribution and sufficient magnetic anisotropy to hold its spin in the direction which existed when the specimen was cooled. COO has these properties. Therefore, Roth&apos; has suggested that cobaltous oxide may form with a {111} plane parallel and coherent with the basal plane of the hexagonal cobalt metal and proposed that, as consequence of the antiferromagnetic interaction between COO and the underlying cobalt, the magnetization direction in oxidized fine Co particles may be ro-lated from the easy c-direction.2 Such a relationship, he proposed, would explain the observed magnetic effect in oxidized cobalt powder. The main purpose of this study was to determine the orientation relationship, if any, between the basal plane of cobalt and the oxide which forms upon it. Attempts to produce a film of COO which was strong enough to be handled were not successful. However we did succeed in making a film of CoCo2O, which was strong enough to be removed from the cobalt substrate and mounted on an electron-microscope grid. Because of the close structural similarity of COO and CoCo2O, we believe the orientation relation found for CoCo2O, on cobalt probably also holds for COO on cobalt. The epitaxial relationship of NiO and Ni also was investigated. To date no shifted hysteresis loop has been observed with nickel powder. However, the similarity of atomic array in cobalt and nickel leads to the prediction that a {111} of NiO may be parallel with a (111} of nickel. Experimental Method Cobalt-Cobalt Oxide—A coarse-grained specimen of a (hexagonal) cobalt was prepared by allowing a large crystal of fee cobalt to transform slowly at 400°C. The crystal was then cut to expose a surface parallel with the basal plane. A back reflection Laue X-ray photograph showed that the 00.1 plane was within 2" of the cut surface. The surface was mechanically polished and then electro polished in 85 pet orthophosphoric acid after which the crystal was held for 30 min at 400°C in air. During this heat-treatment the surface darkened slightly due to the oxide film which formed. The film was not thick enough to give an X-ray diffraction pattern, even by a glancing-angle technique. Glancing-angle electron diffraction was not possible either, owing to the interference of the electron beam with the high magnetic fields which exist at the 00.1 surface of the cobalt. However, it was possible to make an electron-diffraction photograph of the oxide film by stripping it from the cobalt substrate using the method described later. By maintaining reference marks carefully, it was possible to preserve the orientation relationship between the stripped film and the substrate on which it was formed. The oxidized surface of the crystal was first covered with a 1 pet solution of collodion (cellulose nitrate) in amyl acetate. When the film was dry, small rectangles were scored on the surface with a needle point. The specimen was then repolished

Jan 1, 1959

By Emerson F. Heald

A systematic study was made of new and old data on chemical activities in Fe-Pt alloys at elevated ternperatuves. Experimental results may be expressed in terms of the excess free energy using Least-squares analysis of the data gave the following values for the constants: for the temperature range 1130° to 1350°C, and tentatively to 155O°C, B = -3.326564 and C = 0.221051; for the temperature range 650" to 850°C, B = -2.555690, C = 1.762735, and D =0.097196. In the study of iron-containing silicate systems, it is sometimes desirable to have a direct experimental measure of the activity of iron in the system. The well-known solubility of iron in platinum, often a headache in experimental work on iron compounds under reducing conditions, can be used to advantage in this respect. If the activity of iron in an Fe-Pt alloy in equilibrium with the silicate is known as a function of the composition of the alloy, chemical analysis of the alloy will give a knowledge of the activity of iron in all of the phases in the system. The present study was undertaken in order to elucidate the characteristics of Fe-Pt alloys as iron activity indicators. This work is intended to tie together some previous work, which may be summarized as follows. Larson and Chipman1 determined the activity of iron in Fe-Pt alloys at 1550°C by equilibrating platinum metal with calcium oxide-iron oxide-silica melts of known iron activity. Compositions of the resulting alloys were determined by chemical analysis. A similar study was carried out by Taylor and ~uan,&apos; who worked at 1300°C. They brought the Fe-Pt alloys into equilibrium with iron oxide under conditions of known partial pressure of oxygen, and thus, from the work of Darken and ~urr~,~ conditions of known iron activity. Compositions were determined indirectly, by following the change in weight of the sample. Sundaresen et el* used the electromotive force of cells in which the alloy formed one electrode in order to measure the activity of iron in the alloy at 650" and 850°C. These temperatures were chosen to be above and below the first-order phase transition which takes place upon the ordering of Fe3Pt and the second-order transition which occurs upon the ordering of FePt3. EXPERIMENTAL 1) High Temperatures. The starting materials used were thin platinum foil, about 0.002 mm thick, and Fisher Certified reagent ferric oxide, Fez03, which had been heated for 24 hr at 1000°C. An intimate mixture of 80-mesh Fez03 and platinum platelets was placed in a thin platinum foil envelope. The latter was suspended from thin platinum wires in the hot zone of a vertical-tube, platinum-wound furnace of the type described by Muan and ~sborn.~ A capillary gas mixer similar to that used by Darken and Gurry3 was used to prepare a precisely known mixture of carbon dioxide and hydrogen, which was allowed to flow upward through the furnace tube. The partial pressure of oxygen in contact with the sample was thereby fixed at a value which was calculated from the charts prepared by porter? Temperatures were measured with a Pt-10 pct Rh-in-platinum thermocouple, which was calibrated using the melting points of gold (1062 .@C) and diopside, CaMgSizOB (1391.5"C). Temperature control was maintained to within i3"C with a Geophysical Laboratory proportional controller, using the furnace resistance as the sensing element. Samples were quenched by passing a small current through the platinum suspension wires, allowing the sample to drop into a bath of dibutyl phthalate at the bottom of the furnace tube. Prior to chemical analysis the samples were washed with acetone and dried. 2) Chemical Analysis. It proved possible, in almost all cases, to separate the Pt-Fe platelets physically from particles of iron oxide. The platelets were dissolved in a small volume of aqua regia, evaporated to dryness, and redissolved to 0.1 M HC1. In order to determine iron in the platinum alloy potentiometrically, it is necessary first to remove the platinum. A 10-cm column of Amberlite IR-120 cation exchange resin in the hydrogen form provided separation quickly and quantitatively: The mixture of iron and platinum in 0.1 M HC1 was added to the top of the column, and washed with about 100 ml of 0.1 M HC1. Under these conditions, the iron, principally in the form of cations such as FeC1" and FeCl;, is held quantitatively in the uppermost centimeter of the column. The platinum, in the form of anions such as PtC&- , is washed through without being adsorbed. After a qualitative test with stannous chloride indicated all of the platinum was removed, the iron was

Jan 1, 1968

By R. M. N. Pelloux

This study of aluminum-copper alloys had two aims: 1) To determine the effect of the amount and distribution of a second phase, CuAl2, on the creep-rupture strength, ductility, and fracture characteristics of the alloys. The work of Gemmell and grant&apos; on single-phase aluminum-copper alloys provided a basis for the present study. 2) To attempt to provide a relation between the microstructure and strength of two-phase alloys deformed at comparatively high temperatures (greater than 0.45 T,). The creep behavior of the two-phase magnesium alloys Mg;Ce and Mg-A1, has been treated by Roberts. 9 In work reported by Sully and Hardy4 and by Underwood, Marsh, and Manning~ on A1-Cu two-phase alloys, and also in a portion of the work of Gemmell and grant,&apos; aging took place during the creep tests. &apos;In the present work, overaged aluminum-copper alloys were subjected to creep deformation at two temperatures, 500 and 700OF. The overaging treatments to produce the precipitate dispersions were performed prior to the test, at temperatures which were the same as the ultimate respective test temperatures. Hence, the present study is concerned with the creep behavior of stable (in contrast to "underaged") two-phase A1-Cu alloys, i.e., alloys in which no depletion of the solid-solution matrix occurred during: creep. MATERIALS AND PROCEDURE A 2 and a 3 pct Cu alloy, supplied by Alcoa, were studied; Table I shows the compositions. The grain size of the specimens, 0.9 to 1.0 mm, was achieved by annealing the machined test bars for 2 hr at 1000°F, followed by furnace-cooling to 900°F, and homogenizing for 4 hr. Both compositions are in the single-phase condition at 900°F. Two types of overaged dispersions, termed "over-aged I" and "overaged 11" were prepared after the grain size and the homogenization anneals. The overaged I alloys were prepared by quenching to room temperature after homogenization, aging at either 500&apos; or 700°F for 72 hr, and air-cooling to room temperature. The overaged I1 alloys were prepared by furnace-cooling after homogenization, to either 500&apos; or 700°F, and holding at the necessary temperature for 72 hr. The times for furnace-cooling from 900°F to 700° and 500°F were, respectively, about 2 and 3.5 hr. The structures of the A1-2 pct Cu alloys are shown in Fig. 1 (overaged I) and in Fig. 2 (overaged 11). A comparison of Figs. 1 and 2 shows that in the overaged I alloys the precipitate particles along the grain boundaries are much more closely spaced and the grain boundaries are straighter than they are in the overaged I1 alloys. Further, the particle size and spacing in the grains are smaller than in the overaged I1 alloys. Microstructures of the A1-3 pct Cu alloys have not been shown; the structures of these alloys differed from those of the A1-2 pct Cu alloys only in that the distribution density of the particles was greater. The creep specimens had a gage section 1 in. long with a diameter of 0.155 to 0.195 in. Before testing, specimens were electropolished in a solution of 100 cc glacial acetic acid, and 30 cc of 60 pct perchloric acid, at 40&apos; to 50°F, at 24 v. The specimens were kept refrigerated prior to testing to avoid precipitation at room temperature. Constant stress creep tests were conducted; the load was applied 75 min after heating of the specimens to test temperature had begun. EXPERIMENTAL RESULTS Creep-Rupture—Fig. 3 shows the log stress-log rupture life plotfor the overaged alloys studied, and also for the underaged A1-2 pct Cu alloy.&apos; At 700°F, the stress-rupture life relationship is nearly independent of both the amount of precipitate (i.e., composition) and of the type of overaging treatment (i.e., precipitate dispersion). At 500°F it is apparent that the increase in rupture life that was ex-

Jan 1, 1960

By O. D. Gaither

This paper is an analytical study of the flow of fluids through small vertical conduits. Small conduits are defined as 11/4-in. nominal diameter tubing size and smaller, and approximately twice this area for annular conduits (i.e., 1- X 21/2-in. annulus and smaller). Experimental data are presented for the 1-X2-in. and 11/4- X 2%-in. annuli, and the I-in. and 11/4-in. tubing, since these represent the small conduit sizes and configurations generally encountered in oilfield applications. Data have been gathered for these conduits for single-phase water, single-phase gas and two-phase water-gas mixtures, with particular emphasis on high gas-liquid ratios. Water rates in excess of 2,000 BID and gas rates in excess of 2.5 MMcf/D, and two-phase flow ratios in between these two, represent the scope of the data gathered. Existing equations have been applied to predict flowing pressures and compared with experimental data. New correlations have been developed. INTRODUCTION The increased economic pressure on the domestic oil industry in the United States has constantly required the use of new techniques and equipment designed to reduce the cost of finding and producing oil and gas. Since tangible items are most readily apparent in economic analysis, the advent of lower-cost well completions was inevitable. One of the methods used to reduce costs which has received widespread attention is the slim-hole completion technique where tubing is used as the well casing and in which small conduits are used for tubing if necessary. Small conduits, defined by Kirkpatrick1 as "11/4-in. diameter nominal tubing and smaller for tubing flow and less than twice the 11/4-in. diameter nominal tubing internal flow area for annulus flow", have also found widespread usage as siphon strings for de-watering gas wells and as "kill" strings in deep high-pressure oil and gas wells. The growing use of small-diameter tubing has resulted in an increased need for development of improved methods to measure or predict flowing bottom-hole pressures since the physical dimensions generally preclude the use of subsurface-recording pressure gauges. Even in the cases where small bombs are available, the relatively high velocities encountered at nominal flow rates make it necessary to use excessive weight bars or special hold-down devices. Attempts to use recognized correlations to accurately predict flowing or gas-lift performance in wells equipped with small conduits have been generally unsuccessful. Insufficient field data were available to allow the development of a correlation on this basis, and an experimental approach was applied in an attempt to obtain a workable relation. The experimental approach used to obtain the data presented in this paper was actually a compromise between a field installation and a laboratory study. A test well 1,000 ft in length was used to obtain flow data on single-phase liquid, single-phase gas and two-phase water-gas flowing mixtures. Liquid rates up to 2,200 B/D and gas rates up to 3 MMcf/D were used in the single-phase flow studies. Two-phase flow rates from 100 to 600 B/D with gas-liquid ratios from 500 to 8,000 cu ft/bbl were recorded. Experimental data were obtained for single- and two-phase flow through 1-in and 11/4-in. nominal tubing, and through the annuli between 1- and 2-in. and 11/4- and 2%-in. nominal tubing strings. Experimental results for the two-phase flow are compared to the Poettmann-Carpenter correlation2 which is widely used as a comparative standard for development of multiphase flow predictions in flowing and gas-lift wells. Correlations developed by Tek,3 Baxendell and Thomas" were also investigated. The experimental data recorded herein fell in between the two flow regimes as defined by Ros," and this correlation also failed to yield satisfactory results. The fact that existing correlations failed to confirm the experimental data led to the need for development of a new correlation. Although a two-phase flow study was the primary objective of this investigation, data were also recorded for single-phase flow of water and gas, and constants were developed relating to pipe roughness and equivalent diameters for annular flow. These single-phase studies assisted materially in the development of certain of the two-phase flow results. Considerable previous work has been published which presented relationship of surface measurements to bottom-hole condition. The works of Buthod and Whiteley,6 Jones,&apos; Poettmannb and the Texas Railroad Commission" are classic examples of the successful use of mathematical relationships which allow acceptable predictions of subsurface pressures, when gas is the flowing fluid. Darcy and others have derived relationships which may be used with minor modifications to predict subsurface flowing conditions in injection and water-supply wells. As previously stated, the application of the single-phase flow relationships

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By L. A. Rapoport, W. J. Leas

As a preliminary, consicleration is given to the conventional definition of relative permeability and to the conditions governing the simultaneous flow of oil and gas through porous media. For the conditions of flow prevailing throughout most of a gas drive reservoir, the oil and gas can reasonably be supposed to be in capillary equilibrium with each other. Under these conditions, and these conditions only, the relative permeability to liquid can be expressed as a function of saturation. The relative permeability to liquid in that case is dependent upon the distribution of fluids which itself is shown to be related to the capillary pressure, and, in turn, to the saturation. As a consequence, relative permeability to liquid can be expressed in terms of the volume and surface area of a network of liquid channels bounded by the rock and the gas phase. While the volume of this network can be evaluated accurately, the surface area cannot. However, for any such volume, maximum and minimum values of the corresponding surface area can be calculated from capillary pressure data. It is then possible to establish for any saturation the limits within which the value of the relative permeability to liquid must lie. As a consequence of the theoretical development, the validity of an experimental method for measuring relative permeability to liquid which utilizes a stationary gas phase is demonstrated. In this method capillary barriers are cemented to the ends of the core sample to permit the maintenance of capillary equilibrium between the two phases. At the same time, this procedure eliminates undesirable secondary phenomena such as end effects, fissure effects, etc., the presence of which adversely affect the results of other laboratory methods. The results obtained by theoretical calculations, and experimentally, are discussed. In view of the overall precision that can presently be obtained in reservoir calculations, the agreement between the calculated and measured relative permeability to liquid data can be considered satisfactory. In conclusion, for reasons of economy and simplicity, the procedure of calculating limiting relative permeability to liquid curves from capillary pressure data is indicated for general engineering purposes. It is shown that the above procedure can easily be extended to the cases where connate water is present. Its use for reservoir studies is particularly recommended in conjunction with the method for measuring relative permeability to gas&apos; which simultaneously yields the capillary pressure data necessary for the calculations. THEORETICAL Definition of Relative Permeabilities — Basic Equations for Heterogeneous Flow The equations by which the relative permeability concept is defined and upon which the formulation of all of the gas-oil flow problems rests at the present time are expressed as: V, = — Grad PL = — Grad PL ....(la) PL µL, k k KG VG = —K - Grad PG k- Grad Pc .... (lb) Mo where ,. and G refer to liquid and gas; V is the volumetric rate of flow per unit gross area. µ the viscosity, Grad P the potential gradient. and k the specific permeability of the porous medium.* (For horizontal flow, Grad P becomes the pressure gradient; otherwise, gravity must be included.) According to these expressions, each of the constituent phases is considered similar to a homogeneous system where the volumetric rate of flow is proportional to the pressure gradient, and for each of which the constants of proportionality, k, and kG, are termed effective permeabilities. by analogy to the specific permeability as defined by Darcy&apos;s law in its original form. In order to obtain a convenient basis of comparison, the effective permeabilities are referred to the specific permeability, k, of the considered porous medium, with the help of the relations: k1. = KI. k kQ = Kn k..........(2) K,. and Kr. are defined as the relative permeabilities to the liquid and to the gas phase respectively, and frequently expressed in per cent of specific permeability. It may be seen that Equations (1); which appear to be a direct generalization of Darcy&apos;s law, correspond to the assignment at any given time of a set of "local" permeabilities to each point of the porous medium, and represent in a differential form the two fluid flow system as a simple superposition of the individual single flow systems. The above interpretation implies that the effective or relative permeabilities are independent of pressure or rate of flow,

Jan 1, 1951

By N. J. Grant, B. C. Giessen, Paul Predecki

ALLOYS of gold, silver, and aluminum with elements of the groups BII, BIII, BIV, and BV were prepared by a rapid quenching technique (splat) and were examined by X-ray diffraction. Five new intermediate phases were found and will be described briefly herein. For the gold and silver systems, the concentration ranges having an electron/atom ratio e/a of 1.4 to 1.5 ("3/2 Hume-Rothery phases") were studied primarily. Master alloys were prepared from high-purity metals (99.9+ pct or better) by melting either in evacuated fused silica capsules or by nonconsum-able-electrode arc melting in an argon atmosphere. Small pieces, 20 to 50 mg, of each alloy were blast-atomized to form a splat, by a technique similar to that described by Duwez and Willens.1 The technique used for this study is described in detail in Ref. 2; it utilizes a resistance-heated graphite crucible with a small hole at the bottom, directed toward a metal substrate or quenching plate. The prepared alloy rests over the fine hole, through which it is expelled by an explosion shock wave in the form of fine droplets (1 to 50 µ) of molten metal onto a copper or silver substrate, which is maintained at about -190°C. The resulting very high cooling rates (see Ref. 2 for quantitative measurements) can prevent the process of nuclea-tion and growth in many instances, resulting in the formation of metastable phases. The splat particles were transferred to a GE-XRD5 diffractometer and maintained at -190°C, where they were examined with CuKa radiation. The samples were then allowed to warm to room temperature or were heated to higher temperatures until the equilibrium structures formed. Of fifteen alloy systems considered, nonequi-librium structures were encountered in six; these are described below and summarized in Table I. In the system Au-Sb a metastable £ phase (A3 type, hcp, a = 2.898 + 0.002A; c = 4.731 * 0.004A; c/a = 1.633) was found in the concentration range Au + 13 to 15 at. pct Sb. This phase is isomorphous with the stable phases in the systems Au-Cd, Au-In, and Au-Sn, all at an average e/a ratio of 1.4 to 1.5. The concentration range of one-phase metastable was deduced from the small amounts of supersaturated gold solid-solution phase present in the splat product. It was found that ? could also be retained by splatting onto a substrate held at room temperature: however, decomposed into the equilibrium phases Au + AuSb2 after heating to 200°C for 1/2 hr, or on holding the powdered splatted alloy at 20°C for several months. Calorimetric measurements will be made in an attempt to decide the question whether ? is metastable at all temperatures or whether it is a stable phase at low temperatures. There is evidence that another phase, possibly also close-packed but with a different stacking sequence, can be obtained by rapid quenching of alloys with a different antimony content. Klement, Willens, and Duwez3 reported the existence of an amorphous phase on quenching Au-Si alloys (25 at. pct Si) to - 196°C. They found that on heating to room temperature another phase of unknown crystal structure was formed. This was confirmed (see Table I); however, the new crystalline phase, designated as ?, could also be formed simply by rapid quenching to room temperature, and even was found to exist already in the as-cast Au + 20 at. pct Si alloy. It was found that ? decomposed into Au + Si on the specimen surface at room temperature. This behavior, and the question whether or not there is an equilibrium-temperature region for ?, have not yet been resolved. It is probable that ? (Au + 20 to 21 at. pct Si) is cubic of the -brass type (D81-3) with a = 9.60, + 0.01A and N = 52 atoms per cell [compare 6 (CU-Sn)4]. Except for two very weak lines, the powder pattern of about thirty lines could be indexed on this basis; however, a determination of the atom positions has not yet been attempted. For Au-Ge the C phase was observed at about 21 at. pct Ge as reported by Luo et at.5 Lattice parameters a = 2.876A, c = 4.73,A, c/a = 1.64 were found. In the Au-Pb system, formation of a ? phase was not observed, but in the lead-rich region at 75 at. pct Pb, broad peaks belonging to an amorphous phase were found. The maximum diffracted intensity occurred at 28 = 32.4 deg which is about 1 deg larger than the position of the (111) line of lead (Cuka). For Ag-Pb, an amorphous phase analogous to the one found in the Au-Pb system was observed; this metastable phase exists probably at about 75 at. pct Pb. Since no lead-rich alloys were tested, all alloys consisted of silver + amorphous phase at -190°C. In A1-Ge alloys, line-rich and complex powder patterns were obtained at about 30 at. pct Ge; they bear similarities to those of aluminum and germanium, but are of lower symmetry; the existence of more than one intermediate phase is possible. The authors are grateful to the Kennecott Copper Corp. for Fellowship support, and ARPA (Contract

Jan 1, 1965

By J. E. Lawver, J. D. Raulerson, Charles C. Cook

The agriculture industry has made great strides during the past decade to increase agriculture yields through increased use of fertilizers. Increased use of fertilizers may prevent, or at least delay, mass starvation due to the alarming increase in world population. Phosphate was added to soil as a plant nutrient in the form of calcined bones at least 2000 years ago (Anon., 1964), and man has used phosphate minerals as a source of fertilization in one form or another for at least 100 years. During 1977 the world produced about 116 Mt of phosphate rock, with about 86% used for fertilizers and another 4% for animal feed supplements. More than three-fourths of the total production comes from the United States, Morocco, and the Soviet Union. From a mineral beneficiation point of view, the major sources of phosphate rock and the methods of beneficiation can be classified as follows: marine deposits not containing appreciable carbonate minerals, marine deposits requiring a francolite carbonate mineral separation, igneous deposits not containing appreciable carbonate minerals, and igneous deposits requiring apatite carbonate mineral separation. [ ] Guano, mostly from Chile and Peru, accounts for 0.1% of the total world production, and the calcium phosphates from Ocean, Nauru, and Christmas Islands and the aluminum and iron phosphates from Brazil and Aruba account for less than 4% of the world production and are thus not considered in this classification (Lawver, et al.). At present, marine phosphorite deposits account for about 75% of the world&apos;s production; the igneous deposits account for 20%. The igneous deposits low in carbonate minerals are easily concentrated by crushing, grinding, and apatite flotation. The most important igneous deposits are those of the Kola Peninsula, USSR (Woodrooffe, 1972). The igneous deposits high in carbonate materials are of corn appreciably more difficult to beneficiate, but they have been concentrated by froth flotation for a number of years. An interesting but rather complicated flowsheet of this type is at Phalabonva, in the Republic of South Africa (Lovell, 1976). The Phalaborwa deposit is an igneous complex of pyroxenite with a central core of carbonatite surrounded by a serpentine- magnetite-apatite rock called phoscorite. The phoscorite containing about 10% P2O5, 35% magnetite, and 35% calcium magnesium carbonate is currently being processed. The process involves comminuting the material for fiberation and subjecting it to a copper float using a potassium amyl xanthate as collector and triethoxybutane as a frother followed by a magnetic separation of the tailings to produce a feed for phosphate flotation. This process produces a phosphate concentrate containing greater than 36% P2O5 at a P2O5 recovery ranging from 75 to 80%. Considerable success has been claimed for recovering apatite from carbonate-bearing ores at the Jacupiranga Mine of Serrana S/A (Silva and Andery, 1972). The carbonatite currently being mined contains an average of only 5% P205 and is concentrated using a unique flotation process (Andery, 1968) to yield 96% P205 concentrates. The ore contains about 12% apatite, 5% magnetite, 80% calcite plus dolomite, and minor amounts of phlogopite, olivine, zircon, ilmenite, and pyrochlore. Feed preparation consists of crushing to -31.75 mm (-1 M in.), rod milling in closed circuit with hydrocyclones to about 92% (-50 mesh), and two-stage cyclone desliming of the -50 mesh sands at 20 m. Weight recovery in the deslimed feed is normally 85 to 88% and the corresponding P2O5 recovery is usually about 90%. The deslimed feed is conditioned at 60 to 70% solids for 15 min at pH = 8-10 with 0.6 kg/t of causticized starch for iron oxide and calcite-dolomite depression. The conditioned slurry is diluted to 20 to 30% solids, about 0.2 kg/t of fatty acid or soap collector is added to the conditioner discharge, and the reagentized ore is subjected to rougher-scavenger flotation with additional fatty acid added to the scavenger float. The scavenger concentrate is returned to rougher circuit distributor, and the rougher concentrate froth is subjected to two stages of cleaner flotation to yield a final apatite concentrate analyzing 36 to 38% P205. Flotation recovery of P205 is, in general, above 90% when treating fresh carbonatite. The high-carbonate flotation tails normally analyze 1 % P2O5 or less and are suitable for portland cement production. The marine deposits. Types 1 and 2 of central Florida are representative of enormous reserves of phosphate rock that will undoubtedly account for much of the world&apos;s production in the near future. Until very recently the sedimentary deposits high in carbonate minerals (Type 2) have not been considered reserves due to the difficulty in making a francolite-carbonate separation. Although no commercial plant has yet been built to beneficiate Type 2 ore, laboratory and pilot plant data indicate the process is viable. If so, the reserves of Florida and similar deposits throughout the world will be substantially increased. A discussion of the beneficiation of these two types of sedimentary deposits and the relation of the resulting concentrates to the fertilizer industry of the United States is the subject of this paper.

Jan 1, 1981

By J. D. Miller, J. E. Sepulveda

Tar sand deposits in the state of Utah contain more than 25 billion bbl of in-place bitumen. Although 30 times smaller than the well-known Athabasca tar sands, Utah tar sands do represent a significant domestic energy resource comparable to the national crude oil reserves (31.3 billion bbl). Based upon a detailed analysis of the physical and chemical properties of both the bitumen and the sand, a hot-water separation process for Utah tar sands is currently being developed in our laboratories at the University of Utah. This process involves intense agitation of the tar sand in a hot caustic solution and subsequent separation of the bitumen by a modified froth flotation technique. Experimental results with an Asphalt Ridge, Utah, tar sand sample indicated that percent solids and caustic concentration were the two most important variables controlling the performance of the digestion stage. These variables were identified by means of an experimental factorial design, in which coefficients of separation greater than 0.90 were realized. Although preliminary in nature, the experimental evidence&apos; gathered in this investigation seems to indicate that a hot-water separation process for Utah tar sands would allow for the efficient utilization of this important energy resource. The projected increase in the ever-widening gap between the domestic energy demand and the domestic energy supply for the next few years has motivated renewed interest in energy sources other than petroleum, such as tar sands, oil shale and coal. Although a number of research programs on the exploitation of national coal and oil shale resources have already been completed, very few programs have been initiated on the processing of tar sand resources in the United States. In recognition of their significance as a domestic energy resource, investigators at the University of Utah have designed an extensive research program on Utah tar sands. An important phase of this program, and the main subject of this publication, is the development of a hot-water process for the recovery of bitumen from Utah tar sands, as a preliminary step toward the production of synthetic fuels and petrochemicals. The term "tar sand" refers to a consolidated mixture of bitumen (tar) and sand. The sand in tar sand is mostly a-quartz as determined from X-ray diffraction patterns. Alternate names for "tar sands" are "oil sands" and "bituminous sands." The latter is technically correct and in that sense provides an adequate description. Tar sand deposits occur throughout the world, often in the same geographical areas as petroleum deposits. Significantly large tar sand deposits have been identified and mapped in Canada, Venezuela and, the United States. By far, the largest deposit is the Athabasca tar sands in the Province of Alberta, Canada. According to the Alberta Energy Resources Conservation Board (AERCB),2,3 proved reserves of crude in-place bitumen in the Athabasca region amount to almost 900 billion bbl. To date, this is the only tar sand deposit in the world being mined and processed for the recovery of petroleum products. Great Canadian Oil Sands, Ltd. (GCOS) produces 20 million bbl of synthetic crude oil per year. Another plant being constructed by Syncrude Canada, Ltd. is expected to produce in excess of 40 million bbl of synthetic crude oil per year. According to the Utah Geological and Mineral Survey (UGMS), tar sand deposits in the state of Utah contain more than 25 billion bbl of bitumen in place, which represent almost 95% of the total mapped resources in the United States.4 The extent of Utah tar sand reserves seems small compared to the enormous potential of Canadian tar sands. Nevertheless, Utah tar sand reserves do represent a significant energy resource comparable to the United States crude oil proved reserves of 31.3 billion bbl in 1976.5 Tar sands in Utah occur in 51 deposits along the eastern side of the state.4 However, only six out of these 51 deposits are worthy of any practical consideration (Fig. 1). As indicated in Table 1, Tar Sand Triangle is the largest deposit in the state and contains about half of the total mapped resources. Information regarding the grade or bitumen content of Utah deposits is still very limited. The bitumen content varies significantly from deposit to deposit, as well as within a given deposit. In any event, the information available6-8 seems to indicate that Utah deposits are not as rich in bitumen as the vast Canadian deposits which average 12 to 13% by weight.9 Although many occurrences of bitumen saturation up to 17% by weight have been detected in the northeastern part of the state (Asphalt Ridge and P. R. Spring), the average for reserves in Utah may well be less than 10% by weight. Separation Technology As in any other mining problem, there are two basic approaches to the recovery of bitumen from tar sands. In one

Jan 1, 1979

By J. G. Patchett, R. W. Rausch

Fresh waters and the presence of clay in many Rocky Mountain and West Coast sands require special methods of log analysis. Archie&apos;s saturation equation requires addition of a shale correction term, and the SP equation must also be modified to account for clays. Suitable equations were developed several years ago, but have not been widely used due to the algebraic complexity. A computer-oriented method has now been developed to overcome this problem. The basic shaly sand equations are rearranged in four different ways to permit solution for various sets of available input data. Essential to application of the method is the correction of observed SP values to those that would be observed if the resistivity of the formation waters were exactly interchangeable with the activity. A graphic method for doing this is given. Where conditions require consideration of the effect of clay in the sands, the method presented has been found to improve the accuracy of water-saturation determinations. INTRODUCTION Log interpretation in many Rocky Mountain and West Coast basins is complicated by rapid vertical and lateral changes in water resistivity. Calculation of formation water resistivity from the SP curve becomes difficult in zones that contain clay, since changes in SP deflection may be due to changes in either clay content or water salinity. In hydrocarbon-producing reservoirs, the problem is further complicated because hydrocarbon saturation also reduces the SP.1 A log interpretation system using computers has been developed to provide a solution to this problem, based on equations proposed by de Witte.2 Four different simultaneous solutions of de Witte&apos;s equations have been made. Each solution method uses a different set of input data as independent variables. Thus, a choice of solution method is possible, depending upon the logs run and the availability of other data. Two of the solutions do not require a knowledge of water resistivity. This system is intended to be used primarily in multiple sandstone-shale sequences of low and moderate resistivities where the principal contaminant in the sandstones is clay. However, where sufficient regional data are available, interpretation in single-zone sandstone reservoirs can also be improved by using the method. THEORY AND HISTORY OF SHALY SAND ANALYSIS The log interpretation formula originally proposed by Archie3 in 1941 is applicable only to rock-fluid systems wherein the rock has negligible electrical conductivity. In 1949, Patnode and Wyllie4 showed that if the rock itself can be considered conductive due to the presence of clay, a different calculation approach is necessary. During the following years, this problem was investigated at great length, as was the related problem of the effect of rock conductivity on the SP.5-11 These investigations established functional relationships between SP, resistivity, water saturation and water resistivity for such a formation. Refs. 2 and 12 provide summaries of these studies. Unfortunately, practical use of these relationships required that water resistivity be known independently from the SP. Although log interpretation methods for rock systems containing clay were proposed at that time,&apos; they were not generally accepted for routine use. There are three principal reasons for this. First, in many field situations involving high-salinity water, rock conductivity may be neglected (even if present) without introducing appreciable error. This may be seen by considering the following expression for waier-saturated rock.&apos; 1/R2=1/R1+1/FRn....(1) where 1/R, is conductivity due to clay. As Rw becomes small, I/FRw becomes much greater than 1/R, which may be neglected. Where 1/R, may be neglected, the sandstone is called clean. If the term may not be neglected, the sandstone is termed dirty or shaly. For resistivity purposes, the classification between clean and shaly sands then depends not only upon the conductivity due to shale in the sand, but also upon the resistivity of the associated water (shale is used here to mean surface condition due to disseminated clay). A sand of given conductivity might safely be treated as clean in association with high-salinity water, but would require shaly sand methods if associated with fresher waters. Shaly sand methods are not required in many areas having saline waters; but in Rocky Mountain and West Coast sands having relatively fresh waters (often more than 0.3 ohm-m resistivity at formation conditions), the shaly sand methods are needed. Errors Rw calculations from the SP due to the presence of shale are likewise related to water salinity. In saline water formations drilled with fresh mud, the ratio of mud filtrate resistivity to water resistivity is high, the SP is large and the presence of shale can introduce large errors in water resistivity calculated by the conventional method. When the resistivity ratio is low, the errors are smaller. At zero SP, no error would result from shale. Thus, from the SP viewpoint, a given rock could be shaly if associated with a saline water, and clean in association with a fresh water, which is the opposite of the resistivity-oriented definition above.

By William Rawles

THIS report is the first of a series planned by The Mineral Inquiry, organized by the American Institute of Mining and Metallurgical Engineers to make factual studies of the world&apos;s mineral resources in their political and international relations. The efforts of The Inquiry are directed toward filling certain notable gaps in available information, without duplication of many other investigations touching this field which are being conducted by other organizations, public and private. Much of its work is indirect-suggesting topics for study by other individuals and organizations and aiding in their planning and correlation. The present report is designed to show the nationality of the commercial control of mineral resources, in addition to the geographic distribution of production by nation-in other words, the overlapping of political and commercial control. Detailed accuracy cannot be claimed for the figures presented. For some companies division of ownership by nationality has been very difficult to ascertain, as many companies are reluctant to dis-close the extent of their holdings in other corporations, and the impossi-bility of tracing all stock ownership is obvious. Also there are frequent changes in the distribution of stock control. The percentages given rep¬resent only the participation that has been identified definitely, even though it is known in some cases with fair certainty that larger interests exist, as, for instance, American participation in foreign chromite companies. In bringing together these figures we have had the cordial co-operation of the United States Bureau of Mines, the United States Geological Survey, the Mineral Division of the Bureau of Foreign and Domestic Commerce, the Brookings Institution, and of many specialists connected with the mineral industries covered. The report does not cover the entire list of minerals. Coal, although it is, one of the most important minerals, was omitted because not more than one or two per cent of the world&apos;s production comes from mines that are owned by companies foreign to the countries where the coal is produced. Others, such as gold, are omitted because sufficient data have not yet been accumulated. It is hoped that some of these may be covered in subsequent reports. It is not our purpose here to discuss the significance of the figures sub-mitted, but only to point out the general nature of the problems on which they have a bearing.

Jan 1, 1933