Search Documents

By Robert B. Hall, Charles W. Bauer

Alunite, KA13(SO4)2(OH)6, was used from the 15th century until early in the present century as a source of potash alum and aluminum sulfate. Since early in this century its use for this purpose has declined to virtually nil with acid treatment of bauxite or clay replacing alunite as a source of alum. During the First World War alunite served as an emergency source of potassium sulfate fertilizer in the United States and Australia. More recently, alunite has been investigated as a nonbauxite ore of aluminum, with potassium sulfate and sulfuric acid as recoverable byproducts. The Soviet Union established a commercial-scale plant in the mid-1960s producing cell-grade alumina from alunitic ore and recovering byproduct potassium sulfate and sulfuric acid. Interest in nonbauxite ores, including alunite, is expected to continue in countries that presently are dependent on imports of bauxite or alumina to support their aluminum industries. The US has the largest aluminum producing and fabricating industry in the world but imports more than 90% of its aluminum-bearing raw materials. This chapter reviews the status of alunite as a potential source of aluminum with byproduct fertilizer, and the resources available to support this industry, whenever this becomes economically feasible. Historical Background Alunite has received little attention until recently, yet for centuries it held a higher place than many of the industrial minerals considered essential today. The deposits of alunite at Tolfa, northwest of Rome, were mined almost continuously from 1462 until about the middle of the present century (DeLaunay, 1907; Lombardi, 1977). The potash alum made from the alunite, called Roman alum, was a valued commodity in international trade, being used in the making of paper, in tanning leather, and as a mordant in textile dyeing. The industry at Tolfa is said to have been controlled by the Vatican and was a source of revenue for the Papal State until fairly recent times (Lombardi et al., 1977). Since early in the present century, the preferred method for making alum has been by acid treatment of bauxite, clay, or other aluminous material, so that alunite has fallen into near-total disuse as a source of alum chemicals. Appreciable quantities of mixed alunite-kaolin rock still are quarried in the Tolfa district for use by the Italian cement industry (Lombardi and Mattias, 1979), but the role of this material in the making of Italian cement is not clear. A large deposit of alunite-bearing rock in Chekiang province, mainland China, is said to have yielded a "great quantity" of alum over a period of several centuries (Yih, 1931). During the 19th century, alunite deposits in France, Spain, Australia, and other countries were exploited as raw material for making potash alum and aluminum sulfate. Interest in alunite was revived during the First World War, not for alum, but as an emergency source of potassium sulfate fertilizer. Large veins of nearly pure alunite had just been discovered on Alunite Ridge, 11 km southwest of Marysvale in Piute County, Utah (Butler and Gale, 1912), and these were quickly brought into production for making potash fertilizer to replace supplies cut off from the famous deposits at Stassfurt, Germany. An estimated 225 kt of Marysvale alunite was mined and processed at that time (Callaghan, 1973). Similarly, Australian alunite deposits were mined as a source of potash fertilizer during both world wars (Hall, 1978).

Jan 1, 1983

By Frank M. Brooks

Endeavoring to meet the requirements of more oil created by the war effort, the oil industry of Kansas in 1944 maintained a level of exploratory and development work well above the average of previous years. During the year, 1866 test wells were drilled, of which 838 were completed as new oil wells, 89 as new gas wells and 939 as dry holes. These figures compare with 895 new oil wells, 32 new gas wells and 834 dry holes for a total of 1761 tests in 1943. Thus, there was a marked decrease in the number of new oil wells, and an increase in the new gas wells and dry holes. The total daily initial potential of production discovered by the successful wells amounted to 494,689 bbl. of oil as compared with 736,297 bbl. in 1943. The increase in unsuccessful tests and decrease in the daily initial potential of production has caused considerable comment and speculation among the oil operators of the state. This tends to confirm the fact that each year it is becoming harder to find the remaining oil fields. The total number of wildcats drilled in Kansas was 651 in 1944, including the so-called "extension wildcats" (more than 1/2 to I mile from known production). Of this number, 410 are classified as "rank wildcats," being more than one mile from production. The total of 651 compares with the total of 577 for 1943, showing an increase in wildcat drilling of 74 test wells: Of these tests, 81 were successful in finding new oil fields or substafitially extending known areas of production, and 20 either found new gas pools or extended old ones. The ratio of dry holes to producing wells among the "rank" wildcats was approximately 10 to I (372 dry and 38 successful). This compares with a ratio of 10 to I in 1943, 11 to I in 1942, and 6 to I in 1941. A compilation of the results obtained by 34 active companies, including all of the major companies, reveals that their efforts accounted for 286 rank wildcat wells, of which-37 found oil or gas, making a dry-hole percentage of 87. Approximately 55 per cent of all, tests (pool and wildcat) were failures. Thc Kansas Nomenclature Committee of the Kansas Geological Society named and classified 43 new oil fields and 8 new gas fields during the year. The discovery wells of the new oil fields had a total initial potential of 18,358 bbl. for a per well average or' 427 bbl. The average initial potential per well in 1942 and 1943 was 894 and 419 bbl., respectively. With nearly all wells producing at their economic capacity, the total production for Kansas during 1944 totaled 101,896,704 bbl., which compares with total runs of 107,709,679 bbl. in 1943. Production of natural gas during 1944 reached a new high record of 135 billion cubic feet. This compares with 122 billion cubic feet for 1943, and 98 1/2 billion cubic feet for 1942. The large Hugoton field of Southwestern Kansas accounted for about half of the total gas produced. Highlights OF the Year Although 51 new fields were found in the state during 1944, it is apparent that the

Jan 1, 1945

By J. W. Axelson, J. N. S. Kwong, J. T. Adams, E. L. Pire, J. F. Johnson

CRUSHING has always been a major operation in the chemical and metallurgical industries, yet little is known about the theory of crushing, and today, the design of crushers is still based almost entirely on empirical knowledge and accumulated practical experience. In view of the increasing national need for the economic working of poorer grades of ores, the lack of a fundamental understanding of this unit operation hardly presents a satisfactory situation. Basic investigations of crushing have been concerned mainly with three phases of the problem, 1—the mechanism of the fracture process itself, 2— the particle size distribution of the crushed product, and 3—the relationship between the energy input and the amount of new surface produced. Probably adequate information on these phases will be required for a comprehensive understanding of the process of crushing. A bibliography covering these phases will be given as well as a review of the work done in recent years at the University of Minnesota on the relationship between the energy input for crushing and the amount of new surface produced. Fracture of Solids This phase of the crushing problem is the most fundamental because it is concerned with the actual mechanism of fracture. The main problems in a study of the fracture process are concerned with the questions of why and how fracture occurs and why there is such a discrepancy between the actual and theoretical energies needed in a fracture process. In considering the how of fracture, Poncelet40 carried Out experiments On the crushing of glass under compression. Actual photographs were taken of the glass plates in various stages of fracture-From these experiments Poncelet was able to postulate a probable mechanism of fracture. Probably because of the quantitative nature of the problem, considerable effort has been expended to find an explanation for the low tensile strength of materials and the conversely high energy for crushing. Tensile strengths are often only 1/500 of the theoretical while the energy for crushing is usually at least 500 times the theoretical. A review of the work to explain the low strength of glass is given by Weyl,52 in which the work of Griffith," Joffe,32 and Powell and Preston " are discussed. In a recent paper Seitz46 holds that the presence of flaws and fissures, as originally expounded by Griffith, is the most probable reason for the low strength of metals. Since a mass of data has been accumulated on the fracture strength of materials, several investigators have attempted to correlate this data in the form of an equation. Glathart and Preston," Taylor,40 Ponce-let,40 Machlin and Norwich," and Frederickson and Eyring,20 have all considered this problem, and the final relationship has always been that the stress for fracture is proportional to the logarithm of time. Murgatroyd30 arrived at a very similar relationship by a different method. The study of the fracture process has used many methods. Among some of these are the thermody-namic approach used by Furth21 which was based on the work of BornxL and later used by Saibel,43 and the statistical method applied by Fisher and Hollo-mon.10 In recent years metallurgists have concerned themselves with the problem of brittle .fracture as opposed to the usual plastic failure of metals. The great amount of interest shown is reflected by the number of contributors to the recent symposium on the fracture of metals.' Repeated failures of welded ships during World War II have accelerated these Studies.31,50,34 The shifting of thought from an emphasis on the

Jan 1, 1952

By H. H. Rachford

This work presents a first-order analysis of the instability underlying viscous fingering in adverse viscosity-ratio water floods. It extends previous analyses of frontal instabilities, which were carried out with equations for parallel plate models, by including effects of the saturation transition zone observed behind the front in water floods in water wet systems. This zone tends to insulate incipient fingers from the high-mobility watet;. thus conditions for the onset of fingering differ from those in the parallel plate theory. Finite-difference solutions of the two-dimensional equations of displacement in porous media exhibited the predicted stability characteristics in six hypothetical field- and laboratory-scale floods in rectangular reservoirs. In contrast to results with parallel plate systems, this paper concludes that for waterwet reservoirs, laboratory models scaled by the usual criteria are also correctly scaled for frontal instability. Further, fingering in the systems studied can occur in any saturation range behind the front, and may occur at an intermediate saturation even though stability obtains both at the saturation corresponding to the Buckley-Leverett front and near residual oil saturation. Other points of contrast are that the likelihood of occurrence of fingering may not increase as flow rate or viscosity difference increases, but my be sensitive to changes in he relative permealdlity and capillary pressure functions. INTRODUCTION The recovery of oil by water flooding frequently involves displacing the oil by water of a lower viscosity. Displacement of a fluid by a less viscous one may lead to gross channeling or fingering like that observed in solvent floods, in which it severely lowers recovery efficiency.1 In addition to adverse effects on recovery, it has been suggested that the unstable move- ment causing fingering may interfere with interpretation of scaled model studies of proposed water floods, since the instability in the model might not be faithfully scaled to that in the reservoir prototype. In view of the serious implications of this possible breakdown of widely used scaled model techniques4-6 it is the purpose of this paper to examine the question further. Instabilities in the solutions of systems of differential equations imply a loss of smooth dependence on initial and boundary conditions. Thus, the possibility exists that in using models whose scaling is based on the differential system there may arise size- and rate-dependent factors which are not properly scaled. This possibility was examined in detail by Chouke et al., who analyzed the instability of frontal advance in a related problem, water-oil displacement in parallel plate models,2 in which a moving interface separates two regions of constant, unequal mobilities. First-order perturbation theory predicts the existence of a critical wave length for the growth of perturbations: and a wave length of maximum instability of A, = The interpretation is that wave lengths in a perturbation which are longer than A, will grow. Thus, if the width of a two-dimensional channel is greater than A,, fingers will grow, and the spacing of the fingers which grow at the maximum rate will be approximately It is important that the higher the velocity and/or the difference in flow resistance, the lower is A,, and thus the greater the number of fingers that can grow in a given model. In applying these conclusions to porous media, a pseudo-interfacial tension, , was assumed for the invasion front. Since this would not necessarily be equal to the liquid-liquid interfacial tension, an unknown constant was substituted for the foregoing expression for A, This was determined

Jan 1, 1965

By A. J. McEvily, R. C. Boettner, C. Laird

The processes leading to fatigue failure in the low-cycle range were studied to obtain an understanding of the basis of Coffin's law. Particular attention was paid to the manner of mack nucleation and to the determination of the portion of the total lifetime spent in crack propagation. It was observed that cracks are nucleated at surface notches meated during reversed plastic deformation, Principally at grain boundaries, and that crack growth occupies greater than 75 pct of the fatigue life. When inclusions are present, they are often the primary sites of crack initiation. From consideration of the incremental distances advanced by a crack in each cycle and of the lifetimes of different materials, it appears that the universality of Coffin's law arises from the similarity of crack-tip strain for a given applied strain amplitude, irrespecti7:e of material. STUDIES of low-cycle fatigue have been particularly rewarding in that a simple correlation between fatigue lifetime and plastic-strain range has been shown to exist. This empirical relationship takes the forml,2 Nner=C [1] where N is the number of cycles to failure, Er is the plastic-strain range defined as the total plastic strain per cycle (the sum of the tensile plastic strain plus the compressive plastic strain), n, the exponent, is a constant independent of the material to a first-order approximation and often has a value of about 0.5, and C is a constant which for a number of ductile and metallurgically stable materials is approximately 0.65 for fully reversed strain.3 This relationship has been shown to be valid for a wide variety of metals and alloys and is often referred to as Coffin's law. However, it is yet not clear why so many different metals and alloys exhibit the same over-all behavior in the low-cycle range. Less is known about the mechanisms of fatigue failure in the low-cycle range than in the high-cycle range, but it is expected that the processes would have much in common. For example, an important characteristic common to both high-and low-cycle ranges is that fatigue crack propagation can occupy a considerable portion of the total lifetime. This has been shown by Laird and smith4 by a count of the number of ripples present on the fracture surface of copper specimens, each ripple corresponding to one cycle of load. On the other hand, some marked differences have been noted. For example, in copper,5 nickel, and aluminum4 there is a greater tendency for grain boundary cracking to occur as the strain amplitude is increased. The purpose of the present work was to study further the processes leading to fatigue failure in the low-cycle range and to obtain an understanding of the basis for the validity of Coffin's law. Particular attention has been paid to the manner of crack initiation and to the determination of the portion of total lifetime spent in crack propagation. MATERIALS A number of metals and alloys were tested, although more attention was focused on a particular few, once general trends had become apparent. Alloys were selected so as to investigate the effects of crystal structure, number of phases, stacking-fault energy, and impurity content. The materials in this grouping were: high-purity copper OFHC copper 1100 aluminum (commercially pure) 70-30 brass (Cu-30 wt pct Zn)

Jan 1, 1965

By W. J. Helfrich, C. L. Sollenberger

Many present direct reduction processes utilize iron ore concentrates for the production of sponge iron and the sponge iron is usually preferred as an agglomerate. Pelletizing a high grade iron oxide concentrate prior to reduction is a simple approach to the high capacity production of a uniform and agglomerated sponge iron product. The present paper describes the effect of different atmospheres and physical characteristics on the reduction of pelletized iron oxide concentrates. The kinetics of iron oxide reduction have been I studied by numerous investigators1-4 and a rate equation for the hydrogen reduction of hematite and magnetite has been proposed. Mc Kewan' obtained kinetic data for the reduction of hematite spheres sintered at 13'70°C (2498'F'). The kinetics of magnetite reduction was studied by Quets. He compacted the oxide at a pressure of 1 ton per sq cm and sintered the specimens to densities of 95 to 96 pct of theoretical. These studies demonstrated that the reduction of iron oxide is a topochemical process and that the rate of reduction is proportional to the receding surface area of remaining iron oxide. Thus, the measurement of the kinetics of iron oxide reduction depends upon the development of a measurable and well defined interface between metal and reacting oxide. It is evident that the reduction rate will be constant only when the rate of formation of the reaction product layer is constant. It was recognized by Joseph' and Edstrom' that the physical character of the oxide will determine the configuration of the interface between oxide and metal. The reduction of oxides which exhibit few pores or no interconnected pores should proceed by topochemical processes. However, if a massive or dense crystalline ore cracks during the reduction of the oxide, thus allowing the reduction to Proceed away from a well defined oxide/metal interface, then the observed reduction rate will vary or change during reduction. From a process standpoint, it would be difficult to apply absolute reduction rates determined for crystalline or nonporous iron oxide to the direct reduction of a porous pellet of iron ore concentrate. Most present direct reduction processes are attempts to utilize iron ore concentraates in the production of sponge iron. The sponge iron is usually preferred in an agglomerated form for subsequent melting practices. Pelletizing a high grade iron oxide concentrate prior to reduction is a simple approach to the high capacity production of a uniform and agglomerated sponge iron product. One process for preparing such a product has already been de-scribed. This paper describes the effect of different atmospheres and physical characteristics on the reduction of pelletized iron oxide concentrates. ORE CONCENTRATES STUDIED A magnetite concentrate from Reserve Mining Co. and a hematite concentrate from U.S. Steel Corp. were used in preparing pellets for these experiments. A head analysis of these concentrates is given in Table I. To compare the reduction rates of pelletized concentrates, it was necessary that each pellet contain equal amounts of oxygen which could be removed by reduction, regardless of the oxidation state, i.e., Fez% or Fe304. Also, because of the difficulty in preparing spherical pellets of uniform size, weight, density, and porosity for experimental purposes, pellets of a cylindrical shape were formed in a press.

Jan 1, 1961

By B. H. Caudle, A. R. Khan

The oil production performance of thin-oil-column type reservoirs producing entirely by natural water drive are discussed in this paper. The dimensionless production histories were obtained in a series of scaled porous-medium laboratory-model tests under varying effects of rate of production, mobility ratio across moving boundary, ratio of oil- and water-sand thickness and mobility ratio across fixed boundary including the effects of different residual oil saturations in the water-invaded regions. This scaled-model study of thin oil columns indicated that greater recovery of oil results for lower values of the following: rate of production, mobility ratio across the moving boundary, water-sand thickness relative to the oil-sand thickness, and mobility ratio across the fixed boundary or residual oil saturation in the water-invaded region of oil sand The results of the study are given as dimensionless production history plots which may be used to predict the performance of the field they represent. INTRODUCTION Oil fields producing under natural water drive can be classified into two major groups: (1) an "edge-water drive" reservoir group in which the water is encroaching from the reservoir periphery, and (2) the group of reservoirs in which the oil-sand is entirely underlain by an aquifer. Reservoirs in the second group are further subclassified as "bottom-water drive" reservoirs and "thin-oil-column" reservoirs. The latter reservoirs are characterized by their relatively equal oil- and water-sand thicknesses and by having both the oil and water zone relatively thin. When oil is produced from a water-drive reservoir by a partially penetrating well, water breaks through sooner or later before all of the displace-able oil is recovered. Many analyses of this general situation have been made with varying degrees of realism in their approaches, both theoretical and experimental5.6 in nature. Pertaining to the scaled laboratory model studies of natural water drives using porous media, detailed investigations have been made by Henley, et al.7 on bottom-water drive and by Caudle8 on edge-water drive reservoirs. No laboratory model studies, however, have been made to investigate the performance characteristics of thin-oil-column reservoirs. This paper presents the results of a study on scaled models of thin-oil-column reservoirs, operated under varying conditions corresponding to the actual field practices. SCALING PARAMETERS The object of scaling is to produce a one-to-one relationship between physical and dynamic parameters of the model -and prototype. The most comprehensive treatment of scaling techniques for oil reservoirs is given by Geertsma, et a1.,9 which covers two methods for arriving at dimensionless groups for scaling. Table 1 lists the scaling parameters suitable for creating kinematic and dynamic similarities between the model and the prototype in the variables considered for the thin-oil-column reservoirs. In order to achieve similarity it is essential that the dimensionless groups listed in Table 1 have equal values in the model and the protorype simulated. Since it has been established that the role of interfacial forces is neglible in gross fluid movement in a reservoir no attempt was made to scale capillary forces- Although there are reports describing Scaling criteria for these forces,10,ll consideration of the necessary but restrictive conditions for

Jan 1, 1970

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 C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By C. H. Stewart, R. A. Fagin

A two-dimensional, three-phasereservoir simulator was programmed for a large memory digital computer. It was designed to provide a practical solution to describing the complex physical relation between the natural forces and the physical properties of a heterogeneous reservoir when subjected to a specific set of conditions. A reservoir study is briefly described to illustrate application of the model. A full volumetric account of three phases (oil, gas and water) is performed simultaneously throughout an integration net representing the reservoir. Absence of one or two of the phases is treated as a special case of the more general situation. Expansion (or contraction) of all phases, including rock expansion, is performed so that the pressure calculation is the general unsteady-state case. To account for the large variations of subsurface elevation encountered in some reseruoirs, and to allow for segregation of the various phases, a gravity head term is included in the basic drive potential. Appropriate fluid and rock properties are used in polynominal surface form (functions of pressure and/or depth) or they can be entered as space variables at each position of the integration net. An unsteady-state water influx calculation, based on the method of van Everdingen and Hurst, l was connected to the boundary of the matrix to simulate aquifers of various sizes. In addition to reservoir calculations, three-phase fluid flow from the producing depth to the wellhead, including provisions for gas lift, was incorporated in the simulator. A workover routine was also built which can automatically switch to a different set of production relations when a gas-oil ratio or water fraction reaches a limit; or it can shut-in the well if prescribed. INTRODUCTION This paper describes a reservoir engineering mathematical simulator used to represent the complex interaction of natural forces and physical properties of a reservoir during natural depletion or with various injection schemes. The simulator, which was programmed for a large memory digital computer, is a two-dimensional calculation which handles three mobile fluid phases simultaneously (oil, gas and water). Basic requisites for the method are individual well production and pressure data, hydrocarbon fluid properties, geological data (producing depth and net sand), capillary pressure data, relative permeability data and permeability and porosity information. Matching the past performance of a combination drive reservoir often has yielded information concerning continuity and the validity of basic data. Detailed predictions of future performance can be made for continuation of current depletion methods (natural depletion) as well as for various types of recovery by gas or water injection. Combination injection cases and pattern studies can also be performed. Workover programs, gas lift and different types of artificial lift programs have been investigated using a technique similar to that described by Kern and Nicholson2 except that conditions of pressure and saturation at the block within which the well is located are used rather than average reservoir conditions. Drilling additional wells to optimize profit was explored, both as to number and location, by placing wells at different spots within the reservoir matrix. Special depletion processes can be examined, such as upstructure drainage and lateral (or strike) waterfloods in thin oil columns. In one case the mathematics of the simulator were modified to calculate the displacement in the vertical plane rather than in the horizontal plane. In this manner specific reservoir problems can be studied, such as coning of gas and/or water around production points, fingering along permeable stringers or, more generally, frontal advances in a heterogeneous section.

Jan 1, 1967

By G. M. Gay, R. L. Perrine

This paper describes a method of numerical computation for three-dimensional, unstable, miscible displacement behavior useful for heterogeneous systems, as well as for more ideal conditions. In the method, flow equations are first linearized by a perturbation approach. The basic flow process is separated and a solution for its behavior readily obtained. The remaining problem of deviations from the basic flow caused by non-ideal conditions is then subjected to numerical analysis. Results obtained from use of the method are also presented. Although conditions assumed in the test cal-culations were severe, results show the type of dispersing flow expected and appear quite satisfactory. The method has eliminated or reduced in importance problems of oscillating values near steep fronts, excessive computer smoothing, etc. A unique advantage of the method is that the source of variations from ideal behavior can be observed. ,The one serious drawback results from the algebraic complexity of the perturbation approach, and the need for second-order terms to be retained in calculations of interest. Fewer array points are available and more computer time is required than would be desired. However, these difficulties are also experienced with other approaches to the solution of three-dimensional dispIacement problems. INTRODUCTION Prediction of behavior of the miscible displacement process within a porous medium for any system of engineering importance is plagued by a number of difficulties. With typical fluid properties the displacing fluid has the lower viscosity, and there is a natural tendency toward flow instability. The problem of predicting instability is compounded by the fact that every real system is heterogeneous. Permeability will vary from point to point — not entirely systematically, and yet not in a random fashion leading to a readily defined average. Furthermore, permeability is unlikely to be isotropic. These permeability properties, which are characteristic of any real porous medium, accentuate the effects of flow instability. Another factor to be considered is that flow dispersion accompanies the displacement process in a porous medium. Mechanistically, dispersion is due to the fact that flow between any two points in the medium follows multiple tortuous paths, each characterized by slightly different flow properties. The fact that the coefficient characterizing dispersive properties of the real medium JS a tensor with variable elements complicates matters. Further difficulties arise because the parabolic influence, while not negligible, is small. Thus, there is a tendency toward steep and uneven displacement fronts along which dispersion smoothing must be represented accurately. Yet frequently used numerical analysis schemes may tend toward either instability, oscillation or excessive smoothing, none of which gives the desired accurate picture of flow behavior. Thus, while many experimental and analytical studies of the process have been made, predictions of actual performance are still subject to considerable uncertainty and possible improvement. This paper reports on part of a study which attempts to develop improved methods for solution of the flow equations describing the miscible displacement process. Of necessity, the calculations were performed on a large digital computer. A three-dimensional system is represented, and the coefficients defining dispersion and permeability can be varied in a manner representative of a real system. Because the differential equations are nonlinear, a first step is to linearize them using a perturbation approach. Objective of this approach has been to separate the basic flow which would occur in an ideal system and for which predictions are readily made. Then the numerical analysis, which is necessarily complex, is directed toward solution of the perturbation equations. These equations state how the basic flow is modified

Jan 1, 1967

By D. R. Mitchell

MEETING for the thirteenth time in New York as part of the five-ring circus known as the Annual Meeting A.I.M.E., the Coal Division experienced a wartime boom in attendance. Technical sessions were crowded and whereas 65 seats were set up for the annual luncheon meeting of the Division over 100 members showed up, many of whom had to be fed elsewhere. Possibly the Coal Division should adopt "13" as its lucky number. Discussions at technical sessions, formal and informal, were lively. Those papers containing material directly concerned with the war, or containing information that might be usable for war industries, were enthusiastically received. Some criticism was made in regard to the present method of handling discussions and it was suggested that better success would be had in getting discussions in the record by the use of hired reporters. This is being considered by the executive committee of the Coal Division.

Jan 1, 1943

By AIME AIME

EUGENE AUSTIN STEPHENSON, better and affectionately known as "Steve," has long been an active and enthusiastic member of the Petroleum Division of the Institute. As its Secretary-Treasurer for several years, he was the spark plug that so successfully carried the Division over the trying period of the early thirties. "Steve" was born at Franklin, Mich., July 27, 1883; graduated from Adrian College in 1905, and later received his Doctor's degree from the University of Chicago where he was an instructor for two years. From 1917 to 1925 he was chief geologist for the South Penn Oil Co. and from 1925 to 1930 he was associated with Ralph Davis, Inc., of Pittsburgh, consulting engineers. In 1930 the Missouri School of Mines called him to Rolla to head a department of petroleum and natural gas engineering. and in 1938 he was invited to build up a similar department at the University of Kansas.

Jan 1, 1941

SEC. 1. This Institute is incorporated under the Membership Corporations law of the State of New York; its corporate name is American Institute of Mining and Metallurgical Engineers, Inc.; and its abbreviated title A. I. M. E. SEC. 2. The objects of this Institute, as set forth in its certificate of incorporation, are: "To promote the arts and sciences connected with the economic production of the useful minerals and metals, and the welfare of those employed in these industries by all lawful means; to hold meetings for social intercourse and the reading and discussion of professional papers, and to circulate by means of publications among its members the information the obtained, and to establish and maintain a place for meeting of its members, and a hall for the reading of papers and delivery of address, and a library of books relating to subjects cognate to the sciences and arts of mining and metallurgy."

Jan 1, 1943

State of Louisiana, Department of Conservation, Division of Minerals, 126 New Orleans Court bldg , New Orleans, La. Dr. J. A Shaw, Director The Bureau of Scientific Research of the Minerals Division has recently issued the first of a proposed series of Geological Bulletins on the geology of Louisiana, by parishes The bulletin may be had for 25 cents postage by writing to the above address: Geological Bulletin 1, Geology of Iberia Parish (1931) Other bulletins of a more general nature have been published in the past. Bulletin 13, A brief survey of the natural resources of Louisiana (1930) (revised 1931), 14, The clays of the Shreveport area (1927), 16, The clays of Monroe area (1928), 19, The clays of Alexandria area (1929). The latest (Tenth) Biennial Report of the Department of Conservation for 1930-1931, containing a review of the natural resources of the state in the report of the Division of Minerals, may be obtained on payment of 25 cents for postage

Jan 1, 1933

Through the courtesy of Mr. J. H. Batcheller, we are able to give the following extracts from a letter written by Mr. Caetani: It has been a nice idea of yours to collect for me all those scraps of news I received a few days ago and to send me such a long interesting letter. In four years war I have had practically no news. Burch wrote several times but several of his letters were lost and so mine. I believe I even wrote to you. Much has changed since our last days in Telluride; here and there many have disappeared. Of myself I will tell you in a few words: for a year I was in Cadore(?), where I had the good luck to do to the others what they wished to do to me, and had good luck on my side. Then I went to Macedonia for 5 months, did little fighting but lost my faithful fellow officer and companion.

Jan 6, 1919

Discussion of the paper of E. W. Parker, presented at the New York meeting, February, 1913, and printed in Bulletin No. 75, March, 1913, pp. 443 to 451. A MEMBER :-I would like to call attention to the fact that the speaker has failed to indicate where is the center of gravity of the membership of the Institute, while he is dealing with this map. MR. PARKER :-I think you will find that a little further east. CHAIRMAN JAMES F. KEMP :-In connection with Mr. Parker's paper, I would like to inject a little matter of history into our program here, as regards the founding of the Institute. The credit for the founding of the Institute really applies to miners who were engaged in the anthracite coal region. The three men who signed the circular for the call for the meeting which brought the engineers to Wilkes-Barre in May, 1871, were Eckley B. Coxe, Richard P. Rothwell, and Martin Coryell.

Jan 5, 1913

Sec. 1. This Institute is incorporated under the Membership Corporations law of the State of New York; its corporate name is American Institute of Mining and Metallurgical Engineers. Incorporated; and its abbreviated title A. I. M. E. Sec. 2. The objects of this Institute, as set forth in its certificate of incorporation, are: "To promote the arts and sciences connected with the economic production of the useful minerals and metals, and the welfare of those employed in these industries by all lawful means; to hold meetings for social intercourse and the reading and discussion of professional papers, and to circulate by means of publications among its members the information thus obtained, and to establish and maintain a place for meeting of its members, and a hall for the reading of papers and delivery of addresses, and a library of books relating to subjects cognate to the sciences and arts of mining and metallurgy."

Jan 1, 1940

SEC. 1. This Institute is incorporated under the Membership Corporations law of the State of New York; its corporate name is American Institute of Mining and Metallurgical Engineers, Inc.; and its abbreviated title A. I. M. E. SEC. 2. The objects of this Institute, as set forth in its certificate of incorporation, are: "To promote the arts and sciences connected with the economic production of the useful minerals and metals, and the welfare of those employed in these industries by all lawful means; to hold meetings for social intercourse and the reading and discussion of professional papers, and to circulate by means of publications among its members the information the obtained, and to establish and maintain a place for meeting of its members, and a hall for the reading of papers and delivery of address, and a library of books relating to subjects cognate to the sciences and arts of mining and metallurgy."

Jan 1, 1942

WE the undersigned, being all persons of full age and citizens of the United States, and a majority residents of the State of New York, desiring to form a corporation pursuant to the provisions of the Membership Corporations Law for the purpose of incorporating, as provided in Section 5 of Article I. of said law, the existing unincorporated association known as American Institute of Mining Engineers, do hereby make, acknowledge and file this Certificate for that purpose, and do Certify as follows I. That the American Institute of Mining Engineers is an unincorporated association organized and existing with the object of promoting the arts and sciences connected with the economic production of the useful minerals and metals and the welfare of those employed in these industries by means of meetings for social intercourse and the reading and discussion of professional papers, and to circulate by means of publications among its members and associates the information thus obtained.

Jan 1, 1923