Search Documents

By R. B. Burt, James A. Barr, I. S. Tillotson

DISCUSSION Howard Howie (Knoxville, Term.)—The authors are to be congratulated on the presentation of a paper containing so much valuable information on the pipeline transportation of phosphate, as there is very little literature on the subject. The writer is especially interested in the paper, as he conceived the arrangement of the Akin and Godwin plants and was in charge of the design work and the engineering incident to their construction. The Akin and other phosphate deposits in the Tennessee phosphate area lie on beds of limestone that are very irregular. The limestone beds, after the phosphate matrix has been removed, are similar in appearance to land severely eroded by the action of water and denuded of top soil. Depressions in the limestone, called cutters, are irregular in depth with vertical or overhanging walls, having the general appearance of cracks in dried clay. They change abruptly in direction, width, and depth, and vary on the Akin tract from 1 to 25 ft in depth and from 1 to 50 ft in width. Pinnacles of limestone commonly occur in the cutters which appear, when exposed, like small clifflike islands in a river. Limestone floats also occur. The phosphate matrix fills the cutters and covers the uncuttered areas, the thickness of the cover varying continually and sometimes rather abruptly. It occurs generally as stratifications of phosphate rock and clay of varying thickness. The phosphate rock in the matrix varies in hardness and in percents of silica, lime, iron oxides, and fluorine, and the clay varies in toughness. In some places the deposit consists of narrow strata of rock almost devoid of clay streaks. In other nearby locations the clay will predominate. When it is excavated, the phosphate rock breaks into thin irregular lumps, locally known as plate rock. Limestone lumps are also excavated with the matrix. Akin plate rock is generally much softer than that occurring in other deposits in the area. Because of the above described physical and chemical variations of the excavated material, the resultant slurry varies in size distribution, specific gravity, and percent of slimes. When the rock is soft or when there is an increase of clay, the slime fraction is greatly increased as it passes through pumps and pipelines, resulting in reduced pipe friction. It is obvious that the longer the pipeline the greater the reduction of coarse fractions into fines, causing a decrease in pipe friction that cannot be accurately evaluated. The matrix is mined with a dragline that drops it into a hopper with a grid composed of 9-in. parallel bar spacings located above the hammer mill. The matrix on, and passing through the grid, is subject to the action of powerful sprays which wash it down to the hammer mill, together with any limestone lumps that are not removed before passing through the grid. The hammer mill reduces the feed to lumps of plate rock and clay, most of which will pass through the 8-in. pump suction. The mixture discharges into a pool containing the pump suction pipe. Water from a hydraulic nozzle moves the mixture to the pump suction intake. The pump, driven by a variable speed motor, is the same size as the pumps mentioned on p. 279. Provision is made to remove any lumps that lodge in a bend in the suction pipe in a manner similar to that used in the Florida phosphate fields. The hammer mill and pump units are mounted on wide steel skids so that they can be moved as the mining operation progresses. The discharge from the pump flows through an abrasion-resistant spiral welded steel pipe 8 1/4 in. actual inside diam, 8 in. nominal diam, for a maximum distance of 2200 ft, which is the limiting pumping distance for one pump. This pipeline, referred to hereafter as pipeline A, discharges into a ball mill without balls, which in turn discharges into a rotary screen attached to it that separates the slurried matrix into 11/4-in. oversize and undersize fractions. The oversize is returned to the mill for further reduction; the undersize is pumped to a Dorrco washer and then flows into a hydroseparator 160 ft in diam. In spite of the size reduction in the hammer mill and the blunging and washing of the slurry in its passage through the pump, pipeline, mill, and washer, the discharge to the hydroseparator frequently contains mud balls almost perfectly spherical. Sometimes the discharge from the 16,000-ft pipeline at Godwin contains mud balls the size of bird shot and smaller. This pipeline will be referred to subsequently as line B. Liquid caustic is added to the slurry at the Akin plant before its passage through the hydro-separator, which decreases the size of particles in the overflow by dispersion. In passing through pumps 1, 2, and 3 and pipeline B, the slime fraction in the underflow is increased by abrasion and blunging and also by continuing dispersive action of the caustic. The matrix for use in the experimental tests referred to on p. 279 was obtained from three small surface openings on the Akin tract that were made previous to the purchase of the tract by the Authority. Matrix used in the 2 and 4-in. experimental pipeline tests was taken from the three openings and proportioned to obtain a sufficient quantity that would be fairly representative of the average in the Akin deposits. Prospecting samples of matrix had been obtained from drill holes which showed no small variation in physical and chemical properties. Some of the physical variation is evident from the size distribution of solids in samples taken during the tests covering line B flows so thoroughly made by the Authority under the direction of Mr. Burt, see Table V, p. 280. Hydraulic gradients for a pipe of 8-in. diam were derived from the 2 and 4-in. pipeline tests using the so-called representative matrix as above described, and plotted on the profile of pipeline B. Gradients of other materials in slurry form passing through pipelines that bore some similiarity to the Akin matrix slurry were also plotted. After a study had been made of the hydraulic gradients plotted on the profile and the varying slurry flow that would probably occur during actual operation, three pumps were ordered, referred to as No. 1, 2, and 3 on p. 279. No. 1 and 2 pumps were installed and the third kept in reserve should the operation of 1 and 2 pumps prove satisfactory, since installation of the third pump would require an attendant, as well as the laying of 5400 ft of pipe to supply it with seal-water from the Akin plant. Subsequently, it was found desirable to install the third pump to maintain capacity when the slime fraction was low and the coarse fractions were large. Reference to Fig. 6 will show that the flow in B line goes upgrade in three locations. At the outlet at Godwin, the slurry flows between two 45" bends for an approximate distance of 18 ft to rise above the ground a sufficient height to discharge into a launder feeding the first classifier. This condition requires extra energy, which is taken care of by keeping the hydraulic gradient a sufficient distance above the high points. Although there are rather heavy upgrades in the line, the choking condition that might occur at the

Jan 1, 1953

By W. E. Brigham, P. W. Reed, J. N. Dew

By I. W. Alcorn

The cost of drilling in the past few years of proration and ever decreasing allowables has received increasing thought and study. It seems to parallel the strides made with respect to pumping problems in 1932 and 1933. Low profit margins stimulate the search for economies. Any discussion of slim-hole drilling involves a multitude of abstract values, and the problem can be attacked only by a study of records and rational thinking. In a measure, traditional methods must be studied with a questioning attitude. The words "slim hole" may be somewhat misleading, as the term is entirely relative. For that reason I wish to define it as used in this discussion. On the Gulf Coast, I would define a normal hole as one in which 10 3/4-in. surface casing has been set and below which a 9 7/8-in. hole is dug. With this as a datum, holes of larger diameter below surface would be called "large holes," with 121/4 in. as the "extra large" class and holes of smaller diameter classed in the slim-hole range, 6 1/4 in. representing about the minimum size. The optimum size has not been definitely determined, but the preference is for a program that would call for setting 8 5/8 -in. surface casing and digging a 7 3/4 - in. hole. Different sections of the country might very properly deviate from this practice because of differences in depth, formations, regulations, or other factors. This paper will be confined to a discussion of the drilling of exploratory wells, as the drilling of this class of wells is quite a different problem from that of drilling a proven field, although many of the same considerations apply in both cases, as will be apparent. Although for many years shallow slim-hole and core-test holes of small diameter had been dug, it was not, I believe, until the fall of 1937 that serious consideration was given to its application in drilling wells 7000 ft., 8000 ft. and 9000 ft. deep on the Gulf Coast. This does not mean that nearly all such wells drilled on the Gulf Coast now are of the slim-hole variety, but only that efforts have been aimed in that direction with very encouraging results and there has been a very definite trend to reduce the size of holes, a trend that probably will continue. As to exploration policy, there appear to be two distinct and opposite views. The one usually followed in the past was more or less designed to "shoot the works" on one well. Generally on the Gulf Coast, the depth at which production might be obtained is limited by the drilling equipment and economics. With this thought in mind, wildcat wells are often planned with the idea of carrying the depth to 10,000 It., or even 15,000 ft. All of the possible hazards are anticipated and the general drilling and casing programs are so predicated. Such wells contemplate the drilling of large-diameter holes; the setting of at least one protective string of casing and an extensive coring policy, all of which finally adds up to an expensive well. AS opposed to this general policy is the slim hole. This, on the Gulf Coast, recognizes that new discoveries represent gam-

Jan 1, 1941

By I. W. Alcorn

The cost of drilling in the past few years of proration and ever decreasing allowables has received increasing thought and study. It seems to parallel the strides made with respect to pumping problems in 1932 and 1933. Low profit margins stimulate the search for economies. Any discussion of slim-hole drilling involves a multitude of abstract values, and the problem can be attacked only by a study of records and rational thinking. In a measure, traditional methods must be studied with a questioning attitude. The words "slim hole" may be somewhat misleading, as the term is entirely relative. For that reason I wish to define it as used in this discussion. On the Gulf Coast, I would define a normal hole as one in which 10 3/4-in. surface casing has been set and below which a 9 7/8-in. hole is dug. With this as a datum, holes of larger diameter below surface would be called "large holes," with 121/4 in. as the "extra large" class and holes of smaller diameter classed in the slim-hole range, 6 1/4 in. representing about the minimum size. The optimum size has not been definitely determined, but the preference is for a program that would call for setting 8 5/8 -in. surface casing and digging a 7 3/4 - in. hole. Different sections of the country might very properly deviate from this practice because of differences in depth, formations, regulations, or other factors. This paper will be confined to a discussion of the drilling of exploratory wells, as the drilling of this class of wells is quite a different problem from that of drilling a proven field, although many of the same considerations apply in both cases, as will be apparent. Although for many years shallow slim-hole and core-test holes of small diameter had been dug, it was not, I believe, until the fall of 1937 that serious consideration was given to its application in drilling wells 7000 ft., 8000 ft. and 9000 ft. deep on the Gulf Coast. This does not mean that nearly all such wells drilled on the Gulf Coast now are of the slim-hole variety, but only that efforts have been aimed in that direction with very encouraging results and there has been a very definite trend to reduce the size of holes, a trend that probably will continue. As to exploration policy, there appear to be two distinct and opposite views. The one usually followed in the past was more or less designed to "shoot the works" on one well. Generally on the Gulf Coast, the depth at which production might be obtained is limited by the drilling equipment and economics. With this thought in mind, wildcat wells are often planned with the idea of carrying the depth to 10,000 It., or even 15,000 ft. All of the possible hazards are anticipated and the general drilling and casing programs are so predicated. Such wells contemplate the drilling of large-diameter holes; the setting of at least one protective string of casing and an extensive coring policy, all of which finally adds up to an expensive well. AS opposed to this general policy is the slim hole. This, on the Gulf Coast, recognizes that new discoveries represent gam-

Jan 1, 1941

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By G. W. Doepel, W. P. Sibley

A three-dimensional, multilayer, mathematical model has been developed for predicting performance by mis-c.ible displacement. Areal and vertical coverage, total gas-oil ratio, propane-oil ratio and oil-recovery factor for the multilayer system are computed as functions of hydrocarbon pore volumes of solvent injected. The technique has been programed for use on the IBM-7090 digital computer and was used in calculating perfortnance for continuous miscible displacement of a five-spot, 50-layer .system. Results reported are based on the use of published laboratory data for a single, uniform five-spot layer. However, the technique is sufficiently generalized to calculate performance of any system of layers of any pattern configuration and mobility ratio, provided data are available to describe perfortilance of a single layer. The technique has also been adapted for predicting oil recovery resulting from various gas-driven slug sizes. The results presented for both continuous and slug-type injection show, that permeability variation has 0 significant effect on oil recovery but that solvent composition, within the range studied, has a lesser effect on oil recovery. INTRODUCTION Experience from many of the miscible-phase displacement projects currently in operation indicates that the injected solvent flows preferentially through the more permeable zones in the reservoir. This preferential solvent flow in the high-permeability streaks results in a reduction of vertical coverage (C,) and an attendant decrease in oil recovery. The multilayer technique presented in this paper uses areal coverage (C.,) data for a single layer from the literature and from laboratory measurements of displacement efficiency (herein referred to as microscopic coverage, C,,) to calculate areal coverage (C,) and vertical coverage (C,) for the multilayer system. Oil recovery for the multilayer system is obtained as the product of (Ca) (Cm) (Cv). It is assumed that C, in a single layer is 100 per cent; therefore, (C,) (Cm) is a measure of oil recovery from that layer. Data presented in the literature for a single layer generally deal only with areal coverage and/or microscopic coverage —but seldom specifically with oil recovery. Nevertheless, this technique is flexible in that either Ca and C,, of an individual layer or a direct measurement of oil recovery for the layer may be used. In developing this technique, the problem arose concerning a method of determining the relative amounts of solvent that will enter each layer during any increment of solvent injected. To solve the problem, a dimensionless equation was derived for expressing solvent injected into each layer as a fraction of total injection. A finite-difference method was applied in integrating the equation to determine oil recovery for the system of layers as a function of pore volumes injected. The technique subsequently was programed for solution on a digital computer. The method was applied to predict performance of a deep California reservoir under continuous miscible displacement and has been extended for calculating oil recovery as a function of slug size for gas-driven slug injection. The technique may be applied to any petroleum reservoir that can be adequately represented as a layered system, by considering permeability variation,1-3 permeability range and average permeability. Each layer is then assigned its appropriate permeability, porosity and hydrocarbon saturation. RESERVOIR VOLUME CONTACTED BY SOLVENT It has been recognized for some time that high microscopic coverage (C,) can be obtained by injecting a miscible agent, such as propane, into an oil reservoir. These high values of C,,, have been demonstrated in the laboratory and are normally between 80 and 100 per cent, depending on solvent and reservoir oil composition as well as on reservoir temperature and pressure. In miscible-phase displacement, because there is no interface between solvent and oil, effective permeability to solvent is the same as that to oil. (This assumes that oil saturation has not changed and that water is immobile.) Consequently, mobility ratio (M) simply reduces to that of viscosity ratio . Perrine4 and Habermann5 have shown that large

R. H. Sweetser, Sault Ste. Marie, Ont. (communication to the Secretary*):—The recent work of furnace No. 1 of The Algoma Steel .Co., at Sault Ste. Marie, Ontario, using for fuel, first all charcoal, then a mixture of charcoal and coke, and finally all coke, upsets the estimates given by Mr. Hall in the last two paragraphs of his interesting paper. Furnace No. 1 was built for a charcoal-furnace, and it was run on charcoal fuel from March 6 until July 9, 1905; then, for a week, a mixture of charcoal and coke was used, and, since then, coke alone has been the only fuel. It is to be regretted that at no time, during the run on charcoal, was there enough charcoal on hand and " in sight" to warrant blowing the furnace as hard as it should have been blown. During April and May the supply of charcoal limited the amount of blast to about 8,000 cu. ft. of air per min.; in June, to about 8,900 cubic feet. For a few days in July about 10,000 cu. ft. were blown, and then the maximum daily output, 173 tons (2,240 lb.), was reached. At all times the quality of the charcoal was poor, having been made from wood cut from three to four years before. Doubtless this furnace could have produced 200 tons of iron per day, even with this quality of charcoal, provided the supply had been sufficient; and with

Jan 1, 1906

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By B. H. Caudle, A. B. Dyes

In a recent publication' it was shown that wells with a free surface in a homogeneous gravity-drainage reservoir have a hyperbolic decline with index n '. This paper reports efforts to extend this theory to practical field cases. As a first step in making the extension, a large number of lease production curves from gravity-drainage fields in the stripper stage were fitted with the general hyperbolic decline equation, q/q, = 1/(1 + y) y and n in each case being determined by statistical methods. It was found that this equation is reliable for decline curve extrapolation for periods up to at least 25 years. The resrilts obtained were disturbing in the sense that values for the index, n, were often smaller than one, which fact leads to the prediction that a well would produce all infinite amount of oil in infinite time. A theoretical investigation based on the results of the previous paper provides a suitable explanation for these low values of n. It suggests that these are obtained if production occurs from two or more layers of diflerent permeabilitv and thicktress or two lnyers having different skin effects. It is predicted that as the more permeable layers are exhausted, the index of decline, n, will rise from a value smaller than 2, to n = 2. A practical method is given for fitting the general hyperbolic decline Eq. 6 to production decline data. When this method is used, the index of decline observed in gravity-drainage reservoirs in the stripper stage when wing only data over early years 01 decline usually agrees quite closely with that determined using data over 25 years or more. This indicates that the hyperbolic decline equation is reliable for predicting behavior of such reservoirs. A method is also suggested for taking into account any changes in the value of the index INTRODUCTION As a result of both experimental work with models and theoretical deductions, it was recently concluded' that wells having a free surface in a homogeneous reservoir producing by pure gravity drainage should have a hyperbolic decline with n = 2. The present report details efforts to extend this work on ideal homogeneous reservoirs to actual field cases. It should be remembered that both the previous work and the present extension are based on the assumption that the reservoir is at very low pressure, and hence, nearly all the potential available for driving fluid to the well is the gravity head. Attempts to apply conclusions of this study to reservoirs which have considerable pressures supplementing the gravity drive may lead to difficulty. WELLS WITH NO FREE SURFACE In a dipping gravity-drainage reservoir such as that shown in Fig. 1, some of the down-dip wells have no free surface. The previous study' showed that the up-dip wells, which have a free surface, should follow Eq. 1. The following approximate analysis, of a type which was justified experimentally for up-dip wells in the work just mentioned, leads to results for the down-dip wells.

By Edward W. Parker

This plant was built in 1907 at the Bankhead mines to manufacture briquettes by the Zwoyer process under license from the Zwoyer Fuel Go., of New York, N. Y. The building was constructed to contain two units, each of a capacity of 10 tons per hr.; one unit was installed at the time the building was constructed, and the second unit is now being built. The coal used is an anthracite, which is more friable than that of Pennsylvania, and, as a consequepce, a larger percentage of dust or waste is produced, which was formerly thrown upon the slack-pile as waste, but is now passed to the " dust-bin " in the breaker, from which it is conveyed to the briquette-plant. An average analysis of the coal, published by Lewis Stockett and 13. R. Warden,' is as follows: Per Cent. Moisture,.........0.50 Volatile,..8.00 Fixed carbon,... 83..50 Ash,.......8.00 Total,.100.00 Sulphur,.........0.40 per cent. Specific gravity,...... 1.40 Color of ash,.white Heat value,.......14,000 B.t.u. The plant is run under two 12-hr. shifts, the force consisting of 1 superintendent for both shifts, 1 engineer for each shift, 1 briquetter for each shift, 1 helper for each shift, 1 laborer cleaning up for each shift. In the pitch-melting house, 2 men are on duty for one shift. The steam for the engine is furnished from a central boiler-plant, which also supplies steam for the breaker, the electric light and power-plant, machinery-shop, and boiler-shop. During March, April, May, and half of June, in 1907, the

Jan 1, 1909

ROBERT VAN ARSDALE NORRIS, who died at his home at Wilkes-Barre, Pa., on April 20, 1928, was born at Newark, N. J., on May 2, 1864, of parents who were direct descendants of such sterling New Eng-landers as Governor William Bradford, John Alden and Timothy Dwight. He went to the Collegiate School in New York and then to the Colum-bia School of Mines where he graduated with the class of 1885, having among class-mates A. S. Dwight, E. N. En-gelhardt, H. L. Hollis, G. B. Lee, Richard Moldenke, Wil-bur Sanders, Joseph Struth-ers, E. A. Wiltsee, and I. H. Woolson. His marked ability even among so many men of unusual ability caused him to be chosen as assistant to Pro- fessor Munroe in the summer school of mining and survey-ing immediately after gradua-tion, and he shortly after was also inspector of dredging operations on the Maurice River improvement work, Millville, N. J. The follow- ing year he worked in the chemical laboratory of Herman Behr in Brooklyn, and toward the end of 1886 he joined the staff of the Pennsylvania R. R. Coal Co. and the Susquehanna Coal Co. as assistant engineer. Six years later he was made principal assistant en- gineer and in 1900 he was made chief engineer of these companies. In 1904, he withdrew from them and estab- lished an office as a consulting engineer at Wilkes- Barre, later taking his son, Robert V., Jr., into partner- ship with him. They acted as consulting engineers for the coal department of the Delaware, Lackawanna & Western R. R. Co., the Lehigh Valley Coal Co., Coke Bros. & Go., Inc., the coal department of the Pennsyl- vania R. R. Co., the Wilkes-Barre & Hazelton R. R. Co., the Sterling Salt Co., and others, specializing in ex- amination and appraisal of coal properties, and also matters of Federal taxation.

Jan 5, 1928

By Charles Jones

IN the winter of 1902, while occupying the position of mining and examining engineer for the Mountain Copper Co., Ltd., of Keswick, Cal., I had occasion to discuss with the General Manager, among other things, the question of the enormous waste of sulphur in smelting cupriferous, pyrite ore, and to cast about for means to utilize this waste or to render the fumes innocuous to vegetation. The decreasing rate of production of California farms, the consequent need of fertilizers, and the opportunity thus offered to utilize large quantities of sulphuric acid in the manufacture of superphosphates naturally brought the subject to the practical and commercial side of phosphates. I was informed that the company had made efforts to find phosphates in the West, without success. I advised as a start placing inquiries for phosphate deposits on the Pacific coast in the current mining ,journals. This was done, with the result that in January, 1903, from an advertisement in the Mining and Scientific Press we got track of a deposit in Rich county, Utah, through T.J . Wilson, of the Southern Pacific railroad in San Francisco. As soon as snow in the mountains would permit, in May, 1903, I undertook an examination of the property, with the result of recognizing the importance of the find, in spite of its lack of facilities for transportation at a commercial figure. On June 15, 1903, I concluded my report by saying, "This may be one of the world's deposits of phosphates." On Nov. 17, 1903, I concluded my report for the season in the following words : " I have now no hesitation in saying that I believe the Carboniferous formation so largely developed in the Wasatch mountains of Utah and the mountains of the Great Basin will prove to be one of the greatest storehouses of phosphoric acid in the world."

Jan 10, 1913

By Robert Glass Cleland

DURING the closing month of 1929, Walter Douglas found his health impaired by the strain of many difficult years of alternating prosperity and depression, and in April 1930 resigned the presidency of Phelps Dodge. Douglas was succeeded by Louis S. Cates. The choice broke a tradition that had been held inviolate for almost a hundred years. For the first time in its long history, the direction of the company's policies and operations was entrusted to a stranger without previous Phelps Dodge background, training, experience, or tradition. But the choice was wisely made. For a time at least, the path of Phelps Dodge was to run through rough and dubious ground, and no man in the United States was better qualified to blaze that trail than Louis Cates. After graduating from the Massachusetts Institute of Technology, Cates had served a rough and financially un- profitable apprenticeship at one of the century-old mining camps of Chihuahua, Mexico. Following his return to the United States and a brief interlude as assistant to the president of the National Steel and Wire Company, he had resumed his mining career, first as timekeeper, then as

Jan 1, 1952

By N. C. J. Ros

The paper deals with pressure gradients occurring in flowing and gas-lift wells, a knowledge of which can be applied to the determination of optimum flow-string dimensions and to the design of gas-lift installations. The study is based on a pressure-balance equation for the pressure gradient. It appears that a pressure-gradient correlation of general validity must essentially consist of two parts-—one part being a correlation for liquid hold-up and the other part being one for wall friction. Dimensional analysis indicates that both liquid hold-up and wall friction are related to nine dimensionless groups. It is shown that in the field of interest only four groups are really important. On the basis of these four groups a restricted experimental program could be selected that nevertheless covered practically all conditions encountered in oil wells. This experimental program has been carried out in a laboratory installation. Three essentially different flow regimes were found. The pressure gradients in these regious are presented in the form of a set of correlations. Comparison of these correlations with a few available oilfield data showed excellent agreement. INTRODUCTION Prediction of the pressure drop in the flow string of a well is a widely known problem in oilfield practice. Accurate data on the pressure gradient of a simultaneous flow of gas and liquid in a vertical pipe are especially useful for the determination of optimum flow-string dimensions. It is well known that with moderate gas and liquid flows such a vertical string acts as a "negative restriction". The pressure drop decreases (1) when the throughput through a given pipe increases, and (2) when at a given throughput the cross-sectional area is decreased. The reason is that, with increasing velocities, the flow becomes more agitated so that the gas slips relatively more slowly through the liquid. With the resulting increase in gas content in the string, the static head decreases. When the area becomes very small, however, the high velocities entail great wall friction, which causes an increase in pressure drop. For a given flow, therefore, minimal pressure drop is obtained by using a certain cross section. This means that, in principle, each well can be provided with an optimum flow string for minimum pressure drop and, hence, maximum possible production rate. The procedure for the selection of the optimum string has been discussed by Gilbert.' A necessary tool in the procedure, however, is accurate knowledge of the pressure gradient to be expected for various values of the governing variables. Another application of pressure-gradient data lies in the field of gas-lift practice: they provide a means of determining the optimum gas-injection rate, optimum injection pressure and optimum injection depth. Much work has already been done in the study of the pressure gradient of vertical gas-liquid flow. Poett-mann and Carpenter2 presented a pressure-gradient correlation based on measurements in wells. This correlation has been found to provide accurate predictions in high-pressure wells and in high-production wells for flow through both tubing and annuli.2-5 However, when their method is checked on low pressure-low production wells or on wells with viscous crudes, serious discrepancies are found. As we shall see in the next section, this is due to the fact that their correlation factor, representing all irreversible energy losses, is given as a function of only one correlation group. Some important variables, such as gas-liquid ratio and liquid viscosity, are not incorporated in this group so that their specific effects are not accounted for. To study also the mechanism of vertical gas-liquid flow outside the ranges covered by the Poettmann-Carpenter publication and extensions, a laboratory investigation has been carried out. This study is founded on a pressure-gradient equation that is based on a pressure balance. To reduce the number of test runs required, a dimensional analysis has been carried out, followed by a selection of relevant dimensionless groups. These groups guided a subsequent experimental study, and with their aid the experimental program could be minimized while still covering the majority of the situations encountered in oilfield practice. In this paper the choice of a formula for the pressure gradient is discussed first. This is followed by a brief description of the experimental setup. Subsequently, the dimensional analysis is discussed and the relevant dimensionless groups are selected, resulting in the experimental program required. The general relationships of pressure gradient and liquid hold-up are then described; various flow patterns and a certain flow instability (so-called "heading") are discussed and a set of correlations is presented which shows a good agreement with the measurements and a few available field

By J. R. Bell, W. J. Childs, J. H. Bucher, G. A. Wilber

A technique for determining recrystallization times as short as 0.10 sec was developed utilizing the "Gleeble", a commercially available testing system designed for the study of short-time, high-temperaLure themal and mechanical processes. The procedure consisted of heating a small tensile specimen to a given temperature of hot deformation, loading to a given reduction in area, unloading, delaying various intervals at temperature, and then reloading- to failure. The magnitude of the ultimate load obtained upon reloading decreased with delay lime as recrys-lallization proceeded. The technique was applied to austenite recrystallization in AISI 1010 and AISI 1010 uith 0.02 pct Cb steels. For each steel the reduction in area given the specimen on the first pull was mainlairred at 30 ± 5 pct and recrystallization times deterntined at various temperatures. The results indicaled a significantly slower rate of recrystallization for the columbium-modified composition, suggested the presence of- a recovery stage in the softening process , and indicated a greatly increased softening rate at a temperatuve where significant allotropic transformation to a partially ferritic Structure could occur. In recent years increasing attention has been paid to the fact that the process of recrystallization of austenite deformed at elevated temperatures is far from instantaneous at many practical hot-working temperatures.1-3 This realization has given rise to such terms as hot cold-working1 or warm-working,2 These terms generally describe processes where the recrystallization rate at the temperature of deformation is slow enough to have an appreciable effect on mechanical properties despite a relatively high deformation ternperature. The mechanical properties of interest can be either the properties at the deformation temperature as in hot-workability studies4 or the room-temperature properties after cooling as in the many recent studies of various thermomechanical processes172 where heat treatment and deformation are intentionally combined to give a unique set of room-temperature properties. Because of this interest in processes where the austenite recrystallization kinetics can be an important variable, the development of quantitative methods of following the course of short-time, high-temperature recrystallization has received increasing attention.l,3,5 The experimental methods to date have, in general, relied upon rapidly deforming the austenite, holding at temperature for various brief intervals, quenching as G.A.WILBER and W. J. CHILDS, Members AIME,are Research-Fellow and Professor, respectively, Rensselaer Polytechnic Institute, Troy, N. Y. J. R. BELL and J. H. BUCHER, Member AIME, are Research Engineer and Research Supervisor, respectively, Graham Research Laboratory, Jones & Laughlin Steel Co., Pittsburgh, Pa. Manuscript submitted March 13, 1968. IMD. rapidly as possible, and then using room-temperature measurements to follow the recrystallization process. Although such methods can be successfully applied to certain alloy steels, the existence of the allotropic transformation during cooling of plain-carbon or low-alloy steels tends to obscure the results. Thus, such room-temperature measurements as hardness and X-ray line widths do not correlate well with the extent of austenite recrystallization before quenching,5 and results based on room-temperature microstruc-tural observations are dependent upon the success in correlating the observed structure with the prior aus-tenitic grain structure.1,3,5 The purpose of the present work was to develop a quantitative method for the determination of short-time, high-temperature recrystallization rates, based on measurements made at the temperature of deformation. EXPERIMENTAL TECHNIQUE The basic technique consisted of heating a small tensile specimen to a given temperature of hot deformation, loading to a given reduction in area, unloading, delaying various intervals at temperature, and then reloading to failure. The data were obtained in the form of traces of load and elongation as a function of time. Due to the high deformation temperature, the strain hardening introduced during initial loading was progressively annealed out with holding time after unloading and the loads obtained upon reloading decreased as this softening proceeded. Although the value of the second load at any Consistent point On the load-elongation curve could have been used as a measure of the degree of softening, the most convenient to use was the ultimate load. The softening indicated by the decrease in the second ultimate load with time is essentially a process of annealing of cold-worked material at a high deformation temperature. Although some recovery grain growth may contribute to such a softening process, it is generally considered that the major softening which must take place to achieve complete removal of substantial Strain hardening will occur by the formation of new, stress-free grains. As the results of this work indicate that essentially complete removal of strain hardening did in fact occur. the primary softening process will be attributed to recrystallization, and specific reference made where it appears that other mechanisms may be contributing to the total observed softening. It would, of course, be of interest to attempt to correlate the results of this work with the actual austenite fraction recrystallized as determined by other techniques. This was not attempted in the present work because it would have required running a large number of additional specimens and, as discussed previously, there is limited assurance that the results would accurately reflect the prior austenite fraction recrys-

Jan 1, 1969

By K. Hodgson, N. C. Joughin

The detailed mechanism of rockbursts and other forms of damage due to rock failure in underground mining is as yet unknown. Studies of the energy changes which necessarily occur when an excavation is made underground seem to provide the best understanding of the phenomena leading to rock failure and damage (see Cook 1). It is common knowledge that the damage and rockburst hazard increases with depth and is greatest at remnants and in regions mined in an irregular manner. The rate at which energy is released as an excavation is enlarged provides an objective parameter to describe mining geometry in terms which are significant to studies of rock failure. In this paper, it is shown empirically that the rate of energy release per unit area increase in the size of an excavation is closely related to seismicity and incidence of damage in the mine. Consequently, it is suggested that the energy release rate is an important parameter in the design of underground mine layouts. Observations were made at two mines: the East Rand Proprietary Mines, Ltd. (E.R.P.M.), at a depth of about 9000 ft, and at Harmony Gold Mining Co., Ltd., at a depth of about 5000 ft. At both of these mines continuous recording seismic networks were in operation. (See Cook and Joughin.3) DESCRIPTION At E.R.P.M. the gold-bearing reef dips at approximately 35 degrees and is presently being mined at depths from 8000 ft to 11,000 ft below surface. The rock surrounding the excavation consists of a very homogeneous quartzite extending far into the hangingwall and footwall. The stoping width is approximately 40 in.

Jan 1, 1967

By Charles Catlett

AN examination of the census reports' for 1880, which contain a large number of complete analyses of typical American iron-ores, indicates that the existence of barium sulphate in intimate association with iron-ore is rather unusual. Of 94 analyses reported there, but four show any barium sulphate, and the maximum content amounts to 3.04 per cent. The note of the occurrence of considerable quantities of barite in connection with certain iron-ores found in the Pinar del Rio Province of Cuba may be of interest. These deposits are found near the little village of Francisco, about 25 miles west of the city of Pinar del Rio on the Western Railroad of Havana. The principal ridges are anticlines, and consist of very thick deposits of sandy shales and slates containing a large number of little veins of white quartz which, by the gradual disintegration of the mass of the material, is left covering the surface. This quartz occasionally shows sulphide, but is generally barren of any such indication. At times lenticular masses of a more or less pure limestone were observed in these slates. The shales, which are of light color, are occasionally talcose. The formation would tentatively be classed as the Lower Cambrian, and the heavy ridges of comparatively horizontal limestone found in the anticlinal valleys as the Upper Cambrian. Along the crest of these slaty ridges deposits of iron-ore are found which are apparently in place, and were probably formed by replacement of masses of limestone similar to those noted elsewhere. The ready disintegration of the slates, and the

Jul 1, 1907

By A. D. K. Laird, John A. Putnam

This paper describes the application of x-ray theory to design procedures in connection with fluid saturation determinations during fluid flow experiments with porous media. A reliable and rapid method for calibrating the x-ray apparatuy is described. Extension of the method to fluid saturation determinations in three-fluid systems is described. INTRODUCTION In rerearch on oil production problems a method is required which will give quickly the quantity of each component of a fluid flow system present at any cross-section of a porous medium. The sample of porous medium under investigation is usually referred to as a core. The ratio of the volume of one component to the total fluid volume is defined as the saturation of the porous medium by that component. This ratio is generally given as per cent saturation. Some means of measuring saturation which have received consideration include: electrical conductivity of the fluids;1,2 emissions from radioactive tracers dissolved in the fluids; the radioactivity of silver caused by reflection of neutrons from hydrogen atoms in the fluids;' the attenuation of a microwave beam. the diminution and phase shift of ultrasonic wave trains.4,5 and the reduction in intensity of x-ray beams in passing through the fluids. X-rays have already been used with some success. Since every material has a different power to absorb x-rays, the reduction in intensity of an x-ray beam as it passes through a core depends on the fluids present. The strength of the emergent beam can be found by converting its energy into a measurable form such as heat or ionic current. or by its effect on a photographic plate or fluorescent screen. The beam strengths could be interpreted as quantities of known fluids in the core if, previously, these beam strengths had been identified with a known combination of the same fluids. With some fluid cornbinations it might be desirable to dissolve powerful x-ray absorbing materials in one or more of the fluids, to increase the differences in the beam strengths for various fluid saturations. Boyer, Morgan and Muskat6 have described a method of measuring two component fluid saturation. One component was air or water; the other. minerat seal oil in which was dissolved 25 per cent by weight of iodobenzene to increase its absorbing power. The x-ray source was a tungsten target tube operated at 43 kv potential. The beam emerging from the core was measured as ionic current flowing across an air-filled ionization chamber by means of an amplifying circuit and galvanometer. Another portion of the beam from the x-ray tube was passed through a metal plate and measured in another ionization chamber. This portion, called the monitor beam, was used as an indication of the performance of the x-ray tube. The galvanometer readings were calibrated against air-oil core saturations, gravimetrically determined. The method was apparently established by experimental means. In the present investigation the available theory of x-radia-tion was surveyed with a view to extending the usefulness of the method and to developing design procedures for its application to measurement of fluid saturation in porous media. Application of the theory permits prediction of relative meter readings to be expected for any combination of porous matrix, various saturating fluids and auxiliary filtering media. It is thus possible to calibrate the equipment in terms of fluid saturation by an indirect but rapid technique. The results of calculations based on x-ray theory indicate. and results of the saturation calibration technique confirm. that a valid measurement of the saturation of the core can be made for any two components and in some cases for three components. THEORY The strength of an x-ray beam, after it has passed through a distance. 1, of matter of density, p, and mass absorption coefficient, µ at a given wavelength, A, may be expressed by the absorption formula I = I0 e ...........(1) where I, represents the intensity of the incident x-ray beam and I is the intensity of the emergent beam. The expression e is called the transmission factor of the material. The variation of I,, with wavelength depends upon the materials through which the x-ray beam has previously passed and upon the spectral distribution of energy at the source of the x-radiation. A group of curves. called spectra. which show the variation of intensity with wavelength and x-ray tube voltage are given in Fig. 1. These curves represent the general radiation from a tungsten target tube. When the tube voltage is greater than 69.3 kv, the characteristic radiation of the tungsten is emitted and is superposed on the general radiation. At a given voltage the minimum wavelength A,,,,, at which energy can be emitted by an x-ray tube is given by the formula 12,340 xml. = ——..........(2) volts where A,.,,.. is in Angstrom units. The wavelength at which the spectra have maximum intensity a1so decreases with increasing x-ray tube voltaue. The area under each curve represents to an arbitrarv scale the total energy emerging from the x-ray tube for that voltage. The variation of µ with wavelength has been determined for many substances and may be found in such references as those by Compton and Allison7 and by Hodgman.8 The phenomenon of absorption is composed chiefly of the capture of photons by the atoms of the absorbing material with associated displacement of electrons, and of the scattering, or the deflection, of the photons by the atoms. Curves of these mass absorption coefficients show jump discontinuities. or absorption edges. at wavelengths which are short enough for the photons,

Jan 1, 1951

By John A. Putnam, A. D. K. Laird

This paper describes the application of x-ray theory to design procedures in connection with fluid saturation determinations during fluid flow experiments with porous media. A reliable and rapid method for calibrating the x-ray apparatuy is described. Extension of the method to fluid saturation determinations in three-fluid systems is described. INTRODUCTION In rerearch on oil production problems a method is required which will give quickly the quantity of each component of a fluid flow system present at any cross-section of a porous medium. The sample of porous medium under investigation is usually referred to as a core. The ratio of the volume of one component to the total fluid volume is defined as the saturation of the porous medium by that component. This ratio is generally given as per cent saturation. Some means of measuring saturation which have received consideration include: electrical conductivity of the fluids;1,2 emissions from radioactive tracers dissolved in the fluids; the radioactivity of silver caused by reflection of neutrons from hydrogen atoms in the fluids;' the attenuation of a microwave beam. the diminution and phase shift of ultrasonic wave trains.4,5 and the reduction in intensity of x-ray beams in passing through the fluids. X-rays have already been used with some success. Since every material has a different power to absorb x-rays, the reduction in intensity of an x-ray beam as it passes through a core depends on the fluids present. The strength of the emergent beam can be found by converting its energy into a measurable form such as heat or ionic current. or by its effect on a photographic plate or fluorescent screen. The beam strengths could be interpreted as quantities of known fluids in the core if, previously, these beam strengths had been identified with a known combination of the same fluids. With some fluid cornbinations it might be desirable to dissolve powerful x-ray absorbing materials in one or more of the fluids, to increase the differences in the beam strengths for various fluid saturations. Boyer, Morgan and Muskat6 have described a method of measuring two component fluid saturation. One component was air or water; the other. minerat seal oil in which was dissolved 25 per cent by weight of iodobenzene to increase its absorbing power. The x-ray source was a tungsten target tube operated at 43 kv potential. The beam emerging from the core was measured as ionic current flowing across an air-filled ionization chamber by means of an amplifying circuit and galvanometer. Another portion of the beam from the x-ray tube was passed through a metal plate and measured in another ionization chamber. This portion, called the monitor beam, was used as an indication of the performance of the x-ray tube. The galvanometer readings were calibrated against air-oil core saturations, gravimetrically determined. The method was apparently established by experimental means. In the present investigation the available theory of x-radia-tion was surveyed with a view to extending the usefulness of the method and to developing design procedures for its application to measurement of fluid saturation in porous media. Application of the theory permits prediction of relative meter readings to be expected for any combination of porous matrix, various saturating fluids and auxiliary filtering media. It is thus possible to calibrate the equipment in terms of fluid saturation by an indirect but rapid technique. The results of calculations based on x-ray theory indicate. and results of the saturation calibration technique confirm. that a valid measurement of the saturation of the core can be made for any two components and in some cases for three components. THEORY The strength of an x-ray beam, after it has passed through a distance. 1, of matter of density, p, and mass absorption coefficient, µ at a given wavelength, A, may be expressed by the absorption formula I = I0 e ...........(1) where I, represents the intensity of the incident x-ray beam and I is the intensity of the emergent beam. The expression e is called the transmission factor of the material. The variation of I,, with wavelength depends upon the materials through which the x-ray beam has previously passed and upon the spectral distribution of energy at the source of the x-radiation. A group of curves. called spectra. which show the variation of intensity with wavelength and x-ray tube voltage are given in Fig. 1. These curves represent the general radiation from a tungsten target tube. When the tube voltage is greater than 69.3 kv, the characteristic radiation of the tungsten is emitted and is superposed on the general radiation. At a given voltage the minimum wavelength A,,,,, at which energy can be emitted by an x-ray tube is given by the formula 12,340 xml. = ——..........(2) volts where A,.,,.. is in Angstrom units. The wavelength at which the spectra have maximum intensity a1so decreases with increasing x-ray tube voltaue. The area under each curve represents to an arbitrarv scale the total energy emerging from the x-ray tube for that voltage. The variation of µ with wavelength has been determined for many substances and may be found in such references as those by Compton and Allison7 and by Hodgman.8 The phenomenon of absorption is composed chiefly of the capture of photons by the atoms of the absorbing material with associated displacement of electrons, and of the scattering, or the deflection, of the photons by the atoms. Curves of these mass absorption coefficients show jump discontinuities. or absorption edges. at wavelengths which are short enough for the photons,

Jan 1, 1951

By J. E. Kirby, L. B. Schnitz, H. E. Stamm III

The production history of a gas-cap-drive reservoir was reproduced by calculations, and predictions were made for operations under primary depletion, pressure maintenance by gas injection, and pressure maintenance by water injection. Although the general equations used in this analysis have been available to the industry for a number of years, the procedure has not had extensive use because of the unusually large time requirements for developing numerical results. The use of digital computing equipment for processing the calculation work generally reduced the time usually required for an extensive analysis of this type. Calculation techniques were developed that were equally applicable for solving the gar-cap-drive problem, whether digital computing equipment or desk calculators are used. These techniques are outlined and the procedures are illustrated by the field problem. INTRODUCTION The analysis of a gas-cap-drive res-ervoir requires the simultaneous solution of the volumetric balance1 and displacement2 equations. Equations for the analysis of combination drive reservoirs were presented by Wooddy and Moscrip in a previous paper." These equations permit the evaluation of possible operating processes in order that the optimum method of pro- ducing the reservoir may be determined. Although these equations represent the basic technique known for the detailed analysis of gas-capdrive reservoirs, the calculation processes are quite complex and time consuming; hence, the procedure has not had extensive use. Recently, these equations were satisfactorily applied with the aid of a digital computer to investigate the completion history of a gas-cap-drive reservoir under the following operating methods: primary depletion, depletion under full pressure maintenance by gas injection, and depletion under full pressure maintenance by water injection. The use of digital computing equipment substantially reduced the time necessary for making the computations and permitted a more concise investigation of the reservoir's depletion history than was possible when the desk calculator was used. This more rigorous investigation also permitted a detailed analysis of the calculation procedures necessary in the solution of the gas-cap-drive problem. This paper outlines the procedure of analyzing a gas-cap-drive reservoir and illustrates the process with an actual field problem. CHARACTERISTICS OF THE RESERVOIR ANALYZED Production in the field under consideration is obtained from a Frio formation of Oligocene age at a depth of approximately 5,740 ft sub-sea. The accumulation, which is bounded on the flanks by strand lines, is semirectangular in shape with gentle narrowing toward the base of the structure as shown by the sand limits and completion map in Fig. 1. The reservoir development indicates an oil column of 319 ft that is defined by gas-oil and oil-water contacts. Basic data relative to the characteristics of the formation and the formation fluids are shown in Table 1. The formation exhibits considerable shaliaess and thinning near the oil-water contact indicating that little or no benefit will be received from natural water influx. This is substantiated by the fact that wells completed near the original oil-water contact have not produced water. A sizeable gas cap is present; however, considerable faulting is in evidence, and it is believed that only a limited portion of the gas cap will aid in oil displacement and pressure maintenance. Although the reservoir was in an initial stage of depletion, it was believed that as soon as the equilibrium gas saturation in the oil zone was exceeded and a number of upstructure wells were invaded by the expanding gas cap, the resulting high gas-oil ratio production would rapidly dissipate the natural drive and lead to limited oil recoveries and inefficient reservoir operation. This analysis evaluated the effectiveness of the primary displacement mechanism and

Jan 1, 1958