Search Documents

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 J. C. Deppe

A method is presented for calculating approximate injection rates in secondary recovery operations. The method can be applied to cases of unequal fluid mobilities, irregular well patterns and boundary patterns. The steady-state pressure distributions for the four flood patterns reported by Muskat and for five additional patterns reveal that most of the difference in pressure between the injection and producing wells occurs in regions around the wells which can adequately be described as regions of radial flow. This leads to a method of calculating injectivity by approximating the flood pattern with radial flow elements (or a combination of radial- and linear-flow elements for some patterns such as the direct line drive). Irregular and boundary patterns can also be approximated by radial and linear elements. Each of these elements can be described by radial-and linear-flow equations and the results combined as series flow resistances to give an approximate equation for the initial injection rate. The mobility ratio does not affect the initial rate; therefore, if the well pattern is one of those regular patterns for which theoretical rate equations have been derived for unity mobility ratio, the approximate initial rate equation can be improved by adjusting it to match the theoretical equation. The available theoretical rate equations are listed, including five new cases. As the flood progresses, the injectivity changes because in general the flood front will divide the pattern into areas of different fluid mobilities. Simple shapes can be assumed for the flood front so that both areas can be divided into radial- and linear-flow elements. Radial- and linear-flow equations are applied to these elements to account for the change in flow resistance behind the front. To calculate injection rates after breakthrough, it is necessary to know the sweep efficiency at breakthrough. Areal sweep data are available in the literature for a number of patterns, and sufficiently accurate breakthrough sweep efficiency can be estimated from these data if it is not otherwise available. Again, simple shapes for the flood front can be assumed after breakthrough, and injection rates can be calculated for the remainder of the flood. When areal sweep data are available, these data can be used to cbeck. the calculated injection rates after breakthrough. Injection rates calculated by this approximate method compare favorably with available model data. INTRODUCTION When the fluid mobilities in the swept and un swept regions are equal, the injectivity will not change as the flood front advances; and for regular patterns it can be calculated by mathematical formu1as.l When the fluid mobilities in the swept and unswept regions are not equal, the injectivity will increase or decrease as the flood front advances. In this case, the injectivity has not been calculated by analytical means for any practical well pattern; and, furthermore, scale model and analog4,6 results have been published only for the five-spot pattern. The main object of this paper is to present an approximace method for calculating injectivity for the case of unequal mobilities. The method can be applied to regular, irregular and boundary patterns. Before the approximate method is discussed, the analytical formulas for equal mobilities will be summarized, and the analytical solution for radial flow with unequal mobilities will be used to show how the injectivity changes as the flood progresses. EQUAL MOBILITIES, REGULAR PATTERNS In a flooding operation the injectivity, the rate at which fluid can be injected per unit difference in pressure between injection and producing wells, depends on basic physical properties and, in addition, on three variables to be treated here—the area of the swept region, the fluid mobilities in the swept and unswept regions, and the well geometry (pattern,

By K. E. Brown, G. H. Fancher

An 8,000-ft experimental field well was utilized to conduct flowing pressure gradient tests under conditions of continuous, multiphase flow through 2 3/8-in. OD tubing. The well was equipped with 10 gas-lift valves and 10 Maihak electronic pressure recorders, as well as instruments to accurately measure the surface pressure, temperature, volume of injected gas and fluid production. These tests were conducted for flow rates ranging from 75 to 936 B/D at various gas-liquid ratios from 105 to 9,433 scf/bbl. An expanding-orifice gas-lift valve allowed each flow rate to be produced with a range of controlled gas-liquid ratios. From these data an accurate pressure traverse has been constructed for various flow rates and for various gas-liquid ratios. A comparison of these tests to Poettmann and Carpenter's correlation indicates that deviations occur for certain ranges of flow rates and gas-liquid ratios. Numerous curves are presented illustrating the comparison of this correlation with the field data. Poettmann and Carpenter's correlation deviates some for low flow rates and, in particular, for gas-liquid ratios in excess of 3,000 scf/bbl. These deviations me believed to be mainly due to the friction-factor correlation. However, Poettmann and Carpenter's correlation gives excellent agreement in those ranges of higher density. This was as expected and predicted by Poettmann. He pointed out that their method was not intended to be extended to those ranges of low densities whereby an extreme reversal in curvature occurs.l As a result of these experimental tests, correlations using Poettmann and Carpenter's method were established between the friction factors and mass flow rates which are applicable for all gas-liquid ratios and flow rates. Definite changing flow patterns do not allow any one correlation to be accurate for all ranges of flow. INTRODUCTION The ability to analytically predict the pressure at any point in a flow string is essential in determining optimum production string dimensions and in the design of gas-lift installations. This information is also invaluable in predicting bottom-hole pressures in flowing wells. Although this problem is not new to industry, it has by no means been solved completely for all types of flow conditions. Versluys,2 Uren, et al, 3 Gosline,4 May,5 and Moore, et al,6,7 were all early investigators of multiphase flow through vertical conduits. However, all of these investigations and proposed methods were very limited as to their range of application. Likewise, many are extremely complicated and therefore not very useful in the field. Only in the last decade have any significant methods been proposed which are generally applicable. The most widely accepted procedure in industry at the present time is a semi-empirical method developed from an energy balance, proposed by Poettmann and Carpenter8 in 1952. Their correlation is based on actual pressure measurements from field wells. Accurate predictions from this correlation are limited to high flow rates and low gas-liquid ratios. Although this method will be discussed in detail later, it should be pointed out that two important parameters, namely the gas-liquid ratio and the viscosity, were omitted in their correlation. The viscosity was justifiably omitted since their data was in the highly turbulent flow region for both phases, and most wells fall in this category. The gas-liquid ratio was incorporated to some extent in the gas-density term. In 1954, Gilbert9 presented numerous pressure gradient curves obtained from field data for various flow rates and gas-liquid ratios for the determination of optimum flow strings. However, no method is presented for predicting pressure gradients except by comparison to these curves. Baxendell, et al,10 proposed a method based on Poettmann and Carpenter's procedure

By A. J. Conelius, D. P. Sobocinski

This paper presents a correlation for predicting the behavior of a water cone as it builds from the static water-oil contact to breakthrough conditions. The correlation is partly empirical and involves dimensionless groups of reservoir and fluid properties and of production and well characteristics. The groups were deduced from the scaling criteria for the immiscible displacement of oil by water. The correlation is based on a limited amount of experimenfal data from a laboratory sand-packed model and on results from a computer program for a two-dimensional, incompressible system. Because the correlating groups are dimensionless, they can be used to estimate the performance of water coning cases not specifically considered in the correlation. However, despite its dimensionless nature, the correlation is not completely general and will not provide meaningful estimates of cone behavior in many situations. INTRODUCTION AND BACKGROUND The production of water from oil wells is a common occurrence which increases the cost of producing operations and may reduce the efficiency of the depletion mechanism and the recovery of reserves. We will deal with one cause of this water production, namely, coning. The coning of water into a producing well is caused by pressure gradients established around the wellbore by the production of fluids from the well. These pressure gradients can raise the water-oil contact near the well where gradients are most severe. Gravity forces that arise from fluid density differences counterbalance the flowing pressure gradients and tend to keep the water out of the oil zone. Therefore, at any given time, there is a balance between gravitational and viscous forces at points on and away from the completion interval. When the dynamic forces at the wellbore exceed gravitational forces, a cone of water will ultimately break into the well to produce water along with the oil. We can expand on this basic visualization of coning by introducing the concepts of stable cone, unstable cone and critical production rate. For instance, if a well is produced at a constant rate and the pressure gradients in the drainage system have become constant, a steady-state condition is reached. If, at this condition, the dynamic forces at the well are less than the gravity forces, then the water or gas cone that has formed will not extend to the well. Moreover, the cone will neither advance nor recede, thus establishing what is known as a stable cone. Conversely, if the pressure in the system is in an unsteady-state condition, then an unstable cone will continue to advance until steady-state conditions prevail. If the flowing pressure drop at the well is suffcient to overcome the gravity forces, the unstable cone will grow and ultimately break into the well. It is important to note that in a realistic sense, stable cones may only be "pseudostable" because the drainage system and pressure distribution generally change. For example, with reservoir depletion, the water-oil contact may advance toward the completion interval, thereby increasing chances for coning. As another example, reduced productivity due to well damage requires a corresponding increase in the flowing pressure drop LO maintain a given production rate. This increase in pressure drop may force an otherwise stable cone into a well. The critical production rate, well known in the literature, is the rate above which the flowing pressure gradient at the well causes water (or gas) to cone into the well. It is, therefore, the maximum rate of oil production without concurrent production of the displacing phase by coning. At the critical rate, a built-up cone is stable but is at a position of incipient breakthrough. Numerous papers have been published on critical rates. Some of the better known of these include the work of (1) Muskat and Wyckoff,' who first dealt with the coning problem; (2) Chaney, et al,' who developed expressions similar to those of Muskat but who presented results in a conven ient-to-use graphical form (the "Sun" method); and (3) Meyer and Garder,b hose analysis is based on radial-flow formulas. One assumption in critical production rate analyses is that the cone has built-up to just before its breakthrough into the well. But, these analyses reveal nothing directly about the time it takes for the cone to build up to this incipient breakthrough position. Thus, water-free oil can be produced from a well for prolonged periods at rates above the critical rate before the well reaches the condition to which the critical rate applies. The published literature contains little on the rate of growth of a cone. Experimentally, Meyer and Searcy studied the rate of rise of a cone in a Hele-Shaw model.' Additional related work on water breakthrough and produced water-oil ratios in water driven reservoirs was reported by Muskat," Hutchinson and Kemp," Henley, et al,' and Stevens, et a1.B Theoretically, the basic coning equations for a water-oil system can be developed by applying the conservation of mass to each of the phases, relating flow velocities with pressure by Darcy's law, and relating pressures across water-oil interfaces by capillary pressure. With the usual boundaries at the well and reservoir limits, the solution of the resulting equations for the time behavior of a water-oil interface constitutes a free-surface, boundary-value problem. The shape of the boundary depends on the internal potential distribution, while this

Jan 1, 1966

By M. P. Tixier, F. Segesman

Explanations bused on theory and laboratory tests tire offered for certain unusual features of SP logs which may befog their interpretation, namely the "sawtooth" SP and the reduction or increase of SP with time. The sawtooth SP is frequently seen in SP recordings of thick, very permeable watcr. sands of the Louisiana Gulf Coast. Using SP data from a resistor network analog, this anomaly is explained on the basis of a thin shale streak within the sand. under which the invading filtrate accumulates by gravity segregcitiorr. The electromotive force then appearing across the shale streak gives rise to the anomalous SP. Reduction of the SP on successive logging runs Ira often been observed in very permeable sands when as indicated by the resistivity curves, the invasion depth has been reduced to practically zero. Laboratory results verify that this phenomenon may be due, to the presence of air electromotive force, across the mod rake when the fluids on either side differ in activity. Also, there is, in some cases; an increase of SP deflection on successive logging runs, because the emf initially present across the mud cake disappears as the somewhat saline water just behind the mud cake is progressively flushed out by the invading filtrate.. Field example and laboratory data are presented. I N T R ODUCT1ON Some unusual features are frequently observed on certain SP logs which make interpretation difficult. The purpose of this paper is to offer explanations of several cases of interest. All of these arc related to the conditions of invasion. The cases to be discussed are as follows. 1. The sawtooth SP anomaly which is frequently seen within the thick, highly permeable saltwater sands of the Louisiuana Gulf Coast 1. Reduction of the SP deflection from the usual sand line in very permeable sections when the depth of invasion is practically nil. As will be explained, such a reduction is attributable to a membrane electromotive force across the mud cake. 3. An increase of the SP with time just after the iormation is drilled as a result of progressive flushing of the invaded zone. FUNDAMENTALS OF SP PHENOMENA An explanation of SP phenomena must begin with the following basic ideas. 1. The SP log derives its character from ohmic potential differences in the hole produced by currents flowing through the mud.' These currents are driven by electromotive forces which are primarily electrochemical in origin. Other possible types of SP-producing electromotive forces may be disregarded for the purposes of this paper. 2. The electrochemical emf's associated with permeable formations are now understood to consist of liquid-junction emf's and membrane emf's.',' As schematically shown by thc battery symbols in Fig. 1, the liquid-junction emf, E,,,, exists at the contact between waters of different chemical activities, a,,,, (mud filtrate) and tr,,. (formation water); the shale membrane emf, lS.sn. exists across the adjacent shales whcn they separate waters of different activities. The emf's just cited are additivc in producing SP currents (arrowed circles of Fig. I). Thus, neglecting for the moment any possible membrane emf across thc mud cake, the total electrochemical emf, E,, will be given by their sum, where K,., -! K,q,, = 12 + 59, or 71 at 75 °F. In the log interpretations the SP of a thick sand (i.e.. the static SP) is generally related to the activities of the filtrate and formation water through Eq. 1. by setting the static SP equal to - E,.

By K. Harms, T. D. Jones

TO unload large tonnages of lead bullion cast in 100-lb. bars is a problem which has confronted the lead refineries for many years. The bars, on arrival, must be restacked for unloading by truck or hand-charged to charging machines. Selected crews of labor, capable of working at sustained speeds for five- or six-hour periods, are required for this exhausting task. Due to the fatigue involved crew personnel changes rapidly. During the war it became increasingly difficult to maintain crews, especially in critical labor areas, so that the development of casting and shipping ten-ton lead blocks from the smelters to the refineries became a must in our operating program. For many years lead refineries, when making interdepartment or intradepartment transfers, have cast lead blocks weighing five or ten tons which are easily handled by crane by merely attaching the crane hook to a steel anchor embedded in the frozen lead blocks. In previous years, the practice was to use regular smelter slag pots which, due to their rounded shape, produced what were called "buttons" by the refinery personnel. The button-type blocks were unsafe as they rolled and tipped easily and were unsatisfactory for stacking since no piling pattern could be followed. To remedy this situation, flat-bottom molds came into general use and later the bottoms were made concave to prevent slipping.

Jan 1, 1946

By D. H. Henley, F. F. Craig, W. W. Owens

The oil recovery performance of systems producing entirely by bottom-water encroachment has been experimentally determined in a series of scaled laboratory-model tests. The effects of well spacing, fluid mobilities, rate of production, capillary and gravity forces, well penetration and well completion techniques on the oil recovery performance have been investigated. The laboratory tests were performed using two uniform, un-consolidated sand-pack models. The models have ratios of the interwell distance to the formation thickness of 12 and 2, respectively. Tests at constant total fluid production rate were performed simulating a range of uniform reservoir characteristics and operating conditions encountered in field operations. The performance was determined by material balance and by observation of the encroachment of dyed fluids into the models. The results of the model tests agreed with those obtained mathematically when the conditions previously considered in theoretical studies were simulated, that is, when the oil and water are of equal density and no capillary forces exist. The model study of bottom-water drive indicated that certain variables can aoect the oil recovery performance to a greater degree than can be predicted by present analytical methods. In one comparison, the oil recovery at a water-oil ratio of 20 (obtained at a wide well spacing) varied as much as threefold, depending upon the system's properties and the production rate. Lesser effect of mobility ratio and no eflect of capillary forces over the range studied were observed. The test 'results also showed that the deeper the well penetration into the oil column, the greater the total water production to a producing WOR of 20. However, the ultimate sweep efficiency, and so the oil recovery to this level of WOR, did not vary significantly with well penetration. Horizontal fractures at the top of the formation did not significantly change the sweep characteristics of the reservoir models when values of radius and fracture capacity encountered in actual reservoirs were used. Impermeable pancakes at the bottom of the well moderately increased the oil recovery efficiency both at water breakthrough and at high water-oil ratios. A method is outlined by which the oil recovery performance of other uniform bottom-water drive systems can be estimated from the information obtained in these model tests. INTRODUCTION When oil is produced from a well which partially penetrates an oil zone completely underlain by water, the water rises directly beneath the well in a symmetrical cone when the system is uniform. Two different flow mechanisms can cause the water cone to form—coning and bottom-water drive. In coning, the aquifer is relatively inactive and the cone is formed beneath the well by the pressure gradients associated with the oil flow to the well. The oil can be produced by a solution-gas drive, an edge-water drive or other driving forces in the interwell area. In a bottom-water drive, the driving force for oil production comes from an upward encroachment of the underlying active aquifer. Two papers have analyzed the theoretical performance characteristics of bottom-water drive reservoirs. In the initial mathematical investigation, Muskat' established the equations which determine the pressure distribution in this type of reservoir and solved these equations for certain conditions. Specifically, it was assumed that the water and oil had equal mobilities and equal densities, there were no capillary forces, the pressure throughout the oil zone remained above the bubble-point pressure, a constant pressure existed at the initial water-oil contact and the oil was completely displaced by the encroaching water. These assumptions were used in obtaining analytical solutions. In general, Muskat found that the sweep efficiency to initial water breakthrough to the well was larger for the thicker oil zones, the closer well spacings, the lower ratios of vertical to horizontal permeabilities, the smaller the penetration of the well into the oil zone and the smaller the bore size of the well. The production history after water breakthrough was expressed as a volumetric sweep efficiency at a given producing water-oil ratio. The results indicated that cumulative oil production at producing water-oil ratios of 10 is less affected by the well spacing than is the water-free production history. Muskat studied well spacing which today would be regarded as close. The maximum value of his dimensionless well spacing (ratio of interwell distance to formation thickness) was 4.3. This would require the development of a 50-ft-thick oil sand on less than 10-acre spacing if the vertical and horizontal permeabilities were equal, with

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

By R. C. Hepworth

Heaving in roads built on the Mancos Shale is a serious problem in east-central Utah. Heaves are expressed as local bumps, general rises and severe cracking of the pavement. The heaves usually occur in shale road cuts and generally show a seasonal cycle, which reaches a maximum development in the winter. The presence of certain salts, absorption of moisture and physical properties of the subgrade materials were considered, in Utah, to be the main factors in the heaving mechanism. Work done on similar problems elsewhere and studies conducted by the author indicate that the latter two factors are the most important. Results of this investigation indicate that heaving areas are associated with slaking shale or expansive clays, with and without salts present. Laboratory studies show that the initial moisture content of remolded samples and the desiccation history of un- disturbed samples have a major effect on swelling. Swell tests on compacted samples at various temperatures indicate that water has more influence on swell than temperature, except for samples containing sodium sulphate, which swelled with decreasing temperature. Chloride salts inhibited swelling. Tests on shale cores indicate that the stability of slaking shale can be preserved by maintaining its original moisture content. It is concluded that shales and expansive clays, when exposed in road cuts, become desiccated and are subject to expansion when moisture is subsequently available. The seasonal cycle of heaving may possibly be explained by the annual moisture excess and deficiency. To minimize the problem, freshly exposed shale should be kept wet, and the thickness of the subbase and base courses should be increased. The Utah State Department of Highways has been plagued for many years with a distressed condition that develops in roads constructed on the Mancos Shale. Locally these distressed areas are known as heaves, frost boils and alkali boils. These names refer to the surface expression, time of year

Jan 1, 1965

By D. M. Milne, R. F. Stoner, W. Bakiewicz

The Lower Indus Basin deposits form a substantial unconfined alluvial aquifer of excellent quality. This aquifer contains basically highly saline groundwater, overlain in places by considerable amounts of fresh water whose origin can be traced to the present and historical courses of the River Indus. In the past, the main source of groundwater recharge has been direct river seepage and water tables were generally deep. Since the establishment of a perennial irrigation regime however, deep percolation from the increased irrigation supplies has predominated and the water table has risen to within a few feet of the surface. The river system itself is being progressively devel-oped for irrigation and very shortly a situation will arise where the entire low winter flow will be utilized. Fresh groundwater will then be exploited to supplement the surface supplies, and at the same time a drainage system, mainly based on wells will be used to stabilize the water table at safe depths. This paper discusses some of the methods used to ensure optimum economic development of the fresh groundwater and describes the general design philosophy evolved by the authors' firms for both irrigation and drainage wells, with respect to both well configuration and operation. Water quality standards are discussed insofar as they affect subsequent use of the groundwater. A few particular aspects are considered, including the values of the relevant aquifer constants, the construction of wells, and the determination of screen losses. A major problem solved during the Lower Indus Project investigations has been the design of wells to produce fresh water from layers overlying saline water. This solution is also discussed. The recent interest in the groundwater of the Lower Indus Plains of West Pakistan has two main aspects. Firstly, soil waterlogging and the associated salinity hazards caused by rising water tables are a serious problem in agriculture, and the ease with which drainage requirements can be met by aquifer pumping makes this an obvious solution. The second interest stems from the rapid intensification of irrigation agriculture and the growing shortage of water supplies. In the areas where the groundwater is of suitable quality it can be used immediately for irrigation at relatively small cost and with large financial returns. In addition, the aquifer forms a potential storage reservoir which could be recharged from high river flows and depleted by pumping at times of scarcity. SOURCE OF GROUNDWATER The Lower Indus plain is a relatively long and narrow alluvial basin on a north-south axis, some 23,000 sq miles in area, bounded by higher land to the east and west. One hypothesis to explain the distribution and quality of groundwater in the region1 presupposes a river system considerably larger than, but very similar to the present river. This is characterized by a single mainstream, with subsidiary overspill channels in times of flood, running along the crest of a self-deposited ridge and with an extensive advancing delta. Sea water is contained in the alluvial pores of the advancing delta, forming the native groundwater. This has been partially displaced from two sources, namely direct seepage of fresh water from the river channel in meander flood plain areas and percolation from flood spillage in cover flood plain areas. Percolation has been subject to much re-evaporation and is thus slightly saline even as it reaches the groundwater. The whole complex has been further affected by two main influences: evapotranspiration from the land surface that tends to salinize the groundwater and down-valley flow that tends to bring fresher water

Jan 1, 1970

By AIME AIME

BORN a Chicagoan, on Aug;. 31, I879. Oscar Johnson lived in that city until young manhood. As a youth delivering newspapers he set his sight on the l university of Chicago and at sixteen years of age he began working summer vacations in the foundry of the Pullman Company. Misfortune came to his family during the summer of 1898 when tuberculosis took his only sister and brother. His mother, apprehensive of Oscar's health, urged him to try the West at least for awhile. So late summer found Oscar in Salt Lake City, employed by the Trent Engineering & Machinery Co. At that time it was busy building equipment for use in the war with Spain. Though he did not know it then, Oscar war was laying the groundwork for his present membership on the Mining Equipment Industry Advisory Committee of the W.P.B., a post which keeps him commuting between Denver and Washington.

Jan 1, 1944

By Harmon F. Fisher

In the course of a comprehensive investigation for the development of the electrical process for the dehydration of oil-field emulsions, the writer has had an unusual opportunity to direct and execute a number of related, parallel studies concerning crude oil emulsions and extending over a period of nearly six years. Among these was an investigation utilizing the microscope for studying both the general characteristics of crude oil emulsions, the nature of their electrical charge and the action of an alternating current field in causing coalescence of the finely dispersed, more or less permanently suspended, emulsion particles into larger droplets of such mass that subsequent dehydration would follow by the ordinary process of settling, further agglomeration of the water droplets and, finally, separation into commercially dry, pipe line oil and a clean water discharge. It is not the purpose of this paper to enter into a detailed discussion of the highly specialized subject of microscope technique, nor into a theoretical discussion concerning the origin of the electrical charge on the emulsion particle or concerning the mechanism of the electrical phenomena occurring during the electrical treatment of oil-field emulsions. In a short paper of this character one can hope to touch on only a very limited discussion of the theoretical aspects of crude oil emulsions. The subject of emulsions and colloids is comprised in a comparatively young science in which, however, very rapid progress is being made and consequently it is undergoing many changes; a tentative conclusion valid today may be overthrown tomorrow. Nor would it be a simple science even if it were no longer undergoing this rapid progress; a thorough study of emulsions is based upon a knowledge of such fundamental subjects as molecular structure, ionization and other electrical phenomena, surface chemistry, surface and interfacial tension phenomena, viscosity and other factors too numerous to discuss here and the fundamental theories of which are still undergoing more or less change. Therefore this paper will confine itself largely to a generally descriptive discussion, a sufficiently detailed outline of the microscope technique and a brief comment on

Jan 1, 1931

By Walter Miller

Crude oil stocks were some 10,000,000 bbl higher on June I, 1947, than at any time during 1946 but the extremely heavy refinery runs the last half of 1947 cut crude inventories to approximately the 1946 closing level of 225,000,000 bbl. Many refiners were seriously short of crude oil during the year. This shortage was reflected in numerous advances in the posted price of crude, justified in part as an incentive for more production but more generally due to competitive bidding. Some purchasers of crude, endeavoring to keep their plants running, resorted to the payment, either directly or indirectly, of premiums above the posted price. Another common practice employed by some suppliers of crude oil was to demand refined products in exchange or trade for their raw material. The refined products picture has darkened in spite of record throughputs. Stocks of gasoline and light and heavy fuel oils were uncomfortably low at the beginning of the year and were in general still lower at the close. The rationing of some products, principally gasoline, to dealers and stations? initiated during the year, gained momentum as stocks continued to diminish. The emphasis placed by the manufacturers of oil-burning equipment on the convenience and economy of heating with oil, together with the tight situation in coal, has developed a tremendous demand for burning oils. Consumption of liquefied petroleum gas for domestic and household fuel had more than doubled in 1946 over 1943, and the demand seems to be limited only by the ability of the equipment builders to furnish necessary storage. transportation, and burning equipment. The supply situation is likely to become more acute in 1948. Demand for petroleum products is being pushed to new heights not only by the unprecedented domestic market but by a sudden and drastic upswing in estimates of military and governmental requirements. Additionally, oil exports are running far below the volume requested, amounting to less than forty percent in a recent quarter, and export licenses are reportedly being approved on the basis of urgent need *only. AX recent White House release sees little likelihood that estimated requirements of the sixteen nations participating in the proposed Marshall Plan can be met in full over its projected life—1948 to 1951. Present demand dictates maintenance of refinery operations at the highest possible level and the Bureau of Mines has forecast daily average crude runs in 1948 of 5,265,000 bbl. The Economics Advisory Committee of the Interstate Oil Compact Commission recently estimated the 1948 demand for petroleum products at nearly 6 pct above 1947. They further estimate daily average crude runs of 5,375,000 bbl in 1948, some 2 pct higher than forecast by the Bureau of Mines. Projected demand exceeds supply to the extent that Capitol Hill is reportedly investigating ways of enabling oil companies to act in concert without danger of antitrust prosecution, such as elimination of crosshauls, pooling of products, and lowering of octane.

Jan 1, 1948

By Walter Miller

Crude oil stocks were some 10,000,000 bbl higher on June I, 1947, than at any time during 1946 but the extremely heavy refinery runs the last half of 1947 cut crude inventories to approximately the 1946 closing level of 225,000,000 bbl. Many refiners were seriously short of crude oil during the year. This shortage was reflected in numerous advances in the posted price of crude, justified in part as an incentive for more production but more generally due to competitive bidding. Some purchasers of crude, endeavoring to keep their plants running, resorted to the payment, either directly or indirectly, of premiums above the posted price. Another common practice employed by some suppliers of crude oil was to demand refined products in exchange or trade for their raw material. The refined products picture has darkened in spite of record throughputs. Stocks of gasoline and light and heavy fuel oils were uncomfortably low at the beginning of the year and were in general still lower at the close. The rationing of some products, principally gasoline, to dealers and stations? initiated during the year, gained momentum as stocks continued to diminish. The emphasis placed by the manufacturers of oil-burning equipment on the convenience and economy of heating with oil, together with the tight situation in coal, has developed a tremendous demand for burning oils. Consumption of liquefied petroleum gas for domestic and household fuel had more than doubled in 1946 over 1943, and the demand seems to be limited only by the ability of the equipment builders to furnish necessary storage. transportation, and burning equipment. The supply situation is likely to become more acute in 1948. Demand for petroleum products is being pushed to new heights not only by the unprecedented domestic market but by a sudden and drastic upswing in estimates of military and governmental requirements. Additionally, oil exports are running far below the volume requested, amounting to less than forty percent in a recent quarter, and export licenses are reportedly being approved on the basis of urgent need *only. AX recent White House release sees little likelihood that estimated requirements of the sixteen nations participating in the proposed Marshall Plan can be met in full over its projected life—1948 to 1951. Present demand dictates maintenance of refinery operations at the highest possible level and the Bureau of Mines has forecast daily average crude runs in 1948 of 5,265,000 bbl. The Economics Advisory Committee of the Interstate Oil Compact Commission recently estimated the 1948 demand for petroleum products at nearly 6 pct above 1947. They further estimate daily average crude runs of 5,375,000 bbl in 1948, some 2 pct higher than forecast by the Bureau of Mines. Projected demand exceeds supply to the extent that Capitol Hill is reportedly investigating ways of enabling oil companies to act in concert without danger of antitrust prosecution, such as elimination of crosshauls, pooling of products, and lowering of octane.

Jan 1, 1948

By AIME AIME

THE annual meeting of the Auxiliary to the A. I. M. E. was marked by the most delightful cordiality and warm spirit of welcome on the part of the members of the New York Section and an equally charming spirit of cooperation and appreciation on the part of the visitors. On Monday Mrs. Ernest H. Wolfes and Mrs. Charles A. Bohn had charge of registration and most efficiently .provided everyone with the necessary tickets and badges. The reception hall and regisfration booth were gay with beautiful flowers that reflected the welcoming spirit one found there. The smooth efficiency with which the whole program ran was due largely to the tireless efforts of Mrs. Donald M. Liddell and Mrs. Louis D. Huntoon, joint chairmen of the Ladies' Committee, and Miss Margaret Callow, chairman of the Junior Committee. After registering, the ladies proceeded to- the fifth floor, where they had luncheon with the members of the A. I. M. E. On that day Mrs. W. H. Bassett and Mrs. Lucius W. Mayer were the hostesses.

Jan 1, 1929

By Walter C. Smith

THE Cerro de Pasco Copper Corp. began to produce bismuth at the Oroya smelter in 1929, at which time the only important consumption of that metal was in the manufacture, of pharmaceutical compounds, and the normal supply of metal exceeded the demand. With a new and major producer entering the field, additional uses for bismuth had to be developed. A survey indicated that the most promising opportunity lay in fusible alloys containing bismuth. A large number of bismuth alloys were known and a few of them were being used in the manufacture of automatic sprinkler heads and safety plugs. The literature cover in,- the- alloys was found to be sketchy and contradictory. Research, development, and educational work carrried on by the Corporation during the past fifteen years has resulted in greatly enlarged uses for the metal. In general the low-melting bismuth alloys are of high specific gravity, relatively soft, of low tensile strength, brittle to shock, yet subject to plastic

Jan 1, 1945

By L. W. Emery

Undergrorlnd combustion operations were initiated in a 60-acre Bartlesville sand "shoe-string" reservoir in Allen Connty, Kans., in 1956. Tests in separate patterns were conducted using various co~nbinations of air and recycle gas to propagate combustion fronts from the injection toward the producing wells. These patterns were made up of 6 injection and 20 prodrrcing wells Gas and liquid prorluctiorz from each pattern was measured on an individual-well basis, and comparisons were made between the three patterns to ascertain the relative effects of injected gas composition on production behavior. Breakthrough of the combustion front at a well was characterized by an increase in water production from the well followed by an increase in bottomhole temperatrrre to approximately 250" F. After burning fronts had broken through at five producing wells, operations were terminated in 1960. From the total project approximately 79,000 bbl of oil were produced during thermal operations at a cumulative produced GOR of 23 Mcf/bbl. No appreciable change in the character of the produced crude was observed. Combustion in the reservoir was maintained with injected gas compositions ranging fronz 6 per cent oxygen in recycle gas to 100 per cent air. lnjectiotz of large quantities of recycle gas resulted in higher producing GOR's from offset wells than were measured from a pattern into ~vhich straight air ~vas injected. The air required to move the combustion front through I acre-ft of reservoir was computed to be 20 MMscf. This valrre was found to be relatively independent of the quantities of recycle gas injected. The recovery efficiency from the swept area was esti~izated to be about 59 per cent. Areas swept were similar in shape to tlzose obtained with a laboratory potentiometric model. Samples of sund taken from behind the burning front by coring indicated almost total oil removal from the sand. Petrographic analysis of the core samples indicated that the sand had been heated to peuk temperature of rlbout 1,200" F. No rignificant difference in peak temperature was forrnd in two areas where compositions of injected gas were quite different. Compression costs for thermal recovery were estimated to be $1,20/bhl of produced oil. INTRODUCTION The use of the "forward combustion" process as an oil recovery method has received a great deal of attention. This method involves ignition of the formation in an injection well, followed by propagation of a combustion front through the reservoir. Combustion is maintained by the injection of an oxygen-containing gas to react with reservoir hydrocarbons. As the flame front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase and recondensed in the cooler part of the formation. In turn, the condensed fluids push oil into the producing wellbores. Completed field tests of the process were first reported by Kuhn and Koch,' and by Grant and Szasz.' Results from other tests have since been reported by Walter,3 by Moss, White and McNeil,' and by Gates and Ramey." ach of these tests essentially utilized a single injection well surrounded by four or more producing wells. Sinclair Research, Inc., elected to do field experimental work using a number of test patterns in a single field in order that comparisons between various operating schemes could be made. The site selected and purchased in 1955 for this experimental work was a 60-acre Bartlesville sand reservoir located in Allen County, Kans. Combustion operations were initiated in mid-1956. Between that time and termination of the project in mid-1960, combustion fronts were propagated from injection wells to producers in three separate well patterns, using different mixtures of air and recycle gas. The test was terminated before sweep of the three patterns was complete so that information about the effect of combustion on the swept areas could be obtained by coring. Results from the test in the form of injection and producing well performance have been carefully recorded, and these form the general basis for this paper. DESCRIPTION OF RESERVOIR The reservoir in which the combustion tests were conducted is a Bartlesville sand "shoe-string", typical of a number of small reservoirs in Southeastern Kansas. Average reservoir characteristics are shown in Table 1. Fig. 1 is an isopachous map of the producing sand showing the reservoir to be approximately 400-ft wide and 2,500-ft long. Maximum net productive sand thickness is 21 ft. Fig. 2 shows a typical core analysis obtained by coring with water-base mud. The reservoir has no appreciable dip and is closed on the sides by degradation of sand into shale. The main body of sand is heavily laminated with shale stringers, which are not continuous between wells. The main reservoir is overlain by 30 to 40 ft of laminated low-permeability sand and shale streaks. No information is available on the original properties of

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By A. R. Gregory, M. R. J. Wyllie

The literature reveals that scant attention has been paid to the systematic experimental determination of the formation factors of unconsolidated porous media. No experiments appear to have been made on the effect of increasing cementation on the formation factor of an initially unconsolidated porous medium. Measurements have accordingly been made of formation factors as a function of porosity for aggregates of spheres in the porobity range 12-56 per cent and of cubes, cylinders, discs and triangular prisms in the porosity range 30-45 per cent. The results are examined in the light of the theoretical equations of Clerk Maxwell. Fricke and Slawinski. Aggregates of unconsolidated spheres and beach sands have been artificially cemented with silica and the formation factor-porosity relationships determined. A theory is outlined which seeks to explain the results obtained and which postulates that formation factor and porosity data for consolidated sandatones may be used to determine the original porosity of the unconsolidated sands from which the consolidated sandstones were derived. It follows also that for consolidated sandstones F' = C where F is formation factor, fractional porosity and C and k are constants. INTRODUCTION All schemes of electric log interpretation implicitly or explicitly require the use of formation resistivity factors or? as they are now generally called. formation factors. 'The originator both of the concept of formation factor and of the use of formation factor data in electric log interpretation was Archie.' More recently it has been suggested that formation factors are of value in describing fluid flow in porous media. A full discussion of this problem together with adequate reference.: has been given by Wyllie and Spangler.' In spite of the theoretical and practical importance of the subject, there appear to be remarkably few basic formation factor data in the petroleum literature. A wider search reveals that this situation is nut confined to the literature of the to industry. at various times over the last century, attempts have been made to compute the theoretical relationship between the formation factor and porosity of simple systems of unconsolidated particle.. However. the formulae derived are not all identical and the experimental evidence adduced in support of each i; not wholly convincitng. Most of the experimental work has been carried out on system. of unconsolidated spheres, although some attention has been given to systems of spheroids and cylinder. Other simple shapes have received neither theoretical nor experimental investigation. The writers feel that a better understanding of the formation factors of complex systems of practical importance requires a firm foundation of reliable experimental data applicable to simple systems. Accordingly, they have reviewed the relevant literature and have determined experimentally the relationship between the formation factor and porosity of unconsolidated aggregates of particles of regular geometrical shapes. Some experiments have also been carried out to investigate the effect of cementation on the formation factor-porosity relationship. PREVIOUS WORK Although the term formation factor was not used. the relationship between the resistivity of a solution containing dielectric solids and the resistivity of a similar solution not containing solids is an old and important problem. The formation factor of a system of dispersed spheres as a function of porosity was considered by Clerk Maxwell. Maxwell arrived theoretically at the following expression: where F is formation factor and fractional porosity. Lord Rayleigh' also examined the problem of the formation factor of systems of spheres and cylinders and concluded that for spheres Clerk Maxwell'.; expression gave a value of F whirl]. at low porosities particularly, was too small. However,

Jan 1, 1953