Search Documents

By R. L. Slobod, C. Gatlin

This .study defines the basic mechanism of the mis-cible displacerrzent of oil and writer from porous Medici by various water-driven alcohol .slugs. Three distinct alcohol slug processes were .studied. Considerable data concerning the quanity of alcohol required for oil recovery were also obtained. All data were obtained in a 1-in. diameter, 100-ft long, unconsolidated core. The porosity of this system was 35 per cent. and the permeability was approxitrlately 4 darcies. Total core pore volume wax 5,716 cc. All displacements were conducted at a constant injection rate of 5 to 6 cc//min, which corresponds to a frontal advance of 5 ro 6 ft/hr. The first portion of this paper is concerned with 111 use of one alcohol—isopropyl—a\. the slug material. Isopropyl alcohol (IPA) is completely miscible with both oil rind water; however, miscibility of the threc-component systems, oil-water-IPA, requires a relutively. high concentration of IPA. Hence, the displocement is not of the misoible type unless the IPA concentration is maintained above some critical value. A slug of IPA equal to only 13.5 per cent of the pore volume was jound to be sufficient to obtain complete recover?, of residua1 naphtha. In later studies two distinct process variations were developed. The firsr of these utilized methyl alcohol (MA) and IPA as slug materials. It was shown that methyl alcohol may he substituted for IPA at the front and rear of the slug with no loss of oil recovery. A .slug of 4 per cent MA-4 per cent IPA-4 per cent MA was sufficient for complere oil recovery. Because MA is considerably cheaper than IPA, this represents an important step toward economic application. A .second process variation used normal butyl alcohol (nBA) and MA as the composite slug, the nBA segnzent being injected first. This technique requires the smallest total slug size (approximately 10 per cent) of all processes studied. The high cost of nBA, however, precludes commercial application. It is possible that this basic process, subject to changes of alcohol type, may lead to a commrercicil process. INTRODUCTION Within the last 10 years, considerable study has been devoted to the general mechanics of miscible phase dis- placements in porous media. These studies have. in general, dealt with two-component displacernent. All reported field trials to date have been modified gas drives of one type or another. These utilize either lean gas at high pressure, rich gas at lower pressure or an LPG slug. The irreducible water content of the reservoir remains immobile during these processes. Although these techniques hold great promise. certain drawbacks do exist. Poor areal sweep efficiency, inherent in any displacement having a highly unfavorable mobility ratio, is one disadvantage. Furthermore, the miscibility of all processes which use gas is dependent on pressure. These pressure limitations may prohibit application to shallow reservoirs. There are also many areas where large quantities of LPG and natural gas are not readily available. These factors emphasize the need for improved techniques. The alcohol slug process is also a miscible-phase displacement process. It differs from the previously mentioned two-component techniques, however, in that both the reservoir oil and water are displaced by the slug. This behavior is a consequence of the miscibility of certain oil-water-alcohol systems. In the simplest case, a relatively small volume (slug) of an alcohol (such as isopropyl) is injected into the system. Water is then used to drive the slug through the medium. Thus. three components—-alcohol. oil and connate water-—exist at the front side of the slug. Miscibility is obtained at some alcohol concentration, this being dependent on the solubility of the particular system. Miscibility is maintained within the slug until the alcohol concentration falls below that value required to maintain miscibility. When miscibility is lost, the process reverts to a water flood. Certain inherent benefits of this technique are apparent. High pressures are not necessary to obtain miscibility in these liquid systems. Furthermore, water is a more desirable driving agent than gas because of the improved mobility ratio and the attendant improvement in areal sweep efficiency. The main disadvantage is, of course, the relatively high cost of alcohols. The use of an alcohol slug to recover reservoir oil is not a new idea." ' However, this process has received only limited study, presumably due to the lack of industrial interest caused by the seemingly prohibitive cost of alcohols. The main purpose of this study was to investigate the mechanism of the displacement of oil and water by various alcohol slugs and to develop process modifications aimed toward possible commercial application.

By Charles R. Stewart, William W. Owens

Reservoir performance predictions based on laboratory core test data assume that fluid flow is laminar for the laboratory test. A study has been made to determine the validity of this assumption for laboratory tests on various types of porosity found in producing limestone formation. Data are presented which show that turbulence and slippage can occur during laboratory tests on hetero-geneous-type porosity limestones, thus causing serious errors in measured single-phase permeabilities and two-phase relative permeability characteristics. In single-phase flow tests it is possible to eliminate turbulence and correct for slippage or to eliminate both factors by controlling test conditions. It is not always possible to control test conditions and thereby eliminate turbulence and slippage in two-phase .flow tests. A correction method is presented which can be used to calculate the true two-phase laminar flow relative permeability characteristic even though furbulence and slippage exist. .INTRODUCTION It is customary to make use of Darcy's law and modifications of this law, together with laboratory data on formation core samples to predict the performance of producing reservoirs. Such predictions are based on an assumption that fluid flow is in the laminar or streamline region for the laboratory test. It was the purpose of this inves- tigation to determine the extent to which turbulent flow may occur in laboratory fluid flow tests on hetero-geneous porosity limestones. Considering that turbulent flow conditions might exist in some laboratory fluid flow tests, additional emphasis was placed on the development of a method to correct for turbulence when laminar flow conditions could not be attained. FLUID FLOW CONCEPTS FOR POROUS MEDIA The Influence of Pore Geometry on Fluid Flow One of the more important factors influencing fluid flow in porous media is the geometry of the pore space which includes such characteristics of the pores as size, shape, distribution, roughness, uniformity, etc. In general, oil- and gas-producing formations can be divided into two broad types on the basis of pore geometry. One has been called sandstone-type porosity media, which is characterized by a small range in pore size, uniformity in shape of the pores, smooth pore surfaces and a regular and uniform distribution of pores. The other type has been called heterogeneous porosity media and is usually limited to the dolomites and limestones. This type is characterized by a wide variation in the size, shape, and distribution of the pores and rough, irregular pore surfaces. It is therefore apparent that conditions are much more favorable for turbulent flow* in heterogeneous-type porosity media than in sandstone-type porosity media. Interrelationship Between Turbulence and Gas Slippage In studying the problem of turbulent flow in laboratory tests on porous media, it is necessary to be aware of the interrelationship between slippage and turbulence for gas flow. As a result of slippage or the Klinken-berg effecta, apparent perrneabilities to gas are greater than the true value because there is no stationary layer of gas in contact with the walls of the flow channels. Gas slippage decreases as the mean free path of the gas molecules decreases. Since the mean free path of any gas decreases with increasing density, increases in static pressure result in lower apparent gas permeabilities. However, a reduction in gas permeability can also be due to turbulence. Therefore, in studying only turbulent flow in porous media, it is necessary to hold gas density, and slippage, constant or to reduce slippage to a negligible value by operating at high static pressures. Presentation of Laminar and Turbulent Flow Data A graphical relationship between permeability and a pseudo-Reynolds number, N,, will be used to show the two types of fluid flow, i.e., laminar and turbulent. The usual graphical method for such a description has been the use of friction factor-Reynolds number charts4. On such a logarithmic diagram, the laminar region appears as a straight line having a slope of 45 degrees. As the friction factor decreases and the Reynolds number increases, the turbulent region is reached and appears as a deviation from the 45-degree slope line. However, in petroleum engineering literature resistance of por-

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study' that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,' flow in packed columns,' fracture flow capacity' and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By Carl Gatlin, S. Srini-Vasan

It is shown thar k k ratios cornputed on the as.sutnption that the average sandstone reservoir has a pore-size distribution characrerized by a linear 1/P saturatioN plo are in good agreewient with correspotlding k,,/k,, ratios obfairzed by Wahl, Mul-/ins and Elfrink by correlatiotz of field measurements. Methods of calculating relative permeability data have recently been reviewed and extended.' It was shown that for sandstone reservoirs where ' is approximately a linear function of effective saturation (as suggested first by Corey') it is possible to express relative permeabilities in very simple terms. For systems wetted by connate water and producing oil and gas (gas saturation growing at the expense of oil) the expression for the The accuracy of Eq. I has been frequently confirmed in our laboratories; because of its simplicity we make routine use of it in the analytical solution of proble~ns in reservoir engineering. Nevertheless, comparisons of the k,,,/k,, ratios yielded by Eq. 1 with appropriate reservoir data have been few. The publica- tion by Wahl, Mullins and ElfrinkV of an empirical equation for k,,/k,, ratios is thus of interest since this equation is specificalljl stated to be Vased on the correlation of field measurements made on sandstone reservoirs. Not surprisingly these authors found that it was necessary to

By G. R. Countryman, G. H. Thomas, I. Fatt

Procedures are described for scaling up turbulent friction pressure drops of inelastic fluids and laminar friction pressure drops of inelastic/elastic fluids in circular ducts. The laminar scale-up calculation method does not require the use of arbitrarily defined flow models, and can be performed directly from laminar capillary viscometry data, or by means of a numerical solution of the integral equation for Poiseuille flow, inputing true shear rate-shearing stress data from other simple shearing experiments. e.g. Couette and plane Poiseuille flows. The turbulent scale up procedure, based on the Dodge-Metzner correlation for inelastic fluids, requires evaluating the characteristic rheological parameters at the existing wall stress under laminar flow conditions. Since the wall stress is the quantity sought in a pipe flow scale-up problem, a trial-and-err or solution is indicated, and therefore a calculation method has been developed and programmed for a medium-size magnetic-drum-memory digital computer. Laminar flow data from other somces may also be adapted to this turbulent scale-upmetbod for inelastic fluids. Comparisons between predicted and experimentally observed friction losses for generalized Reynolds numbers ranging from 3 — 32.000 indicate an error band width of less than 10 per cent in circular tubes ranging in size from 0.5 to 1.89 in. in diameter. INTRODUCTION Considerable effort has been devoted to the problem of establishing an accurate procedure for predicting the laminar and turbulent flow behavior of non-Newtonian materials in pipes from data obtained with small scale laboratory or pilot plant equipment, as evidenced by the variety of reviews on the subject.1-3 Although there does not appear to be a universally accepted scale-up method for laminar flow, the various methods currently employed can be divided into two principal categories depending upon whether a functional relationship between shear rate and shearing stress is assumed, or whether a generalized correlation is adopted. The methods of Hedstrom4 and of Weltmann are typical of the first category, while the generalized correlation of Metzner and Reed6 illustrates the second category. Under appropriate conditions either approach yields a reasonably satisfactory description of laminar flow behavior. In contrast, the prediction of non-Newtonian transitional and turbulent friction losses is much more complex. For example, certain types of materials have been observed to sustain laminar flow at much higher flow rates than Newtonian liquids. When turbulence does develop, the level of agitation is damped to such an extent that the friction factor is lower than that of the original low viscosity Newtonian medium. Extremely broad transition regions and "diameter effects" have also been noted in the turbulent flow of certain chemically complex materials.' Aside from the latter complications, which are attributed in part to viscoelasticity, the existence of a variety of scale-up procedures for predicting transitional and turbulent friction losses is due principally to the chioce of viscosity parameter used in the definition of the Reynolds number. In contrast, Dodge and Metzner have demonstrated that the friction factorReynolds number relationship for inelastic fluids can be satisfactorily generalized provided the fluid property parameters are evaluated at the existing wall shearing stress.* However, the quantity usually sought in a pipe flow scale-up problem is the wall stress, and hence a trial-and-error procedure is required when dealing with non-power law materials. In this report we discuss two calculation methods: (1) for determining the wall stress, and hence the friction pressure drop, for inelastic fluids under conditions of turbulent flow by means of an iterative procedure, and (2) for predicting laminar friction pressure drops without recourse to arbitrarily defined flow models.

By K. A. Slagle, D. K. Smith

weight obtained. Additives used in conjunction with salt in these slurries have included silica flour, calcium ligno-sulfonate and cellulose retarders, granular lost-circulation materials, bentonite and selected low-water-loss additives that are not significantly deteriorated by the presence of the chloride ion. In some other areas of South Texas, the salt-saturated slurries have been used quite extensively for improvement of the flow properties of slurries and attainment of better circulation characteristics at lower displacement rates. Concurrently with this property, the protection of shales and shaly sands is also realized, as well as useful retardation of the slurry. Resultant superior cementing jobs have been indicated by both communication tests and acoustic logs for bonding to pipe and formation.' In one section of Louisiana, a major oil company has been using salt-saturated API Class A cement with calcium lignosulfonate retarder for cementing through the Miocene at 9,000 to 10,000 ft. This is another of the situations where interbedding of sands and shales exists, creating difficulty in maintaining formation competency when a fresh-water slurry contacts the clay minerals of the formation. Further work has also been done in the shaly Miocene formation at 13,000 to 15,000 ft where fairly close water or gas contacts are encountered. Indications thus far are that better segregation of these various fluids is obtained by use of the saturated salt slurries because of their improved formation bonding characteristics. In addition to the properties of salt in this situation, attainment of turbulent flow at minimum displacement rates has also been accomplished by use of an additive to help provide exceptional dispersion and viscosity reduction of the slurry. Another oil company was encountering considerable expense in completing wells in Southwestern Louisiana due to extensive block squeeze requirements for effective separation of zones. A very effective mud program was being used to minimize washout in the shale sections and, apparently, a nearly gauge-size hole was being obtained. However, primary cementing results with fresh-water slurries were generally poor. On a few occasions when slurry was actually circulated to the surface, large pieces of shale formation were brought out of the hole with the slurry, indicating a severely water-sensitive, sloughing formation. Inhibition of shale heaving was being accomvlished in the drilling program, and immediately indone ipon circulation of the fresh-water slurry even though it contained a low-fluid-loss additive to reduce filtrate damage in the sands. The subsequent change to salt-saturated slurry yielded 11 successful primary cement jobs out of 12, compared to the previous success ratio of practically zero. Since these were deep, high-temperature wells in the range of 13,000 ft with high-pressure zones necessitating 17.5-lb/gal fluid densities, the slurry used was API Class E cement, silica flour, weight material, retarder, salt saturation (which also reduced the amount of weight material) and maintenance of low fluid loss by use of a salt-compatible additive. Other salt slurries have been used to a limited extent in this same general area for similar problems at depths ranging from 5,000 to 17,000 ft. In the shallower wells, the cement has usually been API Class A where the salt functions as a retarder, and in the deeper wells API Class E cement is used with the additional salt advantage being its increased slurry weight and inability to dissolve salt stringers. A considerable number of squeeze jobs have also been done on older wells using the salt-saturated slurries with very good results. MID-CONTINENT In North Texas, salt-water slurries have been used for cementing the Woodbine sand, Strawn sand, KMA sand and Pettit lime. Shales surrounding these formations have created the same difficulty in obtaining separation of producing zones that has been the problem in other areas. Depths in this area range from 3,400 ft for the Strawn to 7,400 ft for the Woodbine, and concentration of salt has varied accordingly. In the deeper wells, where retardation is desired, saturated salt-water cement is used; for the shallower wells, in order to provide shorter waiting-on-cement times, the amount of salt has been 18 per cent by weight of the mixing water. Results have been excellent with no reported failures on any of the salt cement jobs; where acoustic bond logs have indicated indifferent bonds previously, they are indicating very good bonding for the salt slurries. In Oklahoma various shales of Pennsylvanian age exhibit a high degree of sloughing in the presence of fresh water, causing severe washouts above and below sand formations which it is desired to isolate. This situation exists to some degree in practically all parts of Oklahoma and includes formations of other ages such as the Wood-ford shale. For the past few years, salt-saturated cements and displacement rates as high as practical have been used as a remedy for this problem, with very good results. The Layton and Bartlesville formations are two examples of shaly sands where saturated slurries have been helpful. In one area where five wells were drilled through this type of problem shale without obtaining a satisfactory primary cement job, a change was made to salt-saturated cement preceded by a suitable chemical wash for the drilling mud involved. Acoustic bond logging indicated excellent bonding, and final completion bore out this result by being trouble-free. This type of slurry has also been used extensively on squeeze jobs where shales have been heaving around the producing formations. Predominantly, the basic slurry has been either API Class A cement or a pozzolan cement—although, as deeper wells are being drilled, the use of salt in Class E cement is also increasing. Salt cement in West Texas has been used primarily to help prevent channeling through the shales around the Delaware sand, Queens dolomite and Hope lime. Many of the shaly and dirty sands of this area are sensitive to the filtrate from a fresh-water cement. Salt at 16 to 18 per cent by weight of mixing water has been added to centent, and has been effective in controlling formation damage and communication between zones in these formations. Use of these lower salt concentrations is dictated in this area by the relatively low formation temperatures where retardation of the slurry would create unduly long waiting-on-cement times. Also, quite a bit of cementing has been done in this area using salt concentrations in the accelerating range— that is, 2 to 5 per cent by weight of water. Specifically, these concentrations have been used in coiljunction with high percentages of gel to overcome the retarding effect of the calcium lignosulfonate dispersant, although there are probably several shales where these concentrations could provide some degree of formation stability. On several occasions, the salt has also been used to lower the critical velocity or rate for turbulence with the slurry, particularly in the pozzolan cements. On wells in the Hope lime, it has usually been necessary to squeeze the shale above and below the lime to get a water-free completion. Use of salt-saturated slurries has largely eliminated this

By J. B. Campbell, V. V. Valleroy, B. T. Willman, L. W. Powers

A steam drive of heavy oil was field tested in a shallow, low oil-saturation formation near Deerfield, Mo. The pilot was conducted in the Warner formation, a sandstone containing an 18' API oil having 1,000-cp viscosity at the 60F origind reservoir temperature. The formation war at a depth of 160 ft. Steam was injected into nine input wells arranged in an array of inverted five-spot patterns. In the completely confined center pattern, 14 temperature observation wells were installed to obtain thermal data and observe test progress. Late in the test, slugs of ammonia were injected to trace the flow paths of injected fluids. From the test area about 7,000 bbl of oil were produced. Data were obtained on areal and vertical temperature distribution, steam front advance, reservoir fluid movement and terminal saturations. This field test of a steam drive (I) demonstrated the feasibility of the method, (2) confirmed that the low residual oil saturations observed in the laboratory are obtained in the steam-swept region in the field and (3) provided recovery and conformance data for one set of field conditions. INTRODUCTION The Deerfield steam drive pilot test was conducted in a shallow sandstone containing 1,000-cp oil. The venture was undertaken cooperatively by the research and production departments of Carter Oil Co., which organizations have since been consolidated into Esso Production Research Co. and Humble Oil & Refining Co.. respectively. The production department was interested in steam injection at Deerfield because it appeared to be the most promising method of commercially producing this heavy oil deposit. The research department was interested in applying the new recovery method and in evaluating its performance in the field. At the time the test was begun, the initial oil saturation was not well known. Subsequent air coring and early pilot results confirmed that there was too little oil in place for profitable commercial exploitation by steam. Pilot termination at that time, however, would have been premature for evaluating field performance of the process, and the tert was continued to obtain additional data on steam injection as a recovery method. The test was located in Vernon County, Mo., about 10 miles north of the town of Deerfield and only a few miles from the Kansas border. The pilot site was selected as typical of the area. The location represented neither the highest nor the lowest oil saturation region in the acreage under lease in 1954. The steam drive was conducted in the Warner sandstone of Lower Pennsylvanian age. At the test site the top of the Warner occurs at about 160 ft subsurface and the formation is a fine- to medium-grained micaceous sandstone that dips gently to the northwest at the rate of 12 to 15 ft/ mile. A cross-section and permeability profile of the test location are shown in Fig. 1. At the pilot location the average total thickness of the Warner formation is about 43 ft, but the effective thickness for steam drive is 26 ft. Two distinct types of hydrocarbon saturation are apparent. The lower portion of the total sand, averaging about 17 ft thick, contains a very heavy asphaltic material that will not flow under the influence of a steam drive. This bottom interval, referred to as a dead oil residue, was not considered as part of the net sand undergoing steam exploitation. The initial formation and fluid properties of the upper 26 ft in the test area are summarized in Table 1, and variation of oil viscosity with temperature is shown in Fig. 2. Imbibition tests on preserved core samples taken at the end of the pilot test showed that the Warner sandstone was then neutral or slightly water-wet. Initially, the reservoir may have been strongly water-wet as indicated by low relative permeability to water during both water injection testing and early steam injection. PRIOR HISTORY Initial production tests of wells at the pilot site produced water with only a faint show of oil. No gas was produced except at Well 7-W in the pilot area and at another well about 1/3 mile northeast of the pilot. Prior to the start of the steam drive, a two-well water injection test and a two-well air injection test were conducted. No oil was produced by either. Water was pumped into Well I-W in the northeast corner of the pilot area with simultaneous production from Well 1 (Fig. 3). The air-injection tat was run at input Well 9-W and its offset, Well 2, in the southwest corner. Air and water injectivities were about the same when corrected for viscosity and pressure differences.

By J. G. Richardson, K. H. Coats

During the initial growth of a gas bubble in an aquifer storage reservoir the injected gas tends to override the water. The resulting low displacement efficiency and high rate of gas travel down-structure make it difficult to maintain water-free production. This paper illustrates the use of two-dimensional, two-phase calculations to simulate this gas-water displacement. Calculations were performed for a stratified aquifer and for several homogeneous sands differing in permeability, porosity and thickness. Results are discussed and compared with performance predictions obtained from the simpler Buckley-Leuerett and Dietz formulas. This comparison indicates the necessity for the two-dimensional calculations in realistically simulating the gravity override of the water by injected gas. INTRODUCTION In recent years natural gas has been stored near markets in aquifers where insufficient storage capacity is available in depleted fields. Operators of these aquifer storage reservoirs have encountered technical problems relating to the gas-water displacement accompanying initial growth of the gas bubble. Injected gas tends to override the water, with a resultant low displacement efficiency and high rate of gas travel down-structure toward spill points. Low displacement efficiency makes if difficult to sustain water-free gas production. Sometimes the fingering of gas down-structure may be so pronounced that the injection rate must be severely curtailed, which excessively lengthens the time required for bubble growth. This paper is concerned with estimating — for given aquifer characteristics and fluid properties — the displacement efficiency, rate of gas movement down-structure and rate of gravity drainage of water behind the gas front. Several published articles relate to simulative capability necessary to handle this problem. The Dietz formula1 and Buckley-Leverett2 method have some applicability to the problem; Woods and comer3 reported one-dimensional, radial calculations which accounted for the two-phase flow of gas and water. Most recent research directed toward understanding aquifer behavior in relation to gas storage has concentrated on single-phase flow in aquifer.4 Douglas, Peaceman and Rachford5 presented a method for calculating multidimensional two-phase flow in reservoirs; their method has been applied in work reported by Nielsen and Tek,6 Blair and peaceman7 and Goddin.8 Refs. 5 and 7 show comparisons between experimental data and calculations similar to those employed here. One purpose of this paper is to illustrate the use of two-dimensional, two-phase flow calculations in simulating the difficult problem of water displacement by gas in a vertical cross-section. The pronounced fluid density and viscosity differences combine with the disparity in horizontal and vertical dimensions of the cross-section to pose one of the more difficult problems in reservoir simulation. The calculations described account for capillary and gravity forces, relative permeability and reservoir heterogeneity. An example reservoir is described and calculated results are presented for a variety of injection rates and values of permeability, reservoir thickness and dip angle. A second purpose is to compare the two-dimensional calculations with results obtained from the Buckley-Leverett and Dietz formulas. For this reason most calculations were performed for a simplified example reservoir to which the latter techniques might reasonably be applied. However, the computer program employed applies equally well to cases involving arbitrary spatial variations of permeability, porosity, reservoir thickness and dip angle.

By E. G. Woods, A. G. Comer

A mathematical model is presented and solved for determinution of the saturation distribution and pressure in a radial gas-storage reservoir. The model consists basically of two parts: (I) the growing gas-bubble core, and (2) the surrounding aquifer. Since the total pressure at the injection well is a function of the two-phase flow in the gas bubble and the unsteady single-phase flow in the aquifer, the resistance to flow in both zones was taken into consideration. The assumptions involved for both the radial equivalent of the Buckley-Leverett two-phase-flow equation and the injection-pressure equation are as follows: (I) the geometry is radial, (2) the gas bubble is free to expand or contract, (3) compression or expansion of the gas within the bubble may occur at the beginning of a time step, (4) the fluids are immiscible, (5) water is incompressible within the gas-storage region whereas it is compressible outside of this region, (6) a stable gas-water interface exists and (7) gas injection occurs at a constant rate or a series of constant rates. The mathematical model was solved numerically using an IBM 650 computer. A comparison is presented between the predicted results of the model, the results assuming steady-state flow and the actual initial injection-pressure history of an operating reservoir. Using the initial field pressure for a basis, the average deviation between the predicted pressure and the actual field pressure was less than 4.3 per cent. INTRODUCTION Underground storage of natural gas in abandoned oil fields, in abandoned coal mines, in caverns and in aquifers has had varying degrees of success. The lack of these first three facilities in the vicinity of most major gas marketing areas is leading more and more to storage in virgin aquifers. For the most part, this storage in aquifers has led to the study upon which this paper is based. Several excellent papers have been published discussing the problems involved in the underground storage of natural gas; however, to these writers' knowledge, none of them considered the two-phase flow of fluids in the gas bubble and the fact that the gas bubble will grow or shrink in size depending upon the injection and withdrawal history of the reservoir. Not only will both of these factors influence the required injection pressure, but also they will have a definite bearing upon the amount of water production occurring upon withdrawal from the storage area. To solve this problem on an IBM 650 computer, several simplifying assumptions were made. Only the drainage portion of the relative-permeability curve was considered, thereby neglecting any hysteresis effects introduced during an imbibition cycle. The quantity of gas withdrawn during the period of the field study was only a small percentage of the total gas in place; therefore, this assumption is justified. For the problem studied in this case, the combined effects of capillarity and gravity were assumed to be negligible. This assumption becomes less valid as the formation thickness or the average pore diameter increases markedly. The justification of these assumptions should be studied in each case. THEORY This mathematical model for the prediction of saturation distribution and pressure for a radial gas-storage aquifer is based upon the equations for radial two-phase fluid flow'.' and for radial, unsteady-state, single-phase fluid flow.' The two-phase flow is considered to take place in a "core" whose radius is equal to the maximum radius that the gas zone will attain, as indicated in Fig. 1. The injection-withdrawal history is approximated by a series of constant flow rates. Within the "core", gas is assumed to behave as a semi-compressible fluid; that is, the gas is assumed to have a constant density based on the average gas-zone pressure for the flow period. Between each constant-rate time increment, the gas density is allowed to change. This involves an iterative solution of the pressure equation. The liquid phase in this zone is considered to be incompressible, whereas it behaves as a compressible liquid outside of this region. The equations of incompressible flow in the "core" may be written for the gas phase

By Edward Hart

The specimen of glazy iron used for analysis was highly characteristic in appearance. It was made at one of the furnaces of the Glendon Iron Works, working a light burden of ore with a highly silicious cinder. There are very few records of analyses of this variety of iron. I have been able to find but the two given by Bell, and I add these to my analysis for comparison. Bell. 1. 2. Iron,......89.47* 88.18 90.70 Carbon, combined,.... 0.00 0.79 0.71 Graphite,......2.88 2.59 2 68 Silicon,......5.47 5 13 5.13 Manganese,.....1.54 0.77 0.56 Sulphur,......0.04 0.17 0.23 Phosphorus,.....0.60 1.12 1.12 Titanium,.....- 0.26 0.18 Calcium,......0.00 0.22 0.20 Magnesium,.....0.00 0.06 0.03 100.00 9929 101.54

By H. S. Price, R. L. Reed, J. E. Warren

The behavior of the reverse-combustion process in a linear adiabatic system is theoretically investigated by means of an idealized physical model. T his model is described by a pair of non-linear equations involving heat and mass transfer which are coupled by a concentration-dependent reaction-rate function of the Ar-rhenius type. The differential equations governing the quasi steady-state temperature and concentration distributions are approximately solved by using either physical or mathematical simplifications. It is shown that the reverse -combustion process can be mechanistically described by simple physical models whose behavior equations can be solved formally. It is further demonstrated that the derived reverse-combustion equations can be solved numerically within the error limits for the experimental data. The theoretically predicted peak temperature-flux-velocity relationship exhibits reasonable agreement with experimental data obtained from Athabasca tar sand. The principal contribution of this work lies in its usefulness as a starting point for developing a more comprehensive theory. A description of the vaporization-coking process must be obtained; then, the theory must be extended to three dimensions to include geometrical effects and heat losses. INTRODUCTION In another paper,'" he physical processes visualized as occurring during reverse combustion have been described. In this paper, the extent to which an elementary theory can account for experimental observations is determined. Theoretical results that can be used to estimate some of the behavior characteristics of importance to this process are also presented. Since it is not possible to include all the features of the combustion mechanism of which we are currently aware, the mathematical formulation refers to the highly idealized system developed herein. This simplified physical model is based on the simultaneous transfer of mass and heat accompanied by a chemical reaction: it is described mathematically by a pair of coupled, nonlinear differential equations. Similar sets of equations have been studied extensively in the field of flame propagation; an excellent review of these investigations has been published by Evans.1 FORMULATION OF THE PROBLEM It will be recalled that the experiments are performed with a thin-walled cylindrical tube that is uniformly packed with sand and oil. The system is maintained substantially adiabatic by means of heat supplied externally in such a way that the radial temperature gradient in any plane normal to the axis of the tube approaches zero. Initially, the tube is at ambient temperature, except at one end where it is rapidly heated to a predetermined temperature. When the prescribed temperature is achieved, air is caused to flow into the cold end of the tube. As the oxygen in the air stream contacts the hot oil, a localized exothermic reaction takes place. The generated heat is conducted and convected away from the reaction zone so that definite temperature and concentration profiles are rapidly established and move uniformly in a direction opposite to that of the air flow (Fig. 1). As a consequence of this mode of oxidation, all of the oxygen in the air is consumed; then, the remaining nitrogen and the gaseous combustion products sweep out all vaporized hydrocarbons and leave an immobile hydrocarbon-coke on the sand.

By J. H. Moran, E. E. Finklea

The pressure build-up technique is a recognized method of determining permeability from conventional drillstem tests. In this paper an effort is made to extend such techniques to the interpretation of data obtained from the wireline formation tester. Such a study is necessary because of the differences, for this case, in the magnitude of the flow parameters (rate of flow, amount of recovered fluids) and in the flow geometry (flow through a perforation vs flow across the face of the wellbore, etc.) involved in the solution of the equations of flow for compressible fluids. The perforation is replaced by a spherical hole, and the effect of the borehole is neglected, so that the flow can be considered to be radial in a spherical co-ordinate system. Arguments are presented to justify this idealization. Assuming single-phase flow, general relations between pressure and flow rate are developed for a homogeneous medium. The study is then extended to permeable beds of finite thickness. It is shown that the early stages of pressure build-up tend towards spherical flow, while the later stages tend towards cylindrical flow. The thinner the bed, the more quickly flow approaches the cylindrical model. The prevalence of thin beds in practical work makes this analysis quite important. Cases involving permeability anisotropy are treated. INTRODUCTION From wireline formation tester operation, two types of data are obtained: (1) the nature and amount of recovered fluids, and (2) the pressure history recorded during the test. A number of papers have been written dealing with the interpretation of formation production on the basis of the recovered fluids.'.' In general, the methods described have been quite accurate for both high- and low-permeability formations. The present paper will deal with an analysis of the pressures observed. An analysis of the pressure build-up curves obtained in hard-rock country has already been attempted on the basis of the formula proposed by Hor-ner. Although this approach has met with success in many instances, some questions have been raised as to its validity. It is the aim of the present study to place the analysis of pressure build-up in the formation tester on a firmer basis, from which more detailed methods of interpretation can evolve. Because of the great differences between the operation of the wireline formation tester and the conventional drillstem test, modifications are necessary in the interpretation. The major difference relates to the flow geometry. Once the flow geometry has been established other features such as multiphase flow, skin effect, afterflow, etc., well described in the literature, can be introduced. It will be assumed that the mechanical operation of the formation tester is already known to the reader.6 t will suffice here merely to state that the tester provides the means for taking a relatively small sample of the fluid immediately adjacent to the borehole, and for recording the subsequent pressure response. In comparison with conventional drillstem tests, the time required for a satisfactory pressure build-up response is much shorter, because of the relatively small quantity of fluid withdrawn by the wireline tester. This feature is highly desirable in the case of low-permeability formations. For an analysis of the pressure response within the formation, three simple flow geometries are considered— linear, cylindrical and spherical. The spherical and cylindrical flow geometries are most pertinent to the formation tester; therefore, they will receive the major emphasis. Since the configuration of the borehole and the perforation made by the tester complicate the flow geometry, it is necessary to allow for them in the drawdown response. However, because of the volume of formations contributing to the pressure-response, the details of the perforation shape are unimportant in the build-up period. Since relatively small amounts of fluid are withdrawn from the formation, in contrast to a conventional drill-stem test, a study of the "depth of investigation" and the significance of drawdown as well as build-up data will be included. Because the "depth of investigation" will be shown to be rather large, the effect on the build-up curves of the finite thickness of the permeable bed is considered. It is this consideration that leads to the importance of cylindrical flow geometry. Also included is a discussion of permeability anisotropy and its effect on the interpretation of the tester results. The pressure curves recorded by the formation tester will follow two general patterns, depending upon whether the formation is of high or low permeability. Fig. I (a and b) schematically illustrates these two responses. In Fig. 1(a), the high pressure recorded during fill-up of the tool is essentially the pressure differential across the choke in the system. In Fig. l(b), the flow rate is

By G. C. Broadhead

The lead district of Southeast Missouri covers an area of over 3000 square miles, including Maries County on the west, Jefferson on the east, Franklin on the north, and part of Madison on the south, or parts of ten counties. A general section of the rocks of the southeast part of this region would be about as follows, numbering from the top: 1. Sandstone (the second of Missouri geologists),.... 20 feet. 2. Chert beds—beds of passage below,...... 126 " 3. Magnesian limestone, chert, and quartzite,. .. 100 to 300 " 4. Lower magnesian limestone,......100 to 150 " 5. Gritstone and lingula beds,....... 50 " 6. Ozark marble beds,........5 to 20 " 7. Sandstone and conglomerate,. .. .. 5 to 90 " 8. Porphyry, Archæan 9. Granite, The thickness is approximate, for it is not possible to obtain the correct thickness of the various beds. No. 1 may be seen as detached outliers or loose masses of a hard sandstone or quartzite, and is found on the hilltops in the southern part of Madison County, in like topographical position in Reynolds County and the western part of Washington County. No. 2 consists of alternations of chert, clay, and quartzite. It is the formation which contains most of the limonite deposits of Central and Southern Missouri. In Reynolds County shafts have penetrated it 75 feet. Outcrops there and in Madison would indicate it to be at least 125 feet thick. Fossils obtained in various places prove it clearly to be of the age of the calciferous saudrock of the New York system, and Nos. 1 and 2 may also be referred to the age of the second sandstone of Missouri. No. 3 consists chiefly of thick beds of ash, drab and flesh-colored magnesian limestone, both coarse and fine-grained. These beds, where they occur, in Washington, Jefferson, and the southern part of Madison Counties, contain a good deal of quartz in drusy cavities. In Washington County the limestones are beautifully ramified by a system of connected drusy cavities generally lined with minute quartz crystals arranged in botryoidal and mammillary forms, and called " mineral blossom " by the miners. This is the lead-bearing rock of Washington County, but in the southern part of Madison it does not appear to be galeniferous. It is also undoubtedly the equivalent of

By C. P. Lawhon, J. P. Simpson, W. M. Evans

Data obtained under controlled test conditions using a microbit drilling machine showed that oil emulsified in water muds may either increase or decrease the drilling rate, depending upon drilling conditions. A low-viscosity oil such as diesel fuel can give drilling rates in limestone almost equal to that of water. Data obtained for water emulsified in oil muds showed little decrease in the drilling rate in water-saturated cores as the water percentage of the mud was increased above the 5- to 10-percent range. Changes in drilling rate were found to be dependent upon the oil or water concentration of the mud and upon the type of formation drilled. Changes in static filtration on paper (API filtrate) did not correlate with filtration while the mud was circulated across rock. INTRODUCTION Oil additions to water muds have been reported to increase drilling rates, provide hole stability and improve filtration control. Eckel' showed that water-base emulsion muds used in the West Texas area increased drilling rate with increasing oil concentration up to 15 percent oil by volume, but drilling rate decreased at a concentration of 20 percent by volume. Based on laboratory tests using water muds to drill shale, Cunningham and Goins' reported that drilling rates increased and tendency for the bit to ballup decreased with the addition of oil. Percentage increase in drilling rate varied with the particular formation. They showed oil additions to improve drilling rates ap proximately 75 percent in Vicksburg shale and as much as 150 percent in Miocene shale. Each investigation showed an optimum oil content for the particular formation. Most data that indicated improved filtration control due to oil additions were based on static API Eltrates through paper rather than dynamic filtration through permeable rocks. Some types of dynamic test give a better representation of filtration down-hole while drilling and might be more likely to show some correlation with drilling rate. Static filtration would be important, of course, in relation to hole stability and formation damage. This laboratory's drilling tests, conducted on water-raturated Berea sandstone, indicated that improvements in drilling rate were not evident with increasing oil concentration in water-base muds. Investigation also showed similarity between oil-emulsion (water-in-oil) muds and water-emulsion (oil-in-water) muds while drilling these formations. In Lueders limestone high concentration of water-in-oil muds and high concentration of oil-in-water muds provided the same relative drilling rates. In Berea sandstone there was a large reduction in relative drilling rate with both the oil and water muds that contained low percentages of emulsified fluid. Dynamic filtration rates of water muds on rock did not always decrease with increasing oil percentages even though the static API filtration rates on paper did decrease. Data observed in laboratory drilling of limestone and sandstone indicate that improvements in field drilling operations when water- or oil-emulsion muds are wed may not be the result of increased drilling rates but of improved hole conditions. In some cases, actual drilling rates might be slower but improved hole conditions will result in less total time on the hole. DEFINITIONS Mud pressure—Pressure of drilling fluid as measured after leaving the drilling chamber. This is considered as the approximate mud pressure just past the bit and at the face of the formation. Terrastatic pressure—Pressure representing weight of overburden. Formation pressure—Pressure of formation fluid as measured at outlet of drilling chamber. This is considered as approximate pressure of fluid in the pores of the formation. Differential pressure—Difference between the mud pressure and formation pressure. Relative drilling rate, percent—Drilling rate with experimental fluid divided by drilling rate with water times 100 equals percent. LABORATORY EQUiPMENT AND TESTING PROCEDURES The drilling equipment has been described in previous publications The microbit drill is a closed system (capacity, approximately 7 gal) that can be pressurized to 15.000 psi and heated to 500F. Main components are a drilling chamber, filter-heater, rotary-drive and variable-speed cir-

By R. P. Rothwell

A Reference to the geological map of Alabama shows the Goalmeasures of that State to form three distinct fields. The Coosu, or most easterly, contains about one hundred square miles ; the Cahaba, or middle field, which is also the most southern true coal in the United States, contains about 230 square miles ; and the Warrior field, which contains in the State of Alabama some 5000 square miles, is the southern extremity of the great carboniferous deposit, which extends through Pennsylvania, West Virginia, Kentucky, Tennessee, and Georgia. But very little has yet been done towards developing those coalfields, partly owing to the absence of all commercial manufacturing enterprise in the South under slavery, and partly owing to the want of capital and the disturbed condition of the South since the war. During the past three or four years I have devoted a large part of my time to the examination of the coal and iron ores of this range,

By John Henry Harden

It appears that in the earlier days of anthracite coal mining, 1824-25, the Lehigh region mined 76 per cent. of all the coal sent to market. During the same period Wyoming sent 12 and 5 per cent. respectively, whilst the Schuylkill sent 11 and 18 per cent., the number of tons marketed being in those days very small indeed. In the year 1828 the Schuylkill region mined over one-half of all the coal sent to market, and continued to do so for thirty years (with the exception of the year 1831), until 1858 ; since that date the production has declined very steadily in comparison with the other regions, until, in 1876, it was only 34 per cent. Wyoming, in the year 1820, was mining 40 per cent., but the production of all the regions was very small; in six years it declined to 5 per cent., increased to 34 per cent. in 1831, declined to 10 per cent. in 1838. Since which time the increase has been gradual, until 1873-75 Wyoming produced over 50 per cent. of all the anthracite coal sent to market. In 1867 the production exceeded that of Schuylkill, and has continued to do so since that time. In the year 1876, during the existence of the "coal combination," Wyoming lost 9 per cent.; during the same time Lehigh gained 7 and Schuylkill 2 per cent. over their relative production of the preceding year, the decrease of Wyoming being two and a half million tons, the increase of the Lehigh one million, the decrease of the Schuylkill nearly three-quarters of a million. The extraordinary production of the Lehigh region at a time of decreased demand and in opposition to an agreement with the other regions to limit the supply to the supposed demand mas the cause of the disagreement which ultimately broke up the combination of the great coal-producing and carrying companies. The Lehigh region produced very steadily, from 1829 to 1854, an average of about 22 per cent. of the whole production. From 1854