Search Documents

By M. A. Selim, W. H. Somerton

Thermal expansions and heats of reaction of three typical sandstones were measured in the temperature range of 25°to 1,000°C. The significance of these data in subsurface heat-transfer calculations is discussed. Linear thermal expansions were measured both parallel and perpendicular to the bedding planes. Volume expansions are reported as the perpendicular expansion plus two times the parallel expansion. Expansion behavior of the sandstones was found to be controlled by the expansion characteristics of the quartz content. Differential expansion of the quartz grains and other minerals included in the sandstones caused permanent deformation of the heated samples after they were cooled to room temperature. Structural damage resulting from heating is probably an important cause of the reduction of thermal conductivity of heated samples. Measurements of heats of reaction were based on the known heat required for a- quartz inversion. Thermal reactions, which probably include deby-droxylation of clay minerals and decomposition of carbonate minerals contained within the samples, were found to require more than one-fourth of the amount of heat necessary to raise the temperature of the rock alone. In shales and limestones, the reaction heat could be substantially greater than that required from specific heat considerations alone. INTRODUCTION In earlier work, l-3 methods of measuring thermal properties were developed and thermal data for several sedimentary-rock types were reported. Data included thermal conductivity, specific heat and thermal diffusivity. In addition, calculations of heat transfer in subsurface formations may require data on the thermal-expansion characteristics of rock and the heats of reaction of mineral constituents. Thermal expansion of rocks is relatively small in magnitucle and, from the standpoint of change in bulk density, has only minor effects on heat-transfer cliaracteristics. Thermal-expansion behavior, however, may have significant effects on the structure of rocks. Differential expansion of different minerals and along different crystal-lographic axes may result in structural damage which could effectively reduce the thermal conductivity of the rock. Many mineral constituents of sedimentary rocks undergo phase or related changes when heated to sufficiently high temperatures. Although some of these reactions (such as the a-ß quartz inversion) are reversible and the heat absorbed is returned to the system upon cooling, many of the major reactions are irreversible. Through certain temperature ranges, the additional heat energy required to complete the reactions may be nearly as great as that necessary to raise the temperature of the rock in the absence of thermal reactions. If the heat requirements for these reactions are not considered, serious errors in heat-transfer calculations may result. Thermal-tzxpansion and thermal-reaction data have been obtained for three typical sandstones. Methods of measurement and results obtained are reported and discussed in the following. THERMAL EXPANSION OF SANDSTONES The linear expansions of three outcrop sandstones (Bandera, Berea and Boise) have been measured in the temptbratue range of 25° to 1,000° C. Expansion measurements were made in the directions parallel and perpendicular to the bedding. The differential expansion apparatus used is that described by Mitoff and Pask.4 The test sample is heated in an electric furnace with a temperature rise of 6°C/minute. The lengthening of the samples upon heating is compared with the small and known expansion of a fused silica rod. The change in length is transmitted to an X-Y recorder by means of a Statham transducer. A maximum error of ± 1.5 per cent has been obtained with this apparatus

By A. S. Michaels, R. S. Timmins

A study of the effect upon oil recovery from an un-consolidated porous medium of chromatographic transport of selected reverse-wetting additives during water displacement is described. Flow tests were performed on hydrocarbon oil-saturated sand columns (containing connate water) into which small volumes of short-chain aliphatic amine solutions were introduced during water flooding. Concurrently, meusurements of oil-silica contact angles, oil-water interfacial tensions and oil-water distribution coefficients were measured us a function of additive concentration. Certain treatments were found to cause a slight reduction in oil recovery, while others produced a significant increase (up to 50 per cent of the normtzl irreducible rninimum oil saturation). Recovery eficiency was correlated most satisfactorily with the maximum degree of reverse wetting produced during additive transport; conditions which resulted in only mild oil wettability of the sand reduced recovery, while those favoring strong oil wetting increased recovery. A mechanism of reverse-wetting additive action during chromatographic transport is developed which ascribes the results to transient changes in oil-solid contact angle during flow, with consequent release of oil trapped by capillary forces. Implications of these observations with regard to improved waterflood efficiencies are discussed. INTRODUCTION The need for finding improved methods of oil production is amply illustrated by the large volume of oil considered economically unrecoverable by existing production practice. Water flooding has long been an excellent method of secondary recovery. However, even in those areas of a formation which are well swept by a water flood, 20 per cent or more of the oil remains trapped within the formation and is not recovered. This trapping of oil within the pores of the rock structure results primarily from interfacial forces. If these interfacial forces were altered, it is possible that the trapped oil could be released and recovered. The addition of surface-active agents to a water flood to alter the interfacial properties of the oil-water-solid system has received considerable attention as a method of increasing oil recovery. Studies of surfactant treatment have yielded varied and, in some cases, apparently conflicting results. Increased oil recovery using surface-active agents has been reported by a number of investigators. Others have found that use of surfactants in a water flood either has no effect upon oil recovery or decreases it. Since surface-active agents tend to alter the formation wettability and oil-water interfacial tension simultaneously, several investigators have tried to determine the role these two factors play in oil recovery. Data have been presented which indicate that the most desirable condition for the solid surface is strongly water-wet. Other results indicate that, although breakthrough recovery may be less, the greatest ultimate recovery is obtained when neutral wettability prevails.", " The desirability of low interfacial tension seems to depend upon the wettability of the formation. Reports of higher recoveries, when the interfacial tension is high in a water-wet system and low in an oil-wet system, are at odds with reports of increased recovery in water-wet systems when surfactants which decrease the interfacial tension are added to the water flood. A review of previous work studying the effects of using surface-active agents in a water flood shows that these agents produce widely different effects and that opinions as to the role of the various interfacial properties are diverse. The consensus, however, is that the interfacial properties are quite important in the displacement of oil from a porous medium by water. Even in those investigations where appreciable increases in oil recovery resulted from treatment with a surface-active agent, economic considerations generally made practical application of the treatment unattractive because surfactants tend to be adsorbed within the formation. This adsorption would necessitate using large quantities of surfactant if all of the flood water were treated to contain a given surfactant concentration, even if this concentration were quite low. Preston and Calhoun suggested that the principles of chromatography could be applied to the use of surface-active agents as a means of increasing oil recovery. A small volume of surface-active material is injec-

By L. L. Handy

Two equations are compared for calculating imbibi-tion rates of water into porous media. A linear relation was found to apply when the square of the volume of water imbibed was plotted against time for water imbibing in natural rock materials. The effect of adsorbed organic materials on surface wettability was determined from measurements of rate of imbibition, water saturation and water permeability. INTRODUCTION An understanding of the mechanisms by which water displaces oil from porous media requires an understanding of the role of capillary forces in the displacement process. The capillary forces are determined by properties of the water-oil-solid surfaces. Consequently, some measure of the surface properties of reservoir rocks is necessary if displacement behavior in oil reservoirs is to be predicted from laboratory measurements. Several experimental procedures have been proposed to assign quantitative wettability indexes to reservoir rock surfaces. The most recent of these proposals are those of Bobek, Mattox and Denekas' and Amott.' These indexes are designed to show a continuous variation from the preferential oil-wet to the preferential water-wet systems. They require measuring some property of the rock which is a function of surface wettability. The quantities are measured on unaltered core material and compared with values obtained for known oil-wet and water-wet extremes on the same material. These methods are useful but are semi-empirical in nature. They have the disadvantage that the measured quantities may be functions of other variables in addition to surface wettability. Some understanding of the behavior of water-oil-rock systems may be obtained from studies of the simpler water-air-rock systems. Surfaces of a given wettability are more easily maintained for the latter. Also, certain approximations can be made for the flow of water and air which are not justified for water-oil flow. The theoretical and experimental work discussed in this paper is part of a study of the role of capillarity in the dynamic displacement of air by water. The expec- tation is that observations for the water-air system will lead to generalizations which are applicable to water-oil systems. THEORY Imbibition is defined as that process by which one fluid displaces another from a porous medium as a result of capillary forces only. Differential equations have been derived to describe the displacement of oil by water during imbibition, but these equations have not been solved.".' With certain assumptions, however, solutions are possible for water-air imbibition.

By L. H. Robinson

Triaxial compression tests on Indiana limestone indicate that chemical additives in the pore fluids can increase or decrease the strength in plastic failure but not in brittle failure. Tests using a series of the sodium salts of dicarboxylic acids show that odd-number carbon-atom chains decrease the yield strength as much as 3.0 x 10 3 psi, and even-number carbon-atom chains increase the yield strength as much as 9.0 x 102 psi. ,The most effective chemical (sodium citrate) decreased the yield strength of Indiana limestone from 16.0 x 103 to 12.0 x 103 psi. The chemical additives investigated have no influence on microbit drilling rate when drilling with zero differential pressure. However, when drilling at a 2.0 x 103-psi pressure differeniial and using an odd-number carbon chain compound as the circulating fluid, the roller bit drilling rates were slower than with water whereas the drag bit drilling rates were faster. One even-number carbon chain compound added to the circulating fluid produces the opposite effect: the roller bit drilling rate is slightly faster than with water and the drag bit slower than with water. INTRODUCTION A Russian book by Rehbinderl published in 1944 generated much interest in the possibility of aiding and accelerating the destruction of solids by the action of a liquid in which the solids are deformed. This book describes a process in which microcracks are generated during elastic deformation. To describe this process, the rock is divided into three zones: (1) the zone in which the destruction takes place, (2) the predestruc-tion zone (an affected area immediately below the destroyed zone) and (3) the undisturbed specimen. Consider a load applied normally to a surface of an infinite slab of rock. Immediately below the point of application of the load, the rock is destroyed by compression or shear forces. Below this area of destruction, because of a great number of discontinuities ranging from defects to weak points in crystals of discontinuities in crystal latices, very small cracks are created in or between the crystals (often referred to as microcracks). Upon removing the load, the broken rock in the destroyed zone could be removed. In the lower predestruction zone, however, the microcracks in and between the crystals could disappear by the process of rehealing, i.e., the molecular attraction from one side of the crack to the other side being sufficient to completely eliminate the crack. Thus, a large number of these microcracks possibly exist only for a relatively short period of time. The concept of small cracks being generated within the crystals of a rock is not unreasonable. If the rock matrix is deformed below the yield point, some of the defects present within the crystalline structure could start small cracks through the interior of the grains. Upon removal of the load, these microcracks reheal (or the newly generated surfaces come back together) so that all evidence of cracks completely disappears. Under various types of loading, and in various types of fluids, the rehealing of these cracks would be time dependent. If a liquid adhered to the surface of a rock, it could penetrate into the microcracks and thus delay or prevent the rehealing process. This suggests that the prevention of healing of the microcracks may be dependent upon the wetting characteristics of the solid. Thus, rocks which are water-wet should be more easily destroyed in water than in nonpolar hydrocarbon liquid. Adsorption of a liquid could also reduce the free surface energy which could explain the decrease in strengths of materials under nonrepetitive-type loading. Rehbinder's work prompted Boozer, Hiller and Serdengecti 2 to examine the effects of pore fluids on the deformation of rocks subjected to triaxial compression over a range of confining pressures, temperatures and strain rates. Fluids which strongly adsorbed on the rock grain surfaces decreased the strength of sandstone and limestone, which deformed plastically. For certain strain

By V. M. Dobrynin

Experimental data demonstrate that physical properties of porous rocks change under pressure. In this paper an assumption is made and proved that under pressure the changes of physical properties such as porosity, density, permeability, resistivity and velocity of elastic waves are controlled to a large extent by the pore compressibility of rocks. It is also shown that the pore compressibility of rocks can be determined, within the range of pressures from 0 to 20,000 psi, by knowing the maximum pore compressibility and the magnitude of the pressure. Mathematical equations were developed which permit one to define changes in physical properties of porous rocks under pressure. These equations were verified by experimental data obtained from the study of sandstones. INTRODUCTION In studying the behavior of porous rocks under pressure in the field of petroleum technology, the most interesting aspect is the observation of those properties which characterize the rocks as possible reservoirs — for example, porosity, permeability, resistivity, density and the velocity of elastic waves. The literature dealing with this problem mainly contains data concerning the study of only one or at most two of these parameters, but not of the group as a whole. An attempt is made in this paper to find general equations involving each of these parameters, which will permit the study of the behavior of rocks under pressure. All experimental data used here were obtained from the investigation of consolidated sandstones. EXPERIMENTAL In addition to the use of published experimental data, an experiment was carried out which studied the main physical properties of sandstones under pressure. Two homogeneous quartz sandstones were chosen for this purpose: (1) the Torpedo sandstone from Kansas, and (2) the Medina sandstone from Ohio. The porosity of the Torpedo sandstone was 20.2 per cent, and that of the Medina 8.7 per cent. Permeabilities were 45 md and less than 1 md, respectively. Each sandstone contained about 5 per cent clay minerals, consisting mostly of kaolinite and chlorite, which were distributed quite evenly throughout the samples. One cylindrical sample 2 in. in diameter and 5 in. in length was cut from each sandstone and then saturated in a vacuum with a 3N solution of NaC1. This high concentration was used in order to obtain true formation factors and to decrease the swelling of the clay minerals. The methods of mounting the samples and measuring the changes in porosity and resistivity were practically the same as those described by Fatt and Mann.1,2 Changes of resistivity under pressure were studied for sandstones with 100 per cent water saturation, and for sandstones with the irreducible water saturation. The irreducible saturation was obtained by enclosing the saturated rock samples in relatively fine silicate powder so as to remove the water by capillary action. This procedure is described by Orkin and Kuchinski.3 Changes of permeability with pressure were determined at room temperature using nitrogen as the flowing medium. In studying the effects of pressure, one series of measurements was made using an internal pore pressure Pi equal to the atmospheric pressure, while the overburden pressure Pe ranged from 0 to 20,000 psi. A second series of measurements was used over the same range of overburden pressure, but with an internal pore pressure of 1,800 psi. When the results were compared on the basis of net overburden pressure (Pe — 0.85 Pi ), there was practically no difference for these two sandstones. The origin of the factor 0.85 in the expression for net overburden pressure is given by Brandt,4 Fatt5 and Geertsma.6 The relatively large bulk volumes of the sample (about 15 cu in.) made it possible to observe a

By C. C. Miller, A. B. Dyes

The cost of finding and developing new reserves is continually rising. We must meet these rising costs with more economical operations. This can he accomplished if we revise our ideas of proper well spacing and well allowable to consider the concept of optimum well spacing. According to this concept, the optimum spacing is the one which leads to the maximum present worth for a reservoir when ail factors affecting total cost and total revenue are considered and when the wells are produced in the most eficient manner. Application of this principle efficiently utilizes available well potential and properly considers the recovery eficiency in addition to fixing the spacing on the basis of the amount and value of the oil to he recovered. This Study presents an analysis of one producing zone containing low gravity crude to illustrate the effect of these factors on the present worth and on the optimum economic spacing under two production drives—-evolved gas and water drive. The maximum present worth occurs when the optimum number of wells for open-flow operation is employed. Frequently, this optimum development cal1s for very wide spacing and the ideal field rates are not unreasonable. Under other circumstance. where proration is necessary, an optimum combination of well spacing and well allowable exists which permits production at relatively high rates. The optimum well density in a field depends on the recovery efficiency and the valule of the oil. In solution gas-driven reservoirs this optimum spacing for operation at high producing rates can vary from extremely wide spacing to handle viscous low gravity oil in thin formations to relatively close spacing in thick sands where good recoveries are expected. Because of the better recovery from water-driven fields, the optimum spacing in these fields is closer than in solution gas-driven fields. Also, the water encroachment pattern is dependent upon the well spacing, and an adequate numher of wells is needed to assure a good sweep eficiency. The economic optimum well density in a water-driven field is high enough for this purpose. INTRODUCTION From year to year domestic oil becomes more difficult to find. The fields we do locate are frequently smaller and deeper than older fields and the costs for men, material and equipment are continually rising. To replenish our reserves we must continue to search for and develop new fields, but experience has shown we cannot expect prices to rise in proportion to costs. Consequently, we must meet increasing costs by more economical operations. A vigorous effort is being made within the industry to improve exploration methods, to cope with the problems of deeper drilling, and to obtain a secondary yield from older fields. Many significant contributions have resulted from these efforts. The economic operation of new fields can be further improved by developing these fields on optimum spacing and producing the wells at higher rates. This would avoid the drilling of unnecessary wells and provide additional capital for seeking new oil. While the benefits of these practices are obvious, the problem of defining the optimum development of a field for natural depletion can become very complex. A study of the effect of well spacing and several reservoir variables on economic worth of a specific field is reported here to illustrate the problem and to show the magnitude of the benefits to be realized. This study is necessarily limited to the field conditions selected and is not intended as a general solution of the well spacing problem. It does, however, indicate factors to be considered, the trends to be expected, and the direction in which we should proceed in developing new fields. NO consideration is given to land and legal considerations which might arise. METHOD OF ANALYSIS The optimum method of developing and producing a field is to use the combination of spacing and prora-tion which gives the maximum return. In addition, we do not want to lose recovery. These considerations of maximum return and maximum recovery present no serious conflict. The value of each method of operation is conveniently expressed in terms of its present worth. According to the present worth method of evaluation,' an acceptable annual percentage return is assigned to the operation, and all incomes, capital costs, and expenses are discounted at this rate to the start of the operation. This net value is the present worth. If we apply the same discount rate to several alternate methods of operation, the one yielding the greatest present worth is the best method. If we express net income in terms of price per barrel of oil produced, drilling and equipment costs on a well basis, and expenses as cost per well year, this evalu-

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

It is deemed feasible to store atomic reactor fuel wastes in salt dome cavities when the differential pressure acting on the cavity does not exceed 3,000 psi and the temperature is less than 400°F. Tests at pressure increments of 1,000 psi were conducted on a 2-in, cylindrical cavity contained in a 6-in. long X 6-in. diameter salt core. The cavity exhibited stability under pressures up to 3,000 psi and temperatures to 300°F. At temperature ranges of 100 to 400°F and with pressures of 5,000 psi and above continuous deformation of the cavity resulted. Initial movement of the salt was observed at all pressares. This was evidenced by vertical deformation and cavity size reduction. INTRODUCTION The development of storage techniques for safe disposal of atomic reactor fuel wastes is a goal of the Atomic Energy Commission. Salt domes have been used successfully for storage of liquefied petroleum gases and hence should prove to be suitable for liquid waste products of atomic reactors. The purpose of the following work is to evaluate the practicability of utilizing salt dome cavities for such wastes and to de-termine particularly the effect of pressure and temperature on a salt cavity. This work represents a part of a comprehensive study covering the over-all problem of storing re- actor fuel wastes in salt domes and/ or salt beds. The problem of accurately determining the stresses which exist in rocks in the earth's crust is indeed a complex one. Theoretical studies have differed widely in many of the basic assumptions about the physical properties of the rock itself. Some solutions of problems in underground stress analysis assume that rock is elastic, homogeneous and isotropic in character; others assume that rock possesses plastic, viscous, elasticovis-cous properties, or a combination thereof. Experimental evidence has shown that salt possesses plastic, viscous and elasticoviscous properties.' Recent work has shown that salt becomes more plastic with temperature.' APPARATUS AND PROCEDURE For experimental study a cylindrical cavity was chosen for analysis since it provided maximum volume per unit of depth with constant diameter. This particular cavity shape, allowing for irregularities, can be easily leached out of salt domes in field operations. For the purpose of approaching underground storage conditions as nearly as possible, the following factors were considered: (1) control of pressure and temperatures, (2) measurements of temperature gradients throughout salt mass, and (3) measurement of structural stability of salt cavities. In order to obtain triaxial pressure on the salt cavity without resorting to the application of fluid pressure directly on the salt mass, a system was designed whereby the same effect could be produced through axial loading only. By cutting a core to close dimensions to fit a steel test cylinder and applying axial pressure, a lateral pressure was obtained. This equipment is shown in Figs 1 and 2. A steel cylinder was made from

By J. S. Osoba, J. W. Graham, P. H. Monaghan

Of the many factors which affect the productivity of hydraudically fractured wells, the wettability of the propping sand has received little attention in the pas/. This paper shows that the wettability of the propping sand is an important factor which shorcltl not he overlooked. An analysis of relative permeability data shows that fractures packed with water-wet sand should be more permeable to oil than are fracturer packed with oil-wet sand. Laboratory results verify this conclusion. Cnlculations show that higher permeability to oil in the fracture should provide higher well productivity. A method is presented which ensures that the propping sand in an induced fracture is water-we1 even though oil is used as the fracturing fluid. The rnethod has been used in more than 150 field jobs: the field resrrlts are discussed briefly. The fluidized nature of water-wet sand in oil has giver1 unexpected benefits. This characteristic of water-wet sand and how it has been used advantageously in field operations are discussed. INTRODUCTION Hydraulic fracturing of subsurface strata to increase the productivity of oil and gas wells has been practiced for about 10 years. Today, fracturing is used more often than any other method of well stimulation. Basically, the process consists of breaking down (fracturing) a subsurface formation with pressure applied by means of a carrier fluid, usually oil or water, and propping the induced fractures apart by means of an agent, usually sand. suspended in the carrier fluid. Currently about 68 per cent of all fracturing is done with a simple mixture of refined oil or lease crude oil and sand'. The research of other investigators has suggested that dry sand is rendered oil-wet when it is contracted with certain types of crude or refined oils.' Hence, it may be expected that some fracturing operations employing an oil as the carrier for sand leave fractures packed with oil-wet sand. Conversely, those operations which employ water as the carrier for the sand probably leave the fractures packed with water-wet sand. Since the permeability of a sand-packed fracture to a mixture of oil and water is dependent upon the wettability of the fracturing (frac) sand, research was undertaken to determine to what extent the wettability condition (water-wet or oil-wet) controls fracture permeability which, in turn, affects well productivity. Additional research was directed toward the development of a method for obtaining the desired fracture wettability. In this report both theoretical and experimental evidence is presented to show the desirability of employing water-wet sand in fracturing operations. A method to obtain water-wet sand is discussed, and field results are reviewed briefly. DEPENDENCE OF WELL PRODUCTIVITY ON FRACTURE WETTABILITY To determine the dependence of well productivity on the wettability of the frac sand, it is necessary to know, first, the effect of fracture wettability on fracture permeability and, second, the effect of fracture permeability on well productivity. The relationship of these variables is considered separately in the order indicated. Effect of Wettability on Permeability Fig. 1 shows relative permeability relations for 20-30 mesh Ottawa sand. Inspection of these curves reveals that the relative permeability to oil when a small amount of water accompanies the flow is considerably higher when water rather than oil is the wetting phase. For

By P. D. White, J. T. Moss, J. S. McNiel

Results of experimental work on the in situ combustion process were first published in this country in 1953' when Kuhn and Koch described results of a three-well test in Jefferson County, Okla. Shortly there-after Grant and Szasz2 described field studies in Nowata County, Okla., in which several modifications of the basic combustion process were investigated by injection of various oxidizing gas mixtures other than air. In the past five years speculation concerning this method of oil recovery has continued, with a number of interesting papers on various aspects of in situ combustion appearing in the literature3,4,5,6,7. However, only one of these' provided actual engineering information on pattern-type field experiments. The primary objectives of the Oklahoma three-well test were to check operation of specially designed electrical ignition equipment and evaluate problems associated with initiation of combustion in a natural oil sand reservoir. Immediately following successful completion of this test in the Pontotoc sand, plans were activated for a second field experiment to be conducted in the same sand a short distance away. The objective of this second field test1 as to establish a relationship between field observations and information derived from laboratory experimentation and computations. Secondary objectives of the test were to develop operating tech- niques and define operating problems, and to obtain information that would assist in the economic and engineering appraisals of the process. Careful planning was required since this test was designed with the hope of providing basic information on the characteristics and operation of a combustion drive. It was necessary to include test facilities which would not only permit precise control and measurement of the process variables but would also provide sufficient flexibility of control to insure that variables in the combustion process could be studied during the course of the test under a wide range of conditions. In order to accomplish these objectives elaborate instrumentation and several control wells were included in the design of the test installation. This type installation was felt to be well suited for an experimental study, but was in no respect intended as a counterpart of a commercial-scale operation. PRELIMINARY LABORATORY WORK An important objective of the field experiment, as stated previously, has been to obtain reliable data on operating variables which might be compared with similar data obtained in laboratory combustion tests. The value of successful correlations between quantities that can be measured in the laboratory and those obtained in a full-scale commercial venture has been well established for other industrial processes. The difficulty of following a process underground should make such correlations even more essential here. Also. approaching the field test phase of development of in situ combustion, it was considered desirable to precede design of test equipment with laboratory combustion studies under a variety of conditions to establish likely ranges of variables. These included temperature, air injection rate, frontal advance rate and produced gas composition.

By W. F. Hower, W. Brown

Large-scale sand consolidation tests were conducted in an effort to determine the reasons for the successes and failures of this method of sand control. Several different consolidating materials were used in treating both clean and bentonitic sands that were packed in a chamber having a capacity of 3.3 cu ft. The results were essentially the same for all of the different consolidating materials, The data show that low-viscosity consolidating materials pumped at a relatively slow rate gave the best results. Where the formation has produced sand, the treating fluids can compress the formation, thus permitting the channeling of fluids to another horizon. Pressure-packing these zones before attempting to consolidate is recommended. Sands containing more than 4 per cent of water-swelling clays are not good candidates for consolidation. It is indicated that loose sand, particularly when it is bentonitic, can be fractured during the placement of the treating fluids. INTRODUCTION Sand production in oil and gas wells has plagued the industry for many years, and numerous cures for this problem have been suggested. Most methods have been successful to a certain degree, but the great variety of well conditions that exist in the different areas has magni- fied the problem and limited the successful use of the various systems. Four review papers1-4 present a wealth of information concerning the degrees of success that have been obtained by the different sand-control methods. The bridging of sand grains by the use of gravel packs and screens has been quite successful. However, these methods do not leave the casing clear for all types of multiple completions, and the cure does not last for the production life of the well in some instance:;. The control of loose sands by sand consolidation with resins has never been as successful as desired. It has always been hoped that such a treatment would eliminate all sand problems for the life of the well, but. initial applications, starting in the middle 1940's, were only moderately successful. Lott, et a1,3 reported a success ratio of approximately 50 per cent and made the following conclusions. The highest percentage of successes were obtained where: a. Consolidation of a zone was made at the time of initial completion or prior to the production of sand. b. The interval treated was less than 12 ft in length. c. Between 30 and 50 gill plastic/ft of producing interval was displaced through the perforations. REASONS FOR SAND CONSOLIDATION FAILURES Our own experiences in the field of sand consolidation point toward the following conditions as the major reasons for the failure of sand consolidation attempts. 1. Mud-plugged perforations and mud invasion of the formation. 2. Sand in the casing covering all or part of the perforations. This sand could be either formation sand or one of the coarser sands used as propping agents in hydraulic fracturing. 3. Holes in the casing. 4. Channels behind the casing. 5. Attempting to treat too long a perforated section. 6. Too high a percentage of water-swelling clays in the formation. 7. Formations that have produced sand. Recent attempts were made to treat perforated sections ranging from 10 to 30 ft, in wells that have produced sand, by using a straddle packer that was raised and lowered through the perforations as the consolidating material was being pumped. In most instances, the pressure required to pump fluid into the formation varied considerably as the tool was raised and lowered. This suggested the possibility that significant differences in permeability were present or that only part of the formation had produced sand. There were times when a sudden break in pressure indicated that a fracture was being formed. Research conducted several years ago concerning the problem of the control of water in air and gas drilling indicated that shale sections could be fractured quite easily. In addition, it was determined that it was easier to pump fluids into shale bodies by fracturing the shale itself, or the interface between the shale and sand, than to pump into a fluid-saturated formation. Formations that produce sand are usually adjacent to shale bodies and frequently have shale streaks of various thicknesses inter-bedded in the sand. Therefore, where shale is exposed to fluid pressure it

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By V. J. Sikora, T. S. Hutchinson

This paper presents a method for making a water-rlrive ana1gsis without prior knowledge of aquifer geometry and uniformity using a standard desk calculator. Although it is necessary to know the initial oil in place to use this method, this is a minor limitation to the scope of (he method since in most reservoirs it is possible to obtain reasonably accurate volumetric oil-in-place estimates. An equation is developed which relates pressure at the water-oil contact to the water influx rate as a function of time by a factor called the resistance function. The resistance function introduces the composite effect of the aquifer geometry and flow resistance disrribution. It is pointed out that the characteristic shape of this function makes it possible to start with on approximation of the function and successively improve the approximation until the correct resistance function curve is obtained. In this fashion water-drive estimates can be made without the limitation of assuming simplifed aquifer shapes and flow distributions. This is the novel feature of this development. Methods are given for extrapolating the final curve to calculate future aquifer behavior. Equations are developed for adjusting the pressures and water influx rates where it appears that possible errors in these quantities make it difficult or impossible to obtain a useuble resistance function curve without this adjustment. Application of the pressure build-up analysis techiiique to estimate. some of the aquifer properties is also presented INTRODUCTION A great many oil pools are the result of oil accumulation in some type of trap in an otherwise large and continuous porous stratum. The void space of this stratum outside of the oil pool itself is filled with water or brine. In analyzing performance of the oil pool surrounded by this water aquifer, it is quite necessary. in most cases, to include behavior of the water. When the pressure at any point in a fluid system is lowered, such as by opening a well in an oil sand, fluids in the immediate neighborhood of this point will begin flowing towards this lower pressure sink. As pressure in this area drops due to flow towards the sink, fluids from farther out will start to flow towards the lower pressure. As more and more fluid is removed from the system, the distance from the sink or well within which flow is occurring will continually increase; that is, thc region of disturbance will grow. If some rigid boundary such as a fault is reached by the disturbance, this area will cease to grow; but if some movable boundary is reached, such as a water-oil contact, the area will grow on out into the water, although rate of growth may and almost always does change. The relation between amount of pressure drop and amount of fluid flowing at any point and at any time in an aquifer depends on such factors as compressibility and viscosity of the fluids, the porosity and permeability of the rock, geometry of the whole system and withdrawal rate or pressure drop. With these factors known, it is theoretically possible to calculate the pressure-flow behavior of the system, However, in practice true solutions to the problem are next to impossible due to complexity of reservoir systems. A number of approximate methods of solution have been developed based on various simplifying assumptions. One frequent assumption is that the water-oil contact can be located and equations defining oil reservoir and aquifer solved separately by assigning values to various parameters so as to match past pressure and production history. Usually the properties of the aquifer are not known, since few, if any, wells are drilled through the aquifer; but the water influx and the pressure at the reservoir boundary are known over some time. If it can now be assumed that the aquifer is circular, pie-shaped, or linear, and that it has uniform properties, it is possible to fit a theoretical dimensionless curve to the past aquifer performance history and therefore calculate future water drive. These theoretical curves are available in the literature,1,2,3 Of course, the assumption of uniform aquifer properties is almost always somewhat in error. Fortunately, however, moderate variations from the average have little effect on behavior of the system; and, hence, the best fit of a theoretical curve is frequently satisfactory. In other cases variations in aquifer properties and geometry are large enough that none of the available theoretical curves will give an acceptable fit of the data. In these cases, methods of fitting the data with various electrical analyzers have been developed.1,2' It is the purpose of this paper to present an approach

By B. Wilhelmi, W. H. Somerton

Simultaneous measurement of pore and elastic properties of socks under a wide range of triaxial stress conditions may be made with equipment and methods developed in this paper. Tests were run on three outcrop sandstones at several confining pressures and axial stresses up to approximately 80 percent of failure stress. The greater reduction in porosities and permeabilities occurred in the hydrostatic loading portion of the stress cycle, but substantial reductions in both properties did occur upon application of deviator (triaxial) stress. Triaxial elastic moduli were calculated from strain gauge measurements. Both Young's moduli and shear moduli were found to increase with increasing confining pressure. Poisson's ratios increased with axial stress but showed no consistent change with confining pressure. Only very general relations were observed between pore and elastic properties. The lesser decrease in porosity and permeability at higher stress values is related to increase in rigidity of rocks at higher stresses. INTRODUCTION Many investigators have studied the effects of hydrostatic stress on the properties and behavior of rocks. A few have studied the effects of triaxial stress fields on individual physical properties of rocks. To our knowledge no one has measured simultaneously pore and elastic properties of rocks under a wide range of triaxial stress conditions. This is the purpose of this work. Rocks exist in the subsurface under what Jaegerl refers to as "polyaxial" stress conditions (all principal stresses unequal). When two of the three principal stresses are equal, this is referred to as triaxial loading. The triaxial method of loading is most commonly used in modern testing work because experimental difficulties in measuring physical properties of rocks are greatly reduced over polyaxial testing. Although triaxial loading may not reproduce subsurface conditions precisely, subsurface stress conditions are never known precisely anyway. Triaxial testing would seem to be more realistic than testing under hydrostatic loading (all principal stresses equal) for application of results to subsurface problems. The objectives of this work were threefold: (1) to develop methods of measuring simultaneously pore and elastic properties of rocks under triaxial stress conditions, (2) to determine whether changes in these properties are of sufficient magnitude to be considered in subsurface calculations and (3) to attempt to relate changes in pore properties and elastic properties since they will have been measured at the same time under identical loading conditions. REVIEW OF EARLIER WORK PORE PROPERTIES Several investigators have studied the effects of stress on the physical properties of rocks. Fatt and Davis2 showed the reduction of permeability of rock with hydrostatic loading. McLatchie et al.3 and Knutson and Bohor4 related permeability reduction to compessibility of reservoir rock, again under hydrostatic loading. Wble,5 and later Redmond,6 evaluated the effects of applied radial pressure on conductivity, porosity and permeability of sandstones. In this latter work, no axial stress was applied. Dobrynin7 reviewed much of the previous work done . on the effects of overburden pressure on physical properties of rocks. He concluded that changes in physical properties are controlled to a large extent by pore compressibility of the rock. Several sets of curves developed from empirical relations are presented. These curves permit estimation of change of properties with net overburden pressure, knowing the pore compressibility and porosity of the rock. Net overburden pressure is defined as (Pe - 0.85 Pi ) where Pe is the

By T. K. Perkins, L. E. Bartlett

Values of rock surface energy (i. e. energy required to form a unit area of new surface) are useful for interpretation of drilling and fracturing phenonema This paper describes the adaptation of a cleavage technique (splitting of large rocks under controlled conditions) for measuring surface energies. Values of surface energy obtained at room temperature and at low crack extension rates are reported for a variety of sandstones and limestones. Elastic moduli, tensile strengths and compressive strengths of the samples are also reported. The values of Young's moduli (measured in bending) are relatively low, probably because of the low stress levels imposed. Measured values of surface energy me generally higher than anticipated and probably reflect inelastic yielding near the extending tip of the crack. INTRODUCTION The technology of petroleum production includes many operations such as drilling and fracturing, which involve rock breakage. One of the fundamental properties of a rock needed to fully explain breaking and fracturing phenomena is the surface energy (i. e. the energy required to form a unit area of new surface). In this paper we describe an adaptation of the cleavage method (splitting of the rock under controlled conditions) to utilize large rock samples. Application of the theory of elastic cantilever beams and the Griffith concept of mechanical stability permits the calculation of surface energies from simple measurements. Values of surface energies obtained at room temperature on dry samples and at low crack extension rates are reported for a variety of sandstones and limestones. Measured values of Young's Moduli, compressive strengths and tensile strengths are also reported for the rocks. THE EXPERIMENTAL METHOD Equipment designed to cleave large rock samples is shown on Fig. 1. This equipment consists essentially of two steel blocks cemented to a rock sample and forced apart with a steel rod. A ball joint and a pivot prevent bending stresses in the rod. A differential thread arrangement permits the rod to advance only 0.012 in./revolution of the driving nut. The axial force in the rod is measured with type FA-50-12 Baldwin Lima Hamilton strain gauges and a BLH type 20 strain gauge reader. Separation of the blocks is measured with a Starrett 656617 dial gauge. The massive nature and design of the equipment minimize the amount of elastic energy stored in the steel parts as the stress in the rod is increased. This will permit the rock to split under controlled conditions as will become clear later. Fig. 2 shows a sketch of the rock sample during cleavage. Guide slots are cut longitudinally along the rock sample as suggested by Berry.1 An initiating slot is also cut across the top of the block to a depth of 1 to 2 in. During cleavage, elastic energy is stored in the flexed beams of rock of length C. If the length of steel rod L is held constant and the crack extends in length, then the elastic potential energy stored in the rock decreases. However, surface energy must be supplied as the crack extends in length. Griffith6 has proposed the following criteria for mechanical stability. The crack will extend as long as the potential energy available is greater than the surface energy required. However, the crack will stop extending when the change in potential energy with change in crack length is just equal to the change in surface energy with change in crack length. See Eq. 1.

By J. A. Putnam, A. D. K. Laird

In laboratory research on the behavior of oil, gas and water in porous materials, no direct method has been devised to measure saturation without disturbing the flow. Indirect methods involving various forms of radiation have been developed for measuring one or two components. The X-ray absorption method for two components has been described previously1-3. The possibility of three component measurement has also been mentioned1. Since two component measurement still requires considerable care to insure accuracy, the fact that three components can be measured satisfactorily is not generally appreciated. The fraction of the energy of an X-ray beam of one wavelength transmitted through a material is given by the transmission factor. in which 1 is the beam's path length through the fluid. p is the fluid density, and 1,. is its mass absorption coefficient. Eq. 1 implies that all parts of the beam pass through equal amounts of the fluid. Practical applications are, therefore, limited to systems in which the fluid is distributed so that I varies negligibly over the cross-section of the beam. Typical variations of the transmission factor with X-ray tube potential arc shown in Fig. 2 from which experimental points have been omitted for clarity. The X-ray absorption of low pressure gases in a core is practically zero so it cannot he distinguished from that a vacuum by instruments that will differentiate between liquids. The transmission factor of a gas in a core is. therefore, unity. Consequently, the saturation of only one gaseous component can be found. If two fluids completely fill a core, their two saturations add to unity and, therefore, only one value of the transmission factor, T, is needed to determine both of them. When two liquids do not completely fill the core, their saturations do not add to unity and transmission factors at two wavelengths must be measured. At one of these wavelengths, 7 for one liquid must be different from T or the other liquid. The change of T for one liquid us the wavelength is changed between the two values must be different from the corresponding change of T in the other liquid. The two values of fluid saturation will also give that of the gas because the three saturations must add to unity. In Fig. 2 it can be seen that there are many large differences of T between various oil solutions and water solutions at any one potential, and between potentials for any one solution. Consequently, conditions for optimum accuracy of measurement of the saturations of a gas and two liquids can be chosen. Once the operating potentials have been selected, however. it is more convenient to cross plot the data from Fig. 2 on transmission factor-solution concentration coordinates for the chosen potentials, as shown in Fig. 3. Fig. 3 was used to plan the experimental procedure for Core N-38-1. Since the voltage control at 45 kv was good, the slopes of the curves in Fig. 2 at this potential would cause little error. At low potentials, however, the voltage control was not good, so 33 kv was chosen because the small slopes of the curves made voltage control less critical. The lowest strength beam that could be measured reasonably at 33 kv corresponded to a minimum T factor of 0.385. Fig. 3 shows that both 4.80 weight per cent of cadmium chloride in water and 24.0 per cent of iodobenzene in crystal oil had this T value, and, consequently, were indistinguishable to the X-rays at 33 kv. Thus, at 33 kv, the liquid saturation was given by the same curve regardless of the relative amounts of oil and water solutions present. This simple calibration curve from which the gas and liquid saturation can be read, is shown in Fig. 4. Fig. 3 also shows that at 45 kv, T for the aqueous solution is 0.610. and that T for the oil solution is 0.220. These values were transferred to Fig. 4 to give the triangular calibration plot from which the brine and oil saturations can be found, after the gas saturation has been determined at 33 kv.

By R. R. Harvey

Recent evidence has indicated that the films forming at crude oil-water interfaces are the result of naturally occurring su,rface-active components in the crude oil rather than oxidation products resulting from exposure of the oil to air. Since the presence of such films may strongly affect the displacement of oil from porous media, a knowledge of their composition and formatiorz is of fundamental importance. In this study, without the assumption of a specific mechanism, film formation is considered to involve the transfer of a weakly surface-active species in solution in the oil to the oil-water interface. The film gas law and the Gibbs adsorption equation may thus he applied to the system. The rigid films which form upon decrease of the oil-water interfacial area are considered a result of the very slow rate of re-solution of the film material into the oil. When the shrinkage of the interface is very rapid, the amount of film material leaving the interface is negligible, and the system may then be treated in the same manner as an insoluble film. The combination of these treatments leads to a new concept of the film factor in which it may he quantitatively related to the concentration of film-forming material at the interface. Suitable experiments can then yield the following information. 1. Original concentration of surface-active species in the oil. 2. The change of concentration with time. 3. Concentration of surface active species at the interface. 4. The change of this concentration with time. 5. Interfacial tension of the system at any time and its variation with time. 6. Equilibriunz concentrations for a given set of conditions. 7. Concentration requirements for solid film formation. 8. Molecular diameter of the adrorbed species. These calculations may be made solely from data taken from measurements on the original oil. While using the pendant drop method of interfacial tension measurement of crude oil-water systems, Bartell and Niederhauser' remarked on the thin, transparent, pliable, solid films or membranes which were often formed at the interface between crude oil and water. As a result of their observations of these films, Bar-tell and Niederhauser concluded: 1. The films are formed by material from the oil adsorbed at the interface rather than by a reaction at the interface. 2. The adsorption of the film-forming material is a rather slow process. 3. There are good indications that the process is reversible; i.e., the film constituents may re-dissolve into the oil. 4. The film will spontaneously spread over any new surface formed. 5. The process of adsorption at the interface is accompanied by a decrease in the interfacial tension of the system. Although they had proved that film formation was not due to a reaction at the interface, Bartell and Niederhauser were unable to demonstrate conclusively whether these film-forming constituents were present in the original reservoir oil or whether they were formed as oxidation products after exposure of the oil to air. Recent evidence seems to point to their presence in the original oil. Although the true nature of the film-forming constituents is still uncertain, these conditions would be fulfilled by the assumption of the presence of a slightly soluble surface-active component in the original oil. This component would migrate to a freshly formed interface because it would thereby lower the free energy of the system. The component at the surface constitutes an expanded monolayer which would become compressed, upon retraction of drop, due to the slow rate of re-solution of the film components into the bulk of the oil. It is possible, by involved extraction processes, to separate or at least concentrate this film-forming fraction. However, using the assumptions made above, the technique of monolayer film study may be applied to reveal much of the nature of the adsorbing material without resorting to an actual separation of the film-forming constituents.

By J. M. Smith, D. Kunii

Effective thermal conductivities of sandstones filled with stagnant fluids were measured using a steady-state technique. Data were obtained for seven sandstone samples, taken from four different locations and ranging in permeability from 18 to 590 md. The measurements with gases (helium, nitrogen, air and carbon dioxide) covered a pressure range from 0.039 psia to 400 psig. Data were taken for four liquids — n-heptane, methyl alcohol, 79.8 weight per cent glycerol-water solution and pure water at atmospheric pressure. The experimental results were used to evaluate the theoretical equations for predicting stagnant conductivities developed earlier 4 The low-pressure measurements permitted evaluation of the consolidation parameter bpDp/ks (necessary to utilize the theory) for the various types of sandstones. Using these characteristic values, the tbeoretical equations correlated well with the experimental conductivity data for the several fluids and rock samples. INTRODUCTION An aspect of heat transfer in solid-fluid systems of considerable current interest is the effective thermal conductivity of porous media. The stimulus for study of the subject arises from the need for sound procedures for designing thermal methods of petroleum production. The general system occurs when there exists a flow of fluid through the pores of the solid material. However, a logical starting point in developing a theory for predicting the effective thermal conductivity in the general system is to attack the special case when the porous solid is filled with stagnant fluid. Since the flow rates anticipated in thermal production processes are very low, such stagnant conductivities k,O are also of practical significance. Kunii and Smith4 recently reviewed the existing data for stagnant conductivities and proposed a theory for heat transfer in porous rocks. This leads to equations for predicting the stagnant conductivity as a function of the properties of the fluid and solid phases. The existing experimental information1-3,5,7,8 covered a range of porosities and agreed well with the theory. However, the available data were insufficient to examine critically the effects of two important properties on koe, — the thermal conductivity of the fluid and the pressure. The objectives of the present study are twofold: (1) to measure the effect on koe of (a) the pressure for gases, and (L) the thermal conductivity of the fluid kg for bloth gases and liquids; and (2) to determine for rocks of different porosities and permeabilities the apparent solid conductivity ks and the consolidation parameter bpDp/ks required in applying the theoletical equations. Measurements extrapolated to zero pressure are helpful in evaluating ks and the consolidation parameter. However, the extrapolations must be based upon data obtained at low pressures because in small pores the molecular conductivity of the fluid may be unusually depressed. This occurs when the mean free path of the molecules is of the same order of magnitude as the pore diameter, that is, in small-diameter pores at low pressures. The phenomenon has been studied recently by Schotte6in beds of unconsolidated fine particles. With particles in the size range of 45 to 200 microns, a significant reduction of kg with pressure was found even at absolute pressures greater than 1 atm. In the sandstones employed in the present investigation, pore diameters lower than those experienced by Schotte are likely, so that the effect would presumably be greater. SCOPE OF MEASUREMENTS Stagnant conductivities were measured for four types of sandstones using a total of seven samples. The characteristics of the samples, including porosity and permeability values, are given in Table

By R. A. Salathiel, A. D. Kaiser, P. H. Monaghan, B. E. Morgan

The productivities of many hydrocarbon-bearing for-rnntions appear to be reduced during drilling, completing and repairing operations. The reason for these othcarwise unexplainable low productivities may be a drastic reduction of formation permeability immediately surrounding the borehole. A possible cause of the reduced permeability is swelling and dispersion of clays present in porous media when contacted by fresh water or mud filtrates. Hence, the invasion of fresh water or mud filtrate into hydrocarbon-bearing formations that contain clays may cause low productivities as a result of formation damage. Present laboratory investigations have yielded results which supplement existing literature on formation damage. An expanded-type clay, montmorillonite, way found in hydrocarbon-bearing formations not previously analyzed. When sand columns and formation cores containing montmorillonite were contacted by fresh water, their oil permeabilities were reduced. This reduction in oil permeability way restored partially but not completely by a repair procedure, i.e., flooding with kerosene containing organic agents or flooding with water containing electrolytes. Hence, where montmorillonite is present in formations, prevention of formation danzage was indicated to be the only complete solution to the problem. To prevernt formation damage, results suggested use of an approximately 1 per cent solution of calcium chloride or stronger solutions of sodium or potassium chloride. In conventional rotary drilling operations the wall of the hole is exposed to drilling mud for days, weeks, and in some instances for as long as several months. During this time, drilling mud filtrate seeps slowly into formations surrounding the hole. The permeability and stability of many formations, especially those containing clayey materials, are impaired by the mud filtrate. Often interaction between the mud filtrate and the formation causes drilling and production difficulties. Serious reduction in productivity of oil and gas formations during drilling, completing and repairing wells has been noted often in recent years. Cases are known of development wells which produced only small volumes of fluid upon completion. However, after a period had elapsed, they began to produce oil or gas in commercial quantities. In other cases, results of drill-stem tests taken during the drilling of deep wells have indicated good production should be obtained from shallow formations. Later difficulty has been experienced in obtaining production from these shallow formations which had remained in contact with drilling muds for an extended time. In yet other cases wells have failed to attain the productivities expected from previous performance after the use of drilling mud as a circulating fluid during workover operations. It is evident from the field behavior cited that the productivities of certain oil- and gas-bearing formations have been damaged by exposure to drilling muds. This reduction in well productivity, called "formation damage", causes serious consequences. In development wells, expensive remedial operations are incurred. In wildcat wells, good formations may appear to be marginal or even non-productive. The importance of formation damage has been recognized by a number of investigators1-s who have considered the effects of fresh water, saline solutions and aqueous mud filtrates on the permeabilities of formations. A general conclusion from such studies is that interaction between water and formation clays reduces permeability of reservoir rock and that the permeability decreases as salinity of the water decreases. Data in these references, however, concern the permeability change with only a single fluid present. A more pertinent and direct consideration in the production of oil and gas from wells is the permeability change when two fluids, oil and water, are present. Recognizing this, Nowak and Kreuger2 and Bertness1" considered several ways in which mud filtrates and saline solutions caused reductions in the oil permeability of formation cores. Investigations discussed here were designed to supplement information in the literature on changes in oil

By R. J. Wygal, J. Naar

An equation for three-phase (water, oil, gas) imbibition oil permeability is developed, assuming the water to be the dominant wetting fluid. Oil isoperms are obtained for consolidated sandstones characterized by 1/pc2 = CS*. The evolution of an oil-gas system imbibing water from Swi to(Sw,imb )max is shown to proceed along a line of constant oil saturation with increasing oil permeability and decreasing gas saturation. When the gas satuation cannot be reduced further, the system evolves along a line of constant Sg with decreasing oil saturation and permeability. The initial gas saturation is shown to reduce markedly the effect of complete wetting by either oil or water on flow performance. INTRODUCTION Imbibition oil isoperms are required for performance prediction when a well is producing water, oil and gas. This situation occurs in multiphase displacements such as underground combustion, steam injection and the water flooding of highly depleted reservoirs. In a recent paper,1 a model was presented for the prediction of two-phase imbibition characteristics. This paper extends the imbibition model to the case of three phases by assuming that the water is the dominant wetting fluid. The following results were obtained from the model: (1) an analytical expression of oil isoperms; (2) oil isoperrns as functions of reduced water, oil and gas saturations, valid for all sandstones having a capillary pressure curve which can be approximated by 1/PC = CS*; and (3) evaluation of the three-phase flow performance as dictated by complete wetting by either oil or water. The agreement between predicted and observed oil recovery in the presence of a gas phase, reported in Ref. 1, is a partial support for the present development. However, experimental data are not available at this time to check fully the model predictions. Perhaps this paper will stimulate the collection of such data. The imbibition model of a porous medium has been described previously, and the reader is referred to the paper of Naar and Henderson1 for details. In brief, the model is formed by the random interconnection of straight capillaries, with a provision for the blocking of the non-wetting phase by the invading wetting fluid. The following expressions were derived from the model. 1. The amount of oil blocked by the invading water is

By A. S. Odeh

ABSTRACT A simplified model was employed to develop mathematically equations that describe the unsteady-state behavior of naturally fractured reservoirs. The analysis resulted in an equation of flow of radial symmetry whose solution, for the infinite case, is identical in form and function to that describing the unsteady-state behavior of homogeneous reservoirs. Accepting the assumed model, for all practical purposes one cannot distinguish between fractured and homogeneous reservoirs fram pressure build-up and/or drawdown plots. INTRODUCTION The bulk of reservoir engineering research and techniques has been directed toward homogeneous reservoirs, whose physical characteristics, such as porosity and permeability, are considered, on the average, to be constant. However, many prolific reservoirs, especially in the Middle East, are naturally fractured. These reservoirs consist of two distinct elements, namely fractures and matrix, each of which contains its characteristic porosity and permeability. Because of this, the extension of conventional methods of reservoir engineering analysis to fractured reservoirs without mathematical justification could lead to results of uncertain value. The early reported work on artificially and naturally fractured reservoirs consists mainly of papers by Pollard,l Freeman and Natanson,2 and Samara.3 The most familiar method is that of Pollard. A more recent paper by Warren and Root showed how the Pollard method could lead to erroneous results,4 Warren and Root analyzed a plausible two-dimensional model of fractured reservoirs. They concluded that a Horner-type pressure build-up plot of a well producing from a fractured reservoir may be characterized by two parallel linear segments. These segments form the early and the late portions of the build-up plot and are connected by a transitional curve. In our analysis of pressure build-up and drawdown data obtained on several wells from various fractured reservoirs, two parallel straight lines were not observed. In fact, the build-up and drawdown plots were similar in shape to those obtained on homogeneous reservoirs. Fractured reservoirs, due to their complexity, could be represented by various mathematical models, none of which may be completely descriptive and satisfactory for all systems. This is so because the fractures and matrix blocks can be diverse in pattern, size, and geometry not only between one reservoir and another but also within a single reservoir. Therefore, one mathematical model may lead to a satisfactory solution in one case and fail in another. TO understand the behavior of the pressure buildup and drawdown data that were studied, and to explain the shape of the resulting plots, a fractured reservoir model was employed and analyzed mathematically. The model is based on the following assumptions: 1. The matrix blocks act like sources which feed the fractures with fluid; 2. The net fluid movement toward the wellbore obtains only in the fractures; and 3. The fractures' flow capacity and the degree of fracturing of the reservoir are uniform. By the degree of fracturing is meant the fractures' bulk volume per unit reservoir bulk volume. Assumption 3 does not stipulate that either the fractures or the matrix blocks should possess certain size, uniformity, geometric pattern, spacing, or direction. Moreover, this assumption of uniform flow capacity and degree of fracturing should be taken in the same general sense as one accepts uniform permeability and porosity assumptions in a homogeneous reservoir when deriving the unsteady -state fluid flow equation. Thus, the assumption may not be unreasonable, especially if one considers the evidence obtained from examining samples of fractured outcrops and reservoirs. Such samples show that the matrix usually consists of numerous blocks, all of which are small compared to the reservoir dimensions and well spacings. Therefore, the model could be described to represent a "homogeneously" fractured reservoir. la this paper, the fundamental equation of flow

Jan 1, 1966