Search Documents

By I. R. Dunlap, R. F. Burdyn, F. D. Patchen

The composition and properties of a new type of cementing mixture are discussed. In these compositions, API Class A or C cements are added to n stable water-in-oil emulsion instead of the crsual water. The cement disperses in the oil phase and remains out of contact with the water under normal conditions of mixing and storage. At higher temperatcrres, such as those encountered down-hole during placement, the emulsion breaks and cement hydration can then proceed. An extremely low filtration rate is obtained; the major part of the filtrate is oil. Thickening times of the mixture may be increased to meet the 18,000-ft casing cementing schedule by varying the amount of emzilsifier. As in conventional Portland cement mixtures, high-temperature strength retrogression can be overcome by addition of silica flour. Results of several field app1ication.r show that mixing presents no diffculty if emulsion blenders are used and the emulsion is prepared beforehand. INTRODUCTION Water-in-oil emulsions., in which oil is the continuous or external phase and in which water exists as discrete droplets, have for some time been used successfully in drilling muds to reduce filtration loss.1,2 Applied to cement slurries, the same results are obtained; i. e., the fluid-loss is greatly reduced. In addition, thickening times may be extended, lower densities and strengths obtained, and water-base borehole fluids more easily displaced without contamination of the cement. A low filtration rate is often a desirable property in oil well cements since loss of fluid increases placement difficulties and may also cause damage to water-sensitive producing formations.' In the casing cementing operation where the cement slurry is pumped into the narrow annulus between open hole and casing, rapid loss of fluid to a porous formation can cause the slurry to thicken and plug the annulus before the desired fill is obtained. Where a casing or a liner is cemented against a productive formation and the well completed through perforations, the fluid lost between the time the cement first contacts the formation and the time a set is obtained may invade the zone past the depth of the perforations. If the sand contains a material which will swell upon contact with cement filtrate, permanent damage may result. Previous methods for reducing the fluid-loss of cements have employed bentonite4 as an additive, usually with a lignosulfonate as a dispersant, or a derivative of some naturally occurring polymeric material such as cellulose ether." Oil emulsion cements have also been proposedv or fluid-loss reduction, but these have been of the oil-in-water type. Bentonite cements have been used most widely and have been successful in reducing incidence of premature setting of cement attributable to loss of water by filtration. It is the purpose of this paper to describe preparation and properties of a water-in-oil emulsion cement and to relate procedures followed and results obtained in several field applications. DESCRIPTION OF SYSTEM AND MECHANISMS A water-in-oil emulsion cement is prepared by adding an API Class A or C cement to a stable water-in-oil emulsion. Fig. 1 is a photomicrograph of this emulsion before addition of cement. Here the droplets are water and the surrounding medium is diesel oil. High emulsion stability is indicated by the small droplet size [1-5 microns], and by their uniformity.' To form a stable water-in-oil emulsion as illustrated, several conditions must be met. A predominantly hydrophobic emulsifier must be present in sufficient concentration to produce a strong interfacial film; initially, the phase volume ratio (the water-to-oil ratio) must be low;

By F. H. Turan, F. W. Jessen

The gel strength and viscosity of two different suspensions of fractionated montmorillonite clay were measured by using a Stormer viscosimeter and Farm V-G meter. The amount of sodium poly-phosphate (sodium bexametaphosphate, sodium tetraphosphate and sodium acid pyrophosphate) adsorbed by the clay particles was determined. A relationship seems to exist between the amount of polyphosphate adsorbed and the nature of the surface, edge or plate, assumed to be the principal sites for selective adsorption. Reduction in viscosity and gel strength was found to be directly related to the adsorption. Both chemical bonding and physical adsorption seem involved in defloccula-tion Some indication of the axial ratios of particles of .l- to .04-micron average size is possible though study of the adsorption mechanism. INTRODUCTION AND PURPOSE A knowledge of surface areas of clay particles in a suspension and a precise determination of phosphate adsorption should provide a basis for a better understanding of reduction in viscosity and gel strength caused by molecularly dehydrated phosphates. The purpose of this investigation was to study the adsorption of sodium polyphosphates on suspensions of fractionated sodium montmorillonite. Although there is considerable uncertainty regarding details of the structure of montmorillonite, a generally accepted one is that suggested by Hofmann, Endell and wi1m.l The structure shown in Fig. 1 is an idealized one, and Marshall2 has postulated that it is possible for several isomorphous replacements to occur. Thus, aluminum may replace silicon, in which case an exchangeable cation would go along with the aluminum to neutralize what would otherwise be a net negative charge on the structure. There is disagreement among various investigators as to the particle size and shape of montmorillonite clays. Marshall2 reported the following results for the mean lateral dimensions of a Putnam clay (montmorillonite)— 3,730 A° for particles having average diameters ranging from 500 to 1,000 A, and 750 A° for those with an average diameter less than 200 A. The latter particles were estimated at 50A°in thickness, which would correspond to a stack of five unit cell layers. On the other hand, Shaw3 has interpreted micrographs as indicating that particles of montmorillonite ire readilv broken down to the individual cell layers by dispersion in water. Brindley4 and Merings indicate the ultimate particles to be of the order of 300 A° in mean lateral dimension. Van Olphen6 suggested the thickness of particles at four to five unit cell layers. It was suggested also that this number tends to become smaller with increasing dilution. Recently, Kahn and Lewis7 determined the size of sodium montmorillonite particles in aqueous suspension from electro-optical bi-refringence decay curves. The

By R. N. Upadhyay, M. Maleki, I. Fatt

Conventional laboratory core analysis tests on samples of two limestone reservoir rocks indicate that about 20 per cent of PV is in dead-end pores. These tests (electric logging formation factor, mercury injection capillary pressure and miscible displacement) were carried out on 1/4-in. diameter test plugs. Test results show a clear dilference between these samples and sandstone or homogeneous limestone reservoir rock. Although the amount of dead-end pore space can be only roughly estimated, the presence of such pore space seems clearly indicated. Pressure transient studies also show presence of deadend PV. Although they do not give quantitative results, pressure transient data yield a reasonable estimate of the size of the neck connecting dead-end pores to the main flow channels. INTRODUCTION Equations conventionally used to describe reservoir flow behavior contain the implicit assumption that all connected pore spaces contributed to both porosity and permeability. Several authors3-7 have pointed out the changes in pressure transient behavior and in electric log interpretation that may result if this assumption is incorrect and, instead, dead-end or cul-de-sac pores are present. There is a need for laboratory tests that can detect presence of dead-end pores in core samples. With such information on hand the petroleum engineer can make more rational use of the mathematical tools now available for analysis of reservoir flow behavior. This paper describes laboratory studies designed to detect and, if possible, give a quantitative measure of dead-end PV in laboratory-size core plugs. Three reservoir rocks were used, two of which were limestones suspected of having deadend pore spaces and a well-known sandstone, used as a comparison standard, in which there is believed to be little or no dead-end pore space. All the studies were designed to measure the natural dead-end PV; i.e., the pore space which is dead-ended because of rock structure. During multiphase flow in a rock without dead-end pores, some parts of one of the phases can become surrounded by the other, thereby giving (for certain flow behavior) an effective dead-end pV8,9. Such behavior will not be described here.

Jan 1, 1967

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 P. E. Baker

A paper by A. S. Odeh' expounds a capillary-flow theory which purportedly provides a mechanism whereby relative permeability depends on oil-water viscosity ratio and may be greater than unity for high viscosity ratios. Core-measurement results are presented to support the theoretical development. This paper presents arguments to indicate that (1) the theoretical treatment is basically incorrect, (2) data reported are insufficient to support the conclusions and (3) literature cited does not support the theory or the conclusions. PROPOSED MECHANISM The essence of Odeh's theory is that oil flows as a filament along the central portion of a straight, uniform capillary of radius r filling the tube to a radius ri. Surrounding the filament of oil is an annulus of water in which exists vortex motion but no net flow. The water is described as being in "microscopic but not macroscopic" motion. The oil filament is aided in its flow by moving upon the water vortexes as if on roller bearings. This process is presented as a steady-state flow condition. The first objection to the existence of vortexes in the presumed water layer is that the theory provides no mechanism whereby they can be established. This fact was mentioned by C. F. Weinaug in a discussion published with Odeh's paper. The capillary contains no irregularities or obstructions, and the water is not presumed to be in turbulent flow. Apparently, cylindrical symmetry is assumed; therefore, the supposed vortexes must be toroidal, the line of zero velocity being a circle concentric with and in a plane perpendicular to the axis of the capillary. This is a rather difficult picture to accept. The author referred to Yuster's statement' that he had observed the "spinning in place" of otherwise immobile water. Yuster's article does not describe the conditions under which this motion was observed or the means used to observe it. However, the article does make clear that the observed motion was not in a straight, uniform capillary but, rather, apparently oc-cured in some kind of dead-ended pore space. Obvious- ly, such a report cannot be considered as support for Odeh's concept of two-phase capillary flow.

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 J. N. Chew, C. A. Connally

A correlation is presented for predicting the viscosity of gas-saturated crude oils under reservoir conditions. It is based on the dead oil viscosity and the solution GOR. The correlation was developed from a study which showed that at a fixed solution GOR, the relation between the gas-saturated oil viscosities and the correrponding dead oil viscosities is a straight line on logarithmic coordinates. Data from 457 crude oil samples from the major producing areas of the U. S., Canada and South America were used. The best straight lines through the data were fitted by the method of least squares with a digital computer. The correlation is presented in the form of an equation and also in a convenient gra-phical form. INTRODUCTION The viscosity of gas-saturated crude oils under reservoir conditions is an important physical property used in reservoir engineering calculations. Many difficulties are inherent in the sampling of gas-saturated oils and in the laboratory measurement of their viscosities. Therefore, it is frequently neither possible nor convenient to obtain measured values, so the viscosity must be estimated from other, more readily available data. This paper presents a correlation for predicting the viscosity of gas-saturated crude oils. The correlation is based on the amount of gas in solution and the viscosity of the dead oil at the temperature of interest. A previous correlation of these variables was presented by Beal. In comparison, the correlation presented here extends to a higher GOR and to a lower dead oil viscosity. Over the range of variables covered by both, this correlation is based on approximately 11 times as many oil samples. DEVELOPMENT OF CORRELATION A study was made of the variables or parameters which could be used to relate the viscosity of a gas-saturated oil to its dead oil viscosity. Preference was given to those variables which are easily and commonly measured in the field or in the laboratory. It was found that either the differential liquid formation volume factor or the solution GOR of the gas-saturated oil could be used as a correlating parameter. The solution GOR was se- lected for the correlation because of its simplicity and availability. In an attempt to refine the correlation by introducing a second parameter, the following additional variables were tested: differential liquid formation volume factor at the reservoir bubble point, the formation volume factor at the saturation pressure for the GOR of interest, reservoir temperature, viscosity of C+ fraction at 100°F and atmospheric pressure, weighted average gravity of gas in solution, oil gravity, reservoir bubble point pressure, saturation pressure at the selected GOR and reservoir temperature, and gas mol fraction. None of these improved the correlation to any significant extent. It was found that at any fixed solution GOR, the relation between the gas-saturated oil viscosities and the corresponding dead oil viscosities is a straight line on logarithmic coordinates, for the range of variables investigated. This straight-line relationship was used for the correlation. Viscosity data on a total of 457 crude oil samples were used. The data were obtained with rolling ball-type viscometers. The samples came from all of the important producing areas of the U. S., as well as from Canada and South America. Their geographical distribution is shown in Table 1. The ranges of properties represented by these samples are as follows: For each crude oil sample the viscosity of the liquid phase from each step of a differential liberation was plotted against its solution GOR on semi-logarithmic paper. A smooth curve was drawn through the points and values of viscosity in centipoises were read from

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 C. Gatlin, N. E. Garner

Results of impulsive wedge penetration tests on two synthetic, plastically deforming rocks are presented. Basic data obtained were force-time, displacement-time, and force-displacement curves for the impacts, plus the crater geometry. Wedge geometry and blow frequency were varied over a considerable range. The synthetic rocks consisted of wax-sand mixtures; two waxes of diflerent ductilities were used to provide variable "rock" characteristics. Conventional triaxial tests showed that these synthetic rocks exhibited force-deformation curves and Mohr envelopes quite similar to real rocks, except that strengths were much lower. Measured forces from static penetration tests agreed closely with theoretical values; however, dynamic force values were much higher than the static. These latter disparities are attributed to the viscous nature of the waxes. Thus the utility of these or similar rock models must depend on the scaling of rock viscosity, which is as yet unknown for impulsive loadings at elevated stress states. It appears, however, that some macroscopic, static phenomena may be studied with wax-sand rock models. INTRODUCTION The resistance of solid materials to indentation or perforation by projectiles or other penetrators has been studied by workers in many areas. Despite these efforts no universally accepted laws or formulas are available for describing experimental observations. In the metals field the force-deformation behavior of impacting bodies is often analyzed by the Hertz law for elastic collisions, the Meyer law if plastic deformations occur, or some combination of both.' The similarities of these expressions to empirical drilling formulas of the oil industry are apparent. Beginning with the basic contributions of Simon and co-workers at Battelle,' a number of experimental papers concerning the reaction of rocks to vertical impact have appeared in the U. S. mining and petroleum literature.'-' Most published data have, to date, been obtained at atmospheric pressure, although some early high pressure information was reported by Payne and Chippendale.8 Maurer" has recently utilized available brittle impact data to develop a drilling rate equation based on the experimentally observed proportionality between crater volume and blow energy. His result agreed with earlier efforts by both Somerton, who used dimensional analysis, and Outmans, who used plasticity theory. It has long been known that rocks exhibit different modes of failure depending on the state of stress. The literature in this area is considerable; however, papers by Bredthauer, Handin and Hager,13 nd Robinson", are adequate to illustrate the point. Since rocks flow plastically at certain triaxial stress conditions, the mathematical theory of plasticity has been used to analyze the rock drilling problem. Cheatham'" has altered the wedge identation solution of Randtl to rocks, and has developed useful equations for penetrator forces under a variety of conditions. Outmans" has utilized Hill's solution in a similar manner to develop a drilling rate equation. Both Cheatham and Outmans used the linear Mohr-Coulomb rule to relate rock strength and confining pressure. The actual stress at the hole bottom is not easily ascertained, although photoelastic studies by Galle and Wil-hoit," plus the analytical treatment of Cheatham and wilhoiti8 provide some insight. Consequently it is not clear to what extent the highly idealized rheological model of a perfectly plastic solid can be realistically applied to the rock drilling problem. This paper is the first report on a long range experimental study of crater formation in rocks at elevated stress states. The data presented here are from the first phase of the project. Data obtained from impulsive wedge impacts on two synthetic, plastically deforming rocks are presented. MODEL ROCKS Geologists have long been faced with modelling the behavior of the earth and, as a consequence, have studied scaling problems in some detail.' In general, their main problem is handling the wide disparity between laboratory and geologic time. In our studies the time effects (blow velocity or rate of loading, blow duration, etc.) were essentiafly the same for both. model and prototype, as were were geometry and tooth penetration. Thus application of available scaling laws suggests that Similarity is obtained if the stress-strain curves of model and prototype are similar." For this reason Hubbert and Willis''

By R. F. Blakely, P. E. Potter

Any stratigraphic section or bedding sequence can be synthesized if there is a transition procedure from one lithology or bedding type to another, and if thickness distributions of the different lithologies are known. Stratigraphic sections of a fluvial sandstone body were synthesized with five bedding types: cross-bedding, massive beds, parting lineation, ripple mark and mudstone. The transition procedure from one bedding type to another used dependent, Markovian random processes which have a memory that extends one step backward in the depositional process. As observed in nature, median grain size and sand wave thickness (cross-bedding and ripple mark) decline upward in the synthesized sections as proportions of the different bedding types change. Grain size and permeability were also incorporated into the sections. By changing the transition procedures, bed thickness distributions, rate of upward decline or sand wave height and length, different types of sections can be synthesized, thus making it possible to model many different sedimentation problems. INTRODUCTION This paper describes a general method for synthesizing stratigraphic sections and bedding sequences of sedimentary, metamorphic or igneous origin. Synthetic generation is of interest for several reasons. Close correspondence between real and synthetic sections suggests that the factors used in the synthesizing model may indeed be the correct ones, thus giving the investigator a check on his assumptions. Rapid, inexpensive simulation of many stratigraphic sections permits one to synthesize a rock body (sandstone or carbonate reservoir) or, on a larger scale, the fill of a sedimentary basin. Harbaughl gives an example of mathematical simulation of a carbonate basin. He simulated the basin in the hope that improved prediction would follow better understanding of the depositional processes. From the petroleum engineer's viewpoint it seems reasonable to believe that the synthetic generation of rock properties and their distribution in a reservoir should be relevant in the study of reservoirs. Any stratigraphic section or bedding sequence can be generated provided there is a transition procedure from one lithology or bedding type to another, and provided the thickness distributions of the different units are known. The transition procedure involves random processes that are either independent or dependent. If the depositional process is independent, previous deposition will have no influence on present deposition. However, if it is dependent, past deposition will influence either present or future deposition. Such a dependent depositional process can be thought of as having a memory that extends backward in time through one or more pulses of deposition. A process with a memory can be described by a Markov process. Because the concept of memory or dependence appears to be in accord with our understanding of many depositional processes, Markov processes were used to synthesize the bedding sequences of this study (see Appendix). The above methods are perfectly general and are appropriate for any stratigraphic section or bedding sequence: bedding types in a beach deposit, an evolving carbonate bank or the changing lithologic fill of a thick geosyncline sequence. We chose to synthesize a vertical profile of a fluvial sandstone body because its characteristics were well documented, much was known about fluvial processes and fluvial-deltaic sandstone bodies constitute an important class of petroleum reservoirs. CHARACTERISTICS AND ORIGIN OF FLUVIAL CYCLE The fluvial cycle has been well documented in recent years by Bersier,2 Allen3-5and Visher.6,7 Deposits from fluvial cycles range from 10 to 150 ft or more in thickness and are characterized by a "fining upwards": coarse sandstones with occasional conglomerates grade upward into medium- to fine-grained sandstone, and hence into siltstone and mudstone. The dominant sedimentary structure is inclined bedding: thick cross-beds in the lower half

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 P. F. Gnirk, J. B. Cheatham

An experimental investigation has been made to study the effects of tooth angle and confining pressure on the force required for indentation of dry rock samples under confining pressure. In these experiments static single-tooth indentations were made on two standstones, a limestone, a schist and a slate under confining pressures of 5,000 to 15,000 psi. Sharp wedge-shaped teeth with included angles ranging from 30°to 105°were used. In all cases the force-displarenzent curves are approxinlately linear for rocks tested under confining pressure, as contrmted with the discontinuous curves obtained for silnilar rocks at atmospheric conditions. Photonzicrographs of thin sections through the craters indicate very little evidence of chip formation during indentation. In general the experimental data lie between calculated values for perfectly rough and perfectly smooth tooth-rock interfaces. INTRODUCTION The variables associated with ctrilling at depth are numerous and, in many instances, are not easity isolated and controlled. Also, many of these variables are interdependent; however, if the variables which are sufficiently independent of each other are held constant and are appropriately introduced, the dependent variables can be more easily studied and evaluated. It is the purpose of the present investigation to study experimentally the simplest case of bit tooth-rock interaction. In essence, the experiments have involved the static loading of sharp wedge-shaped tools on the surfaces of effectively semi-infinite dry rock samples subjected to confining pressures. During this investigation the following quantities have been maintained constant: 1. Pore pressure: atmospheric (i.e., the rock is dry). 2. permeability: although the rocks may be permeable, they are "jacketed" in such a manner that the fluid is not permitted to flow through them. 3. Temperature: 75F. 4. Rate of loading: essentially static (approximately 0.002 in./second). 5. Bit tooth: a sharp wedge-shaped tooth loaded normal to the rock surface. 6. Rock surface: smooth and flat. 7. Drilling fluid: hydraulic oil. In addition to the above conditions, the rock surface was free of a layer of cuttings. There was no indexing to previously indented surfaces, and the influence of corners of the borehole was neglected because the rock samples were effectively semi-infinite. The following quantities have been varied, and the resulting effects have been studied: 1. Rock a) Type (e.g., limestone, sandstone, etc.). b) Mechanical properties ) Necessarily dependent c) Porosity on type of rock. 2. Bit tooth (included angle of wedge-shaped tooth). 3. Depth of penetration. 4. cobfining pressure. On the basis of laboratory experiments which have shown that many rocks exhibit ductile behavior at confining pressures corresponding to depths usually reached during drilling,'.' a theory has previously been advanced which predicts the force required for penetration by a single sharp wedge tooth into a ductile* material.' The material is assumed to be isotropic, homogeneous and rigidly plastic (no elastic deformation) and to obey the Coubomb yield criterion (Figs. 1 and 2). Equations governing the force required for penelration are given for two extreme cases:** 1. A smooth, perfectly lubricated tooth-rock interface (fridionless boundary), in which the solution gives a minimum value for the penetrating force: in which F = force of the tooth, b = length of tooth edge,

By J. H. M. Thomeer

A method is presented that provides a mathematical description of capillary pressure curves and, probably, of differences in pore geometry of samples. The technique is based on the observation that the location and shape of a capillary pressure curve reflect characteristics of the pore structure of the sample. INTRODUCTION It has been noted1-5 that the observed differences in shape and location of capillary pressure curves reflect some basic properties regarding the pore geometry of the samples. For example, according to basic concepts of capillarity, the following statements can be made regarding the location and shape of a capillary pressure curve.* The location of the curve with respect to the (Vb)pc and PC axes is a measure of the interconnected pore volume and of the cross-sectional area of the pore first entered by mercury, respectively. The shape of the curve depends on the interconnection of the pores and the sorting of the pore sizes. It has been the object of this study to use this information as a starting point in an attempt to find a unique mathematical description of the PC — (Vv)p, relations of different samples. It may be expected that parameters in such a relation can then be used to describe differences in pore geometry of these samples. For a practical description of the PC — (Vb)p, relation, however, it appears desirable to describe each curve (i.e., each sample) by a limited number of parameters so that the curves taken from all samples form a family, each member of which is uniquely determined by one single combination of values of all parameters.

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 F. L. Fayers, J. W. Sheldon

It has been a matter of concern to the petroleum industry to determine what effect the capillary term has on saturation profiles, since these profiles determine ultimate economic oil recovery. In their original solution of the non-capillary two-phase flow problem, Buckley and Leverett&apos; showed that solutions of the two-phase flow equations became multiple-valued in saturation. Since it is physically unrealistic for saturation to have more than one value at a given position, Buckley and Leverett resolved this difficulty by introducing a saturation discontinuity or "shock". They evaluated the strength and position of the shock from material balance considerations. Recently Sheldon, Zondek and Cardwell&apos; have demonstrated that the method of characteristics can be used to solve this and related problems. Owing to capillarity it is physically unlikely that saturation shocks occur in an oil reservoir. Buckley and Leverett suggested that the effects of interfacial tensions should be included in the fluid flow equations, if the solutions are to give continuous single-valued saturation distributions. An analogy to this problem is encountered in the theory of wpersonic compressible flow where the introduction of the viscous term eliminates shock behavior. However, the mathematical form of the capillary term in the porous medium fluid flow equations differs from the viscosity term in the Navier-Stokes equation. This paper demons~rates that capillarity does eliminate the triple-valued behavior and indicates the magnitude of capillary and gravity effects on the saturation profile for a variety of systems. DERIVATION OF LAGRANGIAN AND EULERIAN FORMS OF THE DISPLACEMENT EQUATION The equations representing capillary, two-phase flow are given by applying Darcy&apos;s law and a material balance condition for each phase. The familiar Darcy and material balance equations are as follows. ^U*^0......<3, ^^r=°......c« S is the saturation of the water phase, and B the inclination of the porous medium. Pressures in the oil and water phases differ by the capillary pressure associated with the interfacial tensions, P.-P.=pAS)..........(5) where p, is determined experimentally. The porous medium is said to be oil-wet if P, > Po and water-wet if P, < Po. If the fluids are incompressible and the total volumetric flow rate, q, is constant, then v, f v,, = q...........(6) By eliminating P.. /»,, and v, from Eqs. 1, 2, 5 and 6 it follows that It is convenient to normalize the variables in the problem after the manner of Handy and Hadley% so that Eqs. 3 and 7 become where j,, = —, the fractional flow of water; (F, is a capillary pressure normalizing constant.) Substitution of Eq. 9 into Eq. 8 gives In this partial differential equation for an initial value