Search Documents

By R. W. Raymond

THE question, what is the best proportion among weight, fall, and speed of stamps, is one which has not yet received thorough and systematic examination. In considering the economical application of stamping machinery, we meet, at the beginning, with serious difficulties in obtaining accurate data for comparison. The weight and fall of stamps vary as the shoes and dies wear out; and this may lead to a change of speed also. Moreover, defects in engines, boilers, or machinery for the transmission of power, may occasion serious losses, which cannot fairly be charged to the arrangements of the stamps proper. Again, the capacity of stamp-mills is directly dependent, in some degree, upon the nature and extent of discharge, fineness of screens, and other peculiarities of the battery. Finally, the hardness and tenacity of the rock crushed varies so much that comparisons between different localities cannot be implicitly trusted. The safest experiments are those made in the same mill, by changing first one and then another condition of working; but this is seldom possible for such conditions as weight and lift of stamps, and only within narrow limits for their speed. We may eliminate of friction, transmission, and generation of power, in the case of stamps, by measuring the power actually developed by their fall. Thus, the weight, multiplied into the fall in feet, and the number of drops per minute, gives us exactly the number of foot-pounds exerted by each stamp. Dividing by 33,000, the number of foot-pounds per minute in one-horse power, we have the horse-power per stamp, from which the effective power of the whole mill may be obtained. Dividing the amount of rock crushed daily by the effective horse-power, gives us the daily amount per horse-power ; and this is the best measure that can be obtained for the effectiveness of the stamps. A complete discussion of the subject would require us to determine the exact influence of the discharge, etc., and the exact resistance offered by different classes of rocks, for both of which points the data are wanting. Professor J. D. Hague, in the third volume of the United States Geological Exploration of the Fortieth Parallel, gives a valuable

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 Ellsworth Daggett

The ore smelted in the Winnamuck furnace during the year 1872 consisted, for the most part, of oxidized ores from the Winnamuck mine, only sixty tons of outside ore (from the Spanish mine) having been smelted. The latter, like the principal Winnamuck ore, was oxidized or so-called carbonate ore. There was mixed with these oxidized ores three hundred to four hundred tons, or 7 to 10 per cent. of galena, some of which was mined with the oxidized ore, while a part was mined separately from the lower portion of the mine, and afterwards mixed with the ore, with a view of preventing the formation of deposits of metallic iron in the furnaces. The average assay in silver of all the ore handled was 51.46 oz. per ton, most of it existing as chloride of silver. The lead contents were 34.98 per cent,, all, or nearly all, in the form of carbonate of lead. The relative amount of silver is not at all constant, the best silver ore often being poorest in lead. The predominant gangue was silica, several determinations of which have been made on representative samples, yielding in three such samples 26, 38, and 58 per cent. silica, the latter test being ore containing but little lead. The average contents in silica are about 35 per cent., with 6 to 7 per cent. sesquioxide of iron, and small quantities of alumina and lime. Mechanically, the ore was very fine, and so thoroughly disintegrated that it presented few distinguishing characteristics, rendering sorting or separating of ore from waste difficult, and often impracticable. Experiments on a small scale have been tried, with a view to separate the silica by washing; but these were unsuccessful, as the finest slime, requiring a long time to settle in still water, contained a large amount of silver; and on a careful sizing and washing of the sands and coarser parts, the silver contents were found to be less dependent upon specific gravity

The opening session of the Institute was held at the rooms of the American Society of Civil Engineers, No. 4 East 23d Street, VicePresident R. W. Raymond in the chair. After a short introductory address by the chairman, the following names were presented by the Council for election as members and associates of the Institute. They were unanimously elected : MEMBERS. Abbott, James W.,.....Lake City, Col. Abbott, J. J., Jr.,...... Lake City, Col. Bartlett, J. C.,.....Cambridge, Mass. Billings, G. H.,.....Boston, Mass. Blake, Frank C.,.....Lafayette College, Easton, Pa. Dewey, Fred. P.,.....Lafayette College, Easton, Pa. Douglass, Samuel T.,. .. . Ann Arbor, Mich. Emmons, 9. F.,.....New York City. Engle, George U.,.....Philadelphia. Harrington, Dr. Bernard J.,. .. Montreal, Canada. Hayden, Wm. B.,.....Columbus, O. Hill, Leslie C.,.....Capelton, Quebec, Canada. Land, W. J.,......Atlanta, Ga. May, Wm. A.,.....Scranton, Pa. Mills, James E.,.....St. Louis, Mo. Moore, Charles W.,.... Montezuma, Ind. Page, William Byrd,.... Lexington, Ky. Putnam, Douglas, Jr.,.... Ashland, Ky. Richards, Henry,.....Dover, N. J. Boberts, Percival, Jr.,. .. . Philadelphia. Sackett, S. G.,......Brooklyn, N. Y. Schley, Wm. S.,.....Atlanta, Ga. Setz, Gustav,......Mine La Motte, Mo. Talhutt, John H.,.....Lexington, Ky. Torrey, Herbert G.,.... U. S. Assay Office, New York City. Wernlitz, Fred. W.,..., Cartersville, Bartow Co., Ga. Wiestling, Geo. B.,.....Mont Alto, Franklin Go., Pa.

By H. B. Pulsifer

.ABOut 1.5 varieties, or tnodifications, of the best rnagnesiurn available were prepared and subjected to etching tests, then examined for micro-structure. Of the 30-udd etching reagents that were tried, nearly half, mostly ammonium salts, etched the metal satisfactorily. The surfacing and etching of magnesium is shown to be a very simple and quick operation. The density and hardness of magnesium were determined. The most interesting new observations relate to the finding of the hexagonal etch figures, the crystal laminations, and the fact that the metal is plastic, cold, when restrained or quickly deformed. Under slowly applied pressule cold metal deforms slightly, then shears to fracture without plastic flow. The chief structural features of magnesium arc presented in two reduced photographs and 30 photomicrographs. Materials Tested There are only two sources in the United States from which new metal can be procured: the American Magnesium Corpn. and the Dow Chemical Co. No pronounced structural or property differences in the metals from these two companies were disclosed by the work of this investigation. The following materials were obtained from the manufacturers: 1. Massive crystals of distilled metal (American Magnesium Corpn.). 2. Rods of hot-extruded metal, ½-in. squares and 5/8-in. rounds (American Magnesium Corpn.). 3. Sheet magnesium, 0.005-in. thick (American Magnesium Corpn.). 4. Cast stick metal, 13/8-in. dia. (Dow Chemical Co.). 5. Magnesium-aluminum alloy, hot-rolled plate, ½-in. thick (American Magnesium Corpn.). From these materials, were prepared: 6. Sections from a furnace-cooled ingot made from clistilled crystal. 7. Cold-strained pieces from (6), (4) and (2). 8. Pieces cold-squeezed to fracture from (6), (4) and (2). 9. Hammer-struck pieces from (6), (4) and (2). 10. Steam-hammer struck pieces from (6), (4) and (2). 11. Sections from furnace-cooled ingot of the magnesium-aluminum alloy.

By J. L. Huitt, B. B. McGlothlin, D. K. Lowe

Published correlations for the principal aspects of hydraulic fracturing were combined into a digital computer program to facilitate the study of interrelated variables. The computer program includes individual relationships for fracture width during pumping, fracture area generated, propping agent embedment, flow capacities of propped fractures and transport of propping agents in horizontal fractures. The effects of more than 20 treatment and formation parameters on the predicted results of hydraulic fracturing treatments were studied. The effects of these parameters were determined for (I) fracture width during injection, (2) fracture width after the overburden comes to rest on the propping agents, assumed not to be crushed, (3) generated and propped fracture area, (4) location and concentration of propping agents in the fracture when injection ceases, (5) flow capacities of the various propped sections of the fracture and (6) expected increase in the well productivity. The effects of propping agent, formation and fracturing fluid parameters on well productivity are discussed. The parameters that were found to have the most pronounced effects on hydraulic fracturing treat~nents are injection rate, treatment volume, fracturing fluid coefficient, size and amount of propping agent, spearhead volume, well drainage radius and formation capacity. INTRODUCTION Many correlations have been published for predicting effects of various parameters that are considered in the design of hydraulic fracturing treatments. The Carter equation' can be used to predict generated fracture radius as a function of fracture width, fracturing fluid leakoff and other parameters. Fracture width can be determined by use of the Perkins and Kern correlation' in which the fracture width is related to the fracture radius, fluid injection rate and certain formation and fracturing fluid parameters. Wahl and Lowe et aL4 have reported methods of predicting the location of propping agents in fractures when pumping ceases. The former study is applicable to the case where the ratio of propping agent diameter to fracture width is less than 0.1. The latter is applicable when this ratio is greater than 0.1. These studies showed that the propping agent placement in horizontal-radial fractures depends principally on how the individual particles are transported in the fracture by the carrying fluid. Particle transport in fractures is determined by local fluid velocity in the fracture, fluid and particle properties, and the size of the particle relative to the fracture width. The distribution of propping agents, effective overburden pressure and formation rock strength control the propped fracture width6 by controlling the extent to which the propping agent particles embed into the fracture faces. From the distribution of propping agents and the propped fracture width, fracture flow capacities can be calculated or the various regions of the fracture. The flow capacities and the radial extent of these regions can be combined with reservoir information to predict the productivity increases for fractured wells. In all these studies, the effects of certain treatment and/ or reservoir parameters on one facet of fracturing can be predicted only if other facets which the parameters affect are fixed. For instance, fracture width and radius are interrelated; that is, to calculate the value of one, the value of the other must be known. Also, some parameters influence more than one aspect of fracturing. For example, prop-pant transport is a function of both fracture width and fluid viscosity. but fracture width is itself a function of fluid viscosity. Since these calculations are complex and the parameters interrelated, it is not possible to write an equation with which the over-all effects of treatment parameters can be solved explicitly. For these reasons, the correlations for determining the effects of the parameters which are most significant in hydraulic fracturing treatments have been incorporated into a digital computer program. COMPUTER PROGRAM The program, which was written for an IBM 7094 computer, can be used to predict results of most of the combinations and values for the treatment parameters that are ordinarily considered for fracturing treatments. A spearhead of fracturing fluid and a propping agent-carrying fluid with different fluid properties can be taken into account. Also, the total volumes and relative amounts of the spearhead and carrying fluids can be varied. Two different propping agents (as used in tail-in operations) and a wide range of formation properties and injection rates are considered. The computer program (Fig. 12) consists of several sets of calculations. First, the final flooded fracture radius and average fracture width at the cessation of pumping are calculated. This is done by simultaneous solution of the Perkins and Kern fracture width equation and the Carter equation for flooded fracture radius (equations used in the computer program appear in the Appendix). The next step is to determine local fluid velocity in the fracture as a function of time and radius. Since it is not possible to write this function in closed form expression, a table of velocity values is generated by the program and stored for subsequent use. The time span from the beginning of

By T. Sterry Hunt

The deposit of rock-salt which is known to exist along the eastern shore of Lake Huron, in the province of Ontario, has lately been more completely explored than before, by a boring with a diamond drill, put down by Henry Attrill, Esq., of New York City, and the results obtained are so important in every way that I make no apology for presenting them to the Institute of Mining Engineers. I may be permitted to preface an account of this remarkable exploration, and of its results, with a historical sketch of the discovery and development of this salt-region. It was in December, 1865, that a boring was begun near the town of Goderich, in the hope of finding petroleum. In this the adventurers were disappointed, but, after passing through about 800 feet of limestone, they encountered a series of variegated marls, in which, at a, depth of 964 feet from the surface, a bed of rock-salt, 30 feet in thickness, was met with in May, 1866. The boring was carried to a depth of 1010 feet, ending in hard rock, and yielded, by pumping, a very pure saturated brine when examined by me in August of the same year. In the report of the Geological Survey of Canada for 1863-66, published early in 1867, I described this salt-well, with many geological details, and gave an analysis of the brine. In the next three years a considerable number of wells was sunk in and around Goderich, and numerous trials were made in various other parts of the region. Salt was found at Kincardine, thirty miles north-northeast from Goderich, at a depth of about 900 feet, and also at Clinton, thirteen miles southeast, at 1180 feet. Records of these wells, with analyses of the brine from them, were given by myself in a subsequent report of the Geological Survey, 1866-69, published in 1870 (pp. 211-244), together with accounts of various unsuccessful borings in the neighborhood, analyses of the brines from the various wells (including one analysis by Dr. Goessmann), with many details of the salt-manufacture at Goderich, at Syracuse, New York, and at Saginaw, Michigan. The geological character of the region was there discussed at length, and it was shown that the salt here occurs in the Onondaga or Salina formation, which

By F. L. Fayers, J. W. Sheldon

It has been a matter of concern to the petroleum industry to determine what effect the capillary term has on saturation profiles, since these profiles determine ultimate economic oil recovery. In their original solution of the non-capillary two-phase flow problem, Buckley and Leverett&apos; showed that solutions of the two-phase flow equations became multiple-valued in saturation. Since it is physically unrealistic for saturation to have more than one value at a given position, Buckley and Leverett resolved this difficulty by introducing a saturation discontinuity or "shock". They evaluated the strength and position of the shock from material balance considerations. Recently Sheldon, Zondek and Cardwell&apos; have demonstrated that the method of characteristics can be used to solve this and related problems. Owing to capillarity it is physically unlikely that saturation shocks occur in an oil reservoir. Buckley and Leverett suggested that the effects of interfacial tensions should be included in the fluid flow equations, if the solutions are to give continuous single-valued saturation distributions. An analogy to this problem is encountered in the theory of wpersonic compressible flow where the introduction of the viscous term eliminates shock behavior. However, the mathematical form of the capillary term in the porous medium fluid flow equations differs from the viscosity term in the Navier-Stokes equation. This paper demons~rates that capillarity does eliminate the triple-valued behavior and indicates the magnitude of capillary and gravity effects on the saturation profile for a variety of systems. DERIVATION OF LAGRANGIAN AND EULERIAN FORMS OF THE DISPLACEMENT EQUATION The equations representing capillary, two-phase flow are given by applying Darcy&apos;s law and a material balance condition for each phase. The familiar Darcy and material balance equations are as follows. ^U*^0......<3, ^^r=°......c« S is the saturation of the water phase, and B the inclination of the porous medium. Pressures in the oil and water phases differ by the capillary pressure associated with the interfacial tensions, P.-P.=pAS)..........(5) where p, is determined experimentally. The porous medium is said to be oil-wet if P, > Po and water-wet if P, < Po. If the fluids are incompressible and the total volumetric flow rate, q, is constant, then v, f v,, = q...........(6) By eliminating P.. /»,, and v, from Eqs. 1, 2, 5 and 6 it follows that It is convenient to normalize the variables in the problem after the manner of Handy and Hadley% so that Eqs. 3 and 7 become where j,, = —, the fractional flow of water; (F, is a capillary pressure normalizing constant.) Substitution of Eq. 9 into Eq. 8 gives In this partial differential equation for an initial value

By T. Sterry Hunt

I have in a previous communication called the attention of the Institute to the geognostical relations of the crystalline iron ores belonging to the Eozoic racks of North America, at which time I noticed (1) those associated with the gneisses of the Laurentian series, (2) those of the crystalline schists of the Huronian, as seen in Northern Michigan and in many localities along the Appalachians; (3) those of Southeastern Missouri associated with petrosilex porphyries and schistose rocks, and presumed to be nearly of the same age with those of the Huronian schists; and (4) the highly titaniferous ores of the Norian or Labradorian series from Essex County, N. Y., and from various parts of Canada. It may here be remarked that, although the presence of titanium is most marked in the Norian ores, this objectionable element is sometimes found in the schistose ores of the Huronian, and to a lesser extent in certain Laurentian magnetites. The ores of each of these geological divisions have their individual peculiarities, which serve to distinguish them to such a degree that the magnetites and specular ores of the Laurentian are as dissimilar to those of the Huronian as are the crystalline silicated rocks of the two series. As it happens that these mineralogical dissimilarities correspond to diflerences in composition and in adaptability to different processes of manufacture, which, in the present condition of the iron industry, are objects of special consideration to metallurgists, it will be seen that these geognostical distinctions become of great practical importance. During the past year I was desired by the Geological Commissioners of Pennsylvania to devote a few weeks to certain geological studies in the eastern part of the State, in the course of which I was enabled to make, with regard to its crystalline iron ores, some observations to which I desire to call the attention of the Institute. I may mention, in passing, that the great South Mountain belt, from the Schuylkill to the Delaware, is composed of Laurentian rocks, in which are found the characteristic ores of the series like those of the Highlands of New Jersey and New York, and of the Adiron-

By B. Silliman

THESE ores occur in the Laurentian rocks in the town of Clifton, St. Lawrence County, New York. The Clifton Mining Company have opened these magnetites upon their estate of 23,000 acres, on the waters of the Grasse River, an affluent of the St. Lawrence, which flows for about nine milt",3 across the estate, exposing a section of the formation, and affording ample water-power at numerous points. In approaching this district from De Kalb Junction, upon the Watertown and Ogdensburg Railway, the observer passes from the zone of the specular ores, which abound to a great extent in Jefferson County, across a belt of non-metal-liferous limestones, continuing for some miles intermingled with metamorphic rocks, and gradually giving place to heavy beds of granitic and porphyritic gneiss, the weathered surfaces of which present a gray and rough appearance, from the washing out of the feldspar; leaving the quartz very prominent, and completely disguising the real character of the rock, which, when freshly broken, has a lively reddish appearance from the predominance of flesh-colored and red feldspar. Beds of fine-grained quartzite appear also intercalated in the granitic gneiss. These are the dominating rocks of the magnetite region, but the magnetite is always found, so far as the speaker&apos;s observation extended, associated with beds of calcite, carrying hornblende, black and green, biotite, and brown garnet interspersed with yellow and magnetic pyrites. The strike of these rocks is about northeast, and their dip about 40&apos; to 45&apos; southeast

By E. Topanelian

In recent years considerable progress has been made in the development and application of mathematical techniques for the solution of certain problems involving economic "strategies". Such a problem might involve, for example, the scheduling of shipments of a commodity from a number of sources to a number of destinations. The object would be to schedule the shipments in a manner so as to satisfy the destination requirements and at the same time minimize the transportation costs. The solution to such a problem is not necessarily intuitively obvious. The "obvious" solution is frequently far from optimum. If the shipments are to be made from, for example, only two sources to four destinations, the optimum schedule is readily found. However, if shipments are to be made from, say, 10 sources to several hundred destinations, even a competent and experienced scheduler may spend considerable time in finding a reasonable answer. Even then he is not sure that he has the optimum solution. Furthermore, he has no way of knowing how far from optimum the answer is. Consequently, he does not know whether he should accept this solution or seek a better one. Prior to the advent of large high-speed digital computers, little more could be done with such problems because of their great size and multiplicity of possible solutions. A problem involving 20 sources and 50 destinations would require choosing, from a very large number of possible combinations, the optimum combination of 1,000 variables. The best one could do was to utilize intuition, extrapolation from past experience, and other non-exact approaches. With a high-speed computer, however, such problems can be solved providing a reasonable computational procedure (or algorithm) can be utilized. The purpose of this paper is to describe such a procedure (linear programming) and to apply this procedure to a production scheduling problem. THE RESERVOIR PROBLEM A simple reservoir model is used throughout this analysis to illustrate the use of linear programming. Even the simplest reservoir behavior problem is nonlinear in both space and time, but if the geometry is fixed for a particular study and the time variable is quantized, the resultant system may be described by linear constraints on the variable production rates. Crude oil is available from five separate sources and is delivered to a pipeline. Sources 1, 2, 3 and 4 are ideal reservoirs and therefore are subjected to reservoir flow restrictions. Source 5 is labeled "Outside Source" and crude oil from this source is assumed to be available in unlimited quantities without considering reservoir conditions. The general problem is to determine the schedule of crude oil production from five sources which, over an eight-year period and subject to certain restrictions, will result in maximum profit. "Profit" is defined to encompass all economic factors and expenses involved in producing and selling crude oil to a pipeline facility. Such economic factors are lumped to the extent that the term represents the dollars profit from one barrel of crude oil from any of the five sources entering the pipeline system. Table 1 presents the assumed profit per barrel for oil from any of the five sources over each of the four time periods. In this chart the eight years is broken into four equal intervals of two years. The selection of four time periods of two years each and the values representing the potential profit are completely arbitrary selections for the purpose of illustrating this idealized reservoir problem. In any practical application of this work a careful economic study would have to be made in order to estimate potential unit profits for time periods as far in the future as eight years. In this particular study it is assumed that all of the potential unit profits (C) are known and fixed before the production schedule is attempted. For each reservoir arbitrary values are assigned for

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.&apos; 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&apos;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 W. B. Nelligam, J. Tittman

Refinements in radiation logging techniques during recent years have involved increasing usage of scintillation ditectors. These detectors produce voltage pulses whose heights are related to the energies of the gamma rays which initiate them. Analysis of the gamma-ray spectrum, as indicated by the pulse heights, yields information about the chemical elements composing the formations surveyed. Refined scintillation counter techniques can furnish chemical information concerning earth formations in situ, from a study of the gamma-ray spectra emitted by the formation either naturally or as a result of neutron bombardment."&apos; Accompanying the rising interest in gamma-ray scintillation spectroscopy, there has been increased activity in the development of accelerator-type neutron sources (in contrast to encapsulated chemical-mixture sources). Such neutron generators are attractive for several reasons: (1) they greatly reduce radiation hazards to personnel; (2) there is a great reduction in contamination danger if they are lost in the hole; (3) they can produce larger neutron intensities than can conveniently available encapsulated sources; and (4) they are capable of being pulsed, thus permitting new techniques in logging. In the past, both accelerator and encapsulated neutron sources have been used by others in conjunction with scintillation-detector pulse-height analysis. The results have not been too encouraging, due to the interference among different gamma-ray spectral lines and to the fact that the gamma-ray peaks were not too clearly distinguishable above the large and ill-defined background "noise".&apos;.&apos; This paper is a status report on laboratory studies of a technique using a borehole accelerator as a neutron source, which gives an improved scintillation spectrum, thus permitting more accurate chemical analyses of the formations penetrated. The Schlumberger-accelerator neutron source&apos; is presented; the origins of inelastic and of thermal-neutron, capture gamma rays are dis- cussed, and results are given for some laboratory measurements performed in borehole geometry. THE NEUTRON GENERATOR-TUBE ACCELERATOR In the attempt to develop an accelerator neutron source which would have the desirable properties outlined in the "Introduction" and which would operate satisfactorily under borehole logging conditions, a small-diameter, cylindrical, neutron generator tube has been developed. This tube utilizes the principle of accelerating deuterons (heavy hydrogen nuclei) by a high voltage so that they bombard a tritium (heavy, heavy hydrogen) target. The resulting reactions produce large numbers of neutrons of 14 Mev energy. Furthermore, the tube is permanently sealed and, thus, the use of pumping techniques in the sonde is avoided. The tube consists of a pressure control, an ion source in which the deuterons are stripped of their electrons, an accelerating gap down which the deuterons are sped by the high voltage, a secondary electron suppressor and a tritium-loaded target. Obviously, the tube requires auxiliary circuitry in a sonde for the control of the ion source, tube pressure, high voltage and other operating variables. Neutron yields, both continuous and pulsed, have been produced under simulated field conditions in the range between 1 and 10 times those conventionally used in neutron logging applications. For the experiments discussed in this paper, neutron-pulse repetition rates in the range between 500 and 5,000 pulses/sec are adequate. Furthermore, by making suitable adjustments in operating conditions, one can vary the pulse width from a minimum of several microseconds up to dc. As will be seen later, the present application to well logging does not require that one have available pulses of neutrons which are appreciably shorter than the time it takes for fast neutrons to slow down to thermal energy in water. As a consequence, much of our interest has been directed towards the operating characteristics of the tube with pulses longer than about 10 microseconds in duration. Our experience to date shows that tubes can be operated in this manner with reasonably constant, average neutron outputs.

By H. D. Outmans

In steady vertical flow, the interface of an immiscible liquid-liquid displacement is horizontal for any flow rate below the critical. In nonvertical flow, however, the shape of the interface in the steady state does depend on the flow rate, and the purpose of this paper is to calculate the unsteady interfaces during the transition of one steady state of flow to another. A knowledge of these transient interfaces is of considerable importance in reservoir engineering where the calculation of breakthrough recovery depends on the instant the interface reaches the producing wells and on the shape of the interface at that time. Although the emphasis is put on transient interlaces, which eventually approach stable equilibrium, it is shown that if the displacement exceeds a critical rate no equilibrium is possible. The interface is then unstable and viscous fingers are formed during the displacement. The critical rate and the shape of the transient and equilibrium interfaces are affected by the effective interfacial tension; but since this effective interfacial tension appears in the calculations only in combination with the inverse square of the thickness of the medium, its effect in the reservoir would appear to be negligible compared to its significance in model experiments. INTRODUCTION Stability criteria and the early growth of interfacial disturbances in a plane parallel to the boundaries of a dipping formation in which oil is displaced by an incompressible fluid were described in a previous paper.l This type of instability is significant in thin reservoirs. However, if the reservoir has appreciable thickness, then interfacial stability in vertical planes, normal to the upper and lower boundaries, also becomes important (the displacement is supposed to be parallel to these vertical planes). The difference between the two stability problems is that, in the first case, the intersections of the interface with planes parallel to the boundaries are normal to the direction of the displacement; in the second case, the intersections, this time with vertical planes, are not normal to the displacement. Instead, they are tilted at an angle which depends on the displacement rate. The tilt of steady interfaces was calculated by Dietz2 who also determined the critical rate of displacement for stability in the vertical plane by assuming that this rate would coincide with an interfacial tilt equal to the dip of the formation. The critical rate thus calculated is the same as has been found for thin reservoirs (see Eq. 1.1 of Ref. 1 and of the present paper). Dietz&apos;s calculation of the stable tilt was verified by laboratory experiments and the agreement was found to be fairly good.3 It is doubtful, however, that stable tilts actually exist in the reservoir because a change in production rate is not followed by an instantaneous adjustment of the interface to the new rate but, rather, by a transition period during which the interface changes from one equilibrium tilt to the other. The principal objective of this paper has been to describe these transient interfaces without putting any restrictions on the flow conditions or the shape of the interface, as had been done previously. The second objective was to compute the critical velocity, taking into account capillary effects, and the third was to evaluate, at least qualitatively, the shape of the front at rates above the critical, again without making the simplifying assumptions introduced by previous investigators.2,3 In the following sections two examples are given of the calculation of interfacial motion. The first describes this motion for an initially horizontal interface in a dipping layer, and the second for a vertical interface in a horizontal layer. The mathematical formulation of the problem is nonlinear in the boundary conditions, and this prohibits its solution in closed form. Instead, the solution is obtained in the form of higher-order approximations. 1 Before proceeding to a description of the mathematical model, however, we define two quantities

By P. Barnes

Several interesting and important considerations may be based upon an analysis of the cost of producing Bessemer rails, and the facts thus set forth may be much more clearly emphasized by re ducing each item of cost to its own fractional proportion of the sum total in each branch of the manufacture. The estimate given below is not absolutely complete in minute ness of detail, but it may be fairly taken as a basis for the suggestions to be noted : Perhaps the most striking comments that can be made upon this statement are those arising from a consideration of some of the smaller items, and even the least of these quite loses its insignificance when its proportion of the whole is expressed in actual money, as for the work of a year, or even for a single day. Thus, if the item as be considered, the importance can be clearly

By T. F. Witherbee

The following calculation of heat requirement covers the working of the furnace from January 25th to February 14th, inclusive. A short time previous to the first date the furnace had been working rather badly, so that 300 lbs. of ore and 500 lbs. of limestone were taken off. Within the time mentioned the quality of the iron passed

By M. R. Tek, K. H. Coats, D. L. Katz

In natural gas storage operations, seasonal pressure fluctuations in the gas reservoir cause the water from the surrounding aquifer to flow into and out of the gas sand. The theory of unsteady-state liquid flow through porous media developed by Van Everdingen and Hurst&apos; has been applied to predict the water movement into and out of the gar bubble for .several .postulated pressure cycles. Those cycles with as many pound-days above the original aquifer pressure as pound-days below, may cause the gas reservoir to slowly grow or shrink rather than hold to a constant volume. Applications to an actual field case study give the predicted gas reservoir monthly pressures and volumes as compared with the observed monthly pressures and volumes. INTRODUCTION A large number of natural gas storage reservoirs are bounded by or adjacent to large water-saturated formations called aquifers. The presence of these aquifers is usually evidenced by the production of water from wells delineating the gas-hearing sands. Some gas storage reservoirs in use today have been purposely developed on aquifers. Cyclic pressure variations in a gas storage reservoir cause water influx and efflux from the surrounding aquifer. This, in turn, results in a varying gas reservoir volume. The prediction of the effect of this aquifer fluid movement on the size and size variation of the gas reservoir can provide information valuable in the study of several reservoir engineering problems. The economics of gas storage operations are directly influenced by the influx of aquifer fluid into the reservoir since a shrinking storage reservoir requires increasing pressures for the storage of the same quantity of gas. Material balance and reserve or recovery calculations also obviously require knowledge of reservoir volume. Various other reservoir engineering calculations such as interpretation of well interference data, evaluation of physical characteristics of porous media, water coning, gas injection and pressure maintenance studies are typical problems where variations of volume due to edge or bottom-water encroachment becomes important. The solution of equations describing influx and efflux of water in the gas sand, application to gas storage reservoir behavior and interpretation of results are the main objectives of this paper. While this work was directed toward gas storage reservoirs, it applies equally well to repressuring of partially depleted oil reservoirs with an active water drive. MATHEMATICAL ANALYSIS As just mentioned, the cyclic pressure variations in an aquifer-surrounded gas storage reservoir cause water movement and varying gas reservoir volume. Prediction of this volume variation from the reservoir pressure cycles or the injection-production schedule requires adoption of a flow model, certain assumptions concerning reservoir and fluid properties, and formulation and solution of the governing differential equation. The flow model used in this study is a gas bubble located on a very large aquifer (Fig. la). Gas bubble radius and aquifer sand thickness are denoted by rr and h, respectively. Fig. lb shows the horizontal, radial two-dimensional porous reservoir attended by the unsteady, compressible motion of the aquifer liquid. Volume variation of the gas bubble is determined from aquifer water flow across the cylindrical boundary, r = r0

By N. W. Ratliff, P. J. Closmann

Production of oil by expansion from a cylindrical reservoir composed of two concentric regions of different properties has been determined as a function of time for a reservoir producing at constant terminal pressure. The parameters involved are permeability, porosity, compressibility and oil viscosity. Results agree with those of a generally accepted method for uniform reservoirs when all parameters are taken as uniform. Cumulative production at any given time is reduced below that of the uniform reservoir when a region of considerably lower permeability adjoins the well; cumulative production increases when a high-permeability region adjoins the well. Curves are presented illustrating quantitative effects of these variations. INTRODUCTION Most petroleum reservoirs can be exploited by release of pressure and consequent expansion of underground fluid. During part of the production history of certain reservoirs by this mechanism, fluid compressibility can be considered small and constant. Cumulative fluid produced by this means as a function of time was calculated for a uniform reservoir by van Everdingen and Hurst.l For certain cases where a production well has been damaged or where an acid treatment has been used, a zone of properties different from those of the reservoir may be created around the production well. In such a case, expansion of the reservoir fluid takes place in a system composed of two concentric cylindrica1 regions of different reservoir properties. Such composite systems were studied in connection with heat flow by Jaeger2 who presented a solution for temperature distribution in a radial system with an infinitely large outer radius. Similar heat flow problems have been studied by other authors Oil production from composite reservoirs was studied for constant production rate by Hurst,4 Loucks and Guerrero,5 and Carter.6 THEORY This work considers constant terminal pressure. Cumulative flow is calculated for a composite bounded system in which no flow takes place across the outer boundary (Fig. 1). Since the flow is considered to be purely radial, the problem involves only one space dimension and time. It is convenient to define the following dimen-sionless pressure drops (Fig. 1).