Search Documents

By D. H. Gray, G. Bergamini, I. Fatt

A method developed by soil scientists for measuring permeability anisotropy of soils has been adapted for use in studies of sandstone permeability aniso-tmpy. Studies on Berea and Bandera sandstone show that permeability of these materials decreases upon application of simulated overburden pressure. Permeability reduction of Berea and Grubb sandstone cores subjected to simulated overburden pressure is shown to be a function of the ratio of radial to axial stress. INTRODUCTION The effects of mechanical stress on reservoir rock parameters are now generally accepted as being important in reservoir engineering calculations. In several oil fields the unique behavior of the reservoir has been attributed to the influence of overburden pressure or internal fluid pressure on rock permeability or fluid content. l Although the gross effects of overburden pressure on permeability, porosity, waterflood efficiency, electrical resistance and sonic velocity have now been established,3-9 there are some very interesting, and perhaps important, details of rock behavior under stress that have not heretofore been studied. In this paper we will show that permeability anistropy of several sandstones is a function of overburden pressure. We will also show that permeability reduction of cylindrical samples of two sandstones, when subjected to mechanical stress, is a function of the ratio of radial to axial stress. Most studies of overburden pressure effects made in the past have attempted to duplicate reservoir conditions by applying hydrostatic external pressure on all surfaces of the sample. Assuming perfect elasticity, that is, no viscoelastic flow during geological time, reservoir rock in situ cannot be under hydrostatic stress as is apparent from the value of 0.25 for Poisson's ratio for most rock. In the absence of tectonic forces the horizontal stress component If there are directional permeability effects, as has been reported from field and laboratory studies, l0,11 then one would expect hydrostatic stress as used in the laboratory to give a different permeability reduction from that occurring under reservoir stress conditions. The influence of stress on directional permeability was measured by a modification of the method first described by Maasland and Kirkham for use on soil samples.12 Their method permits calculation of the permeabilities along each of the three mutually perpendicular axes of a cubic sample. The modification reported in this paper extends the method to observing these permeabilities when the sample is under hydrostatic stress. The Hassler sleeve type of core holder is satisfactory for studying overburden pressure effects on permeability if a known force is applied to the end plugs. Wyble6 used a sleeve core holder but did not load the end plugs. In this paper we combine sleeve pressure and axial force to give permeability changes which we believe are close to those occurring in reservoirs when reservoir pressure changes. The combination of sleeve pressure and axial force assumes, however, that the sleeve pressure simulates both vertical stress and horizontal stress in one direction. METHOD PERMEABILITY ANISOTROPY The Maasland and Kirkhaml2 method of measuring permeability anisotropy uses a cube of porous material which is Open On three faces and has discs with a Center hole for an air inlet cemented to the other three mutually perpendicular faces (Fig. 1). The rate at which air pressure falls in a bottle of known volume, and pressure slightly above atmospheric, when the bottle discharges through each disc, is

By L. R. Raymond, G. G. Binder

One primary purpose of hydraulic fracturing as a well stimulation technique is to overcome formation damage. The literature provides ways of designing fracture treatments and evaluating their results but not of incorporating formation damage in vertically fractured wells. Results of an investigation of this problem are presented in this paper. Prediction of stimulation ratios in vertically fractured, damaged wells is accomplished with a mathematical model relating stimulation ratio to relative conductivity of fractures whose lengths are not more than about half the drainage radius of the well. Comparison of results from the new model to results in published predictions verify its utility; these results also show that the range of stimulation ratios that can be predicted for undamaged wells is extended to include relative conductivities of less than 300. This extension is important when using fracturing to stimulate wells with high production rates and high native formation permeabilities. For example, large increases in daily oil production rate can be obtained with stimulation ratio increases as low as 1.25. The importance of complete fracture fill-up (uniform proppant packing) is shown through use of the mathematical model. If, at the mouth of a fracture, only a small fraction (1/2 percent) of the total fracture length is not packed with proppant, nearly all the polential stimulation increase is lost. Proppant crushing, compaction and embedment have been shown in laboratory studies to be responsible for low fracture conductivities in wells producing from highly stressed formations. Equipment and methods for testing the effect of stress (overburden) on conductivity of fructures in consolidated and unconsolidated sands are discussed in this paper. Laboratory tests with simlilated fractures in cores from both types of formations showed that crushing, compaction and embedment seriously affect conductivity. Results indicate that similar laboratory tests should be made when accurate knowledge of fracture conductivity is needed to assure good stimulation results for important wells. The chief factor in stimulation ratio reduction was found to be overburden pressure, but the size and type of proppant and the hardness of the formation have significant effects. Fracture conductivity reductions of up to 50 percent were observed with sand propping fractures in consolidated cores; a reduction of 83 percent was measured for an unconsolidated core. The range of effective overburden pressures for which conductivities were measured was from 100 to 5,000 psi. An example fracture design and evaluation problem indicates the usefulness of considering formation damage in planning well stimulation jobs. When damage exists, but its extent and the degree of permeability reduction are not estimated from diagnostic tests, the formation permeability used in planning the job may lead to under-designing. As the example shows, too low a target stimulation ratio can lead to much lower production rates (by half) than could be attained otherwise. Solutions of equations representing several special cases of the mathematical model are presented in graphical form and details of the derivations of the equations are given in the Appendix. INTRODUCTION Since its inception in 1947, hydraulic fracturing has proven to be an effective and widely accepted stimulation technique. In the past 18 years the ability to execute a successful hydraulic fracturing treatment has been substantially increased. The development of design and evaluation procedures1,2 has been one of the major contributions to this increased skill. However, as the art of hydraulic fracturing has moved closer to a science, new problems concerning the design and evaluation of the optimal hydraulic fracturing treatment have arisen. Three questions are pertinent to these problems. I. How is a fracturing job evaluated in a damaged well? 2. What is the effect on the stimulation ratio of not filling the fracture in the vicinity of the wellbore? 3. What is the effect of overburden pressure on fracture conductivity and, consequently, the stimulation ratio? A primary objective of fracturing a well is to stimulate production by overcoming wellbore damage. Presently. however, there is no rational basis for designing or evaluating the optimal fracturing treatment in a damaged well. All present fracture design and evaluation techniques assume that proppants can be uniformly placed in fractures. This assumption may be in serious error, particularly for the portion of a fracture directly adjacent to the wellbore. In this area, turbulence of the injected fluid can cause the proppant to be swept farther into the fracture. Also, loss of fluid from the fracture to the

By F. Selig, A. S. Odeh

A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well's shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner's' theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner's widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs' analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner's' method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner's method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner's original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner's approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink

By Albert J. Phillips

As early as 1824, Haidinger' described crystals of native copper that were, according to Dana,2 " probably twinned parallel to the octahedral plane and normal to this axis." In 1837, Rose3 very clearly demonstrated that native copper crystals were commonly twinned on the octahedral plane. Lasaulx,4 in 1882, described polysynthetic octahedral twins and stated that enclosed twinning lamellae were observed. Even on crystals formed artificially by electrolysis, twinning was observed as early as 1886, by Brown.5 Many other descriptions of twinning in native copper and in gold and silver are to be found in the literature, and in every case with which we are acquainted the twinning is conformable to the spinel law; i. e., the octahedral plane is both twinning and composition plane. Metallurgists have long been familiar with the straight-sided bars or bands that commonly occur within the irregularly outlined grains of metallic crystalline aggregates. They are certainly found in nearly all metals of the face-centered cubic type after cold working and annealing. If a specimen showing these bands is rotated, keeping the angle of incidence of the illuminating beam constant, they usually exchange color contrast periodically with their respective matrix crystals. Also slip lines, when produced in these specimens, nearly always show a change in direction upon passing through the bands, giving rise to the "stairway" effect often illustrated in metallurgical textbooks. Apparently the only quantitative observations on bands of this sort were made by Osmond and Cartaud6 on austenitic steel. They selected,

By W. B. Gogarty

With the use of polymer solutions in secondary recovery operations, the need has developed to understand the mobility control mechanism. This study investigated mobility control by considering both permeability and rheological effects. Experiments used a high molecular weight, partially hydrolyzed polyacrylamide polymer. Flow studies took place in reservoir and Berea cores having zero oil saturation. Effective size of the polymer flow unit was inferred from Nuclepore filter tests. Clay studies indicated the particle size capable of decreasing the core permeability. Flushed permeabilities measured the approximate core permeabilities with flowing polymer solutions. These permeabilities were considerably lower than original values. With mobility data and the flushed permeability, maximum effective viscosities were determined for polymer solution flow in a core. Effective viscosities showed that rheological properties play an important part in mobility control with polymer solutions. The study showed that permeabilities decrease and stabilize with polymer flow. At the lower permeabilities, high shear rates exist in the cores. Because of the pseudoplastic character of the polymer solution, the high shear rates caused low effective viscosities. This condition pointed to the inefficient use of the potentially high viscosity of the polymer solution at low shear rates. INTRODUCTION In the oil industry, a great deal of interest is being shown in the use of polymer solutions for secondary recovery and a number of polymer floods are being performed in the United States.1 Some of these floods have become commercial while others have been reported as failures. A number of floods are still in progress and remain to be evaluated. With the advent of polymer flooding, the need developed to understand the mobility control mechanism in porous media. Rheologically, the polymer solutions behave as pseudoplastic fluids. Investigators have studied the rheology for this type of non-Newtonian fluid in porous media?-6 Results have also been reported on the ability of polymer solutions to decrease mobility.7-10 Some investigators indicated that polymer solutions show a different rheological behavior in cores than on the viscometer.7,8 This study is representative of an initiatory investigation on mobility control with polymer solutions. Some effects are considered on a qualitative basis. Mobility control is considered from the standpoint of both permeability and rheology and, to the extent possible, permeability and rheological effects are considered separately. Results show that mobility control occurs from a reduced permeability caused by polymer retention and an effective viscosity determined by the average shear rate in the core. In this respect, rheological behavior is as important in mobility control as permeability. PROCEDURES A high molecular weight, partially hydrolyzed, polyacrylamide polymer was used in the investigation. Solutions of various concentrations were made with this polymer in different types of water. Some solutions used distilled water and others used tap water. The tap water contained approximately 500 ppm dissolved solids. As used in this paper, concentrations in ppm refer to parts by weight solute to one million parts by weight of solvent. In the tap water, calcium and magnesium accounted for from 60 to 80 ppm of the solids, the remainder being mostly sodium chloride. All solutions contained a bactericide, and anti-oxidants and antiprecipitants were added to the tap water. Capillary and concentric-cylinder* viscometers determined the rheological properties of solutions. The capillary viscometer is described elsewhere.6 Berea and reservoir sandstone were used as core materials for this investigation. Cores were 1 in. in diameter by approximately 2 to 3 in. long. Prior to mounting, the Berea cores were heat-treated for 3 hours at 850F. This treatment decreased the water-sensitivity of the clays in the cores. The reservoir cores were cleaned with solvent and

By J. M. Smith, G. P. Willhite, J. S. Dranoff

Heat transfer rates were measured in sandstones with flow of gases perpendicular to the direction of energy transfer. Effective thermal conductivities ker ranged from 0.7 to 1.7 Btu/(br)(ft)(°F). The contribution of the solid phase appeared to be the most important in these consolidated materials, although the thermal conductivity of the gas had some effect. The velocity of the gas through the pores of the sandstones had no influence upon ker up to values of 168 lb/(br) (sq ft) in agreement with data obtained for unconsolidated beds of glass beads. The present results indicated that gas mixing, and hence heat transfer by convection in the pores, is less for perpendicular transfer of energy than when fluid flow and energy transfer are in the same direction. HEAT TRANSFER PERPENDICULAR TO FLUID FLOW IN POROUS ROCKS Heat transfer in porous media with pore sizes in the micron range depends upon the fluid in the pores and the geometry of the solid phase. The best characterized system is a bed of solid spherical particles. Heat transfer in this system has been studied extensive1y (Ref. 8 summarizes the literature up to 1959) when the pores contain stagnant fluid. When the fluid is in motion the directions of flow and energy transfer have an effect on the heat transfer rates, as demonstrated by comparing the work of Willhite, et al l3 for perpendicular flow and that of Kunii and Smith9 for parallel flow of energy and fluid. For perpendicular flow, no increase in effective thermal conductivity ke was noted up to mass velocities G of 77 lb/(hr) (sq ft). In contrast, for parallel flow ke, increased with G in the same range of flow rates. These higher values of ke in the direction of flow also have been observed3,4,10 in beds of larger particles, 0.1- to 0.5-in. diameter. For beds of consolidated materials, such as porous rocks, data are not available for these compar- isons, although Adivarahan1 reported results for the parallel case. Hence the primary objective of this work was to measure ke values for perpendicular flow of fluid and energy in porous rocks. Of interest also was the variation in effective conductivity with fluid velocity. APPARATUS AND PROCEDURE In the experimental method, a constant heat flux was applied to the inner surface of an annular section of the porous rock. By cooling the outer wall, a temperature gradient through the annular sample was established and measured with thermocouples placed within the sample at various radial positions, and at three elevations (A,B.C). The location of the 2-in. O.D., 3.75-in. long sample in the apparatus is shown in Fig. 1. Fluid entered the bottom (1) of the 3-in. I.D. (approximate) steel shell, flowed upwards through the sample and out at the top (3). pressure taps (2,4) were used to check the permeability of the sample. The energy flowed radially from the centrally-located electric heater through the sample and was absorbed in the water-cooled jacket. The samples studied were naturally occurring sandstones from different locations with the properties given in Table 1. These materials are identical with those used by Adivarahan 1 for the parallel flow of energy and fluid. Prior to use they were refluxed with toluene to remove hydrocarbons and leached with distilled water to remove soluble salts. Each sample was visually examined and discarded if large nonhomogeneities, such as cracks or stone particles, were noted. Eight copper-constantan thermocouples were inserted in holes drilled radially into the sample with

By J. W. Marx, R. H. Langenheim

The injection of heat-bearing fluids may offer a wider application to secondary and tertiary recovery from conventional oil reservoirs than inderground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. In all cares, conduction heat losses to overburden and un-derburden will impose an economic limit upon the size of the area which can be swept out from any one injection point for any given set of conditions. The present report describes a method for c.rtin~ating thermal invasion rates, cumulative heated area, and theoretical economic limits for sustained hot-fluid injection at a constant rate into an idealized reservoir. Full allow-ance is made for non-productive reservoir heat losses. Required solutions can be obtained directly from the numerical tables given here, as soon as all operating conditions are specified, withou extensive mathematica1 ~nanipulation on the part of the user. The injection of a heat-bearing fluid into a reservoir is frequently proposed as a means of secondary or tertiary oil production. This form of thermal recovery may offer a wider application to conventional oil reservoirs than underground combustion since the process is more easily controlled and the reservoir requirements, in general, are less critical. Because of its large gross heat capacity, steam appears to be the most efficient heat injection medium, although mixtures of steam and other gases, hot water, hot oil and hot non-condensable gases also have been employed. Whatever the nature of the injection medium, economic evaluation of such thermal recovery processes will depend upon thermal invasion rates, or rates of productive reservoir heating, at any given time after the start of heat injection. In all cases, conduction heat losses to overburden and underburden will impose an economic limit upon the size of the area which can be swept out from any one injection point, for a given set of reservoir conditions, at any given heat injection rate. The present report describes a method for estimating thermal invasion rates, cumulative heated area and theoretical economic limits for sustained heat injection at a constant rate into an idealized reservoir, making full allowance for non-productive reservoir heat losses. Required solutions can be obtained directly from the nu- merical tables given here, as soon as all operating conditions are specified, without extensive mathematical manipulation. It should be noted that ordinary water flooding may be used to move an annular-heated zone outward, even though the economic areal limit for straight heat injection has been reached. Thus well spacing need not be confined to the areal limits as determined here. However, the present calculations do not cover such supplementary heat scavenging programs. The operating element consists of a radial flow system, concentric about the point of injection, with the idealized step function temperature profile illustrated by the dashed line in Fig. 1. The solid curve shown in Fig. 1 gives a qualitative picture of the true radial temperature distribution. It will be shown later that the expressions given below, which are derived from the simplifying step idealization (dashed line in Fig. I), also can be interpreted in terms of the more realistic temperature distribution shown in Fig. 2. For this case, the thermal invasion radius is simply redefined as the distance from injection well to the mid-point of the temperature distribution.

By R. E. Jenkins

The Bell Canyon member of the Delaware Mountain group has yielded quite a large number of fields in which completion and production problems have been numerous and complex. Reserves are difficult to estimate due to the problem of evaluating the formation water saturation, and the feasibility of water flooding is questionable. Extensive laboratory investigations were undertaken to determine if normally measured rock and contained fluids properties could account for many of the peculiar characteristics of the formation. A statistical study or routine core-analysis data was made. Capillary-pressure data were obtained by several techniques. Fresh-water and brine permeabilities were measured and relative permeabilities of oil, water and gas were determined. Some qualitative wetta-bility tests were performed, and the petrography of a few thin sections was studied. The results of this work are presented, with typical or average data shown for each test. The permeability -to- water characteristic of the Delaware was found to be abnormally poor. Correlations of laboratory and field performance data indicate high water saturations in much of the formation. Pore geometry is highly uniform. INTRODUCTION The Delaware formation has undergone extensive exploration and development in the past few years. This has resulted in the discovery of quite a number of fields—producing oil for the most part. A large percentage of the productive wells have produced water along with the oil. Often there is no apparent pattern for the percentage of water cut experienced. In several instances vertical displacments of over 200 ft have been established between communicating wells in which the formation appears to be about the same, yet the structurally high wells yield a much higher water cut than the lower wells. Drill-stem tests have contributed little to well completion because they nearly always show small amounts of water (usually mud filtrate), regardless of what is produced later in the life of the well. Other formation evaluation methods have had no more than limited success in indicating ultimate productivity. Consequently, predicting the type of fluid productivity of Delaware wells has been very difficult. Hydraulic fracture treatments are almost universal, and unstimulated Delaware wells have not exhibited the productivity that normally would be expected of wells with the specific permeabilities encountered. Some operators have indicated that their Delaware wells need frequent re-stimulation. The uncertainty of reservoir water saturations reported for the Delaware has made difficult the task of evaluating oil in place in the reservoirs and estimating primary recovery. The unreliability of this data also will cause feasibility studies for water flooding or other secondary-recovery methods to be problematical. Several opinions have been voiced to explain why the Delaware performs as it does. Chronically poor well completions, dynamic water and tilted water tables, capillary inequi-libria, peculiar relative permeability characteristics and oil-wetness are among the suggested reasons. The dual purpose of this paper is (1) to describe the principal rock characteristics, and (2) to present an evaluation of those characteristics so that approaches can be made to explain the irrational behavior of wells completed in the Delaware. STUDY OF ROUTINE CORE-ANALYSIS DATA A statistical study of routine core-analysis data has been made to relate fundamental rock properties and characteristics of the formation to well productivity. For the most part, well production data were obtained from operators and augmented by production data from commercial information services. Only data which were available in detail and which were considered to be accurate were used. The core data used were obtained on samples from the perforated or open interval only and, thus, may not always represent the total productive interval in a well. All the core-analysis data included in this study are from the Bell Canyon member of the Delaware and were obtained by "plug-type" techniques, as opposed to whole-core or full-diameter techniques. Porosity values were obtained by the summation-of-fluids technique, permeability values are air permeabilities corrected for slippage to an equivalent liquid permeability, and the oil and water contents were determined by means of a high-temperature downdraft retort. All the procedures used are discussed at length in Ref. I. The core-analysis data were sorted

By J. B. Cheatham, R. B. Paslay

Rock stresses and steady-state flow rates induced by the pressure gradient associated with the flow of formation fluid into a borehole have been analytically determined for a permeable, elastic material saturated with an incompressible fluid. In this analysis, the material properties and loading are considered to be symmetric about the axis of the borehole and independent of axial position. For Case I the material is assumed to have uniform permeability in the radial direction, whereas for Case 11 the permeability is assumed to have been reduced in a localized region adjacent to the hole by either normal well completion and production operations or deliberate plugging during air drilling. Results of a numerical example indicate that, in the absence of plugging, the rock shear strength must be approximately two-thirds the formation fluid pressure in order to prevent rock failure. The required rock strength is high for small radial zones of plugging and decreases as the region of reduced permeability becomes larger; however, a depth of plugging can be reached beyond which there is no real gain in strength, although the flow rate can be further reduced. INTRODUCTION During normal production of oil from a well, it is often desirable to increase the production rate of the formation fluid by increasing the pressure gradient through the formation adjacent to the borehole. Depending upon the magnitude of this pressure gradient and strength of the rock material, this production-rate increase can cause sloughing of the hole wall. In many cases, the production - rate increase can result in excessive sand production, increased wear of production equipment, lost production time and expensive workover jobs. In addition, the phenomenon of increased rock bit penetration rate with the use of a gaseous instead of a liquid drilling fluid has been observed in oilfield drilling operations and experimentally demonstrated by various investigators for several years. The improvement obtained by employing this technique can be quite significant and offers apromising method for reducing drilling costs. However, air drilling is currently limited to geographical locations where high-capacity water-bearing formations are not encountered. This limitation has prevented widespread adoption of air-drilling techniques, because the water influx into the borehole interferes with efficient removal of the drilling cuttings and usually results in a condition such that the bit becomes "balled-up" or stuck in the hole. In an attempt to remove the water - intrusion limitation from air drilling, various chemical and mechanical water shut-off methods have been proposed. Goodwin and Teplitzl suggested one such proposal whereby the permeability of the water-bearing rock structure was reduced in the vicinity of the borehole. Although the development of a shut-off method based upon this approach would certainly be welcomed by the oil industry, it is conceivable that, under certain conditions of the pressure gradient, strength of the rock material and depth of the modified permeability zone, a stress field can be created that will result in an unstable hole. As part of their study, an analytical solution is given for stresses in an idealized model of a hole and the surrounding rock material. The purpose of the present study is to extend the analysis of Goodwin and Teplitz to gain more insight into the details and consequences of excessive production rates and formation water shut-off. In particular, simplified models of these problems have been analytically examined, which makes possible the evaluation of the type of stress fields that can be anticipated as a result of these production and drilling practices. Both problems solved concern the determination of the steady-state volume flow rate of the formation fluid and the resulting steady-state stress and displacement distribution in a hollow, cylindrical configuration. The cylinder of Case I, corresponding to the production-rate problem, consists of a material with a constant permeability from the inside surface to the outside surface; the cylinder of Case 11, corresponding to the water shut-off problem, consists of a material with a constant permeability from the inside surface to an intermediate concentric cylindrical surface and a second constant permeability from the intermediate surface to the outside surface-

By H. D. Outmans

A method has been developed for determining the relative water wet-tability (fraction of the surface wet by water) of porous media. This method involves the adsorption of methylene blue dye from an aqueous solution onto the solid surfaces of rock contacted by the injected dye solution. Experiments on Berea snndstone cores have shown that water-wet cores adsorb a large amount of dye from solution. In oil-wet cores, however, the injected aqueous dye solution is prevented from contacting the adsorptive surfuces by a film of oil and virtually no dye is adsorbed. This method was tested by experiments on mixtures of water-wet and oil-wet sand. Using the dye adsorption method, it was possible to determine quantitatively the fraction of water-wet and oil-wet sarzd in such mixtures. The dye adsorption test is readily applied to determinitrg the relative water wettability of fresh oilfield cores. The cores are not extracted before testing. Observed values of the wettability of field cores tested to date have varied from 100 to 6 per cent water wet. INTR OD UCTIO N Recent articles1 ' have emphasized the important role that wettability plays in almost all phases of reservoir rock behavior. It has been shown that wettability has a profound effect on the interstitial water saturation, residual oil saturation, capillary pressure, relative permeability, waterflood behavior, and resistivity index of oilfield cores. A quantitative investiga- tion of the effect of wettability on these core properties is severely handicapped by the lack of a simple technique for measuring the wettability of a porous medium. If such a technique were available, it would provide investigators with a means of gaining insight into the distribution of phases within a core. A number of qualitative methods of estimating wettability are known. These include contact angle measurements, flotation tests and studies of displacement pressure: imbibition7 and relative permeability data.' Using these tests, it is usually possible to classify a core as water-wet, oil-wet, or some intermediate degree of wettability. One recently proposed method involving measurement of nuclear magnetic relaxatione offers promise of determining quantitatively the fractional wettability of cores. The equipment is rather elaborate and necessary treatment of the cores before testing probably alters the natural wettability. This latter point also applies to most of the cited qualitative methods. This paper presents a new technique for determining water wettability of cores; it is based on measurements of the dye adsorption capacity of the core. This method has the advantage of requiring relatively simple equipment and involves minimum handling of the core before testing. It is important to distinguish between the various concepts of wettability referred to in the literature. The term "wetting" is defined by Bartell as that phenomenon which occurs when a solid and liquid phase come in contact in any manner so as to form a solid-liquid interface. A preferentially water-wet solid is defined by Jennings in terms of the advancing contact angle in the solid water-oil systems. Brown and Fatt regard fractional wettability as the fraction of the total surface area that is preferentially oil-wet or water-wet. However, none of these definitions adequately describes the type wettability measured by the dye test. This calls for the introduction of a new term, "relative water wettability", which we define as the fraction of the total surface area of a core contacted by injected water. DESCRIPTION OF METHOD The ability of a porous medium (especially the clays) to adsorb large quantities of surfactants, dyes, etc., is well known. It is reasonable to assume that a substance can be adsorbed from solution onto a solid surface only if the solution contacts the surface. In this connection Shapiro" has shown by experiments on silica gel that completely dry gel will adsorb methyl red dye from a benzene solution. On the other hand, if a film of water more than several molecules thick is present on the silica gel, the adsorption capacity of the gel for the dye is drastically reduced. Thus, a study of the dye adsorption capacity of cores should give some clue as to the fluid distribution on the surface of the core. TEST OF METHOD As a qualitative check on this theory, dye adsorption tests were run on similar Berea sandstone cores (Table 1, Cores 1-5). Cores saturated initially with water and flooded with 0.01 per cent aqueous methylene blue adsorbed 0.58 mg dye/gm of core. Similar cores saturated initially with water, driven to irreducible water with oil, and then flooded with the same dye solution also adsorbed 0.58 mg dye/gm of core. Another core, treated with a 5 per cent Dri-Film solution, dried and saturated with oil (Soltrol) adsorbed no dye when the dye solution was injected. These tests confirm

By D. M. Beeson, G. D. Ortloff

Water-propelled banks of carbon dioxide recovered both high- and low-viscosity crude oil substantially in excess of that recovered by water flood in linear flow model experiments. The increase in oil recovery is explained by oil swelling and viscosity reduction. Additional oil was recovered by solution gas drive at the end of the carbon dioxide floods when the model systems were depres-sured. This study was an investigation of the use of a water-driven carbon dioxide bank to recover crude oil. Techniques involving the injection of carbon dioxide have aroused con-siderable interest in recent years among methods considered for increasing oil recoveries. A number of investigators have studied the process in the laboratory1-' and in the field.7,5 Laboratory studies referenced indicate that carbon dioxide in conjunction with water flooding recovers more oil than does water flooding alone. These re- ports diverge widely, however, in stating magnitude of improvement. We believed that a water-driven carbon dioxide bank would increase recovery of crude oils by two mechanisms: viscosity reduction and oil swelling. Due to reduction in the oil viscosity caused by solution of carbon dioxide, displacement efficiency should be improved. In addition, the residual oil should be highly swollen by the carbon dioxide held in solution. Thus, the final residual oil saturation, in terms of stock-tank oil, would be lower than in the absence of carbon dioxide. To test this hypothesis, experi-

By F. H. Deily, D. S. Rowley, R. J. Howe

By D. A. Flanagan, P. B. Crawford

A study has been made of the feasibility of storing liquid meihane at low pressures in undergrohd caverns. Methane liquefies at — 258°F at atmospheric pressure. It is shown that the methane evaporation rates will rapidly decrease and cool the surrounding rock so that at the end of one month they would be between 50 Mcf/hr for a cavern of 25-ft radius and 700 Mcf/hr for a cavern of 100-ft radius. At the end of 10 years, the evaporation rates would be 18 and 100 Mcf/hr, respectively, for caverns of the same radii. The evaporation rates may be reduced by a factor of 2 to 10 by the application of insulation. The cost for mining the caverns is estimated to be $.75 to $1.25/Mcf of storage. This is substantially less than surface storage; it is believed to be safer and to result in lower maintenance, savings in space and savings in strategic materials. INTRODUCTION During the past few years, there has been an increasing interest in the economic feasibility of liquefying methane. Methane liquefaction is being considered for tanker transport; in addition, liquefaction is being reconsidered for shaving peak gas demands. Several articles have described natural gas liquefaction and the progress of the tanker in making trial runs from the United States to Great Britain to determine the feasibility of tanker transport. At the present time, pipelines are not designed to supply peak gas loads during extremely cold periods such as are oiten encountered in the North and Northeast. Gas is being stored in underground reservoirs en-route to its destination, but in many instances satisfactory storage in porous reservoirs has not been practical, especially along the Eastern seaboard where few petroleum reservoirs have been found. In England and other foreign countries, it is unlikely that satisfactory porous structures could be found, and it may be desirable to mine or excavate the rock to obtain storage. By storing the methane near the consumption point, product availability can be increased during periods of need. Methane liquefies at - 258OF at atmospheric pressure, and 1 cu ft of liquid methane will make about 600 cu ft of gaseous methane. Methane liquefaction for shaving peak gas demands was conducted in Cleveland, Ohio in the early 1940's. The Cleveland plant was designed to liquefy 4 MMcf/D of natural gas.' Development of several low-temperature processes has resulted in some improvements in the design of low-temperature storage vessels. Liquefied-methane storage vessels can be designed for various evaporation rates. The desired rate will depend on the reason for liquefaction. For shaving peak gas demands, the desired evaporation rates would be low, possibly between 0.2 and 1.5 per cent/day. For methane tanker transport and supply, the evaporation rates must be high and must equal the total daily rate of supply. It is apparent that the cost of storage facilities will depend on the objective and will vary appreciably. Surface storage of methane may cost between $3,000 and $20,000/MMscf, depending on the size and design. The failure of the surface storage vessels at Cleveland, with the resulting injury and loss of many lives, has caused considerable emphasis to be placed upon developing new, safer and improved methods of storing liquid methane. Storage of liquefied propane and butane in underground salt domes has been very satisfactory. The per-barrel cost of underground storage for any sizeable capacity is very small compared to the cost of surface storage for LPG mixtures. The larger the capacity of the underground storage, the less is the per-barrel cost. Numerous advantages appear to exist if liquefied methane could be satisfactorily stored underground, the principal advantages being increased safety, lower initial cost, lower maintenance cost, savings in space and savings in strategic materials. The purpose of this paper is to determine the feasibility of storing liquid methane underground in mined caverns. If a large spherical cavern is dug underground and is filled instantaneously with liquid methane, the surface of the sphere may approach the liquid methane temperature almost instantaneously. The temperature distribution for a sphere of radius r. in an infinite medium, initially at one temperature with spherical surface kept at T, from time t = 0, is given by Eq. 1.'

By Eben E. Olcott

The Mine.—For some years attention has been drawn to the copper deposits of the Appalachian range of mountains, and especially to those in that portion crossing the corners of Virginia, North Carolina, and Tennessee. Certain it is that in many places in this locality available deposits of copper await the miner, and one after another are being opened and worked. The parts of these three States alluded to are but sparsely settled, and communications are insufficient to meet the demand for transportation. Railroads are wanted, as the existing wagon roads are unfitted for heavy hauling; these circumstances increase very much the difficulties in the way of mining enterprises. Notwithstanding this, the success of some of the mines of this section is an established fact, and Ducktown has for some years held a high rank among the copper-producing localities of the United States. of all these deposits, however, none is more interesting or more deserving of attention than the one now under consideration. The works of the Ore Knob Copper Company have been energetically pushed forward, since the reopening of the mine, in the spring of 1873, by Mr. J. E. Clayt,on, the superintendent and one of the chief owners of the property, and if one thing more than another can insure success it is the combination in this gentleman of rare business capabilities with the deepest personal interest in the prosperity of the undertaking. The copper vein at Ore Knob* is a clearly defined fissure traversing the gneissic rocks of Pre-silurian age in a direction 61' east of north.

By Benjamin Smith Lyman

THE importance of topography to geology is so commonly underrated as to deserve to be pointed out again and again. The relation of topography to the different branches of geology may be seen best by a glance at them all. Geology is well understood to consist of the observation of facts about rocks as materials of the earth (geognosy), and the explanation of those facts (geognosy). The facts to observe are : I. GEOMETRICAL, as (1) structure, including the dip and strike at different points, the size, shape, and position of basins and saddles; (2) vertical section, that is, the thickness of rock groups and of their parts ; (3) horizontal or geographical position. 11. PHYSICAL, such as texture, whether sedimentary or crystalline ; and, if crystalline, whether granitic, porphyritic, amygdaloidal, globular, or concretionary ; whether compact or granular; coarse or fine; pervious or impervious; hard or soft; firm or tender; with cleavage or without ; massive or thin ; the microscopic structure ; also whether magnetic or not ; soluble or insoluble ; fusible or infusible ; and the constituent minerals. 111. CHEMICAL that is, the composition of the rocks and of their minerals and their affinities. IV. PALAENTOLOGICAL, as what kind of fossils are in the rocks,