Search Documents

Discussion of the paper of Max Hebgen, presented at the Butte meeting, August, 1913, and printed in Bulletin No. 80, August, 1913, pp. 1910 to 1933. A MEMBER :-Have there been any definite studies made in the hydrology of this system ? That is to say, what is the total extent and area of the drainage or catchment basin above these plants ? What is the total annual precipitation, and what is the percentage of run off ? That is a very interesting question to engineers and a question engaging a great many, because, unfortunately, there is not the amount of information on that subject which there ought to be.

Jan 11, 1913

By L. D. Mullins, W. R. Bartlett, R. H. Barham, W. L. Wahl

This paper describes an electrical model and its applicution to the analysis of four reservoirs in Saudi Arabia. The model has 2,501 mesh points and represents 35,000 sq miles of the Arab-D member. Details of modeling such as mesh size, control problems and standards of performance in matching reservoir history are discussed. The particular performance match achieved for the Arab-D rnember is presented. Details such as permeability barriers, aquifer depletion and iwerference between oil fields are given. The performance match realized in the Abqaiq pool is presented in detail. INTRODUCTION The resistor-capacitor network and associated control equipment described in this paper comprise an electrical analog of a reservoir system. Similar equipment has been used to study the transient response of reservoirs for many years. The unique feature of the model and application to be described is the extremely large size of the model and reservoir system, and the detail observed in simulating the reservoir with the model. The Arabian American Oil Co. first became interested in analog computers for simulation of oil reservoirs in 1949. Since that time, several models have been developed, each more elaborate and refined so that the reservoir system might be more closely simulated. The current model is the latest in a series designed, built and operated by the Field Research Laboratory of Socony Mobil Oil Co. in collaboration with Aramco. It has been and continues to be used to study the regional performance of the Arab-D member limestone reservoir. The Arab-D member is one of the Middle East's most prolific producing horizons.. THE MODEL The theory of simulating a reservoir system with an electrical system has been presented in the literature.1-5 Therefore, this paper will not discuss the theoretical aspect of the problem except to point out the correspondence between the fluid system and electrical system, as shown in Table 1. In general, the complete model is made up of input devices, output devices, central control and a resistance-capacitance (RC) network. At times, the RC network alone is referred to as the "model". However, it should be evident from the text which meaning is attached to the word "model". A discussion of the equipment follows. THE RESISTANCE-CAPACITANCE NETWORK The RC network consists of 2,501 capacitance decades interconnected through 4,900 resistance decades. The components are arranged to form a rectangular network of 2,501 mesh points in a 41- X 61-mesh array. Imposing the mesh grid system on the continuous reservoir system divides the reservoir into discrete areal segments. These discrete segments may be of various sizes. More precisely, the mesh size need not be uniform throughout the model. The RC network is fabricated in two sections which are connected at the top. An inside view of the "tunnel" formed by the two sections is shown in Fig. 1. The height and width of the tunnel are shown in the figure. Numerals appear along the bottom and along the back opening of the tunnel. These numbers denote the x and y co-ordinate positions of mesh points. Fig. 2 presents a rear view of one-half the model. The length dimensions of the model, as well as a rear view of the capacitor decades, are shown in this figure. The control dials used in adjusting the resistance and capacitance values on the model can be seen in the enlarged portion of the model shown in Fig. 3. The electrical capacity at any mesh point can range from 0 to 1.0 microfarads set to the nearest tenth of a microfarad. The electric resistance connecting any two mesh points can range from 0 to 9,990,000 ohms set to the nearest 1,000 ohms. External capacitors may be added to any or all mesh points if the need arises. The values of electrical resistance and capacitance are adjusted manually by manipulating the two types of decade units. INPUT EQUIPMENT A considerable quantity of equipment is used to control the input to the RC network. These input devices are

By Clyde E. Williams

MORE has been said about the iron and steel situation on the Pacific Coast than has been done .about it; but perhaps as much has been done as conditions have warranted. The production of finished steel in an amount sufficient to take care of the local demand for certain sizes and shapes is an established industry on the .coast. The amount .and the variety of products, so made, changes with the demand or them which is largely dependent upon the growth in population. An increasing demand will not always increase the quantity of steel produced locally because the cost of production may rise so -high as to let in steel from the outside. The reason for this is that all of the steel made in the West comes from steel scrap, the supply of which is .limited. From the standpoint of the metallurgist, it would. be desirable to subititute pig iron for some of ' the steel scrap, but this cannot be done profitably because the pig iron now available is too expensive. In other words, the large future demand for steel and pig iron on the Pacific Coast must be supplied by reduction of ores in the West, unless we are to be content with paying high transportation costs for material produced elsewhere. ."

Jan 1, 1924

By G. Rowan, M. W. Clegg

The basic equations for the flow of gases, compressible liquids and incompressible liquids are derived and the full implications of linearising then discussed. Approximate solutions of these equations are obtained by introducing the concept of a disturbed zone around the well, which expands outwards into the reservoir as fluid is produced. Many important and well-established results are deduced in terms of simple functions rather than the infinite series, or numerical solutions normally associated with these problems. The wide range of application of this approach to transient radial flow problems is illustrated with many examples including: gravity drainage of depletion-type reservoirs; multiple well systems; well interference. INTRODUCTION A large number of problems concerning the flow of fluids in oil reservoirs have been solved by both analytical and numerical methods but in almost all cases these solutions have some disadvantages — the analytical ones usually involve rather complex functions (infinite series or infinite integrals) which are difficult to handle, and the numerical ones tend to mask the physical principles underlying the problem. It would seem appropriate, therefore, to try to find approximate analytical solutions to these problems without introducing any further appreciable errors, so that the physical nature of the problem is retained and solutions of comparable accuracy are obtained. One class of problems will be considered in this paper, namely, transient radial flow problems, and it will be shown that approximate analytical solutions of the equations governing radial flow can be obtained, and that these solutions yield comparable results to those calculated numerically and those obtained from "exact" solutions. It will also be shown that the restrictions imposed upon the dependent variable (pressure) are just those which have to be assumed in deriving the usual diffusion-type equations. The method was originally suggested by Guseinov,1 who (a) postulated a disturbed zone in the reservoir, the radius of which increases with time, and (b) replaced the time derivatives in the basic differential equation by its mean value in the disturbed zone. In this paper it is proposed to review the basic theory leading to the equations governing the flow of homogeneous fluids in porous media and to consider the full implications of the approximation introduced in linearising them. The Guseinov-type approximation will then be applied to these equations and the solutions for the flow of compressible and incompressible fluids, and gases in bounded and infinite reservoirs obtained. As an example of the application of this type of approximation, solutions to such problems as production from stratified reservoirs, radial permeability discontinuities; multiple-well systems, and well interference will be given. These solutions agree with many other published results, and in some cases they may be extended to more complex problems without the computational difficulties experienced by other authors. THEORY In order to review the basic theory from a fairly general standpoint it is proposed to limit the idealising assumptions to the minimum necessary for analytical convenience. The assumptions to be made are the following: 1. That the flow is irrotational. 2. That the formation is of constant thickness. 3. Darcy's Law is valid. 4. The formation is saturated with a single homogeneous fluid.

By E. H. Timmerman, A. F. van Everdingen, J. J. McMahon

The prevent paper contains a method which combines the material balance equation' with the water influx equation' to obtain reliable values for the active oil originally in place and a quantitative evaluation of the cumulative water influx. The method is illustrated by an application to a reservoir without original gas cap. In the absence of an original gas cap. results may be obtained using only field production and pressure data, PVT analyses. and a minimum of subsurface information. CHARACTERISTICS OF RESERVOIR AND RESERVOIR FLUIDS Production in the field under review is obtained from the top of the Wilcox formation of Eocene age, at a depth of approximately 8,100 ft subsea. The structural map, Fig. 1, shows that the accumulation, half elliptical in shape with the long axis in east-west direction, is trapped to the north and on the upthrown side of a normal fault. A number of east-west faults with additional minor faulting in random directions are found in the general area, and seismic information and production performance indicate in this particular case the existence of a second fault a short distance downdip, with a strike more or less parallel to the fault controlling the accumulation. It is estimated that 1,830 acres were originally underlain by oil and that the maximum thickness of the original oil column, about 37 ft, was only slightly greater than the maximum thickness of the oil-bearing sand. Maximum net sand thickness was estimated at approximately 26 ft. From logging and subsurface information the gross sand volume has been placed at 37.400 acre-ft. which reduces to about 27,500 acre-ft after probable non-productive intervals have been deducted. The following table presents laboratory analyse; data for the Wilcox sand cores. Average Average Samples Considered Permeability Porosity All Samples 236 md ..19.9percent Samples with k>l md 256 md 20.5 per cent Samples with k>8 md 275 md 20.9 per cent Interstitial water was estimated at 15 per cent. The initial gas-oil ratio was approximately 900 cu ft/bbl. and the initial shrinkage factor was determined by the laboratory to

Jan 1, 1953

By W. B. Braden

This paper presents a suitable method for predicting gas-free oil viscosities at temperatures up to 500F knowing only the API gravity of the oil at 60F and the viscosity of the oil measured at any relatively low temperature. The API pravity and the one viscosity value are used as parameters to determine the slope of a straight line on the ASTM slanaord viscosity-temperature chart. Then, knowing the slope of the line and one point on the line, the vrscosities at higher temperatures can be determined. The line slope correlations were developed at I00 and 210F since viscosity data are frequently measured at these temperatures. A procedure is given for predicting line slopes from measurements at other tetnperatures. A nomogram is furnished for solving the relationship. The correlation has been evaluated at temperatures up to 5OOF for oils varyzng in gravity from 10 to 33 " API. The correiution is applicable only to Newtonian fluids. Comparison at 500F of true viscosities and those predicted from values at 100F shows an average deviation of 3.0 per cent (maximum deviation of 6.0 per cent). Predictions from the values at 21 0F for the same oils how an average deviation of 1.5 per cent (maximum deviation of 3.4 per cent). INTRODUCTION Correlations have been developed by Beal' and by Chew and Connally' for predicting viscosities of gas-saturated oils at reservoir conditions. Each of these correlations requires a knowledge of the solution gas-oil ratio and the viscosity of the gas-free oil at the reservoir temperature. For temperatures below 350F, measurements of the gas-free oil viscosities can be made easily using commercially available equipment. In thermal recovery processes, however, reservoir temperatures well in excess of 350F are encountered. Viscosity measurements at such conditions are more difficult and time consuming and require modification of existing equipment or the construction of new equipment. Measurements are further complicated by the difficulty of handling highly viscous oils associated with thermal recovery processes. Therefore, it is desirable to have a correlation which allows accurate prediction of viscosities at high temperatures. A commonly used technique for predicting viscosities at high temperatures is to measure the viscosities at two lower temperatures, plot the values on ASTM standard viscosity-temperature charts and extrapolate to the temperatures desired. If either of the values is slightly in error, the extrapolated value can be significantly in error. To justify an extrapolation, three points are actually necessary. This procedure can consume much time, particularly with heavy oils. Considering the cost of viscosity measurements, it would be desirable to eliminate the need for direct measurements by having correlations which would allow viscosity predictions from other physical or chemical properties. Beal1 investigated the possibility of correlating viscosity with oil gravity at temperatures from 100 to 220F. While showing that a general relationship exists, he also found significant deviations. It is possible that correlations will be developed based on oil composition as more information becomes available. While not eliminating the need for viscosity rneasurements, the method presented herein requires that only one viscosity measurement be made. The API gravity must also be known. The theory is based on the fact that the viscosity of paraffins (high gravity) changes less with temperature than does the viscosity of naph-thenes or aromatics (low gravity). The gravity. therefore, is used as a parameter to determine the slope of a straight line on the ASTM standard viscosity-temperature charts. The correlation is applicable only to Newtonian oils, and deviations due to thermal decomposition and nonhomo-geneity cannot be predicted. Oils containing additives have not been evaluated. PROCEDURE Fifteen oils were used in developing the correlation; eight were crudes and seven were processed oils. Oil gravities varied from 9.9" API (naphthene base) to 32.7' API (paraffin base). The temperature range studied was 81 to 516F. Each oil used had a minimum of three viscosity measurements and each plotted essentially as a straight line on the ASTM charts. In all, 91 viscosity measurements were used in the correlation. Saybolt, rolling ball and capillary tube viscometers were used for the measurements. Viscosity data for Samples 1, 2, 4, 7, 10, 11 and 14 were obtained in Texaco, Inc. laboratories. The data for Samples 3, 5, 6, 8, 9, 12 and 15 were from Fortsch and Wilson,3 and data for Sample 13 were from Dean and Lane.' All data points used in the correlation are plotted in Fig. 1. It is seen that some of the viscosity values deviated slightly from the straight-line plots at the higher temperatures. Properties of the oils after exposure to the

Jan 1, 1967

By T. Fraser, R. L. Sutherland

Under a cooperative agreement between United States Bureau of Mines and the Truax-Traer Coal Company, some operating-scale experiments have been made with the cyclone thickener in the preparation plant of the Shamrock No. 1 mine at Kayford, West Virginia. These tests were made with a 14-in. cyclone of the type developed by the Dutch State Mines in Limburg. It was applied experimentally to the thickening of solids in washery water to obtain performance data to show the range of sizes that may be recovered and other operating characteristics that would indicate the field of usefulness of this device in coal preparation practice. This is a progress report only. It covers results obtained in a series of 17 tests of the operation of what has been considered a unit of commercial size operating on suspensions of fine coal mixed with a relatively small amount of freeimpurities. The underflow material, after further dewatering in centrifuges, is a marketable grade of coal.

Jan 1, 1949

By J. J. Arps

n approach to the study of the electrical properlies of drilling muds, their cakes and their filtrates was made by treating the drilling fluid as a porous medium. Lluring this study it was found that the weight pararneter was omitted from previous investigations on electrical properties. Patnode in 1949 made a similar approach to ,his investigation. His work showed that there was a variation of resistivity relationships with per cent Puid present in a given suspension. This present study ex-tends Patnode's investigation to one in terms of the solids fraction of the drilling fluids. Previous investigations were made on low weight muds where the effect of the rolids fraction was obscured. The following equationr were fitted to data taken from 94 field muds and 47 laboratory muds (see Tables I and 2). The parameter, C, is a complex function of the solidr content, or mud weight, of a drilling fluid. It was also found that for all drilling muds examined. A nomograph is presented combining Eqs. I and 2. INTRODUCTION It is generally recognized that the materials present in drilling muds affect the relationship between mud filtrate and mud resistivities. However, most electric log interpretations proceed with the assumption that the mud resistivity and temperature are the only important parameters needed to obtain a filtrate resistivity or a mud cake resistivity. Lamont1 suggests that only average relationships exist between mud cake, mud filtrate and mud resistivities. He indicates that there is limited accuracy when these intercorrelated factors are used since no other measure-able variable is apparent. A previous paper by Patnode indicates that drilling fluids may be considered to be a system of fluids in solid particles which may be correlated electrically in a manner similar to rock solids. 1 he equation for the resistivity factor. F = (per cent fluid) is satisfied for a given clay type in water, but the value of m varies as different materials are used in the water. The resistivity factor is not constant for a given amount of clay. As is suspected it varies as the solids fraction conducts electrical current. Since both weighting material and clay solids are present in commonly used drilling fluids, it is expected that both materials contribute to the variation of the m factor in Eq. 3. An extension of Patnode's work is presented in this paper to show some of the interrelations which exist between the resistivities of mud cake, mud filtrate and mud. Use of these data enable a more accurate calculation of mud filtrate and mud cake resistivities if only the mud resistivity and weight are measured. EXPERIMENTAL The measurements of mud weight, fluid and solids resistivity, alkalinity, and chloride content of mud samples were made in this research. Mud samples representing drilling fluids from 94 Louisiana and Texas wells (Table 1) and 47 laboratory-constructed systems (Table 2) were used in the experimentation. The measurements of the mud properties were made both in the field and in the laboratory. The tests on the properties of weight, alkalinity and chloride ion were standard tests as described by API specifications. The resistivity measuring devices were constructed in the laboratory for this study. For details on the resistivity equipment see Fig. 1. The fluid cell

By AIME AIME

This section contains announcements of general interest to the members of the Institute, but not always of sufficient permanent value to warrant republication in the volumes of the Transactions.

Sep 1, 1906

By A. Werenfels

No PETROLEUM, in commercial quantities, is found in France except the small production of the Pechelbronn field (Fig. 1) in Alsace, which produced 75,000 tons in 1923; therefore the test of the petroleum prospects of the Limagne is most interesting. Under the superintendence of M. Termier, chief engineer of the French mines, a well near Crouelle (Fig. 2) was drilled in 1921 to a depth of 856.3 m. (2808 ft.); at different depths favorable showings were encountered. The well, however, was abandoned because of drilling difficulties. A second one is now being drilled near the same place.

Jan 1, 1925

By E. L. Anders

The Mile Six pool is located on the La Brea-Parinas Cullcession of International Petroleum Co., Ltd., in northwestern Peru on the west coast of South America. The reservoir pressure in this pool has been maintained within 200 psi of its initial value throughout its history, and gravity drainage has played an important role in the production behavior. It has now produced 95 per cent of its estimated ultimate recovery. It is estimated that this interesting oil pool will ultimately produce 67 per cent of the initial oil in place and that the resulting residual oil saturation may be as low as 19 per cent of the pore volume (29 per cent of the hydrocarbon pore volume). An evaluation of reservoir rock and fluid characteristics and ultimate oil recovery is presented. INTRODUCTION This study of Mile Six pool was made to evaluate its performance according to latest available information. The production performance of this pool has been discussed in various articles in the past. and the reported behavior has been used as an example for application of computation procedures for gravity drainage depletion' and as an illustration of field behavior under gravity drainage or expanding gas cap drive. There have been wide variations in reported values of initial oil in place, reservoir oil volume factor. connate-water sauration, volume of effective sand, and ultimate recovery because of the paucity of reliable basic data. These various factors have been determined as accurately as practicable with the latest available information, and this evaluation is presented herein. The production history of Mile Six is an excellent example of gravity drainage depleion with effective pressure maintenance by gar injection. GENERAL Mile Six pool was discovered by cable-tool drilling in November. 1927. when well 1996 was completed in the Parinas sand. After slow development with cable tools and sporadic production. the pool was opened to continuous pro(iuction in November. 1933. and develpment was completed with rotary rigs. Pressure maintenance was started in December, 1933. by returning gas to upstructure wells. Most of the development was 'completed by 1937, but some additional wells were drilled in the period 1939-1947. and several old wells were deepened. A total of 46 oil and gas wells and 4 dry holes were drilled on approximately 7-acre spacing. Of the producers. 21 are now flowing. 2 are pumping. 44 are gas input wells. 3 are abandoned. I is a gas well shut in. and 15 are shut in because of non-commercial production or high gas-oil ratio. The locations of all wells are shown on the map of Fig. I. Total oil production on Dec. 31. 1952, was 30,867,373 bbl: cumulative gas production was 22,023,777 Mcf; and 26,410,946 Mcf of gas had been returned to the reservoir. These figures do not include oil and gas lost ill a blowout in January. 1940. GEOLOGICAL DESCRIPTION Mile Six pool is located on the northern end of a structural spur projecting from the La Brea-Negritos uplift.' The spur is probably a reflection of a basement structure. It plunges gently to the north, i broken into a complex series of fault blocks. and contain.; the Verdun Alto. Section Sixteen. and Mile Six pools. The Parinas handstone (lower Eocene!. which is the producing formation in Mile Six. occurs at an average depth of 2,200 ft in the pool and dips north and east at from 15' to 20". The pool covers an area of approximately 350 acres. Mile Six is downfaulted about 600 ft from Section Sixteen pool to the soutb. and a major fault forms its western boundarv. The north and east boundaries are formed by the intersection of the sand top with the water-oil contact which occurs at approximately 2.440 ft subsea. An original gas-oil contact probably existed at about 1.875 ft subsea. Fig. 1 presents the latest structural interpretation of the pool. and Fig. 2 is an isopach map showing thickness of the total Parinas formation above the original water-oil contact. The heavy lints of Fig. 1 are contours on the sand top, and the fine lines are contours on the fault planes. This type of straight-line structural map was developed 1)) International's geologists to reflect structural conditions where the bedding planes dip and have no curvature. 'The La Brea-Parinas Concession is highly faulted by normal fault.. The beds are flat wherever exposed. The Parinas formation is approximately 635 ft thick. and it is etimated that 62 percent II the formation is effective sand. The original oil zone was about 565 ft thick. Fig. 211 presents an electric log showing typical Parinas sand development ill Mile Six pool. The Parinas band in Mile Six i. a well-sorted. medium- to-coal-se-grained. cross-bedded sand with minor lenses of shale and small lenses and pockets of pebble conglomerate. The sand grains are subangular to rounded and consist chiefly of quartz with feldspars. biotite hornblende, and augite as accessory minerals. Because of faulting of the Parinas- formation to the east and north of the pool. there is probably l possibility of a significant, natural water drive in Mile six. The faults within the pool. as indicated in Figs. 1 and 2. are of smaller disI,laceInent and seem to act only , partial barriers to fluid movement within the reservoir. RESERVOIR CHARACTERISTICS Core analysis data are available from five wells. The data were obtained from three wells (Nos. 3401. 3586 and 3719) at the time of their completion and from well 1996 when the original liner was sidetracked and the well was deepened ill 1946. Data from well 2779 were obtained in 1943 from old cores taken when the well was deepened in 1934. From these core analyses. the average porosity was estimated to be 22.6 per cent. and the average permeability to dry ail. was estimated to be 780 and Measured productility indices varied from 3.1 to 71.4 B/D per psi differential. Specific productivity indices varied approximately from 0.1 to 0.3 B./D per psi per ft of sand.

Jan 1, 1953

By Claus G. Goetzel

SOME of the reactions and procedures upon which modern techniques in the production of metal powders are based were used for 2000 years by the ancients to reduce iron and other metals from their ores. Despite this, however, the art of powder metallurgy is relatively new and probably stems from the introduction of incandescent electric lamp filaments. Attempts to improve upon Edison's carbon filaments led to the production, in order. of metallic filaments of osmium. tantalum.. and finally of tungsten, which ha, replaced the others and is the commercial product of today. These metallic filaments were originally manufactured by the extrusion and sintering of a mixture of the metal or its oxide with a reducing agent which also served as a binder. Further progress in this direction was in the pressing of a small percentage of nickel with the tungsten. The sintering operation was carried out in an atmosphere of hydrogen, at a tem-

Jan 1, 1944

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

The nature and the limits of validity of Darcy's law US applied to the flow of natural gas through reservoirs has been considered in order to resolve some controversial aspects of the effect of turbulence on pressure drops. The equivalence between various concepts and viewpoints advanced in the past by several investigators to explain how and why a gas well does not necessarily perform according to Darcy's law is shown. Starting with generalized equations of flow of fluids through porous media, a partial differential equation has been derived which accurately represents the flow at all rates. This equation has been numerically solved using an IBM 704 digital computer. The results permit plots of unsteady radial pressure distribution curves from which specific isochronal backpressure curves may be constructed. These back-pressure curves show the effect of the ,8 factor on the slope of the back-pressure curve. The calculations further indicate that the drainage radius for a gas well in turbulent flow propagates at a rate dependent upon the rate of production at the wellbore. This is quite different from the case with liquid flow or natural-gas flow in laminar regime. Additionally, the effect of reservoir inhomogeneities and crossflow between layers of different permeability on the back-pressure performance of gas wells has been conridered. In light of the current numerical results the significance and linlitation of the rate of flow function Y proposed by Smith" has been discussed. INTRODUCTION The relationship between the pressure drop and flow rate in problems of fluid flow through porous media is known to be affected by the nature of flow through the porous matrix. It has been observed by many that, for a range of flow rates, the pressure drop remains proportional to the rate of flow. When some flow rate is reached, however, it is usually observed that the pressure drop gradually begins to increase more than proportionally to the flow rate. It is well known that this phenomenon was first observed by Osborne Reynolds in 1901 in experimenting with flow through pipes. In his classical experiments, Reynolds made visual observations on the condition of streamlines evidenced by injecting a dye into water flowing through glass tubes. In these experiments, the abrupt transition between steady, "streamline, laminar" flow and unsteady random turbulent flow was found to be a function of the dimensionless group (Dvp/u), now known as the Reynolds number. During these experiments, in addition to observations on the nature of flow regimes, the proportionality between flow rate and pressure drop in laminar flow was contrasted with the nonlinearity between these variables in turbulent flow. Fancher and Lewis' reported data on various consolidated and unconsolidated sands in 1933. Their conclusions were that ". . . the flow of fluids through these porous materials closely resembles that through pipes; that there is a condition of flow in porous systems which resembles viscous flow, another which corresponds to turbulent; that the change from one type to the other takes place at a definite and reproducible condition for each system". In 1947, Brownell and Katz published a method to predict the laminar and turbulent flow behavior from the particle size, bed porosity and the particle sphericity, employing the friction factor-Reynolds number charts for pipes. Several investigators have verified the work of Fancher, Lewis and Barnes and presented their data as friction factor-vs-Reynolds number plots.' The equation which would represent the pressure gradient over the whole range of velocity must have an added term over that represented by Darcy's law. Accordingly, the pressure gradient necessary to sustain flow at the velocity (v) through a porous medium may be represented by the following equation, suggested by Forscheimer." The nature and the range of validity of Darcy's law has been the subject of studies by many investigators over the past years. While everyone seemed to agree on the need for a quadratic correction term to Darcy's law to make it effective over the range of velocities, the concept of inception of turbulence and the use of the term "turbulent flow" remained controversial. Some fluid dynamicists define turbulence as a flow

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. S. Odeh

This paper analyzes the effect of limited entry to flow at the wellbore on the steady-state productivity of a well. Wells that have been opened to flow along a fraction of their productive interval are termed wells with limited entry. Previous work treated the cases of a partially penetrating well, a well producing from the central portion of the productive interval and a well in which several intervals equally spaced were open to flow. In this paper the open interval can be located anywhere within the productive interval. Thus, in a sense, it generalizes previous work. The finite cosine transform was used to arrive at a solution for steady-state flow of a slightly compressible fluid. The solution was programmed for a CDC 1604 computer. Numerical vaIues for rd = 660 ft, r, = 1/4 ft, and range of sand thickness of 20 to 200 ft are presented in graphical form. The effect of rd and r, values on the result is shown in a table. The correct calculation of skin and damage ratio in the presence of limited entry to flow is explained and illustrated by examples. Moreover, the paper shows how to calculate the net decrease in productivity due to the combined effect of limited entry and perforations. INTRODUCTION In some wells only a fraction of the productive interval is open to flow. Location of this fraction is usually dictated by formation characteristics and reservoir behavior. For instance, if a gas cap exists, the open interval is located away from the gas-oil contact to prevent any possible gas coning. Wells that intentionally have been opened to flow along a fraction of their productive formation are tened wells with limited entry. Obviously, unintentional completions of this type also exist. Limited entry to flow decreases well productivity. Magnitude of the loss depends on the fraction of the formation open to flow, on the thickness of the sand, on the location of the open interval and on the ratio of rd /r, where r , is well radius and rd is the drainage radius of the well. The use of pressure buildup data on producing wells to calculate the condition of the formation around the wellbore is an accepted practice. van Everdingenl and Hurst2 introduced the concept of she skin factor s considered to be due to a thin layer of different permeability immediately around the wellbore. These authors dealt with the case of a well of complete radial geometry, i.e., a well with open-hole completion that completely penetrates the formation. The presence of a low-permeability skin results in a loss of productivity, as does limited entry. Therefore, if pressure buildup data obtained on a well with limited entry are used to establish the presence or absence of skin (i.e., formation damage), and a correction is not made for this loss of productivity, the calculations would result in an erroneous skin value. They might indicate the presence of formation damage when in reality there is none, or they might indicate a value larger than the true value. This could lead to an incorrect basis for planning remedial measures. Muskat3 studied the problem of partially penetrating wells for the case of incompressible flow. He presented equations and figures which allow estimation of loss in productivity. Brons and Marting,4 using equations based on Nisle's work,5 studied the loss of productivity for three cases. The first was for a partially penetrating well; the second was for a well producing from only the central portion of a productive interval; and the third was for a well in which several intervals equally spaced were open to flow. Their work was for steady-state depletion-type reservoirs wherein the well radius of drainage is established and the fluid is considered to be slightly compressible. Considered in this paper is the problem of wells with limited entry in which the open intervals are located anywhere within the productive sand. The finite cosine transform is used to arrive at a

Jan 1, 1969

By L. A. Wilson, N. H. Harrison, R. R. Clay, R. L. Reed. D. W. Reed

Experiments have been performed in a linear near-adiabatic system for the purpose of extending data on reverse and forward combustion from atmospheric pressure to 1,000 psig. Results obtained from reverse combustion appear to conform qualitatively with the existing description of the process when reaction kinetics are suitably modified to account for the increase in pressure, i.e., increasing tbe pressure decreases the peak temperature and increases the combustion zone velocity. Forward combustion appears to be a fuel dominated process wherein peak temperature and combustion zone velocity are not very sensitive to changes in pressure. The moderate effects of pressure that do exist at low flux virtually disappear at high flux providing all oxygen is consumed. With this provision. increasing the pressure decreases the frontal velocity and increases the peak temperature. Results me shown graphicaIIy which demonstrate the effects of presswe on peak temperature, rate of advance, oil recovery, air/oil ratio, carbon oxides produced and temperature distributions. INTRODUCTION Numerous field tests of oil recovery using the technique of underground combustion are now in progress. Operating pressures used are always substantially in excess of atmospheric pressure. Nevertheless, the literature contains no laboratory data pertaining to the effects of pressure except for those of Martin, Alexander and Dew.1 Unfortunately, these data appear to reflect heat losses which may have obscured some of the effects of increased pressure. The underground combustion processes are exceedingly complex, and general concepts relating to them have necessarily been described in relatively simple terms. Particularly with regard to forward combustion, the interaction between multiphase fluid flow, heat transport and various rate mechanisms is so involved that to date a rigorous treatment is not available. In this paper a linear near-adiabatic system is described which was found capable of operating at pressures up to 1,000 psig and temperatures of at least 1,100 F. Using this equipment, experiments on both forward and reverse combustion were accomplished at elevated pressures and the data presented as a function of air flux or temperature. It will be assumed that the reader is familiar with the papers of Wilson, Wygal, Reed, Gergins and Henderson on forward combustion and that of Reed, Reed and Tracht3 on reverse combustion, so that it is not necessary to repeat the details of the two processes. EQUIPMENT AND PROCEDURE It was decided to study the effects of pressure over the interval 0-1,000 psig using a near-adiabatic system. Since satisfactory operation of such a system requires very thin tubing walls, it proved necessary to insert the entire combustion tube and heater assembly into a high pressure jacket in such a way that the internal tube pressure could be counterbalanced by pressurizing the annulus. This meant that the numerous heater and thermocouple wires had to pass through gas-tight seals that required considerable effort to maintain effective. An unexpected difficulty arose in connection with the nature of the annulus pressuring gas. It turned out that the composition of this gas critically affected the functioning of the adiabatic control system. Operation of the equipment at elevated pressures was not particularly difficult. The problems that did arise were a direct consequence of the large numbers of connections and seals involved. EQUIPMENT The equipment consisted of the combustion tube which fitted within the pressure jacket. This jacket was in turn mounted on a wheeled carriage which facilitated movement and provided for rotation of the

By J. S. Lasater

Resu1ts of experinmental measurernents of heat capacities and thermal conductivities of some typical porous rocks are presented. Measured heat capacities agree closely with va1ues calculated front known chenlical . compositions of the rocks. On the basis of this agreement, heat capacitir.7 of fluid-satttrated rocks were calc~ililted. Therrnal conductivities of the rock samples were measured under varioits conditions of fluid saturation. From these data therrnal dif-fvsivities were calculated. A significant variation of thermal diflusivity with trnlperature is indicated. Knowledge of the thermal characteristics of fluid-bearing porous rocks has become increasingly important with the development of thermal oil-recovery processes. Reliable thermal data for petroleum reservoir rocks are not available and approximate values are generally assumed for reservoir calculations. 'The effects of temperature, pressure and fluid saturation on thermal properties of porous rocks have not been carefully studied. The purpose of the present work was to measure and report the thermal characteristics of some typical porous rocks. A group of eight sedimentary rock samples, including limestone, shale, sandstone and silt-stone, was selected for the study. Properties reported include heat capacities, thermal conductivities and thermal diffusivities. The effect of fluid saturation on these properties was also considered. Determinations of thermal properties require precise measurements which are time consuming and difficult to reduce to a routine basis. Methods have not been developed for measurement of thermal proper-ties under conditions of pressure, temperature, and fluid saturation, which may be encountered in thermal recovery processes. As one approach to these problems, some preliminary work on methods of estimating thermal properties from other more readily measurable characteristics of the rock-fluid system is reported. The samples selected for the studies were oilfield cores from within, or closely associated with oil-productive zones. Descriptions of individual samples are given in Table 1. Clay minerals present in the samples were determined by differential thermal analysis (DTA). Quartz content was estimated from DTA cooling curves. For all but the thermal conductivity tests, a representative 200-gm sample of each core was disaggregated, thoroughly extracted in diox-ane and oven dried at 100°C. The cleaned and dried sample was care- fully split into appropriate portions for the several tests. Chemical analyses of the samples were obtained to permit calculations of heat capacities by methods which are outlined later. Resuits of chenlical analyses are given in Table 2. Heat contents of the test samples were measured by a Runsen-type calorimeter. The heat content of the unknown sample is measured relative to the known heat content of platinum. The apparatus is referred to as a relative error type because of the comparative measurements, and has the advantage of canceling systematic errors. Experimental results are reproducible within 0.5 per cent. A temperature base of 298.16 K is taken for heat content measurements. Measurements were made at five temperatures (260, 440, 620, 800 and 980° F). It was necessary to preheal the sample above the maximum experimental temperature to eliminate volatile constituents. Assuming the weight loss on pre-heat, ing was mostly water, a correction to the observed heat content was made to account for the water originally present in the sample. This latter procedure was followed in the determination of calculated heat contents as will be described later. The apparatus and the method of measure-

By G. Thodos, R. J. Smialek

A method has been developed that permits the correlations of Mehra and Thodos8 and Dastur and thodos2 for binary systems to be used for the prediction of the vapor-liquid equilibrium constants of the light and heavy constituents, Kland Kh, of a ternary hydrocarbon mixture. Furthermore, the unique relationship observed in this study, i.e., a plot of log Ke/Ki vs log Kp/K h is linear and passes through the origin, permits the calculation of the K-value for the intermediate component. This approach has been used to establish K-values for the constituents of the three ternary hydrocarbon systems used to develop this method and has also been applied to two additional ternary hydrocarbon systems to yield an average deviation of 7.7 per cent resulting from a total of 294 points for the five ternary systems considered. INTRODUCTION The rigorous calculation of vapor-liquid equilibrium constants for systems containing two or more components is a difficult problem and becomes extremely complex at elevated pressures and for conditions approaching the critical point. Although thermodynamic relationships can be developed to account for the vapor-liquid equilibrium behavior of simple systems, their use is limited because the necessary PVT data for these systems are not generally available. To establish a method for calculating K constants directly, without any recourse to the PVT behavior of a system, Mehra and Thodos8 applied the principle of corresponding states to the vapor-liquid equilibrium constants of binary hydrocarbon systems. They developed a series of charts which made possible the prediction of K-values for all binary hydrocarbon systems not containing methane. This work has been extended by Dastur and Thodos2 to include binary hydrocarbon systems containing methane. The charts developed in these two studies cover the high-pressure region and are applicable to conditions in the vicinity of the critical point. The ability of this approach to predict vapor-liquid equilibrium constants of binary hydrocarbon systems suggests that this method may be extended to the prediction of the K-values of the components of ternary hydrocarbon systems. In this investigation, the charts developed for both the methane and non - methane binary systems,2,8 together with mixture rules developed in this study, have been used to predict the K-values of the light and heavy components of a ternary hydrocarbon system. These calculated values are then used to establish the K-value of the intermediate component through the use of a novel relationship developed in the present study. DEVELOPMENT OF METHOD FOR PREDICTING K-VALUES OF TERNARY SYSTEMS The application of the corresponding state charts developed for binary hydrocarbon systems is restricted to those systems containing two components. The extension of these correlations to three-component systems requires that the ternary system be reduced to an equivalent binary, which can be treated in a manner consistent with the corresponding state development. To establish such an approach, the vapor-liquid equilibrium constant of the light component K is first calculated by assuming that the ternary mixture consists of the light component and a pseudo-heavy component of the intermediate and heavy constituents of the original ternary. To establish the vapor-liquid equilibrium constant of the heavy component Kh, the light and intermediate constituents of the ternary are combined to form a pseudo-light component of another equivalent binary system. The calculated values Kl and Kb are then utilized to calculate the vapor-liquid equilibrium constant of the intermediate component EQUILIBRIUM CONSTANT FOR LIGHT COMPONENT K To calculate the K value of the light (most

Jan 1, 1965

By S. A. Hovanessian

Analytical solutions are obtained for calculating the pressure distribution in rectangular fields due to injection and/or producing wells located anywhere within the field. The field is assumed to be homogeneous with either constant pressure or no-flow boundaries. These solutions are extended to include the calculation of average reservoir pressure and permeability from pressure build-up curves. Numerical examples are given to illustrate the application of equations to practical problems. INTRODUCTION Pressure distribution analyses are of considerable importance in the field of reservoir mechanics. A number of papers dealing with this subject have appeared in the petroleum literature. perrinel in 1956 presented a summary of the methods available for pressure build-up calculations and discussed the applicability of each in some detail. More recently (1958), Hazebroek, et a1,2 discussed a theoretical method for obtaining pressure distribution in a radial field. Also in 1958, Nisle3 gave a theoretical solution of the pressure distribution in a field extending to infinity with a partially penetrating line source at the origin. Mathews, et a1,4 applied the method of images to existing solutions for pressure distribution in radial fields and obtained expressions for calculating pressures in bounded reservoirs. However, a literature survey revealed that a general analytical solution for the pressure distribution in a rectangular field, with injection and/or production wells located anywhere within the field, was not available. Such a solution could be used to approximate the pressure distribution in a field the shape of which is nearly rectangular. Furthermore, it could be applied to pressure calculations pertaining to a single well taken from a system of wells which were drilled in a rectangular pattern. In this case the boundaries of the well drainage area are considered to be the perpendicular bisectors of the lines joining the well and its nearest neighbors. The basic problem considered in this paper consists of a rectangular field with either constant-pressure 01 no-flow boundaries and with either an injection or production well located anywhere within the field. The physical properties of the rock and fluids present are considered to be constant. Analytical expressions are obtained for the pressure distribution within the field during the production or injection periods. These expressions, which are infinite series, can be used to evaluate pressure at any point within the field as a function of time and position. For fields containing more than one well (production and/or injection), the solutions for each well acting independently are superposed. That is, at any point in the field, the pressures due to each injection 01. production well acting alone can be algebraically added to give the pressure resulting from all wells acting simultaneously. Evaluation of the series solutions given in this paper is readily accomplished by means of an electronic digital computer; and, in most cases, sufficient accuracy is obtained by setting the upper limits of summation at 20. These solutions are extended to include the calculation of average reservoir pressure and permeability from build-up data. A number of numerical examples are given to illustrate application of the solutions to actual field problems. ANALYTICAL SOLUTIONS The partial differential equation representing the pressure distribution in a homogeneous rectangular field (see Fig. I.), having a single source or sink of strength (2, can be written*

By C. D. King

A QUARTER century ago this country was producing an extraordinary quantity of iron and steel, with a decisive influence on the outcome of the first World War. Today this country is again demonstrating its industrial capacity with telling results. If this lapse of time has confirmed the opinion that man learns little from experience with respect to human relationships, at least we can derive some consolation from the knowledge that in the same period we have made huge strides industrially, and particularly in the steel industry. The differences represented by 25 years are striking, and have resulted from a multitude of changes that have contributed to the present-day high state of our pig iron and steel manufacture.

Jan 1, 1944