Search Documents

By D. D. Donald

C. J. Barber (U. S. Smelting Refining & Mining Co., Salt Lake City)—In his paper Donald describes how New Jersey Zinc Co. made surveys for a connection between the Ivanhoe and Van Mater shafts at Austin-ville, Va. Except to say that the two faces had to meet "accurately" on line and grade, Donald does not indicate the required precision. Assuming that there were 24 angles in the 11/2-mile traverse and 15 in the one-mile traverse, it can be shown that if the average error in plumbing each shaft was 230" and the average error in measuring each traverse angle was 210". then the average error at the point of -connection would have been about ±1.9 ft normal to the line between the shafts. This calculation assumes that any errors in the triangulation would be negligible compared with the errors in the plumbings and traverses, and it also neglects taping errors. With no constant errors or blunders, the latter would be important only in lines normal to the line between the shafts. To make the average error at the connection less than 1.9 ft would, therefore, require either reducing the error in the plumbings to less than ±30", or that in the traverse angles to less than ±l0", or fewer stations, or a combination of these. Referring briefly to the triangulation, because of the problem of fitting a new triangulation into older surveys of the district the orientation deserved some mention, even though the connection could have been made with an assumed bearing. It would be interesting to know how many triangles were required and what the average summation error was before making any adjustments and without considering the algebraic signs. Perhaps this is referred to indirectly in the statement that the maximum angular error distributed was 2". Turning to the shaft plumbings, it would be helpful to know how many men were employed and how long each shaft was in use. Donald says that the surface positions of W-2 and W-3 were carefully surveyed from the collar position of W-1, without indicating how this was done. The length of the backsight would be particularly important. There must have been some error in setting W-1 vertically below the stations in the headframes. How immovable were the headframes, especially the Van Mater, which appears higher than the Ivanhoe and subject to more vibration because of skip hoisting? Donald does not say whether the plumbing wires had been previously restraightened to minimize spinning (otherwise they behave like weak helical springs). The use of light steel weights is most surprising because there seem to be excellent reasons for using heavier, nonmagnetic weights. Did the shafts contain no steel sets, pipes, power cables, etc., which might attract steel? The plumbing method described by Donald was designed for deep shafts in South Africa but differed from the South African practice in two important respects. As described by Browne,6 in South Africa the line between the wires was made parallel to the long axis of the shaft, whereas in the Ivanhoe shaft the lines between the wires were diagonally across the shaft. The main reason given for the South African practice is to insure that the gravitational attraction between the wires and the rock walls is the same on both wires, and therefore does not affect the bearing of the line between them. It seems probable, however, that the effect of air currents might be minimized in the South African procedure, and might be serious with the wires in the diagonal position at the Ivanhoe shaft. In the South African case cited by Donald the wires were swinging freely (although the plumb bobs were sheltered from air currents) but in the Ivanhoe case they were dampened with the plumb bobs set in water. In the discussion of Browne's paper R. St. J. Rowland said:' It has been the practice for a long time to damp the oscillations by immersing the bobs in oil or water. The time per oscillation is thus increased, thereby extending the time taken to complete the work. The longer the suspended wire the less there is to recommend the practice . . . The theoretical time for one swing of a simple pendulum 1050 ft long is approximately 36 sec, which would be increased by dampening the plumb bob in water. Hence very few complete swings would be observed in the 5 min intervals used at the Ivanhoe shaft. In the two South African cases described by Browne, the length of plumb line in one shaft was 5425 ft, the calculated period of swing was 81.6 sec, the average actual period was 76.6 sec, and 94 complete swings were observed in 2 hr. In the other case the length of plumb line was 3116 ft, the calculated period of swing was 61.8 sec, the average period was 63.5 sec, and 86 complete swings were observed in 1 hr 31 min. Browne concluded that observations of more than 30 swings are not likely to result in sufficient gain in accuracy to be justified. Returning to the Ivanhoe and Van Mater plumbings, an objection to the method used is that all four azimuths were taken from fixed points instead of swinging wires, and that each pair of observations would— barring blunders— check closely, and so perhaps give a false feeling of security. In fact, it seems that only two azimuths were obtained from one plumbing, and not four as stated by Donald. Nevertheless, the tying in of each pair of wires from both sides of the shaft has much to commend it. Donald's description leaves the impression that if each shaft was plumbed only once, the engineers were fortunate indeed if the average error in the underground orientation was as little as 30". Because the survey was done over a period of three years, it seems likely that the plumbings were repeated, perhaps more than once. The underground traverse angles were measured by conventional methods, but because the number of angles in the overlapping traverses was not given, the angular closure given by Donald does not indicate the accuracy with which this was done. Donald's description of a method of taping lines of irregular length is welcome. The literature on taping is usually confined to lines of about one tape length, generally 100 ft. Such lines are rare in metal mining because the time, trouble, and cost of setting points at 100-ft distances underground are not warranted. (Nevertheless civil engineers may go to this expense

Jan 1, 1960

By J. C. Cordell, C. F. Weinaug

A study of the behavior of retrograde condensation from gas mixtures was made in the presence and absence of sand in order to determine if the condensed liquid would re-vaporize in the Presence of sand. Methane-butane and methane-pentane mixtures that would form liquids by retrograde condensation when produced from a constant volume cell were used. The methane-butane mixtures of similar composition were produced by three different methods. The first two were charged into an empty cell and were produced in one case within a period of 9 hr and in the other within a period of 3 days. The third mixture was charged to a sand-packed cell and produced within a period of 9 hr. The curves relating composition of the produced gas to pressure obtained from these experiments show that equilibrium was maintained as long as liquid was condensing. However, during the portion of the pressure decline where the condensed liquid should be revaporizing, equilibrium was maintained only when the mixture was produced from the sand-packed cell. The methane-pentane mixture is produced only in the presence of sand. The data obtained for this system also show that equilibrium is maintained at all times during the pressure decline. These results indicate that revaporization is aided rather than prevented by the fact that the condensate "wets the sand." The significance to the petroleum industry of the behavior of hydrocarbons in the retrograde region is becoming increasingly important because a large percentage of reservoirs being discovered today arc of the gas-condensate type. The important characteristic of a gas-condensate reservoir is the retrograde condensation of a liquid phase throughout the reservoir if the pressure is allowed to decline. In order to prevent the loss of this retrograde liquid in some of these reservoirs, they are "cycled"; that is, the material produced from the reservoir is processed to remove the heavier hydrocarbons and the light fractions are returned to the reservoir to maintain the reservoir pressure. Another method of producing these reservoirs is Pressure depletion. There are other methods but they are not in general use. Although a large number of factors are involved in determining the optimum method of producing a reservoir, only the problem of whether the retrograde condensate resulting from pressure depletion will revaporize from sand at equilibrium conditions is considered in this paper. Also the data are limited to the methane-n-butane and the methane-n-pentane systerns at loo°F and one type of sand. Three types of tests were run in order to compare the revaporization of retrograde liquid formed in a cell with that formed in a sand-packed cell, and to establish the effect of time. The results of tests presented in this paper indicate that both time and the presence of sand promote the revaporization of

Jan 1, 1949

By A. K. Snelgrove

This symposium is designed to marshal a state-of-the-art account of approaches and solutions to problems arising out of waterlogging and salinization in the world's largest irrigation system which was dismembered by legal beheading of the Eastern Rivers under the lndus Waters Treaty of 1960. In this introduction, the geological and historical backgrounds of surface waters and groundwaters in the Indus Basin are presented with a summary of remedial measures undertaken and proposed to combat waterlogging and salinization. At a time when water-resource planning is proceeding on a continental scale, a study devoted solely to the Indus River Basin may seem to require an apologia. None is offered, however, for three major reasons: (1) Although the Indus is a single river basin, it is one of great magnitude and importance. The Indus River, which comes from Himalayan sources to form part of the great Indo-Gangetic Plain drainage system, is vital to agriculture and hydroelectric power in the thickly populated northern subcontinent of Indo-Pakistan. (2) The Indus has become a plexus of extraordinary problems since the dismemberment of the world's largest irrigation system of which it was a part. The break-up of this system occurred with the legal beheading of the Eastern Waters under the Indus Waters Treaty of 1960. (3) Th River Basin has been the target of a Herculean effort to contain the twin encroachments of waterlogging and salinization that now threaten it with a shortage of fresh water. The British irrigation system of barrages and canal on the Indus was a venture admirable in conception. Unfortunately, it is now a legacy burdened with unwanted consequences due to interference with the hydrological regimen, particularly the rise of water table occasioned by leakage. To cope with this situation, there has been a possibly unprecedented concentration of domestic and international engineering and scientific talent on the Indus and its problems (about a score of organizations by recent count). To marshal a state-of-the-art account of their approaches and solutions is the purpose of this Symposium. The coordinated effort now being made on a variety of fronts under the general supervision of West Pakistan Water and Power Authority (WAPDA) is directed toward the modification and control of interactions and relationships between geological materials — soils, water and its salt balance, and air —in what has been termed a huge "underground lake" in an alluvial basin. The problem is an intriguing one in that a highly complicated case history of river-course changes visible in the present surface must be projected downwards through Recent and Pleistocene sediments in order to decipher the environmental controls of deposition in terms of vertical and lateral variations in lithology and patterns. Such studies must be extended to compaction, deformation, and fabrics of the alluvium. Geohydrology is often defined as groundwater hydrology. As employed herein, it includes the distribution, quality, and behavior of surface and subsurface waters, including their evolution, vagaries, and geological, climatological, and botanical milieux." GEOGRAPHY The Indus (Sanskrit, Sindhu = river) rises at nearly 17,000 ft above the sea in the glaciers on the northern slopes of Kailas Parbat, Tibet. It is about 1800 miles long and has a total drainage area estimated at 372,000 sq miles. The average discharge is 196,000 cfs at mouth. Leaving the Plateau of Tibet, the river flows at a level which is miles above the sea, and in places miles below Himalayan peaks, in a structurally controlled northwestward course. Joined in this course by six tributaries and swollen by the melt waters of giant glaciers, it swirls about the base of the famous

Jan 1, 1970

By E. W. Engelmann

A NEW and interesting type of drier has been developed and operated at the Magna plant of the Utah Copper Co. for the past year for the drying of a filtered concentrate in the molybdenum recovery plant of this company. The several units installed have ranged from 12 to 30 in. in width, and from 14 to 30 ft. long. Their economy was immediately apparent after the first one was installed. The machine is not only low in capital cost, but has proved efficient for the drying of fine concentrate at low operating cost. The Utah vibrating drier consists of a steel conveyor trough supported on flexible reeds from a structural steel foundation anchored in concrete. The steel pan, which is also the drying surface, is vibrated by a heavy Utah electromagnetic vibrator fastened either to the head or discharge end of the pan. The vibration thus obtained agitates the material and causes it to be conveyed over the heated surface.

Jan 1, 1938

By J. L. Tomlinson, B. D. Lichter

Electrical resistivities 01 liquid Cd-Bi, Cd-Sn, Cd-Pb, In-Bi, and Sn-Bi alloys were measured using an electrodeless technique. The resistivities ranged from 50 to 160 microhm -cm, temperature dependences were positive, and no sharp peaks in the composition dependence of the resistivity were observed. On the basis of these observations, it was concluded that the alloys are typical metallic liquids. The electron con-cent9,ation was calculated from the measured resis-tizlity and available thermodynamic data using a model which attributes electrical resistivity to scattering by density and composition flzcctuations. A correla-tion was shown between the departure of the electron concentration from a linear combination of the pure component valences and the value of the excess integral molar free energy. Calculation of the temperature dependence of the electrical resistivity showed a need for more detailed thermodynamic data in these systems and led to suggestions for improvement in the concept of residual resistivity in the fluctuation scattering model. THE electrical resistivity of liquid metals provides information regarding interatomic interactions and their effects upon structure. In this experiment an electrodeless technique was used to measure the electrical resistivities of liquid alloys of Cd-Bi, Cd-Sn, Cd-Pb, In-Bi, and Sn-Bi, and the results were used with thermodynamic data to calculate a parameter which reflects the tendency toward localization of electrons due to compositional ordering. It was found that the resistivities of these alloys are generally metallic in magnitude and temperature dependence. The electrical and thermodynamic properties are discussed in terms of the fluctuation scattering model'22 which supposes that the electrical resistivity arises from scattering due to a static average structure and departures from the average due to fluctuations in density and composition. Further, this model is compared with the pseudopotential scattering model of Ziman et al.3-5 EXPERIMENTAL PROCEDURES Alloy samples were prepared from 99.999 pct pure elements obtained from American Smelting and Refining Company (except tin which was obtained from Consolidated Smelting and Refining Company.) J. L. TOMLINSON, Member AIME, formerly Research Assistant Division of Metallurgical Engineering, University of Washington, Seattle, Wash., is now Physicist, Naval Weapons Center, Corona Laboratories, Corona, Calif. 0. D. LICHTER, Member AIME, is Associate Professor of Materials Science, Department of Materials Science and Engineering, Vanderbilt University, Nashville, Tenn. This work is based on a portion of a thesis submitted by J. L. TOMLINSON to the University of Washington in partial fulfillment of the requirements for the Ph.D. in Metallurgy, 1967. Manuscript submitted May 31, 1968. EMD Weighed portions were sealed inside evacuated silica capsules, melted, and homogenized before the resistivity was measured. The resistivity of a liquid alloy was measured by placing the sample inside a solenoid and noting the change in Q. According to the method of Nyburg and ~ur~ess,~ the resistivity of a cylindrical sample may be determined from the change in resistance of a solenoid measured with a Q meter as T7--5--W =R7JT^ ='Kc-lm(Y) [1] where L, R, and Q = wL/R are the inductance, series resistance, and Q of the solenoid. The subscript s refers to the solenoid with the sample inside; the subscript 0 refers to the empty solenoid. Kc is the ratio of the sample volume to coil volume and y = 2 [bei'0(br)-j ber'o(br)~\ br\_bero(br) +j bei0 (br) expressed with Kelvin functions which are the real and imaginary parts of Bessel functions of the first kind with arguments multiplied by (j)3'2. The argument of the function Y is hr where r is the sample radius and b2 = po~/p, i.e., the permeability of free space times 271 times the frequency divided by the resistivity in rationalized MKS units. Since Eq. [I] cannot be solved explicitly for p, values of Kc. lm(Y) were tabulated at increments of 0.1 in the argument by. A measurement of Q, and Q, determined a value of Kc . lm (Y) and the corresponding value of br could be read from the table. From the known r, uo,, and w, the resistivity, p, was determined. The change in Q was measured after letting the encapsulated Sample reach equilibrium inside a copper wire solenoid. The solenoid was contained in an evacuated vycor tube in order to retard oxidation of the copper while operating at high temperatures and heated inside a 5-sec-tion nichrome tube furnace capable of obtaining 900°C. Temperature was determined with two chromel-alumel thermocouples, one in contact with the solenoid 30 mm above the top of the sample and the other inserted in an axial well at the other end of the solenoid and secured with cement so that the junction was 2 mm below the bottom of the sample. Temperature readings were taken with respect to an ice water bath junction, and the voltage could be estimated to the nearest thousandth of a millivolt. The lower thermocouple was calibrated by observing its voltage and the Q of the coil as the temperature passed through the melting points of samples of indium and tellurium. The upper thermocouple reading was systematically different from the lower thermocouple reflecting the temperature difference due to a displacement of 60 mm axially and 6 mm radially. Calculations show that the gradient over the sample was less than 2 deg. Q was measured by reading a voltage related to Q from a Boonton 260A Q meter with a Hewlett Packard

Jan 1, 1970

By R. E. HORE, J. B. Tyrrell

IN the vicinity of Kirkland Lake, northern Ontario, several gold mines are producing, and three are at present being operated very profitably. At these three mines and on the adjoining properties east and west of them there have been important developments during the past year, and additional equipment has been installed to permit a large increase in output. At the end of 1924 the total production amounted to $13,688,600. For the year 1926 the production is estimated at about $19,000,000, and during 1926 the output should be much greater than that of 1925. The mines which are being operated most profitably and are producing a large part of the ore now being milled are the Wright-Hargreaves, Lake Shore and Teck-Hughes. Others which have contributed smaller amounts are the Tough-Oakes-Burnside and the Kirk- land Lake Gold mines.

Jan 1, 1926

By Dudley L. Davis, E. C. Peterson

Introduction The Darwin District is 30 miles east of Olancha which is 220 miles north from Los z4ngeles via U. S. Highway No. 6. The ore deposits occur in the Darwin hills that have been elevated above the Darwin plateau by typical basin and range-type faults. Some of these faults have been active quite recently as shown by the displacemerit of Quaternary (?) basalt flows and the rejuvenation of erosion in Darwin wash. The ore bodies occur as fissure fillings and replacement of favorable beddings in a thick series of Paleozoic limestones, dolomites, shales and quartzites. Rich oxidized lead-silver ores were dis'Overed in the early seventies. BY 1875, Darwin had a population Of about 5000, and by 1880 several small mills and smelters had been built. Early exhaustion of surface ores, the isolated location of the district, and fluctuation Of metal prices caused intermittent operations until World War 11. The district production is estimated at about $3,000,000 prior to 1900 and perhaps $4,000,000 between 1900 and 1945 when Anaconda purchased the major producing mines. 'or the Past two Years, daily production has averaged 75 tons of direct shipping and 150 tons of milling-grade ore. At present, 300 tons of mixed oxide and sulphide ore are treated by flotation, the oxidized lead minerals being activated by addition of sodium sulphide after galena and sphalerite have been recovered from the circuit. Extraction averages 85 pct of the lead, 80 pct of the silver, and about 33 Pct Of the zinc in a bulk concentrate. The mill is designed so that minor changes permit selective flotation of pyritic lead-zinc ores encountered at depth. General Geology The oldest rocks in the Darwin hills are a series of Paleozoic limestones, dolomites, shales and quartzites probably pennsyl-vanian in age. These rocks have been intruded on the Rest by the Coso granite batholith and on the east by a granodiorite stock. Numerous sills and dikes grading from orthoclase granite to gabbro are exposed in the underground workings. The sedimentary rocks have been considerably folded and faulted, the most prominent structure being a northwest pitching anticlinal fold, the crest of which lies just west of the ridge line of the Darwin hills, The east limb of this fold has been intruded by the granodiorite stock and east of the stock the sediments show several closely spaced anticlines and synclines. The minor flexures on the west flank of the major fold are structurally important in localizing the ore in the Defiance and Essex mines of the Darwin group. There are three principal systems of faulting and fissuring, most of which are post-intrusive and premineral in age, with minor amounts of postmineral movement. The Darwin tear fault is the largest, most persistent fault in the district. It strikes N 60 to 75°W, dips steeply south and is traceable for about 10 miles. It is a normal

Jan 1, 1949

By R. J. Burgess, A. R. Ramey, A. R. Adams

During interpretation of pressure buildup tests on gas wells in a tight dolomite gas reservoir, peculiar behavior was noticed. Two straight lines were apparent. Effective permeability to gas taken from either straight line was about the same, and the Miller-Dyes-Hutchinson dimensionless time check for the straight line was proper for both straight lines. Geological data indicated the likelihood of scattered trending fractures in the reservoirs. Since the first straight Iine yielded permeability values close to the geometric mean permeability from core analyses, it was postulated that the reservoir model was that of an acidized well completed in the tight dolomite, but that widely scattered hairline fractures caused the mean permeability of the reservoir distant from the well to be higher than the matrix permeability. Because all other studies of fractured reservoirs to the authors' knowledge assumed that the fracture matrix was dense enough to communicate directly with the well, no interpretative methods were available. The Hurst line-source solution for a radial change in permeability for interference between oil reservoirs was adapted to pressure buildup testing. The result indicated that the first straight line should yield the proper matrix permeability and wellbore skin effect. The second straight line may be extrapolated to obtain static pressure. The time of bend between the straight lines was used to estimate distance to a fracture. Application to field test data is shown. It is believed that the methods developed and the case history presented will add to present tools available for pressure buildup interpretation. Introduction Since the pioneer studies by Miller, Dyes, and Hutchin-son1 and Horner' in 1950 and 1951, well test analysis has become recognized as one of the most powerful tools available to both production and reservoir engineers. Well test analysis serves as a logical basis for well stimulation and completion analysis, and for long-term reservoir engineering. Since the early 19501s, much effort has been placed on the development of well-test analytical methods. Reservoir and well conditions of increasing complexity have been considered systematically to provide the analyst with a catalog of causes and effects. Matthews and Russella state that some 200 papers dealing with this subject have been published in the last 35 years. Developments in well test analysis appear to have originated in one of two ways. Either a physically realistic field condition was anticipated and analytical solutions for the condition achieved, or anomalous field test behavior was recognized and interpretative methods sought for the anomaly. In recent years, it has appeared that the latter has inspired an increasing number of studies. The analyst today finds an increasing number of known cause and effect studies available for well test analysis, the classic of which is that of finding the specific flow problem that generated the answer — the well behavior. Although it may be impossible to achieve this goal uniquely, the analyst often is able to select a useful interpretation that combines all known performance and geologic data — or to show that various logical alternatives would not significantly affect the interpretation. During a recent reservoir study, we observed gas well test behavior that did not appear to fit behavior described previously. Although it cannot be said that we have found a unique interpretation, we shall present in this paper the peculiar behavior observed, and describe the reservoir and interpretative methods developed. Reservoir Description The subject gas reservoir is a 9-mile-long, narrow dolomite reservoir lying within a limestone of Ordovician age. (See Fig. 1.) The dolomitized rock in the field consists of dark brown to buff, dense to coarsely crystalline, vugular dolomite containing numerous hairline fractures, many of which may have been closed in the reservoir and parted when cores were brought to the surface. Larger fractures are also apparent in core, but usually are filled and sealed with euhedral dolomite crystals. Portions of the north flank of the reservoir are known to be cut by a sealing fault downthrown to the north. Gas wells located near the fault have higher open flow potentials than those more distant from the fault. This is believed to be a result of higher permeability near the fault due to more extensive and open fractures. Detailed coring and core analysis have been performed on several of the wells in this reservoir. Fig. 2 presents permeability variation' plots for both horizontal and vertical

Jan 1, 1969

By R. G. Agarwal, Ramey Jr. H. J., Al-Hussainy R.

Although it has long been realized that gas recovery from a water-drive gas reservoir may be poor because of high residual saturations under water drive, it appears that only limited infomlation on the subject has been available until recently. This study was performed to show the qiiantitative potential importance of water influx. Results indicate that gas recovery may be very low in some cases: perhaps as low as 45 per cent of the initial gas in place. Gas recovery under water drive appear to depend in an important was on: (I) the prodirction rate and manner of production; (2) the residual gas saturation; (3) aquifer propertie.); and (4) the volumetric displacement effciency of water invading the gas reservoir. The manner of estimating water-drive gas reservoir recovery can vary considerably. Examples are: the steady-state tnethorl. the Hurst modified steady-state method, and various unsteady-state methods such ac. those of van Ever-dingen-Hurst, Hurst, and Carter-Tracy. The Carter-Tracy water influx expression was used in this study. In certain cases, it appears that gas recovery can be increased significantly by controlling the production rate and manner of production. For this reason, the potential importance of water influx in particular gas reservoirs should he investigated early to permit adequtrtr planning lo optirtize the pay reserves. INTRODUCTION In recent years, the economic importance of natural gas production has become increasingly apparent. This has been evidenced by more intensive exploration efforts aimed at gas production, and exploitation of both deep, as well as low-permeability gas reservoirs. Technical developments such as deep-penetration fracturing have made development of such formations economically feasible. Unfortunately, water influx has forced abandonment of a number of gas reservoirs at extraordinarily high pressures. Although reservoir engineering methods for estimating water influx have long been available, it appears that application of these methods to the water-drive gas reservoir has been sporadic.'a Available methods for estimating water influx which can be applied to the water-drive gas reservoir problem include the steady-state method,1 the Hurst modified steady-state method and various unsteady-state methods such as those of van Everdingen-Hurst. Hurst, and Carter-Tracy. Interesting applications of these solu- tions to gas reservoir and the aquifer gas-storage problems have appeared recently.3,12,14 The experimental study of residual gas saturations under water drive by Geffen et al. in 1952 indicated that residual gas saturations could be extremely high. A value of 35 per cent of pore volume is often used in field practice when specific information is not available. The study of Geffen et al. showed that residual gas saturation might be much higher in some cases. Naar and Henderson concluded that the residual non-wetting phase saturation under imbibition should be about half of the initial non-wetting phase saturation. The Naar and Henderson result that residual gas saturation under water influx should be about half the original gas saturation is recommended as an estimate if laboratory measurements are not available. Thus, it is clear that a considerable portion of the initial gas in place might be trapped in a water-drive gas reservoir as residual gas at high pressure. A full water-drive would result in loss of residual gas trapped at initial reser.voir pressure. Consideration of transient aquifer behavior leads to the conclusion that high-rate production of water-drive gas reservoirs could result in improved gas recovery by reduction of the abandonment pressure. However, there appears to be little quantitative information on this possibility. One of the few advantages of water-drive gas production appears to be improved deliverability through water-drive support of the reservoir pressure. There may also be an advantage in higher condensate recovery caused by pressure maintenance for gas-condensate water-drive reservoirs. In view of the preceding, this study was made to assess the potential importance of water-drive in gas reservoir engineering. The Carter-Tracy approximate water-influx expression was used because this equation offers some advantages in hand-calculation which do not appear to have been generally recognized.' However, calculations were performed in the main with a high-speed digital computer to permit evaluation of the effect of water-drive under a large variety of conditions. CALCULATION METHOD Water-drive gas reservoir performance can be estimated in a manner completely analogous to oil reservoir calculations: a materials balance is written for the reservoir, and a water influx equation is written for the aquifer. Siniltaneous solution provides the cun~ulative water influx and the reservoir pressure. When reservoir performance data (gas produced and reservoir pressures) are available, it is usually possible to match performance data to determine the initial gas in place and the water influx parameters —

Jan 1, 1966

By T. D. Mueller, H. Dykstra

This paper describes a method for obtaining compositions of gas and oil phases in equilibrium with each other at a given reservoir temperature and pressure as lean or enriched gas is injected into a reservoir oil not in equilibrium with the gas. The method is based on calculation of equilibrium constants, K, from the equation log Kp = A + B F, where p is pressure and F is component characterization factor. A and B are constants that vary with a composition parameter. The paper explains how A and B could be determined with a minimum amount of data. The method consists of calculating equilibrium constants as a function of liquid composition and using the K's to calculate changes in compositiort of injection gas as it successively contacts original resetvoir oil at the gas-oil front and the change in composition of reservoir oil as it is successively contacted by the injection gas in the vicinity of an injection well. An iterative process is required. Two examples of the method are given which show how oil and gas compositions and properties are calculated with the aid of a computer. INTRODUCTION The calculation of recovery from lean-gas or enriched-gas injection into a saturated or under-saturated reservoir oil requires a knowledge of how phase properties change in a reservoir between injection well and producing well. To determine the phase compositions or properties experimentally would be difficult and time consuming, especially if one wanted to look at the effect of various injection gas compositions. For some fields an engineer would want to make preliminary calculations on the economics of gas injection before he makes the decision to spend a corrsiderable sum of money for laboratory work to obtain the necessary fluid properties. A method for calculatine compositions and properties based on a minimum of laboratory data is described. The minimum amount of data required is injection gas and reservoir oil composition and injection gas solubility as a function of pressure for an undersaturated oil, or the compositions of mixtures of reservoir fluid with two or more enriched gases, which mixtures have a bubble-point pressure that is the same as the expected average displacement pressure. On the other hand, if the economics appear to be attractive, more laboratory data can be obtained in the form of equilibrium constants for three or more mixtures of the reservoir fluid and injection gases. The method is based on the calculation of K's from the equation, log Kp = A + B F, where K is equilibrium constant, p is pressure and F is component characterization factor. This method of plotting K data by means of the above equation is described by Hoffman, Crump and Hocott.1 A and B are constants that vary with a composition parameter. It is shown how experimental data are used to determine how A and B vary with composition. BASIS OF METHOD EQUILIBRIUM CONSTANTS Equilibrium constant (K) data are required for this method. The K of each component is obtained from an equation that gives the product of K and pressure as a function of component characterization factor. The characterization factor F = b(l/Tb — l/T), where b is a constant characteristic of a given component, Tb is the boiling point in OR at atmospheric pressure of the given component, and T is the temperature of the system in OR. This method of plotting K data is described by Hoffman, Crump and Hocott.l The straight-line equation of the log Kp - F plot is log hp = log Kop +BFi, where Kop is at F =0, i designates the component, and B is the slope of the line. The factor F can be

Jan 1, 1966

By J. Lohrenz, C. R. Clark, B. G. Bray

Procedures to calculate the viscosities of in situ reservoir gases and liquids from their composition have been developed and evaluated. Given a composition expressed in methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide together with the molecular weight and specific gravity of the heptanes-plus fraction, the procedures are capable of calculating the viscosity of the gas or liquid at the desired temperature and pressure. The procedure for reservoir liquids was developed using the residual viscosity concept and the theory of corresponding states, and was evaluated by comparing experimental and calculated results for 260 different reservoir oils ranging from black to highly volatile. The average absolute deviation was 16 per cent. This is the first known procedure for calculating the viscosity of reservoir liquids from their compositions as normally available, i.e., including the heptanes-plus fraction. The procedure for reservoir gases uses a sequence of previously published correlations. Evaluation of the procedure was accomplished by comparison of 300 calculated and experimental viscosities for high-pressure gar mixtures in the literature. The average absolute deviation was 4 per cent. The calculations are useful for (I) determining viscosities in compositional material balance computations and (2) predicting the viscosity decrease which occurs when gases, LPG, or carbon dioxide dissolve in reservoir oils. INTRODUCTION Methods to predict viscosities of reservoir fluids from the normally available field-measured variables have been presented. Beal,1 Standing,2 and Chew and Connally3 orrelated oil viscosities with temperature, pressure, oil gravity and gas-oil ratio. Carr, Kobayashi, and Burrows4 and Katz et al.5 have presented correlations for reservoir gas viscosities as a function of temperature, pressure and gas gravity. Lie all intensive physical properties, viscosity is completely described by the following function: Eq. 1 simply states that viscosity is a function of pressure, temperature and composition. These previous correlations1- hay be viewed as modifications of Eq. 1, wherein one assumes more simple functions may be used. The assumptions are practical, because the composition is frequently not known. Further, the assumptions are sufficiently valid so that these correlations are frequently used for reservoir engineering computations. In compositional material balance"' computations, the compositions of the reservoir gases and oils are known. The calculation of the viscosities of these fluids using this composition information is required for a true and complete compositional material balance. For reservoir gases, Carr, Kobayashi and Burrows4 have presented a suitable compositional correlation. For reservoir oils, no correlation is available, and data from reservoir fluid analyses have been used7-9 for compositional material balance calculations.* From a theoretical point of view, this is entirely invalid. The reservoir fluid analysis, whether flash, differential, or other process, does not duplicate the compositions which occur during the actual reservoir depletion process, therefore the viscosities measured during reservoir fluid analysis are not those which occur in the reservoir. From a practical point of view, the "error" of using viscosities from reservoir fluid analysis is of varying and unknown significance. One can say qualitatively that the error is greatest where compositional effects are greatest, i.e., for volatile oil and gas condensate reservoirs and pressure maintenance operations. The first requirement to obtain a quantitative estimate of the significance of the error is to develop a reliable compositional correlation for the viscmities of reservoir oils. No such correlation has been available. Consistent with this requirement, the objective of this study was to develop a procedure to predict the viscosity of reservoir fluids from their compositions. Normally, the compositions of reservoir fluids are available expressed as mole fractions of hydrogen sulfide, nitrogen, carbon dioxide and the hydrocarbons methane through the heptane-plus fraction, with the average molecular weight and specific gravity of the latter. The final correlation was to use the composition in this form. While the more challenging objective of the study was the development of a correlation for the viscosities of reservoir oils, the viscosities of reservoir gases were also studied. The end result of the study was a procedure to calculate the viscosities of reservoir gases and liquids suitable

Jan 1, 1965

By D. R. Parrish, K. W. Beaver, H. W. Wood, R. W. Rausch

A pilot test of forward combustion in the Shannon pool, Salt Creek field, Wyo., is described. The Shannon sand, 950-ft deep, contains a heavy (25" API), viscous (76 cp) oil. Natural reservoir energy is limited. Primary production, intermittent since 1889, recovered only about 2 per cent of the oil in place. The field is operated by Pan American Petroleum Corp. for the Midwest Oil Corp., the owner. The original pilot was a 1.32-acre five-spot. The expanded pilot has eight producing wells surrounding a roughly triangular area of about five acres with the injection well near the center. A control or comparison well was also recompleted in another part of the field. Operation of the pilot has been little different front an ordinary gas drive. Little special equipment was found to be absolutely necessary. Except for some use of a temperature-resistant cement, all wells were conventionally completed. In spite of poor oxygen consumption, the over-all performance of the pilot has been good. Total oil recovery to June 1, 1961, was 73,971 bbl. The wells of the original pilot alone had produced about 24,000 bbl, equivalent to 50 per cent of the oil in place, when fire breakthrough at the first well occurred. These wells have now produced oil equivalent to more than 74 per cent of the oil in place in the original pilot area and production is continuing. It appears that ultimate recovery will approach theoretical maximums before the wells must be abandoned. Performance of the pilot has been encouraging, and expansion to a fieldwide combustion operation is being investigated. INTRODUCTION The results of both laboratory investigations and field tests of underground combustion have been reported previously.' However, most of the field tests were primarily experimental. More information and experience are needed before forward-combustion operations can be engineered with confidence. The purpose of this paper is to present the results of a successful pilot test of forward combustion. These re- sults should increase confidence in forward combustion as a practical method for commercial oil recovery. HISTORY OF THE SHANNON POOL The Shannon pool is located on the north end of the Salt Creek field in Natrona County, Wyo. It is approximately 50-miles north of the city of Casper. The Shannon pool's3 discovery well was completed in 1889, making this one of the oldest oil fields in the Rocky Mountain region. Three more wells were drilled in 1890. First production was hauled in wooden barrels by horse and wagon to the railroad in Casper. In 1894 a small refinery, the first in Wyoming, was built in Casper to process the Shannon crude. In the following years the field changed hands several times. It appears that each new owner did some development drilling as several wells were completed in each of the years 1895, 1902, 1905 and 1912, with negligible development in the intervening years. Forty-eight wells were drilled in this period, but many were later abandoned. Discovery of the more prolific Salt Creek field proper ultimately forced suspension of operations at the Shannon pool. After 1915 there was only sporadic production, mostly to supply cheap boiler fuel to drilling rigs in the Salt Creek field. But even this was discontinued in 1931. Since then the pool had been dormant until the recent operations, which are the subject of this paper. The Shannon pool is now owned by the Midwest Oil Corp. Field operations are conducted for Midwest by Pan American Petroleum Corp. THE RESERVOIR Fig. 1 is a map showing subsurface contours of the Shannon pool. The reservoir is on a nose of the Salt Creek anticline dipping to the north at about 500 ft/mile. The trap is provided by a shallow fault on the updip side of the productive area. The downdip limits are bounded by water, but this water has not provided an effective source of reservoir energy. The Shannon sand outcrops at many places in the immediate vicinity, providing good surface indications of the Salt Creek anticline. At the Shannon pool the sand has been lowered by faulting and is overlain by about 900 ft of shale and other sands. The Shannon sand consists of two members. The upper, a water sand, is about 40-ft thick and is separated from the lower member by several feet of sandy shale. The oil pay occurs in the lower mem-

By J. C. Savins

The practical analysis of the hydrodynamics of the wellbore has long been a subject of interest to engineers. This paper presents a simplified solution to the problem of computing the pressure drop for the flow of drilling mud in the annulus of the wellbore. This solution is, however, an exact and rigorous solution under the assumptions which have been imposed. These assumptions are that the drilling fluid is a Bingham plastic fluid* and that the annulus is formed by two concentric, stationary, cylindrical pipes. It is further assumed that the fluid is incompressible and that its motion is isothermal and in a steady state. This problem under the same assumptions has been attacked by previous authors. Beck, Nuss and Dunn' proposed that the equation for the flow of a Bingham plastic fluid in a cylindrical pipe could be applied to an annulus if the pipe radius were replaced in the equation by the hydraulic radius. This equation, known as the Buckingham-Reiner equation' (see Appendix 1), was also used in an approximate form. Van Olphen pointed out that even for a simple or Newtonian- fluid the pipe equation (Poiseuille's law) could not be converted to the Lamb equation escriptive of flow in an annulus (see Appendix 1) by using the hydraulic radius. Van Olphen further attempted to give a solution for the annular flow of a Bingham plastic fluid by introducing approximations similar to those which have been used in the case of the Buckingham-Reiner equation. Other attempts to provide approximate or exact solutions have been made by Grodde' and by Mori and Ototakeq. The present authors some years ago in unpublished work derived the correct expressions relating the pressure drop and flow rate for this problem. It was found that the solution consisted of two simultaneous equations, one of which contained a logarithmic term. Thus, obtaining numerical results for any particular case of interest involves very tedious trial-and-error computations. Very recently Laird presented the correct derivation of the two equations which are given in full detail in Appendix 1. In order to reduce the amount of calculation time which would be involved in providing a complete tabular or graphical solution to the problem, a high-speed electronic digital computer has been utilized. For this purpose the two simultaneous equations were transformed into more compact expressions by introducing reduced variables. These expressions are given in the following theoretical section. A similar procedure in this problem has been developed by Fredrickson and Bird1" Their tabular results, however, are very incomplete in the range of practical interest for problems of wellbore hydrodynamics. We have furthermore been able to express our graphical results in terms of convenient and familiar dimensionless groups. THEORETICAL DEVELOPMENT Use of Reduced Variables In terms of reduced variables the two simultaneous equations just discussed take the following form, The reduced variables q, x, a and z are defined in terms of the various measured quantities, where Q is volumetric flow rate, AP/L is pressure gradient, Dl is OD of inner pipe, D, is ID of outer pipe, is plastic viscosity, and is yield point. Thus, we have a dimensionless volume flux, a dimensionless reciprocal pressure gradient, and the ratio of the pipe diameters Before introducing the fourth reduced variable, z, it is of interest to consider the physical significance of the parameter x. As may 'be seen from the velocity profile of Fig. 1 the Bingham plastic fluid has the interesting property that a portion of the stream flows at a uniform velocity without shearing action. This section of the stream is situated approximately in the center of the conduit and is known as the "plug flow" region. Its existence is due to the fact that the shearing stresses within the region do not exceed the yield point, which is one of the two flow properties characterizing the fluid. The parameter x then turns out to be the ratio of the

By R. D. Carnahan, G. A. Stone, R. J. Arsenault

An experimental technique has been developed for carrying out a dynamic tensile stress-strain test in which plastic strain is measured continuously throughout the microstrain region extending through the macroflow region to total deformations of 5 pet. The tests were carried out on niobium, samples having interstitial impurity levels of 160 and -800 ppm at temperatures ranging from room temperature down to 100°K. A band spectrum of activation energies was obtained from calculations based on the measured activation volumes and temperature dependence of flow stress. The inability of a single rate-controlling process to predict such a phenomenon has led to the proposal that several sequential rate-controlling dislocation mechanisms are operative in the preyield micro-strain region. These are thought to be: the motion of geometrical kinks, the formation of double kinks in edges, and finally in the macrostrain region the formation of double kinks in screws. In less pure niobium the effect of interstitial impurities is shown to be dominant in the microstrain region, suggesting that a fourth probable mechanism is overcoming interstitial barriers instead of geometrical kink motion and the formation of double kinks in edge dislocations. THE principal fundamental aim of most mechanical properties studies of crystalline solids is to determine the dislocation mechanisms responsible for the plastic behavior. The extension of such studies to the region of preyield microstrain behavior has received renewed attention in the past few years as the details of macroscopic flow behavior have become better understood. For the most part, microstrain studies have been carried out using one of two variations of a single experimental technique.1 The first method involves loading the sample at a constant extension rate to a predetermined stress level (less than the macroscopic yield stress), instantaneously unloading, and then determining the extent of plastic deformation to a sensitivity of 10"6. The microscopic stress-strain curve is then constructed from the results of a repeated series of such individual measurements, each of which involves loading to an incrementally higher stress level, by relating the accumulated plastic microstrain for each step to the maximum applied stress in that step. The second variation, also carried out by step-wise increases in stress, measures the total strain to a sensitivity of ~10"6 and relates the maximum stress to the area of closed hysteresis loops that are generated under certain conditions (namely, prestrained samples). The foregoing methods have yielded interesting information about the nature of preyield deformation and plastic response to stresses below the yield stress, but they offer certain disadvantages for studying the dynamic behavior of dislocations. In the first method it is physically impossible to specify a plastic strain rate inasmuch as the test sample is reexposed, as it is loaded to its maximum stress in a particular step, to lower stresses at which flow will take place. The strain corresponding to a given stress level on the stress-microstrain curve thus generated will then always be greater than that at the corresponding stress level of a truly dynamic test carried out at the same extension rate. The second method has been used to obtain correlations with dynamic macroscopic flow parameters by applying dislocation damping theory to the hysteresis loop areas.2 The approach has had some success but suffers in that only samples that have been prestrained a few percent can be studied. Although the stresses

Jan 1, 1968

By M. Prats, D. G. Russell

The performance of a well in a bounded, layered reservoir with interlayer crossflow has been investigated mathematically. The system studied comprises a centrally located well in a bounded cylindrical reservoir composed of two layers of contrasting physical properties. The reservoir contains a single fluid of small and constant compressibility. In the case of a well producing at constant pressure, the production rate of the crossflow system declines exponentially after the early transient behavior dies out. For wells which produce with constant rate, the pressure-time relationship is linear following the transient. For very early production time, the production rate of the two-layer crossflow system behaves as that of the two-layer stratified system (the case in which the layers are in communication at the well only). Except for early time, in most practical cases the performance of a two-layer reservoir with cross-flow can be duplicated by that of a single-layer reservoir. The equivalent single-layer reservoir must have the same pore volume and the same drainage and wellbore radii as in the two-layer case. In addition, the equivalent reservoir must have a "kh" (permeability x thickness) product equal to the sum of the "kh" products of the respective layers in the crossflow system. If the porosity and permeability distributions are known, the production and pressure behavior of a well with and without crossflow can be predicted. From a comparison of the observed pressure or production performance of a well and the theoretical behavior, the occurrence of interlayer crossflow can be detected. INTRODUCTION When flow is possible across the bedding plane from a low-permeability layer into an adjacent layer of moderate or high permeability, the earlier depletion of the more permeable layer will cause such flow to occur. If this situation exists, the oil recovery from tight layers can be appreciable; thus, the economic ramifications of interlayer crossflow are apparent. The possible importance and extent of flow between adjacent and connected layers — whether resulting from capillary, gravitational, or viscous forces — have been the subject of considerable speculation for many years. Jacquard1 presented a rigorous analysis of the crossflow due to single-fluid expansion occurring between two adjacent and connected layers when the well produces at constant rate. His presentation was entirely mathematical, however, and no numerical results from evaluation of the mathematical solutions were presented. Our objectives in the preparation and presentation of the present paper are essentially twofold— (1) to analyze the behavior of a well which produces a single compressible fluid from a bounded, layered reservoir with interlayer crossflow, and (2) to present practical numerical results and simplified formulas to facilitate engineering analysis of reservoirs with crossflow. The constant-terminal-pressure case was analyzed rigorously, whereas the constant-terminal-rate case was treated only for semisteady-state conditions. DESCRIPTION OF THE IDEALIZED TWO-LAYER RESERVOIR The reservoir to be studied is assumed to be bounded and cylindrical in shape, with a completely penetrating well at its center. The reservoir is composed of two adjacent horizontal layers between which fluid is free to flow. Each layer is assumed to be homogeneous, isotropic and completely filled with a fluid of small and constant compressibility. The total thickness of the reservoir is small enough

By J. A. Fellows, Earnshaw Cook, Howard S. Avery

HEAT-RESISTANT alloys of the higher nickel and chromium ranges have been empirically developed through the practical experience of the past two decades to a position of significant industrial importance. Few quantitative metallurgical data have been published, however, that could be related directly to engineering specifications and design. Inadequate information led to the inauguration, in 1934, of a research directed toward a better understanding of the properties and performance of these steels at temperatures between 1400° and 2000°F. Investigation was concentrated upon the 26 per cent Cr:12 per cent Ni type, which owes its extensive application to an excellent combination of economy, strength and surface stability at elevated temperatures. The requirements of the oil-refining industry have been prominent in influencing the scope and evolution of the program. Creep strength at the maximum allowable deformation ranges has been employed to evaluate the merits of these alloys. Permissible elongation rates vary from 1.0 per cent in one million hours for turbine blades to 1.o per cent in ten thousand hours for tube supports and heating-furnace equipment. Throughout this paper, the limiting creep stress (L.C.S.) for a given temperature has been defined as the stress that will produce a uniform elongation rate (stage II) of one per cent in ten thousand hours as extrapolated from rates per hour or tests of less than two thousand hours duration. Certain current specifications require a minimum L.C.S. of 1600 lb. per sq. in. at 1800°F., of 3600 lb. per sq. in. at 1400°F., and include a tensile acceptance test at room temperature, after aging for 24 hr. at 1400°F., as a measure of embrittlement in service. An elongation of 4 per cent is considered satisfactory for general use, while a minimum of 9 per cent is prescribed on occasion. Many differences of opinion concerning the properties and performance of these alloys appear to result from their sensitivity to unappreciated variations of chemical analyses which, for example, may change their L.C.S. from 3500 to 350 lb. per sq. in. at 1800°F. Within broad chemical limits, it is possible to balance the six or more important elements present to attain desirable properties with a number of different combinations. One method for adjusting the composition has been described under the sponsorship of the Alloy Casting Institute,1 in which, by maintaining a substantially austenitic* alloy, optimum creep strength and minimum tendency to embrittle after aging are expected. Unless, however, the limits are undesirably dose, it has not been possible to ensure a moderate range of strength and of ductility by any chemical specification alone. There is a definite need

Jan 1, 1942

By E. Gybbon Spilsbury

(Pittsburg meeting, March, 1910.) THE recovery of gold and silver from their ores by means of the cyanide process has been so successful in the last few years that any radical improvement would seem impossible; yet the appliance to which I wish to call attention in this paper is really a radical departure from the methods now in general use. The most modern and approved of these, known as the all-sliming method, depends for its success on the grinding of the. ore so fine that practically 90 per cent. of it will pass through a 200-mesh screen. The slimes thus produced are then agitated and aerated in tanks of various types. The main object of this treatment is to insure such a thorough admixture of the pulp and solution that every particle of the ore is surrounded by a volume of solution sufficient to insure the dissolving of the whole of the gold- and silver-content. In addition, in order to expedite the action of the cyanide, and to oxidize such elements as would, if left in their active state, become cyanicides, air is blown in under pressure. In the Pachuca tank, this air is the active and sole means of agitation. The chief objections to all the methods off agitation used heretofore are the expense of operating and maintaining the mechanical devices for keeping the pulp in suspension, and the length of time required to obtain a fairly complete extraction. of the values. These objections are inherent to any method of agitation effecting a circulation of the pulp and the solution together, by giving the whole mass a circular movement, as in the low-tank system, or a vertical circulation, as in the Pachuca-tank system. In both methods, the particles of ore are kept traveling in the same direction as the solution, and with very little difference of speed; so that, while the whole mass is in

May 1, 1910

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 Roy P. Hudson

A micropyrometer known as the De Graaf apparatus has several advantages over the gas-furnace method for determining ash fusibility. When the De Graaf apparatus is operated by a modified method of procedure, the results closely check those of the gas-furnace method. This is indicated by tests on 27 samples of coal ash, with an extreme variation in fusibility values of 654" F. However, most of the ashes fused at temperatures between 2300º and 2500º. This, presumably, is the temperature interval within which the micropyrometer gives the best results. Additional comparisons are desirable. DISCUSSION (Thomas G. Fear presiding) A. C. FIELDNER,† Washington, D. C. (written discussion‡).-—The Pittsburgh Experiment Station of the Bureau of Mines has made two investigations1 of the micropyrometer for determining the fusibility of coal ash. These studies were undertaken with a view to adopting the micropyrometer method us a substitute or alternate to the present standard gas-furnace method, in view of the obvious advantages of greater speed, smaller ash sample, smaller laboratory space requirements, less equipment and greater comfort to the operator. However, the investigation showed that the micropyrometer type of furnace as then operated in five laboratories having this equipment could not be accepted as an alternate to the standard gas furnace, which is capable of close duplication in different laboratories, as shown in a recent colmparison between the Pittsburgh Experiment Station of the United States Bureau of Mines and the Fuel Research Station of Great Britain.2 The algebraic meau

Jan 1, 1934

By Howard N. Eavenson

WHEN the A.I.M.E. was formed 75 years ago the bituminous coal industry was in its swaddling clothes, although it had been operating for more than a century and coal was being mined in every state now producing it, and its output had long passed that of anthracite. In 18/1, the value of bituminous coal was about seventy per cent of the value of all metals produced and more than one fifth of the value of all products of mines, both metallic and nonmetallic.[1a] While many of the operating mines then shipped by rail, most of the development was along waterways, although the canal era of transportation was nearly at its end. Coal mining is a heavy industry and one of the best measures of its progress is its output throughout the period. Fig. 1 shows this graphically from 1800 to 1945. From 1870 to 1910, the production-which is very nearly the consumption-almost doubled every decade, a most amazing growth. In 1918, in World War I, a peak was reached and a decline began, which ended in 1932. A gradual recovery then started, interrupted in 1938 but culminating in the record output of 1944, when, owing to World War II, a tonnage of 619,576,000 was produced. It is evident that the forces leading to the continuous increase ending in 1918 have spent themselves and that only unusual conditions will require so large a production again. The postwar conditions of World War II will probably

Jan 1, 1947