Search Documents

By E. R. Brownscombe, Francis Collins

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By D. W. Fuerstenau

Recently R. J. Charles1 showed that comminution of brittle or semi-brittle materials in batch operations is described more appropriately by a variable energy relationship than by the specific relationships proposed by Rittinger,2 Kick,3 and Bond.' That is, the value of n in the empirical differential operation which applies to crushing and grinding (Eq. 1) is variable. Independently, Holmes' developed an identical generalized equaton for size reduction. The experimental data of Charles and Holmes show that the value of n is determined by the material and the comminuting machine and is constant for a given combination of these two conditions. The general objective of the present article is to extend the Charles-Holmes concepts to continuous ball milling. The specific objective is to analyze the operation of a vibratory ball mill dry-grinding chromite in a continuous open circuit and to develop an equation for the comminution of chromite in this mill. Basic Principles: In a batch mill the length of time during which the material is ground is a direct measure of the energy expended on the material, except possibly for exceedingly fine grinding. However, in a continuous operation, the fraction of the total expended energy which is useful for grinding will vary according to feed rate and amount of material in the mill. If the amount of material within a continuously operating vibratory mill is low enough so that only part of the expended energy is used for grinding, the extent of size reduction at constant power input must depend upon the length of time the material is retained within the mill, namely, upon retention time. The following analysis pertains to a continuous mill operating under these conditions. The empirical differential equation that applies to the comminution of brittle materials is:" dE = -Cdx/x" [1] where dE is an infinitesimal energy expenditure, C a constant, dx an infinitesimal size change, x the particle size, and n a constant. For continuous grinding in an open circuit at constant power input, P, Eq. 1 must be written as (dtr)P = -C dx/x" [la] where dt, is an infinitesimal time during which a particle of size x is in the mill, and C' is a constant. A ground product never consists of particles of single size but always comprises a distribution of sizes. In comminution, the expenditure of energy causes a lateral shift in the size distribution.', ' " The general objective of any engineering study of comminution is to correlate the magnitude of this size distribution shift with the energy expended, the ultimate aim being to cause the greatest shift with the least amount of expended energy. Although a number of criteria can be used to measure this shift, Schuhmann's size modulus has been found convenient.' Mathematically, the Schuhmann size distribution can be expressed as follows:" y = l00(x/k) [2] or in differential form: dy = [F- xa-' dx ] K, where y is percentage of material finer than size x, a is a constant that determines weight distribution among particles of various sizes, and k the size modulus that denotes the theoretically maximum-sized particle in any size distribution. By considering the energy required to reduce an element of weight of material, dy, from size x, to size x, the element of time during which the material must remain in the mill at at power will be (dt..).= J'-Cax!(dy) [3] Summing up all the elements of weight from zero to 100 for a product that conforms to the Schuhmann size distribution yields C C f Cdx,\( 100«x°-'dx\

Jan 1, 1960

By C. D. Gilpin, Franklin Moeller

The application of electric power for driving mine hoists handling heavy loads at high speeds has recently been extended by the installation of what is probably the largest electrically driven hoist in this country at the Granite Mountain shaft of the North Butte Mining Co., Butte, Mont. While the weight of rock which is hoisted by this machine is not the largest now being handled by electric mine hoists, the speed of hoisting is so, and the combination of the high speed with the heavy load has resulted in a machine of considerable interest to all those actively engaged in mining. Aside from mere considcration of size, this installation deserves attention because of the careful attention given, not only to the question of economy of operation, but to the details of construction and installation of the entire equipment. The conditions to be met in this installation were as follows: Weight of rock to be raised per trip, 7 gross tons........... 15,680 1b. Weight of skip.......................................... 8,000 1b. Weight of cage.......................................... 1,800 1b. Maximum depth of hoisting.............................. 4,000 ft. Diameter of rope.......................................... 15/8 Total weight suspended on one rope from drum............. 42,000 1b. Normal hoisting speed.................................... 2,700 ft. per minute. Maximum hoisting speed.................................. 3,000 ft. per minute. Desired capacity-200 tons per hour from 4,000 ft. depth. The ore is mined in several levels and hoisting is carried on from the 800-ft. level down to the 2,700-ft., which is at present the maximum depth from which the ore is raised regularly. Development is proceeding and a depth of 2,900 ft. has already been reached. The rock is handled underground in cars holding % ton, and about 10 cars are the usual load for a skip. All of the ore is dumped into underground bins, which are of diffcrent capacities varying up to 250 tons. From the bins the ore is fed into the skips without passing into measuring pockets. At the surface the

Jan 1, 1916

By E. J. Haug, K. L. Clifford, K. L. Purdy

St. Joe’s Minerals Corp.’s Southeast Missouri Mining and Milling Division produces lead, zinc and copper sulfide concentrates at three of its operations. The Viburnum and Fletcher mills have been described in previous papers. The new Bushy Creek mill, dedicated in May of 1973, is an updated version of the Fletcher mill incorporating 8.5m3 (300 cu ft) lead and zinc

Jan 1, 1980

Jan 1, 1934

By R. W. Wittkopp, J. D. Stephens

In a device known as the electron-beam microprobe, the mining industry has lately acquired a sophisticated new analytical tool that can dig out answers to a broad range of questions spanning all phases of mining activity. Essentially, the microprobe is a specialized X-ray fluorescence unit that makes it possible to obtain chemical analytical data for volumes as small as 1 cup. With this capability, it can be used both for determining the composition of small inclusions in rocks or metallic products and for monitoring changes in composition that proceed from a reacting surface. At the Research Center of Kennecott's Metal Mining Division, this versatile probe has successfully solved problems in mineral operations ranging from exploration through refining. A few instances of how the probe has been applied at the Center are presented here. In Mineral Exploration-A Good Look Before Leaping Since the objective of mineral exploration is to locate deposits of materials which can be processed at a profit, the mode of occurrence of elements or minerals in the deposits is important.

Jan 1, 1970

By Linn M. Farish

The reported production of petroleum in the Dominion of Canada for 1930 was 1,555,199 bbl., an increase of 455,373 bbl. over 1929. Alberta accounted for nearly all the production with Ontario and New Brunswick, the other two productive provinces, contributing less than 1 per cent, of the total. The completion of about 25 producing wells in Alberta was more than sufficient to offset the decline of the older wells in that province and the small decline in Ontario. The production of about 5000 bbl. in New Brunswick represents a restricted output to meet the market demand. Table 1 shows the production by provinces for the past 10 years. The consumption of natural gas for all purposes in Canada reached a new record of practically 40,000,000,000 cu. ft. in 1930, which was a considerable increase over 1929. The estimated value of the gas output was $10,500,000, against a little over $5,000,000 for the total petroleum production. Gas wells in Alberta accounted for over 75 per cent. of the total, with Ontario second and New Brunswick a poor third. The outlook for 1931 is for very little change in production of either gas or oil in Ontario and New Brunswick but for a substantial increase in both in Alberta.

Jan 1, 1931

By Luiz C. Correa da Silva, Robert F. Mehl

A. D. Le Claire and R. S. Barnes (Atomic Energy Research Establishment, Harwell, Didcot, Berks., England)-—This much awaited paper admirably confirms that the Kirkendall effect is a true diffusion phenomenon and likely to be found in all diffusing systems. The most likely explanation put forward so far to account for the shift of interface is that the two components of the system, diffusing in opposite directions, do so with unequal rates so that there is a net loss of material from that half of the diffusion zone rich in the more rapidly diffusing component. This loss is accompanied, and at least partly compensated for, by a closing up or decrease in volume of this latter half of the diffusion zone and an expansion of the other half so as to tend to keep constant the overall volume per atom. This is manifest as a movement of the original interface. If this is a true picture of what is taking place then there are at least two important matters requiring further study."1—The quantitative relating of the observed shift with the independently measured separate diffusion coefficients of the two diffusing species."2—A detailed study of the mechanism and characteristics of the closing up process."We would like first to make a few remarks pertinent to these topics. In connection with the first of these, it is a pity that more and detailed measurements were not made on the Au-Ag system, the only one for which the separate diffusion coefficients of the two species have been measured.'" One of us (A. D. Le C.) has calculated from Johnson's results the shift to have been expected in samples D.2 and D.3. Darken4 showed that the velocity v with which the markers move is related to the separate diffusion coefficients Dau, Dag, by the eauation:"dN"v = (Dag — Dau "dx"dN"where .is the concentration gradient at the mark-ax"ers. Under certain conditions,4, 20 it can also be shown that:"Xm"v = "2t"where x, is the total marker shift after time, t:"dN". . xm = 2t (Dag — Dau) [9]"dX"Johnson measured Dag and DAu in a chemically homogeneous alloy and so his reported figures represent the true or ideal diffusion coefficients"; the same quantities in eq 9 represent the actual diffusion coefficients as would be measured in the presence of a chemical concentration gradient, so we need to multiply Johnson's reported values of Dag and Dau by the"1 before insert-dlogN "ing them in eq 9. 7 is the activity coefficient of Ag at the concentration prevailing near the markers which we shall assume to be - 50 pct. The values of the thermodynamic factor for N = 0.5 and the temperatures at which D.2 and D.3 were annealed were obtained from thermodynamic data quoted by Darken.' aN/ax was computed from the coefficient of inter-diffusion of Au and Ag given by Johnson on the assumption that in the samples D.2 and D.3, the concentration distance curve could be represented by the standard equation:""No is the original difference in concentration between the two halves of the couple, and + is the error or probability function."aN/ax at x = 0 is then given by:"adN/dx "2/nDt"and we assume that this represents the concentration gradient at the markers, although it is probable that the figure is slightly high."When all the appropriate quantities are inserted in eq 9, we find that the expected shifts are: for D.2, 0.0134 cm and for D.3, 0.0119 cm."The agreement between these results and those

Jan 1, 1952

By Fred C. Bond

Many grindability tests have been made in the Allis-Chalmers laboratory since the tabulated results were last published.' The lists of standard ball-mill and standard rod-mill grindability tests are brought up to date in this paper, and considerable new material is added. This includes a tabulation of open-circuit grinding tests, with impact and compression-crushing tests. The standard ball-mill grindability tests (Table I) were made in a standard laboratory ball mill, 12 by 12 in. inside, running at 70 rpm, with a charge of 285 iron balls ranging from 1½ to 3/4 in. in diameter and weighing 20,125 grams. Seven hundred cubic centimeters of minus 6-mesh stage-crushed dry feed was used, and the circulating load was maintained constant at 250 pct by adjusting the number of mill revolutions for each grinding period. Tests were made at all mesh sizes from 28 to 200, and the grindability is the number of net grams of screen undersize produced per revolution. The table includes 479 tests. It has been calculated that the ball mill does 52.3 joules of useful work on the ore per revolution2 with 93 joules input to the mill. The standard rod-mill grindability tests (Table 2) were made in a laboratory tilting rod mill, 12 in. in diameter by 24 in. long inside, with a wave-type lining, running at 46 rpm. The grinding charge consisted of six I1/4-in. diameter and two 13/4-in. diameter steel rods, all 21 in. long and weighing 33,380 grams. The mill was rotated level for eight revolutions, then tilted up 5' for one revolution, down 5° for one revolution, and returned to level for eight revolutions. It was discharged by tilting downward at 45° for 30 revolutions. Feed consisted of 1250 cc of minus in. stone and the circulating load was maintained constant at 100 pct by adjusting the number of mill revolutions for each grinding period. Tests were made at all mesh sizes from 3 to 65, and the rod-mill grindability is the number of net grams of screen undersize produced per revolution. The table includes 105 tests. The impact and compression-crushing tests (Table 3) are given in a table that is an extension of a previous publication8 in which the method of testing was described in detail. Two equal hammers weighing about 30 lb each were arranged to strike simultaneous blows on opposite sides of the test specimen. Broken pieces of rock passing a 3-in. square and retained on a 2-in. square were set up with the smallest dimension between the hammers, and this distance was measured. The height of the hammers was increased slightly before each blow until the specimen broke, and the foot-pounds per inch of thickness was calculated. Ten pieces

Jan 1, 1949

By A. J. Phillips, C. G. Grant, Wm. B. Price

Admiralty nickel is a new corrosion-resisting and heat-resisting white metal alloy composed of 70 per cent. copper, 29 per cent. nickel and 1 per cent. tin. It has been given the trade name "Adnic." In 1924, Wm. B. Price' described this alloy and its successful use as a diaphragm metal. The physical properties of Adnic at ordinary and elevated temperatures were given in Mining & Metallury.~ A patent3 was granted on Adnic, April 28, 1925. Work was first started on this alloy during December, 1921. An investigation of the patent situation in this country and abroad, and also a search of the literature, did not reveal any work on the ternary diagram, copper-nickel-tin; and the only wrought copper-nickel-tin alloy found was that known as "Newloy," composed of copper, 64 per cent.; nickel, 35 per cent.; tin, 1 per cent. This was being manufactured and sold by Wm. Gallimorr & Sons, Ltd., Sheffield, England, for unplated spoon and fork stock, and first came to our attention in January, 1922. The preliminary work in connection with the development of Adnic was started by casting four binary copper-nickel alloys of the compositions shown in Table 1. Various a~ilounts of tin were added to each of these mixtures, to produce the series of ternary copper-nickel-tin alloys also given in Table 1. A 200-gm. sample of each of these was cast by iiielting in a clay graphite crucible and allowing the melt to slowly cool with the furnace. A longitudinal section from each of thcsc castings was polished and etched for examination. The castings were not homogenrous solid solutions, but showed the usual cored structure associated with cast alloys. Howcvrr, frorn an exarrlination of these sections we roughly determined

By Richard H. Olson, Lawrence A. Roe

Talc, when it can be isolated as a pure mineral, has a composition of 63.36% Si02, 31.89% MgO, and 4.75% H20. However, as an industrial commodity, talc rarely approaches theoretical purity. Nevertheless, the impure talc products of commerce find a multitude of uses and few substitutes in many industrial applications. Minerals commonly associated with and sold together in talcose mixtures are tremolite, chlorite, dolomite, mica, and magnesite. Steatite was originally a mineralogical name applied to pure talc. In today's commerce steatite generally refers to the massive variety of talc suitable for electrical insulator manufacture. Block steatite ore can be machined into various shapes. Impure varieties of massive or block talc are still commonly referred to as soapstone. Soft massive talc, suitable for crayon manufacture, has been referred to as French chalk. It has been common practice to discuss talc, soapstone, and pyrophyllite under the same general heading. In the case of soapstone and talc, this is a natural thing since many different types of platy, soft minerals exhibiting a high degree of lubricity (commonly referred to as "slip") have been grouped together and called soapstone or talc. Further, there is a mineralogical relationship of pyrophyllite to montmorillonite and talc to hectorite. Thus, talc and pyrophyllite are sometimes referred to as clay minerals. When finely divided talc or pyrophyllite are combined with water in proportions to make a slurry, the end product does have the appearance of a clay-water mixture. Soapstone for utensils and ornaments was mined by prehistoric Indians on Santa Catalina Island, CA. In the mid-1800s, soapstone from deposits along the western foothills of the Sierra Nevada mountains was used by white settlers for building and ornamental stone and in the linings and foundations of furnaces (Anon. 1956). Previous to 1916, the annual recorded production in California did not exceed 2000 tons. In the period 1912 to 1918, output rose sharply when the Talc City, Western, and Silver Lake mines were put into operation. From 1916 to 1935, the state's annual output of talc was in the range of 9,000 to 20,000 tons. In the mid-1930s, the use of talc in wall tile grew rapidly. In California, the production of talc grew to 63,000 stpy in 1943. The post-war building boom helped California production to grow to 120,000 tons in 1951. In 1968, production was 165,000 tons, 1969 was 145,000 tons, 1970 was 185,000 tons, and 1972 was 155,000 tons. Talc mining in New York state dates to about 1878 when a Colonel Palmer and associates opened the first commercial talc mine on the Nelson Freeman farm near Talcville, New York. In 1893, this operation was sold to the International Pulp Co. which changed its name to the International Talc Co. in 1944. The company was acquired by R.T. Vanderbilt Co., Inc. in 1974. A new talc operation, the Gouverneur Talc Co. owned by the R.T. Vanderbilt Co., began operations near Balmat, NY, in 1948. The initial capacity of the processing plant was 200 stpd (Gillingham, 1950). Subsequent expansions have increased capacity to over 600 stpd. The Vermont talc industry began just after the start of the 20th century. Talc was discovered in the area of the present Johnson, VT, talc mine in 1902. The American Minerals Co. initiated plant operations at Johnson in 1904. The Magnesia Talc Co. opened a plant at Waterbury, VT, in 1913 (Burmeister, 1963; Trauffer, 1964) and acquired American Minerals Co. in 1923. The Waterbury mine has a long history as a producer of talc crayons. The Eastern Magnesia Talc Co. was formed in 1924 by a merger of the Eastern Talc, Co.

Jan 1, 1983

By K. V. Hoover, D. M. Lapham

In Pennsylvania, 310 sampled sites of clay and shale have been studied with regard to potential use, lithology, physical properties, fired properties, major-element chemical analyses, and quantitative five-component X-ray mineralogy. Approximately 80% of the samples are of some potential use. Binary correlations of test data vs. potential use indicate that certain empirically determined parameters may be critical. These are: lithology, particularly for medium to superduty refractories and lightweight aggregate; kaolinite content for refractories vs. low-temperature fired product sources; quartz analyses with an apparent necessary minimum quantity for brick; quartz content, which is low in high grades of lightweight aggregate; mica content, which is low for refractories and possibly high for lightweight aggregate; total iron, which is low for refractories and shows a 4.0% minimum for lightweight aggregate; carbon dioxide, which is low for refractories, beneficial in lightweight aggregate, and probably restrictive for some tiles; chlorite-vermicalite, which is low for refractories and possibly beneficial in lightweight aggregate. In Pennsylvania, clay and shale materials constitute one of the state's major economic mineral resources. To date, the testing program indicates considerable development potential in brick, tile, and lightweight aggregate industries and perhaps a future potential for aluminum from shales that exceed 20% AID..

Jan 1, 1972

By Jack R. Elenbaas, Donald L. Katz

A radial turbulent flow formula has been developed which permits the computation of the pressure drop for radial flow in gas wells whether the flow is laminar, turbulent, or partially laminar and partially turbulent. Using the formula a complete back-pressure curve has been calculated and analyzed by a comparison with existing back-pressure curves. A procedure is presented for computing the permeability of the porous media when the porosity, sphericity and average particle diameter are known. Introduction Recent studies on the flow of fluids through porous media have provided new methods for computing flow under turbulent conditions. The present study was initiated as an investigation of pressure -drop computations for flow through porous sands in gas wells and as an analysis of present-day back-pressure tests. Although laminar flow exists in the producing formation of gas wells under normal flow rates, turbulent flow does take place adjacent to the well bore. As the flow rate is increased, turbulent flow exists further and further into the producing formation and under open-flow conditions in relatively deep wells a considerable portion of the flow through the sand will be turbulent. In order to calculate a complete backpressure curve for relatively deep wells without resorting to extrapolation of curves drawn from data obtained in back-pressure tests or calculations made in the completely laminar region, it was necessary to have some means of calculating pressure drop for flow through porous sands under turbulent conditions.4.6 As a result, a radial turbulent flow formula was derived from the equations recently presented by Brownell and Katz.1 A rigorous solution of this equation involves a trial and error graphical integration. By making the assumptions of average viscosity, temperature and compressibility factor of the flowing gas, it was possible to graphically integrate the relation for the general case of radial flow of gases through the porous sand of gas wells. A chart has been prepared relating the value of the integral to the radius of the well bore and to the Reynolds number of the flowing gas. The chart actually performs the graphical integration and the solution of the formula is then reduced to a simple trial and error calculation. A typical back-pressure curve is calculated to show the transition from laminar to turbulent flow. This paper will briefly describe the pro-cedurcs for computing flow through porous media, including a method for predicting sand permeability from porosity and grain size and shape. Flow through Porous Media Brownell and Katzl recently reported a correlation for computing flow of fluids through porous media by the use of an enlarged Reynolds number and friction factor in which the porosity of the bed is included as an additional prime variable. The enlargcd Reynolds number com-

Jan 1, 1948

By Donald L. Katz, Jack R. Elenbaas

A radial turbulent flow formula has been developed which permits the computation of the pressure drop for radial flow in gas wells whether the flow is laminar, turbulent, or partially laminar and partially turbulent. Using the formula a complete back-pressure curve has been calculated and analyzed by a comparison with existing back-pressure curves. A procedure is presented for computing the permeability of the porous media when the porosity, sphericity and average particle diameter are known. Introduction Recent studies on the flow of fluids through porous media have provided new methods for computing flow under turbulent conditions. The present study was initiated as an investigation of pressure -drop computations for flow through porous sands in gas wells and as an analysis of present-day back-pressure tests. Although laminar flow exists in the producing formation of gas wells under normal flow rates, turbulent flow does take place adjacent to the well bore. As the flow rate is increased, turbulent flow exists further and further into the producing formation and under open-flow conditions in relatively deep wells a considerable portion of the flow through the sand will be turbulent. In order to calculate a complete backpressure curve for relatively deep wells without resorting to extrapolation of curves drawn from data obtained in back-pressure tests or calculations made in the completely laminar region, it was necessary to have some means of calculating pressure drop for flow through porous sands under turbulent conditions.4.6 As a result, a radial turbulent flow formula was derived from the equations recently presented by Brownell and Katz.1 A rigorous solution of this equation involves a trial and error graphical integration. By making the assumptions of average viscosity, temperature and compressibility factor of the flowing gas, it was possible to graphically integrate the relation for the general case of radial flow of gases through the porous sand of gas wells. A chart has been prepared relating the value of the integral to the radius of the well bore and to the Reynolds number of the flowing gas. The chart actually performs the graphical integration and the solution of the formula is then reduced to a simple trial and error calculation. A typical back-pressure curve is calculated to show the transition from laminar to turbulent flow. This paper will briefly describe the pro-cedurcs for computing flow through porous media, including a method for predicting sand permeability from porosity and grain size and shape. Flow through Porous Media Brownell and Katzl recently reported a correlation for computing flow of fluids through porous media by the use of an enlarged Reynolds number and friction factor in which the porosity of the bed is included as an additional prime variable. The enlargcd Reynolds number com-

Jan 1, 1948

By H. P. Greenwald

THE rapidly increasing use in coal mines of portable auxiliary fans, which are generally "blowers" employed in connection with canvas tubing, raises questions concerning the hazards of such equipment, especially where the practices are opposed to established methods of ventilating face workings. Half a dozen or more coal-mine explosions and some fires have been attributed to the use of auxiliary fans; moreover, competent mining engineers have found many dangerous installations and methods of use-for example, intermittent running of the fans in gassy mines with consequent accumulations of firedamp at the faces. These facts have led the U. S. Bureau of Mines to issue this recommendation: "In the interest of safety, the Bureau of Mines, Department of Commerce, recommends that auxiliary fans or blowers should not be used in coal mines as a substitute for methods of regular and continuous coursing of the air to every face in the mine." The comparatively recent introduction of new types of mining machines, mechanical loaders, face conveyors, and slabbing machines has encouraged the trial of new systems of face mining, all of which have for their object the rapid advance of the faces. Although the latter has certain economic advantages, in gassy mines the rate at which gas is given off is generally found to increase in proportion to the speed of advance into new territory, which in turn increases the difficulties of efficient ventilation of the face workings. Some of the systems of advance mining being tried in connection with new mechanizing methods practically compel the use of blower fans and tubing, as single narrow places are now driven, 300 ft. or more without crosscuts. In some instances a series of parallel headings are driven without crosscuts pre-

Jan 1, 1928

By J. R. Kyte, R. L. Dalton, V. O. Naumann, R. L. Ridings, R. W. Greene

This paper presents a one-dimensional numerical method for calculating dissolved-gas-drive behavior for a single well or a linear reservoir. Pressure and saturation distributions can be calculated, including the effects of capillary and gravity forces. The analysis has been programmed for solution on the IBM 7090 digital computer. A comparison of experimental and calculated behavior for two linear models demonstrates that the method is suitable for calculating dissolved-gas-drive behavior under reservoir conditions of flow. Dissolved-gas-drive tests were run, at different production rates and with oils of different viscosities, on short (6.ft) and long (173ft) flow models. For tests approximating reservoir pressure gradients and flow rates, there was excellent agreement between experimental and calculated results. A series of calculations of well behavior for thin, horizontal, homogeneous, dissolved-gas-drive reservoirs shows that ultimate recovery is independent of production rate or well spacing in such cases. These calculations also show that Muskat's method and similar material-balance analyses are not applicable under some conditions particularly for high pressure drawdowns. INTRODUCTION The problem of calculating dissolved-gas-drive reservoir performance has received a great deal of attention throughout the oil-producing indusfry for many years. Really significant advances, however, have been made only in recent years. The rapid evolution of high-speed computing equipment and techniques has opened up new approaches to the problem and led to realistic, quantitative analyses free of many of the assumptions and limitations of earlier methods. The one-dimensional numerical method described in this paper includes the effects of both gravity and capillary forces. A number of test calculations, based on different rates, spacings, and rock and fluid properties, have shown that the method is adequate for systems that can be reasonably represented by a one-dimensional model (linear or radial) and in which fluid segregation perpendicular to the direction of flow is not important. Along with the development of mathematical methods, an experimental program was conducted and the data were used to test the validity of the theoretical basis of the mathematical method. Linear models were used and all experiments were run using 1,000-psig bubble-point oils. The use of long, high-pressure models enabled us to more closely approach reservoir gradients and flow rates. MATHEMATICAL FORMULATION DIFFERENTIAL EQUATIONS The differential equations describing dissolved-gas-drive behavior in one dimension dare summarized here. They include Darcy's law for flow of oil and gas (Eqs. 1 and 2) and the equations for conservation of mass for oil and gas (Eqs. 3 and 4). The oil and gas pressures are related by capillary pressure (Eq. 51.

By T. A. Rickard

A knowledge of the principles of punctuation is essential to effective and intelligible writing, for the ease and pleasure of the reader, and even his understanding, may depend upon the choice and the placing of punctuation marks. To 'punctuate' is to mark with stops or points; the word comes from the medieval Latin punctuare (punctuatus), which in turn is derived from the Latin punctum, a point. The purpose of punctuation is to show the manner in which the writing is to be understood; to indicate the relations between parts of a sentence and between successive sentences. To some extent "it does for the eye what vocal stress does for the ear", but especially it marks logical relations. Some punctuation marks, such as those indicating admiration, exclamation, and interrogation, are rhetorical as well as grammatical, but most of the common ones are logical, or grammatical. My discussion will be confined to grammatical punctuation in ordinary writing. .'The problem of punctuation in text matter", as George Summey * says, "is to employ words, points, and paragraph breaks in such a way as to achieve at the same time clearness; proper distribution of emphasis, and the desired kind of movement." The criticism of the punctuation of a sentence or of a group of sentences is not final until the context, within the paragraph, is read. A paragraph is a sentence or a group of sentences forming an independent unit of the composition; it represents a unit of thought. A multiplicity of paragraphs produces a choppy effect; 'on the other hand, paragraphs of. excessive length are tiresome, because the belaboring of one unit of thought wearies the reader. The aim should be to

Jan 1, 1931

By S. Isserow

RECENTLY, a description of the wartime work in this laboratory on the U-Si phase diagram was published. This diagram was available earlier in the open literature; as were Zachariasen's crystal structure identifications for the compounds in the system. The uranium end of the phase diagram is shown in Fig. 1. The most recent work on this system has emphasized the peritectoid E phase, of particular interest for a number of reasons. The early work provided evidence of this phase's resistance to corrosion in water and of its slight ductility. Further information on these properties has been obtained here. Work elsewhere4 has provided evidence of stability under irradiation. This body-centered-tetragonal phase is relatively isotropic compared to the more usually applied a or ? uranium phases. Unlike the other compounds in this system, it contains no Si-Si covalent bonds.3 In addition, consideration on the basis of Pauling's concepts leads to the conclusion that in the E phase, uranium has a higher valence than in any of the other phases. Definition of the composition of the ? phase has been of prime concern in the study of this alloy. Evidence is presented below that its silicon content is slightly higher than that of stoichiometric U8Si, the composition frequently assigned this phase on the basis of Zachariasen's identification. This work on the composition required careful attention to the method of dissolving the alloy for analysis, since apparently low silicon contents could be found under certain conditions, see Appendix. Preparation: Melting, Casting, and Heat Treatment—Almost all the samples were from castings with a diameter of 13/4 or 2 in. and about 1 ft long, weighing 10 to 15 lb. The U-Si charge was melted by heating to about 1650°C in an induction heated vacuum furnace. Beryllia crucibles were used at first, but were replaced by zirconia-washed graphite as the technique of keeping the carbon content low enough (about 500 ppm) was mastered. The alloy was bottom-poured at about 1650°C into graphite molds. A few small buttons, weighing no more than 50 g, were prepared in a small arc melting unit.' Since cast material has a preponderance of uranium and U8Si2, heat treatment is necessary to ob-tain the ? phase by the solid-solid peritectoid reaetion. Once obtained, the ? phase is stable to room temperature, obviating the need for speeial attention to the method of cooling after epsilonization, see Fig. 1. It was found desirable to epsilonize at 800°C for a week. For some purposes, in which residual uranium and/or U8Si2 are tolerable, shorter periods may be adequate. The temperature of 800°C is not critical, since essentially the same results were obtained from 785° to 825°C. Further departure from this epsilonization range is inadvisable.

Jan 1, 1958

By GUY C. RIDDELL, Donald M. Liddell

THE electronic arts today constitute the outstanding development in the field of rare metals, if not indeed in the arena of scientific progress at large. The year 1930 may become known as the year in which the electronic tube took hold of the wheels of industry; it is the year in which the tube's emancipation from the art of communication was proclaimed. The little tube which made a $500,000,000 industry out of a vacuum has now been applied to mining, manufacturing, transportation, chemistry, and other industries. Not many appreciate the extensive utility that e1ectroni.t tubes have already found in transportation, where they control the speed and security of express trains, as well as operate traffic lights on important highways-in the case of the Detroit-Windsor bridge counting the traffic and instantly turning the reports into the superintendent's office some distance away. In an Ohio steel mill, ten-ton ingnts 'rash back and forth through the rolls, reversed by their own shadows. The vacuum tube also starts, stops, and levels, fast passenger elevators.

Jan 1, 1931

By Herman Dykstra, Irving Fatt

Relative wetting phase permeabilities calculated from capillary pressure-saturation data are compared with measured relative permeability data. The equation relating relative permeability to capillary pressure-saturation data is derived by assuming that a porous medium is analogous to a bundle of capillary tubes. The equation includes a term to correct for the difference between fluid path length and length of core. Relative oil permeability data in the water-oil-gas system are presented, and it is shown that if the water is considered part of he rock matrix, the relative oil permeability curve is typical of a wetting phase relative permeability curve. Apparatus for measuring two and three-phase relative permeabilities are described. INTRODUCTION The difficulty in measuring relative permeability of cores has made it desirable to have a correlation between relative permeability and some more easily measured property of porous media. Such an easily measured property is the capillary preisure-saturation relation. In the past there has been a tendency to separate the capillary pressure concept from the very complex pore geometry and to consider that the capillary pressure-saturation curve give? only some characteristic distribution of interfacial curvature between two fluids.' Recently, however, it has been realized that if capillary pressure data are to yield information concerning pore size distribution2 and fluid flow in porous media. the pore must be assumed to have a simple shape such as a cylinder or a sphere. Childs and George' showed that relative water permeability could be calculated from the capillary pressure-water saturation curve by assuming that the pores are cylinders in which the fluid flowing obeys Poiseuille's Law, and that the capillary preisure curve indicated the size and number of pores. Purcell5 derived an equation relating the permeability of a porous medium to the capillary pressure curve by assuming the porous medium to he analogous to a bundle of capillary tubes. Gates and Lietz6 calculated relative permeability by an equa-tion derived from Purcell's equation. In this paper an equation is derived giving the relative water permeability as a function of water saturation by as-sulning that the core sample can be represented by a bundle of capillary tubes in which the fluid path length is not the same as the bulk length, and in which the fluid path length varies with saturation. THEORY For a bundle of N capillary tubes the flow through dN tubes will be dQ = qav dN ..........(1) where Q is the total flow rate through all the tubes and q., is the aberage flow rate through the tubes in the interval dN. If the interval dN is made small the average flow rate through a tube in the interval dN is given by Poiseuille's Law, 8 nl where ?P is the pressure drop across the tube of length 1 and radius r, and µ is the viscosity of the fluid. Substituting for qav in Equation (1), dQ =pr4 ? P/8µl dN.........(3) Darcy's Law for linear flow of an incompressible fluid in a porous medium is Q = KA?P/µL...........(4) where K is the permeability and ?P is the pressure drop across a bundle of tubes which has replaced the porous medium of length L and cross-sectional area A. Differentiating Darcy's Law with A, ?P, µ, and L constant, dQ = A?P/µL dK.........(5) Equating (3) and (5) gives dK------------dN.........(6) 8Al If the pores are assumed to be cylinders dV =pr2ldN..........(7) where V is the volume of flowing fluid in the pores. Substituting for dN in Equation (6) gives dK = —dV..........(8) 8AI' By definition the saturation, S, of the core is v v S =V/Vp..........(9) Vv -PAL

Jan 1, 1951