Search Documents

By B. K. Larkin, H. K. van Poollen, H. C. Bixel

A detailed treatment is given of the transient pressure behavior of a well located near a linear discontinuity. On either side of the discontinuity, the values of permeability, viscosity, compressibility, and porosity are uniform, but may be different from those on the other side of the discontinuity. These results generalize previously published works which assumed the permeability ratio for media on different sides of the discontinuity was zero, one, or infinity. One application of this material would be to facilitate detection of facies changes or fluid-fluid contacts by transient testing. Solutions, presented graphically, show pressure decline for constant rate fluid production. The range of variables examined includes dimensionless time from 0.1 to 100, rnohility ratio from 0.001 to 1000 and specific compressibility ratio from 0.001 to 1000. For dimensionless tinze less than 0.4, the pressure behavior is essentially that of a well in an infinite reservoir. Thus, reservoir properties near the well may be estimated in the usual manner. Using an overlay technique to match an experimental curve with one of the theoretical curves, it is possible to estimate the distance to a discontinuity by substituting the actual time t and the corresponding dimensionless time , at which a match occurs into the equation where is the diffusivity near the well, and may be estimated from data taken at early time. Because of both the complexity of the solutions and limited computer facilities, well interference problems were not examined extensively. However, it was found that well interference may be masked by the discontinuity in some cases. By superposition of the results of constant rate fluid production, the problenz of pressure buildup may be treated. Such results show that for early shut-in times correct values for the transnzissibility are obtained from conventional semilog plotting. However, erroneous values of the static reservoir pressure are obtained unless data at large shut-in times are taken. INTRODUCTION Transient well testing has yielded estimates of such reservoir parameters as static reservoir pressure, permeability. porosity, and well bore effects. In most tests, it was assumed that the reservoirs were homogeneous and iso-tropic. Anisotropic media' and radial discontinuities were considered.'.' Various studies have been made on multiple-layer reservoirs." ' Linear discontinuities have been described for problems of one-dimensional flow.' In the present paper, a mathematical treatment is given for a well located near a linear discontinuity in an otherwise infinite medium (see Fig. 1). An actual reservoir analogy would be a well in a relatively thin reservoir near a fluid-fluid contact, or a well located in a reservoir exhibiting a sudden change in rock and formation characteristics, such as thickness, porosity, or permeability. A sudden change in porosity or permeability may be caused by clay fill or other types of facies changes. Sudden thickness changes niay be the result of such geological formations as buried shore lines and channel deposits. In otherwise uniform formations, a linear discontinuity may be caused by fluid-fluid interfaces such as oil-water contacts. 3s-oil contacts, or gas-water contacts.

By A. J. Garnier, N. H. van Lingen

Rock downhole is known to be lesc. drillable than when brought to the surface. This must be ascribed mainly to the presence under downhole conditions of a pressure differential across already made chips, which hinders their being lifted. The pressure differential has partly a static and partly a dynamic origin. Balling-up of bits is another consequence of this pressure differerztial. Reduction in penetration rate owing to an increase in the strength of the rock is governed by the difference between the mud pressure and the pressure of the for/nation pore liquid. Rotationally symmetric geostatic stresses have 170 effect on drillability. Fracture of rock when drilling will be brille in most cases. The above is supported by laboratory drilling experiments with drag bill and roller bits an elevated mud, pore, and confining pressures On rocks differing in strength and permeability. INTRODUCTION In oil well drilling drillability of rock is found to decrease with increasing depth of the hole. Naturally deep rock will be more compacted and, therefore, harder to drill than shallow rock of the same type. However, apart from this the drillability of a sample of deep and compacted rock brought to the surface is generally much higher than in its original location downhole. In view of the economic implications of this reduction in drillability, it seems worthwhile to analyze its causes. The origin clearly has to be sought in the difference of environment. The only conceivable factors would seem to be the presence of mud under pressure, the pressure of the formation pore liquid, and the overburden of the rock. Down the borehole the rock is compressed triaxially by mud pressure and overburden. It is well known that the strength of rock is increased when confined by external pressure1. Various authors have, therefore, ascribed the difference in drillability mainly to the strengthening of rock by triaxial compression. Another factor, mentioned by Bobo and Hoch5, is that forces. including "pressure differential forces", tend to hold a dislodged particle in place. However, the conditions determining their magnitude are not clarified nor is their effect on drilling rate assessed. In the laboratory, drilling experiments on pressurized samples of rock yielded evidence that pressure differential forces holding the chips down are the major factor in reducing rate at depth. This paper describes the experiments and shows how the results enable both a qualitative and a quantitative interpretation of the factors determining the effect of chip hold-down on drilling rate. The implication with respect to balling-up and jet action is also discussed. The greatcr part of the experiments were performed in the pressure vessel shown in Fig. 1. The spacc between the sample of rock and the vesscl is divided by "0" rings into three separate chambers. The rock sample is confined laterally by pressurized oil in thc middle chamber. Penetration of oil into the sample is prevented by partly jacketing it with brass foil. Thc pressure of the drilling fluid in the hole can be adjusted via the upper chamber. Permeable rock specimens arc water saturated before being drilled. With a properly plastering mud as drilling fluid the pressure of the pore water can be adjusted independently via the lower chamber.

By J. W. Graham

For many years the author has been cognizant of the difficulty encountered by some in treating with the water influx formulas for unsteady-state fluid flow as pertain to the material balance equation. This has particularly applied in establishing reservoir performance and identifying reservoir pressure, which to the practicing engineer has entailed a trial-and-error procedure, and for others has necessitated resorting to computing devices and reiteration processes. In retrospect this difficulty stems from the fact that reservoir pressure in the material balance formulas, as well as associated with the water influx equations, is an inexplicit term, and the work reported in the past is irrefutable. However, what will be presented in this paper is another approach to the problem, whereby the entire material balance equation will be treated by the Laplace transformation, and reservoir pressure which hereto has been inexplicit, can now be isolated by mathematical procedure to relate that parameter with all the factors contributing to its change. This is the simplification entailed, that treats first with an undersaturated oil reservoir as an integrated effect from the inception of production. The second phase pertains to saturated oil reservoirs that encompass a survey traverse. Although both methods of approach are necessarily different in aspect, the most interesting fact is that the mathematics so deduced are identical. Both the linear and radial water-drive systems are incorporated. for which an illustrated factual example is offered for the latter, treating with a saturated oil reservoir. INTRODUCTIO N What is performed in this work is the simplification of an involved computation by advanced analysis. Although such may be construed as a contradiction when one treats with higher mathematics; nevertheless, when direction is given to such an undertaking the results car. be most revealing. Likewise, it is to be mentioned that the bases for these mathematics have been developed on the expediency of the occasion. This is not to be inferred as a qualification of this work, but rather the demands frequently placed upon the author in his private prac- tice in meeting a time limit. A situation, instead of being fraught with hazards, often has given emphasis to creative thought. What will be entailed in this work is the simplification of the material balance formulas by the Laplace Transformation., Although this reveals entirely new horizons that will be given expression in a forthcoming tract, it suffices in the present instance to limit our attention to this phase of the development that treats both with an undersaturated and saturated oil reservoir. To orient the reader's thoughts as to what is involved in this simplification is the recognition that reservoir pressure, as such, is an inexplicit term in the material balance equation. This is the independent parameter that defines the total history of performance in the author's' unsteady-state water influx formulas, as well as the basis for the physical dependency of fluid behavior within the formation as prescribed in the Schil-thuis' material balance equation. Therefore, to isolate reservoir pressure, which is the most essential factor in any reservoir study, is rather a cumbersome procedure entailing either a trial-and-error calculation for the engineer; or as some prefer, a reiteration process performed on a computing device. However, once such an equation can be transcribed as a Laplace transformation, this inexplicitness so expressed can be alleviated to identify reservoir pressure as an explicit function of all the factors contributing to its change. This is the simplification encompassed, that will treat first with an undersaturated oil reservoir as an integrated effect from the inception of production, and secondly, with a saturated oil reservoir as a survey traverse. Although the two approaches are necessarily different because of the uhvsics involved. it is an interesting commentary that the mathematics are identical, showing the interdependency of the two methods. In order to acquaint the reader with this development, the simplest case will be treated first; namely, an under-saturated oil reservoir subject to a linear water drive. However, what may be construed for this example as an idealistic case is actually a most practical application in certain parts of the world, where the size of the fields are so large that radial water-drive approaches the configuration of a linear drive. Further, to avoid the repetition of much symbolism, frequent references will be made to the work of the author and an associate on Laplace Transformations3,

By Martin Coryell

Ehrenberg was the first to publish and direct attention to the peculiarities of the diatomaceous deposit of Richmond; but how and by whom he was supplied with the material for his microscopical investigations, is not generally known; but we may presume W. B. Rogers, formerly State Geologist of Virginia, furnished the material, and may have given the geological position or horizon. OF these facts we have no available evidence, nor have we any estimate of the thickness of this stratum of diatomaceous or infusorial earth. At the request of a member of this Institute, I have investigated the subject, and offer the following facts as the result, derived principally from information supplied by Mr. John Ott, of the Southern Fertilizing Company, Dr. William Taylor, State Chemist, C. Peticolas, microscopist, Captain Robert H. Temple, civil engineer, Major C. M. Bolton and Mr. C. M. Smith, civil engineers, and formerly engaged in constructing the tunnel under Church Hill, in Richmond. The city of Richmond is generally described as being built on seven hills. This is not the fact, geologically speaking. Originally, there was an immense plateau of land, about two hundred feet above tide, in which a deep excavation was made from the west to east by the James River cutting down to the granite. Opening southward and meeting the valley of the James are various ravines and gullies, in the principal of which are Shockoe and Gillies Creeks. The formation is tertiary, subdivided into all the classified divisions, but at the present moment no one has been found to make this section a special study. Rogers's early classifications and diagrams, with all their imperfections, are the only published records, and the printed editions are long since exhausted. The city of Richmond is really situated in these ravines of the James, Shockoe, Gillies, and their connections. The profile marked A, Plate V, represents the section of the plateau

By Robert H. Richards

It will hardly be worth while to spend time over the discovery of this mine, how lumps of galena were picked up and brought to town, and how legends were told of an old mine from which Revolutionary bullets were made. This is an entirely new mine, with no evidences of former or ancient workings in its immediate vicinity to guide or attract the prospector. The first pit was opened in May, 1874. In August I made a visit to this pit, which was 6 feet wide in the northwest direction, 12 feet long north 80° east, and 6 feet deep..About 4 tons of a galena carrying gray copper, siderite, and quartz, had been taken from this pit. The lumps were all of them evidently float specimens, and varied from 8 to 12 inches thick, and showed very distinct indica-

By B. Silliman

The experiments detailed by Dr. Hunt,* having demonstrated the fact that the copper contained in the "clay ore" of Jones's Mine, was rendered completely soluble in the bath of ferrous chloride, used in the Hunt and Douglas process, after heating in contact with carbonaceous matter in a close vessel, I set myself to devise a form of muffle furnace adapted to the treatment of large quantities of these and similar ores by a continuous process. The result is shown in the diagrams on Plate X, reduced from the working drawings, after which this furnace was built early in 1876, at Phœnixville, by the Chemical Copper Company on their works at this place. The peculiar character of this ore determined the form, dimensions, and position, with reference to charging and discharging of these muffles. The ore arrived from the mine with from 20 to 25 per cent. of moisture, and when dried at 21o F., or more slowly at lower temperature, it falls to a light incoherent powder, with occasional lumps of undecomposed rock. In this condition it is readily mingled with coal-dust or any like reducing agent, and requires no other preparation for the muffle than the use of the shovel, to mix it well with the reducing agent. As it is a remarkably good non-conductor of heat, it was obvious that the mass, to be heated through in a reasonable time, must not be too thick, while the weight of the mass must be sustained in a way to avoid undue strain upon the walls of the muffle. These walls must be as thin as practicable, to favor the more rapid transmission of heat, and must, therefore, be so constructed as to admit of being stayed on the sides at frequent intervals to resist the lateral thrust of the very mobile mass of pulverulent ore, which, for obvious reasons, must be charged at the top and drawn from the bottom of each chamber. These considerations led to the form adopted, viz., two vertical muffles, standing upon very strong bridge tiles, seen in longitudinal section, on the line E, F, and in plan in Fig. 3, on the same plane. (Plate X.) The vertical section of the muffle is seen in Fig. 4, drawn in the plane G, H, of the longitudinal section, Fig. 6, which is the key to all sections. The walls of the muffles are built of the best firebrick, very care-

By John F. Blandy

IT is not my intention in this article to consider this subject in the light of the geographer or geologist, but rather in that of the mining engineer, and to endeavor to show the necessity and value to him of studying the minor details of the field or district in which he may be engaged. I t has been said by a distinguished geologist that "the natural boundaries of the geographer are rarely described by right lines. Whenever they occur, however, the geologist may look for something remarkable." And again : " Physiognomy is no idle or doubtful science in connection with geology. The physiognomy of a country indicates almost invariably its geological character." Our country furnishes examples enough of the truth of these remarks, and in no part of it is it so beautifully illustrated as in our own State, which has been so ably described by Prof. Lesley in his admirable little work on the "Topography of Coal." As the great mountain lines become guides to the geologist in determining and mapping out the geography of the formations, so do the minor details of topography become guides to the mining engineer in conducting his explorations, discovering of dislocations, and determining the character of the deposits which are the object of his search. I t is therefore requisite that he should, whilst making his surveys, pay the strictest attention to every variation of the landscape and surface of the ground. Until he has mapped out every feature he may be unable to read the geological changes which may have

By Willard P. Ward

Having been engaged during the past year in the manufacture of ferro-manganese in a blast furnace, I have imagined that some further information on this subject might be of interest to that large number of members of the Institute who are engaged in the iron and steel industries. The great question of the economical production of good steel, ingot metal, or homogeneous iron (we scarcely know now what to call it) from good materials, and at reasonable prices, has, thanks to the intelligence and energy of a few American engineers, been successfully solved in this country. Nearly all first-class railroads, and many that would not rank so high, have adopted the Bessemer rail. Martin boiler-plates are rapidly gaining ground, in competition with

By I. Lowthian Bell

THERE has been probably no improvement introduced into the manufacture of iron which created more surprise in the minds of practical smelters and of scientific men than Neilson's discovery of the hot blast. In 1829, Messrs. Dunlop & Co. consumed at the Clyde Works, near Glasgow, nearly 79 tons of coal to make a ton of pig iron, of which about 20 cwts. was employed under the boilers of the blowing engines, leaving nearly 7 tons for the consumption of the furnace itself. In the year 1833, by heating the air to 612' F. (322' C.), they redoccd the 79 tons to 3 tons by the mere burning of 8 cwts. of coal in the apparatus for raising the blast to this temperature. The statenlent just given is strictly true in a commercial sense, but when we come to consider the question as one of physical science, it is necessary to eliminate some of the conditions which conduced to so extraordinary a result. Anterior to Neilson's time the fuel employed in smelting iron was coke, and it was supposed at that period, erroneously, however, that coal in its raw state, when burned with cold air, was totsally unfit for the purpose in question. This mineral, as it was obtained from the Lanarkshire coal-field, contained about 35 per cent. of volatile matter, but the process of charring it was performed in so unskilful a manner that instead of receiving 88 cwts. of coke from 135 cwts. of coal, 60 cwts. only was the product of the operation. Again, the reduced colisumption of carbon in the furnace itself was accompanied by a corresponding diminution of blast, and this was followed by an economy of half a ton in engine coal. After making a proper allowance, however, for all these collateral circumstances, we are within the mark when we admit that in actual coke, which is really the form the fuel has to assume before it is burnt in the blast-furnace, not less than 30 cwts. was saved by heating the blast to 615° F. (322° C.). This unexpected economy, be it remembered, was obtained by the combustion of the 8 cwts. of coal, somewhat wastefully applied in the hot-air apparatus. Numerous have been the opirlions advanced by different authorities in the scientific world to account for the apparent anomaly, bat none,

By T. Sterry Hunt

This remarkable mine, to which attention has lately been drawn, is situated not far from the New River, in Ashe County, North Carolina, on a spur of the Blue Ridge which lies between the main crest of this name and the so-called Iron Mountain, a part of the Unaka or Smoky Mountain range. It is on a true fissure-lode, which cuts the strata of gneiss and mica-schist of the region. These, which have the lithological characters of what I have called the Montalban or White Mountain series, have here a dip of about forty-five degrees to the southeast, while the lode is vertical, with a course north 60' east. Both the country-rock and the lode, as is usual in this region, are decomposed to considerable depths; and the latter exhibits a cap of hydrous peroxide of iron or gossan, blocks of which, left by the denudation of the softer inclosing rock, are scattered over the surface, and first called attention to the locality as a source of iron ore. In sinking upon this lode, it was found that the porous gossan at a certain depth became charged with carbonate and red oxide of copper, and lower down was replaced by sulphuretted ores of this metal of remarkable richness. This deposit was worked for copper, in an irregular manner, before the late civil war, but was afterwards abandoned, and only reopened during the present year, by Messrs. S. S. & J. E. Clayton, of Baltimore. The outcrop has been traced for 1900 feet, and the sulphuretted copper ores have been met with in five shafts in a distance of 661 feet, a drift having moreover been carried for the above distance through the solid ore. The breadth of the lode, so far as opened, varies from six to eight feet to fourteen feet in places, and that of the outcrop of gossan is, in some parts, twenty feet. After comparing this deposit with the copper mines of Ducktown, Polk County, Tennessee, and those of Carroll County, Virginia, it is evident that we have, at Ore Knob, to do with the same kind of black ores as were met with in these mines, and will be noticed further on, Beneath the gossan, at Ore Knob, is found, in some parts, an iron-black, friable, drusy, crystalline, sulphuretted ore, inclosing grains of quartz, garnet, and magnetite, besides a black non-magnetic mineral, not attacked by nitric acid. This ore contains about 36 per cent. of copper, and portions have the mineralogical charac-

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By C. F. Weinaug

In recent years, several methods have been proposed for correlation of physical properties and behavior of hydrocarbon systems. These correlations were found to be unsatis. factory for predicting differential gas liberation and shrinkage data now obtained experimentally. To achieve the desired accuracy, methods were developed to adjust the correlations with two experimental pieces of data. The required data are the saturation pressure at reservoir temperature and the specific volume or density at these conditions. The saturation pressure is used to adjust the vapor-liquid equilibrium-ratio correlations. The density measurement is used to determine the apparent density of liquid methane. With these adjustments, differential gas liberation and shrinkage are predicted within the experimental error. EQUILIBRIUM RATIOS Equilibrium ratios published by the National Gas Assn. of America' were used as a basis for the calculation. The convergence pressure for the sample was determined by using a hypothetical binary system, with methane as one component and the remainder of the system lumped together to form the other component. For the critical temperature of the heavy constituent, the weight average of the critical temperatures of the ethane and heavier components was used. The correlation of pseudo-critical temperature, density and molecular weight of Matthews, Roland and Katzl was used in this calculation for the hexane and heavier fractions. Points for the critical properties of methane and the hypothetical heavy constituent were placed on a chart of the critical loci of hydrocarbon systems. A critical loci for the hypothetical binary system was drawn, and the critical pressure of this system was determined at the reservoir temperature. This critical pressure became the convergence pressure for which equilibrium ratios were read from the NGAA4 charts. The NGAA equilibrium constants for C, through C, were used for the C,; through C, fractions since they were separated on the bases of normal paraffin boiling points. The correlation of Hoffman, Crump and Hocott5 as used to estimate the equilibrium constant for the C fraction. The NGAA eauilibrium-ratio charts are based primarily on paraffin systems and have to be corrected for the presence of naphthenes and aromatics. The correlations in the literature proved too unsatisfactory for determining these corrections. EXperimental data, the saturation pressure, is used with the correlation of Soloman2 to make this adjustment. The correction factors of Soloman were cross-plotted so that the factors for ethane and heavier components were made a function of the factor for methane. a2 = 0.7927 at + 0.1646 .........(1) where a, = correction coefficient for for methane NGAA equilibrium ratio, and a2t = correction coefficient for ethane and heavier NGAA equilibrium ratios. This equation and the equation for the bubble point are solved simultaneously for the correction coefficients. The bubble-point equation is anKnXn = 1.0 . . . . (2) The equilibrium constant for the saturation pressure at reservoir temperature is used. This procedure not only corrects for the relative volatility between components, but also assures that the results will be directly related to the experimental saturation pressure. METHANE DENSITY After determining a set of equilibrium ratios, the apparent density of condensed methane at 60°F and 14.7 psia was estimated. This was found to be necessary because the correlations of Standing3 ere found to be unsatisfactory. The law of additive volumes was employed to make this calculation.

By H. Engelmann

The very extensive development of a brown coal formation, in the region of the Rocky Mountains, is well known to all of yon. The existence of these coals was known years ago, but they were of no practical importance until within a few years, when the mineral wealth of that country began to be developed and railroads were built. As their qualities became better known, their adaptability to domestic manufacturing and metallurgical purposes was more appreciated. They are still much undervalued, although they are mined at some points on a very extensive scale. A large portion of this coal formation is of cretaceous age, while it appears to extend into the tertiary epoch; but whatever difference there may be in the ages of these coals of different localities, they all have some common characteristics which distinguish them from our Eastern stone coals of the

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing

By J. J. Bodmer

In 1855, Franz Uchatius patented, in England, his process of manufacturing cast steel. The first experiments, on a practical scale, were made at the Ebbw Vale Iron Works, Monmouthshire. The charge consisted of a mixture of cast iron with about 20 per cent. of a pure iron ore, and with or without other ingredients. In order to obtain the cast iron in a subdivided condition, Uchatius granulated molten pig-metal by running it into water, which during the operation was kept well agitated, and found that "the finer the iron is granulated, the softer will be the resulting steel." The writer is not aware that granulated iron has been used practically, otherwise than in the above-named process; and until 1866 no other process of subdividing iron than by granulation was known. In April, 1866, an English patent was obtained by the writer for subdividing blast-furnace and other slags and metals in a molten condition, by passing the same through one or more pairs of rolls. In reference to iron the idea was, in following up the direction pointed out by the Uchatius process, to obtain better results in the puddling process by means of a thorough amalgamation (not simply a mixture) of the iron with the oxides and other ingredients in the charge itself. Experiments were made. A quantity of iron direct from the blastfurnace was subdivided (or laminated) by passing it through a pair of plain rolls. The rolls were hollow, and water was made to pass through them to keep them sufficiently cool. Without giving differential speed to the rolls, sheets were obtained about one-sixteenth inches thick, and from about one hundred square inches surface, downwards. With differential speed, the iron falls from the rolls in

By D. Havlena, A. S. Odeh

The material balance equation used by reservoir engineers is arranged algebraically, resulting in an equation of a straight line. The straight line method of analysis imposes an additional necessary condition that a successful solution of the material balance equatiott should meet. In addition, this algebraic arrangement attaches a dynamic ineuning to the otherwise static material balance equation. The straight line method requires the plotting of one variable group vs mother variable group. The sequence of the plotted points as well as the general shape of the resulting plot is of utmost importance. Therefore, one cannot progrm the method entirely on a digital computer ar is usually done in the routine solution of the material balance equation. If this method is applied, then plotting and anaIysis are asential. Only the appropriate equations and the method of analysis and interpremtion with comments and discussion are presented in this paper. Illustrative field examples for the various cases treated are deferred to a subsequent writing. INTRODUCTION One of the fundamental principles utilized in engineering work is the law of conservation of matter. The application of this principle to hydrocarbon reservoirs for the purpose of quantitative deductions and prediction is termed "the material balance method of reservoir analysis". While the construction of the material balance equation (MBE) and the computations that go with its application are not difficult tasks, the criteria that a successful solution of the MBE should fulfill have always been a problem facing the reservoir engineer. True and complete criteria should embody necessary and dcient conditions. The criteria which the reservoir engineer uses possess a few necessary but no sufficient conditions. Because of this, the answers obtained from the MBB are always open to question. However, the degree of their acceptability should increase with the increase in the number of the necessary conditions that they should satisfy. Generally, the necessary conditions commonly used are (1) an unspecified consistency of the results and (2) the agreement between the MBE results and those determined volumetrically. This second criterion is usually overemphasized. Actually, the volumetrically determined results are based on geological and petrophysical data of unknown accuracy. In addition, the oil-in-place obtained by the MBE is that oil which contributes to the pressure-production history,' while the volumetrically calculated oil-in-place refers to the total oil, part of which may not contribute to said history. Because of this difference, the disagreement between the two answers might be of paramount importance, and the concordance between them should not be overemphasized as the measure of correctness of either one. In this paper, a third necessary condition of mathematical as well as physical significance is discussed. It is not subject to any geological or petrophysical interpretation, and as such, it is probably the most important necessary condition. It consists essentially of rearranging the MBE to result in an equation of a straight line. This straight line method of the MBE solution has invalidated a few long time accepted concepts. For instance, it has always been advocated that if a water drive exists, but one neglects to take it into account in the MBE, the calculated oil-in-place should increase with time. The straight line method shows that in some cases, depending on the size of the neglected aquifer, the calculated oil-in-place might decrease with time. The straight line method requires the plotting of a variable group vs another variable group, with the variable group selection depending on the mechanism of production under which the reservoir is producing. The most important aspect of this method of solution is that it attaches a significance to the sequence of the plotted points, the direction in which they plot, and to the shape of the resulting plot. Thus, a dynamic meaning has been introduced into the picture in arriving at the final answer. Since the emphasis of this method is placed on the interpretation of the sequence of the points and the shape of the plot, one cannot completely automate the whole sequence to obtain "the best value" as normally done in the routine application of the MBE. If one uses the straight line method, then plotting and analysis are musts. The straight line method was first recognized by van Everdingen, et al,2 but for some reason it was never fully exploited. The advantages and the elegance of this method can be more appreciated after a few cases are carefully treated and worked out by it. SOLUTION OF THE MATERIAL BALANCE EQUATION SATURATED RESERVOIRS The MBE for saturated reservoirs written in AIME symbols is

By K. E. Brown, F. W. Jessen

By utilizing an 8,000-ft experimental field well equipped with 10 gas-lift valves and 10 Maihak pressure recorders, gas-lift tests were conducted with port sizes ranging from 5/16 through I in. The well was equipped to provide accurate means of measuring surface pressures, temperatures, quantity of injection gas and fluid production. The tests were conducted in 2%-in. OD tubing, and the well was making 95 per cent water. A complete evaluation of gas-lift-valve port sizes shows the relationship of per cent recovery, gas-liquid ratios, minimum pressure created at the operating valve and horsepower requirements for each port. The length of time necessary for the fluid in the tubing to reach equilib.rium conditions after each cycle was recorded. Fallback of fluid at depths of 477, 969, 1,685, 2,493 and 4,290 ft was noted. For each port size, pressure loads of 2.50, 300, 350, 400 and 450 psi were lifted with a valve operating at approximately 550 psi at 6,000 ft. Gas-liquid ratios for each load were varied from excess gas to a gas volume per cycle whereby the load failed to reach the surface. Numerous curves are presented in evaluating the accumulated dara. The results show a 1-in. port to be the most efficient under all conditions. The production of intermittent liquid slugs against different-sized surface chokes was evaluated. These tests were conducted from a 7/16-in. ported valve at 4,072 ft. Tests indicate that, when possible, a %-in. in diameter choke or larger should be used at the surface. In the past few years most of the advancement in gas-lift operations has been made in continuous-flow operations. Yet, it is estimated that at least 70 per cent of the wells on gas lift in the United States are of the intermittent type. Since the term "slug flow" is sometimes used in both intermittent- and continuous-flow operations, it would be well to distinguish between the two types of flow. Continuous-flow gas lift is defined as a method whereby the fluids are produced at a continuous rate at the surface. This generally requires a continuous injection of gas through a surface choke; however, various other control devices sometimes are installed to eliminate freezing, to shut-off gas during natural flow periods, etc. The actual flow of fluids in the tubing may be of the slug type (one of many flow patterns known to exist in continuous flow). Intermittent flow is defined as a method of gas lift whereby the liquid is produced in separate piston-type slugs. Perhaps this type of flow could best be thought of as a ballistic-type flow where the liquid leaves bottom as a piston, propelled by a slug of expanding gas. Gas generally is injected through some type of control at the surface at predetermined intervals. However, the valve may have characteristics whereby gas can be injected through a small choke and still result in a ballistic-type flow. The purpose of the experimental work was to evaluate the most efficient port size to be used on the operating valve for the ballistic type of lift and, in addition, to establish the importance of utilizing a surface choke large enough to allow slugs to be produced without detrimental effects. This work is part of a compre- hensive study of both intermittent-and continuous-flow gas lift, representing a joint project conducted by the Ohio Oil Co., the Sun Oil Co., Otis Engineering Corp, and The U. of Texas. The problem of evaluating port sizes has been given little previous attention. Some work undoubtedly has been done which has not been published to date. Some tests were conducted when the wireline, mechanically-opened valve (Nixon) first came on the market. This valve was capable of utilizing full tubing area as its port size. It is known that this was a very efficient valve, but to the authors' knowledge the results of tests have never been published. EXPERIMENTAL EQUIPMENT These tests were conducted on an actual field well, the Ohio-Sun Unit Well No. 2-E, in the North Markham-North Bay City field, Matagorda County, Tex. The well incorporated 23/8-in. OD tubing and produced 95 per cent water. Since the running of equipment was to be quite elaborate and expensive, a well was selected in which both intermittent- and continuous-flow tests could be conducted. This particular well was capable of producing in excess of 1,000 B/D of liquid (95 per cent salt water), yet with a 3/64-in. in diameter bottom-hole choke, production was controlled to 82 B/D. Most of the intermittent tests were conducted at this low rate. Figs. 1 and 2 show all the surface and down-hole equipment. As can be seen, every attempt was made to insure that ample equipment was available for reliable testing procedures. Fig. 1 shows the surface testing equipment. The input gas was controlled first by a regulator, then

By R. D. Vaughn

The limiting cases of non- Newtonian fluids flowing inside a concentric annular duct are developed without using a model of the fluid behavior. The solutions provide limits with which to test the various models of fluid behavior such as the power law and Bingham plastic models. The results of previous theoretical work are discussed in terms of limiting cases. This limiting case study also shows that experimental work on flow of non-Newtonian fluids in annular ducts should be confined to ducts for which the ratio of the radius of the inner wall to that of the outer wall is less than 0.3 and preferably less than 0.2. INTRODUCTION During the last 10 years the problem of laminar flow of non-Newtonian fluids in concentric annuli has received much attention large1y because of its application to the hydrodynamics of the wellbore. Recently solutions utilizing the power law1,2 and Bingham plasticl,3,4 models have been published. In this paper the method of limiting cases, which has been successfully applied to laminar-flow heat transfer516 will be applied to the problem of flow of non-time dependent, non-Newtonian fluids through annuli. This method permits solutions for the limiting cases to be made without using a model of unknown validity. The solutions, therefore, provide limits with which to test the various models which have been or will be proposed. A pertinent conclusion concerning the region of experimental work is also provided. DEVELOPMENT OF LIMITING CASES The limiting cases for the axial flow of fluids in concentric annuli may be defined with reference to Fig. 1. It is possible to define two limiting cases which pertain to the physical dimensions of the annulus. First, the annulus must degenerate to a circular pipe as the radius of the inner wall decreases or, as k = (k R/R) ? 0. Second, the annulus must approach the limit of parallel plates of infinite extent as the spacing between the inner and outer tubes becomes small in comparison with the radius R of the annulus, or as k?1. It is also possible to ascertain three limiting cases which pertain to fluid behavior. With reference to Fig. 2, as a fluid becomes progressively more pseudoplastic, the shear stress-shear rate relationship progressively approaches the indicated horizontal line more closely. At this limit the shear stress becomes independent of the shear rate. At the other extreme of increasingly dilatant behavior, the vertical asymptote is approached. Intermediate between these two limiting cases lies the case of the Newtonian fluid. Fluids which exhibit a yield shear stress also approach the limiting case of "infinite" pseudoplastic behavior. The remainder of the theoretical development is concerned with the description of and the solutions