Search Documents

Jan 1, 1966

By K. K. Clark, W. H. Somerton

The feasibility of reducing pressures needed to fracture formations by heating the borehole intensively has been investigated on a laboratory scale. Reductions in fracture pressures of heated Bandera, Berea and Boise sandstone cores were observed at various confining pressures, simulating reservoir stress conditions to a depth of 5,000 ft. Both permeable and impermeable borehole boundary conditions were tested. Three heating conditions were used in the tests. One group of samples was heated in a furnace to a temperature of 600C, cooled to room temperature., and then tested for fracture pressures. Two other groups of samples were heated with a borehole heater to maximum borehole wall temperatures of 625 and 792C, then cooled and tested. In all cases, heated samples showed substantial reductions in fracture pressures compared to unheated samples. The higher temperature borehole heated samples showed the greatest reduction, amounting to 70 per cent in the extreme case. In actual application to wells, it was concluded that fracture pressures may be reduced by 20 to 50 per cent by intensive borehole heating before fracturing. INTRODUCTION Hydraulic fracturing of oil and gas wells is an accepted and effective technique for increasing well productivity in many areas, but is notoriously unsuccessful in other areas. The technique becomes less interesting for production stimulation purposes when excessive fracture pressures are necessary. Formations which do fracture but require high pressures, probably form horizontal fractures. There is good evidence to indicate that some formations do not fracture at all, but yield plastically under high pressure. Thus, any technique which will make formations easier to fracture or will reduce fracture pressures is potentially interesting to petroleum producers. Horizontal fractures generally occur in planes of weakness, such as bedding planes, and at pressures equaling or exceeding the effective overburden stress. In the unusual case where planes of weakness do not exist, a slightly higher pressure must be applied to overcome the tensile strength of the rock. In general, high fluid-loss fracturing fluids are conducive to the formation of horizontal fractures.1 Vertical fractures, when they do occur, are initiated at pressures generally lower than overburden pressure and are favored by low fluid-loss fracturing fluids or impermeable boreholes. The pressure necessary to create a vertical fracture in an impermeable, elastic medium can be predicted from the so-called "thick-walled cylinder" equation: Pf = 2 or+ st.............(1) where pf = fracture pressure or = radial stress st = tensile strength of medium. This equation indicates that to fracture it is necessary to overcome a tangential compressive stress, induced by and equal to two times the radial confining stress, plus the tensile strength of the medium. If the above equation held strictly true for rocks, very little could be done to reduce fracture pressure. The tensile strength of the rock might be reduced somewhat, but any change in st would be small compared to the magnitude of the tangential stress at reasonable subsurface depths. Laboratory tests have shown that actual fracture pressures exceed values predicted from Eq. 1 by factors of 1.4 to 28.6.1,2 Eq. 1 thus does not properly describe the

Jan 1, 1966

By J. U. Messenger

A technique is described whereby the resistance of an emudian to breaking can be quantitatively determined. Produced ailfield emulsions are usually the water-in-oil type and, accordingly, do not conduct an electrical current. However, there is a threshold of A-C voltage pressure above which an emulsion will break and current will flow. The more stable an emulsion, the higher the required voltage. A Fann Emulsion Tester, modified so that low voltages (0 to 10 v) can be accurately measured, is suitable. This technique has application in evaluating the effect of a demuksifier on the stability of an emulsion. Emulsions can, in essence, be titrated with demulsifiers by adding a quuntity of demulsifier, stirring, and measuring the voltage required to cause current to flow. Any synergistic effect of two or more materials added simultaneously can be followed accurately. A demulsifier that significantly lowers the threshold voltage (from 100 to 400 v to 0 to 10 v for the emulsions in this study) is effective and can cause the enlulsion to break. A demulsifier that will bring about this drop in the threshold voltage at low concentration ir very desirable. The technique is also well adapted for rapidly screening demulsifiers. INTRODUCTION Stable emulsions in produced reservoir fluids resulting from certain well stimulation and completion procedures are common problems. The use of suitable demulsifiers can often mitigate these difficulties. At the present time, a rapid and efficient method for selecting satisfactory demulsifiers is not available. It is badly needed. Reliance is now placed primarily on trial-and-error procedures. A new test method has been developed which permits a more rapid and precise selection of demulsifiers. It involves measuring the electrical stability potential of an emulsion before and after a demulsifier has been added. This paper describes this method and shows where it should have application in field emulsion problems. NATURE OF OILFIELD EMULSIONS Two immiscible components must be present for an emuhion to form; we are concerned here with crude oil and water. An emulsifier must be present for tin emulsion to be stable. J Emulsifiers can be substances which are soluble in oil and /or mter and which lower interfacial tension. They can be colloidal solids such as bentonite, carbon, graphite, or asphalt which collect at the interface and are preferentially wet by one of these phases. Unrefined crude oils can contain both types of emulsifiers, A popular theory is that, of the two phases in an emulsion, the dispersed phase will be the one contributing most to the interfacial tension.' Usually this phase contains the least amount of emulsifier. The stability of a water-in-oil emulsion is affected by the fol1owing: l) viscosity; (2) particle or droplet size; (3) interfacial tension between the phases; (4) phase-volume ratios; and (5) the difference in density between the phases. A stable emulsion is usually characterized by high-viscosity, small droplets, low interfacial tensions, small differences in density between its phases, and slow separatian of the phases. It also has low conductivity (high electrical stability potential). Water-in-oil and oil-in-water emulsions"' are both common; however, oil field emulsions are predominantly water-in-oil emulsions. The emulsions which commonly occur during oompletion and stimulation operations contain a combination of several of the following: acids, fracturing fluids (oil, water, acid), and formation water and oil. Produced emulsions usually contain formation water and oil. Emulsions form in oil wells because oil and water are mixed together at a high rate of shear in the presence of a naturally occurring or unavoidably produced emulsifier. During the completion and stimulation of productive zones, and while formation fluids are being produced, oil and water are very often commingled. These mixtures are formed into emulsions by agitation which occurs when the fluids are pumped from the surface into the matrix of the formation or produced through the formation to the surface. Restrictions to flow (such as perforations, pumps, and chokes)".'" increase the level of agitation; tight emulsions are more likely to form under these conditions. Often an emulsified droplet is an emulsion itself.'" Therefore, emulsion-breaking problems can be quite complex. The complexity can be even greater if a third phase (gas) is included. Demulsifiers operate by tending to reverse the form of the emulsion. During this process, droplets of water become bigger, viscosity is lowered, color becomes darker, separation of the phases faster and electrical stability potential approaches zero. Any of these effects could be followed as a means of determining emulsion stability. However, electrical stability potential is the most reproducible and most easily measured parameter for following the stability of a water-in-oil emulsion.

Jan 1, 1966

By A. R. Hendrickson, C. F. Smith

Hydrofluoric-hydrochloric acid mixtures have been successfully used to stimulate sandstone reservoirs for a number of years. Hydrofluoric acid (HF) has a specific reactivity with silica which makes it more effective than HCl for use in sandstone. Kinetics of the reactions of HF have been studied to determine the related effects of reservoir composition, temperature, acid concentration and pressure on the spending rate of HF. Secondary effects from by-product formation are noted and described. Predictions are made concerning the improvement in productivity resulting from HF treatment of skin damage. The kinetic order of HF reaction in sandstone was experimentally determined to be first order, i.e., the reaction rate is proportional to concentration. HF reacts faster on calcite than on clay, which, in turn, is faster than the reaction rate of HF on sand. Static conditions retard the HF reaction rate. As HF is forced into cores, there is a temporary reduction as a function of flow rate and acid concentration. Extensive deposition of calcium fluoride in acidized cores was not observed. Although some CaF, was defected, it was not considered a major source of damage in cores containing moderate amounts of carbonate. Other fluosilicates could be potentially more dangerous than CaF, in reducing permeabiliry. INTRODUCTION Hydrofluoric acid has been widely used in stimulation treatments since 1935, when mud acid was introduced to the petroleum industry. Originally, this hydrochloric-hydrofluoric acid mixture was intended to remove mud filter cake, but it has since been successfully applied to many other oilfield problems. Mud acid treatments have been unusually successful in sandstone reservoirs where hydrochloric acid is unreactive due to a lack of enough calcite in the formation. The relatively small amount of hydrofluoric acid present (2.1 per cent) reacts with sand grains, clays and traces of calcite which are generally present in sandstone reservoirs. Since hydrofluoric acid (HF) is the key to mud acid success, this research effort has been dedicated to gaining a more thorough understanding of the basic chemical and physical principles involved as HF reacts. Hydrofluoric acid's reactivity with silica makes it unique in application. Other mineral acids such as hydrochloric, sulfuric or nitric are unreactive with most silicious materials which comprise sandstone formations. A typical sandstone reservoir may contain 50 to 85 per cent silicon dioxide, more commonly called sand or quartz. Hydrofluoric acid reacts as follows: 4HF + SiO2 + SiFO + 2H2O The silicon tetrafluoride (SiF,) is a soluble gas, in some ways similar to CO2, and is capable of undergoing further reaction when held in solution by pressure. These reactions will be considered in detail later. Kinetics of the reactions of HF have been studied to determine the effect of reservoir composition, temperature and pressure on the spending of the acid. Secondary effects from by-product formation have been noted and described. The individual reactions of HF on quartz, glass and clay are reported. Mathematical correlations have been drawn, then applied to studies of HF spending in cores obtained from actual producing sandstone formations. The research reported herein is only the beginning of a continuing approach to better understanding and use of HF in petroleum reservoirs. THEORY AND DEFINITIONS Through the years, a concentrated effort has been made to understand the effects of many variables on hydrochloric acid (HCI) spending in limestone. Hendrickson el al., have given mathematical relationships for HC1 reactions which made possible the engineered approach to acidizing. The same variables—temperature, acid concentration, formation composition, pressure and permeability-porosity relationships—which affect HC1 behavior in limestone also govern HF behavior in sandstone. Insoluble by-products of HF reaction have been isolated and identified. Their effect on fluid flow has been measured under varying conditions in an attempt to evaluate the extent of possible damage and means of eliminating it. In general, HF follows the same reaction paths as HCI. It will react with limestone and dolomite with speed and ease. Thin sections of acidized cores show the reaction of HF with limestone or calcite faster than its reaction with either clay or sand. When HF reacts with calcite (CaCO,), theoretically, calcium fluoride (CaF2) is precipitated, and has been blamed as a major cause of reduced permeability. On the other hand, pH and pressure such as that encountered in an underground formation under acid treatment definitely retard CaF2 formation,' so the whole question of CaF, deposition in wells is a subject for study.

Jan 1, 1966

Jan 1, 1966

By B. H. Caudle, L. H. Silberberg

Reservoir depletion by natural water drive is typified by the movement of water from an aquifer into the adjacent oil-bearing formation. Prior studies of this type 01 water movement have generally neglected the resistance to flow in the aquifer and its effect on the movement of water into the oil bearing zone. A method for designing and operating scaled models of such reservoir systems is presented. Experimental data on a model of an edge-water-drive reservoir are shown and discussed. INTRODUCTION Most petroleum reservoirs derive at least a part of their productive capacity from water influx. This water may be injected from the surface or it may come from an aquifer adjacent to the oil zone. In either case, the reservoir engineer must be able to estimate the advance of the water as a function of either the elapsed time or the fluid produced. This paper describes a type of fluid flow model which includes the effect of viscous and gravitational forces in the reservoir and the surrounding aquifer. Natural water influx can be divided into three general (but widely overlapping) categories according to the direction of flow in the aquifer. These are shown schematically in Fig. 1. Fig. l(a) illustrates an edge-water-drive mechanism. In this case, water advances updip along the stratum but the movement of water is essentially horizontal and very little of the oil is actually underlain by water. The bottom-water drive (Fig. l(b)) is characterized by a thick aquifer underlying the oil zone. The water movement is generally vertical in the aquifer. The third category is illustrated in l(c). Referred to as the "thin oil column", this type of oil accumulation consists of oil over water in a relatively thin stratum. Fluid movement is horizontal in both the oil and water zones, except close to the producing well. One characteristic which all three of these types have in common is that part of the water flow takes place in the water invaded region of the oil reservoir while the remainder of the water flow occurs in the aquifer. The resistance to the flow of water usually will not be the same in these two regions. This occurs because the porous rock normally contains only water in the aquifer, while the water-invaded region always contains microscopic globules of bypassed oil (residual oil) which interfere with the flow of water through the rock. The production history of any type of a water-drive system is a function of two phenomena — where the water goes, and how it displaces the oil in the area invaded by the water. It is usually impractical to study these two phenomena at the same time. In most cases one type of analysis is used to predict the gross movement of water in the reservoir, while a different type of analysis is employed to determine the amount of oil to be recovered from the gross volume contacted by the water. The portion of the oil reservoir which is invaded by water is mainly a function of the resistances to fluid flow in the several parts of the flow system. This gross water influx in the three types of natural water drive typically results in water cusping to the production wells either from the side or from below (water coning). In 1747, Muskat described a greatly simplified mathematical model for predicting the area invaded by water in a bottom-water-drive reservoir.1 This model used the Laplace equation and suitable boundary conditions to describe the isopotentials and streamlines in the flow system. The use of this model assumes that (1) the reservoir rock is homogeneous in nature, (2) the oil and water mobilities are equal, (3) there is no oil flowing in the water invaded region, and (4) external forces — such as gravity — do not affect the flow. This method has been used for both natural water drives and water injection projects. In this, as in most models, the method of images is used to reduce the size of the model necessary to describe the reservoir. The assumption of equal viscosities and permea-

Jan 1, 1966

By K. E. Brown, A. R. Hagedorn

A 1,500-ft experimental well was used to study the pressure gradients occurring during continuous, vertical, two-phase flow through 1-in., 1 1/4-in. and 1 1/3-in. nominal size tubing. The test well was equipped with two gas-lijt valves and four Mathak electronic pressure transmitters as well as with instruments to measure the liquid production rate, air injection rate, temperatures and surface pressures. Tests were conducted for widely varying liquid flow, rates, gas-liquid ratios and liquid viscosities. From these data, an accucrate pressure-depth traverse was construcred for each test in each of the three tubing sizes. From the results of these tests, correlations have beer] developed which allow the accurate prediction of flowing pressure gradients for a wide variety of tubing sizes, flow conditions, and liquid properties. Also. the correlations and equations which are developed satisfy the necessary condition that they reduce to the relationships approprate to single-phase flow when the flow rate of either the gas or the liqutid phase becomes zero. All the correlations involve only dimensionless groups. which is a condition usually Sought for in similarity analysis but not always achieved. The correlations developed in this study have been used to calculate pressure gradients for pipes of larger diameter than those upon which the correlations are based. Comparisons of these calcuclaed gradients with experimetally determined gradients for the same flow conditions obtained from the literature indicate that extrapolation to these larger pipe sizes is possible with a degree of accuracy sufficient for engineering calculations. The extent of this extrapolation can only he determined with additional data from larger pipe diameters. INTRODUCTION The accurate prediction of the pressure drop expected to occur during the multiphase flow of fluids in the flow string of a well is a wideiy recognized problem in the petroleum industry. The problem has been brought even more into prominence with the advent of tubingless or slim-hole completions which use small-diameter tubing. Many of the correlations which give reasonably accurale results in the lager tubing sizes are greatly in error when applied to small-diameter conduits. Small-diameter conduits are defined as 1 1/2-in. nominal size tubing or smaller. The study of the pressure gradients which occur during multiphase flow of fluids in pipes is exceedingly complex because of the large number of variables involved. Further difficulties relate to the possibility of numerous flow regimes of widely varying geometry and mechanism and the instabilities of the fluid interfaces involved. Consequently, a solution to the problem by the approach normally used in classical fluid dynamics based on the formulation and solution of the Navier-Stokes equation has not been forthcoming. This is primarily the result of the nonlincarities involved and the difficulty of adequately describing the boundary conditions. As a result of the foregoing, most investigators have chosen semi-empirical or purely empirical approaches in an efford to obtain a practical solution to the problem. Much of the previous work in this area was done in short-tube models in the laboratory. A number of problems arise, however, when attempts are made to extrapolate these laboratory results to oilfield conditions where a much longer tube is encountered. In those few studies where data taken in long tubes were utilized, the data covered only a limited range of the variables, and as a result inaccuracies are introduced when the correlations are extended outside the range of the original data. Also, as a consequence of the limited amount of data available for these studies, the effects of several important variables were overlooked. The problem of predicting the pressure drop which occurs in multiphase flow differs from that of single-phase flow in that another source of pressure loss is introduced, namely, those pressure losses arising from slippage between the phases. This slippage is a result of the difference between the integrated average linear velocities of the two phases, which in turn is due to the physical properties of the fluids involved. In contrast to single-phase flow. the pressure losses in multiphase flow do not always increase with a decrease in the size of the conduit or an increase in production rate. This is attributed to the presence of the gas phase which tends to slip by the liquid phase without actually contributing to its lift. Many investigators have attempted to correlate both the slippage losses and the friction losses by means of a single energy-loss factor analogous to the one used in the single-phase flow problem. In an approach of this type, hourever, many of the important variables, such as

Jan 1, 1966

Jan 1, 1966

By R. G. Agarwal, Ramey Jr. H. J., Al-Hussainy R.

Although it has long been realized that gas recovery from a water-drive gas reservoir may be poor because of high residual saturations under water drive, it appears that only limited infomlation on the subject has been available until recently. This study was performed to show the qiiantitative potential importance of water influx. Results indicate that gas recovery may be very low in some cases: perhaps as low as 45 per cent of the initial gas in place. Gas recovery under water drive appear to depend in an important was on: (I) the prodirction rate and manner of production; (2) the residual gas saturation; (3) aquifer propertie.); and (4) the volumetric displacement effciency of water invading the gas reservoir. The manner of estimating water-drive gas reservoir recovery can vary considerably. Examples are: the steady-state tnethorl. the Hurst modified steady-state method, and various unsteady-state methods such ac. those of van Ever-dingen-Hurst, Hurst, and Carter-Tracy. The Carter-Tracy water influx expression was used in this study. In certain cases, it appears that gas recovery can be increased significantly by controlling the production rate and manner of production. For this reason, the potential importance of water influx in particular gas reservoirs should he investigated early to permit adequtrtr planning lo optirtize the pay reserves. INTRODUCTION In recent years, the economic importance of natural gas production has become increasingly apparent. This has been evidenced by more intensive exploration efforts aimed at gas production, and exploitation of both deep, as well as low-permeability gas reservoirs. Technical developments such as deep-penetration fracturing have made development of such formations economically feasible. Unfortunately, water influx has forced abandonment of a number of gas reservoirs at extraordinarily high pressures. Although reservoir engineering methods for estimating water influx have long been available, it appears that application of these methods to the water-drive gas reservoir has been sporadic.'a Available methods for estimating water influx which can be applied to the water-drive gas reservoir problem include the steady-state method,1 the Hurst modified steady-state method and various unsteady-state methods such as those of van Everdingen-Hurst. Hurst, and Carter-Tracy. Interesting applications of these solu- tions to gas reservoir and the aquifer gas-storage problems have appeared recently.3,12,14 The experimental study of residual gas saturations under water drive by Geffen et al. in 1952 indicated that residual gas saturations could be extremely high. A value of 35 per cent of pore volume is often used in field practice when specific information is not available. The study of Geffen et al. showed that residual gas saturation might be much higher in some cases. Naar and Henderson concluded that the residual non-wetting phase saturation under imbibition should be about half of the initial non-wetting phase saturation. The Naar and Henderson result that residual gas saturation under water influx should be about half the original gas saturation is recommended as an estimate if laboratory measurements are not available. Thus, it is clear that a considerable portion of the initial gas in place might be trapped in a water-drive gas reservoir as residual gas at high pressure. A full water-drive would result in loss of residual gas trapped at initial reser.voir pressure. Consideration of transient aquifer behavior leads to the conclusion that high-rate production of water-drive gas reservoirs could result in improved gas recovery by reduction of the abandonment pressure. However, there appears to be little quantitative information on this possibility. One of the few advantages of water-drive gas production appears to be improved deliverability through water-drive support of the reservoir pressure. There may also be an advantage in higher condensate recovery caused by pressure maintenance for gas-condensate water-drive reservoirs. In view of the preceding, this study was made to assess the potential importance of water-drive in gas reservoir engineering. The Carter-Tracy approximate water-influx expression was used because this equation offers some advantages in hand-calculation which do not appear to have been generally recognized.' However, calculations were performed in the main with a high-speed digital computer to permit evaluation of the effect of water-drive under a large variety of conditions. CALCULATION METHOD Water-drive gas reservoir performance can be estimated in a manner completely analogous to oil reservoir calculations: a materials balance is written for the reservoir, and a water influx equation is written for the aquifer. Siniltaneous solution provides the cun~ulative water influx and the reservoir pressure. When reservoir performance data (gas produced and reservoir pressures) are available, it is usually possible to match performance data to determine the initial gas in place and the water influx parameters —

Jan 1, 1966

By B. S. Gottfried

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By Harry A. Wahl

Precent fracture design procedures are bared on the total fracture area created. A method to distinguish beI,,.ecn total area and [he propped or effective fracture area has not been available. This paper presents a solution to this problem, applicable to horizontal fractures. The difierences between effective fracture area and torn] area are demonstrared in example calculations. This work is hayed on experimentally determined transport efficiencier of solids in sand-liquid slurries. Newtonian and non-Ne~vtoninn systems are considered. INTRODUCTION Fifteen years after commercial introduction, hydraulic fracturing remains the most successful stimulation technique in the oil field. This success is primarily due to ability of induced fractures to penetrate and alter permeabilities deep within formations. Many fields producing today could not have been developed without the hydraulic fracturing process. Because of wide usage, fracture-treatment design has received a great deal of engineering and research effort. This work, resulting in improved equipment and materials, has increased the benefits from fracture treatments as well as the applicability of the process. A major contribution was the development of fluid-loss additives. Necessarily, the number of parameters to be considered in treatment design has steadily increased, resulting in more complicated design techniques. Almost all present design procedures are based on the precepts set forth by Howard and Fast. Relating the fluid volume lost into the formation, the volume required in extending the fracture, and the total slurry volume injected, they developed an expression for the total fracture area created in terms of pertinent treatment parameters. Fluid loss during treatment was expressed as a function of time for three flow mechanisms. Although modifications of fluid loss equations have been made, the total fracture area concept has remained unaltered. A vast amount of field data indicate that induced fractures must be propped and held open to be effective. A notable exception is the Mesa Verde formation in the San Juan basin. However, analysis of these treatments shows that improved well productivities are obtained when propping agents are incorporated in the treating fluid. Although propped fracture area has been recognized as an important design parameter, a method to distinguish between total area and effective fracture area has not been available. The necessary information on slurry-sand transport in fractures has been lacking. Interest in the propped region of induced fractures is not restricted to areal extent alone. The distribution of sand within fractures is important from the standpoint of fracture flow capacities. Flow capacity affects the increase in well productivity after stimulation. The work of Huitt and Darin4 hows that partial monolayers of sand have large flow capacities compared to thick* dense sand packs. It has been postulated that gelled fluids have the ability to transport sand within the fractures at the deired low concentrations. An early contribution in the area of sand placement in fractures was made by Kern et al.' They studied sand movement in a transparent vertical fracture model. It was observed that the sand tended to settle out in the bottom of the model before moving very far. When the fluid velocity exceeded a certain critical value, all of the sand injected began moving through the crack even though it settled to the bottom. This critical velocity was determined under several flow conditions. Some work on sand movement in horizontal fractures has been reported in Russian publications. Sand movement was studied by Izyumova and Shan'gin' using a transparent "pie-shaped" flow model to simulate a horizontal radial flow system. However, the data were limited, especially in a quantitative sense. Dorozhkin, Zheltov and Zheltow studied the behavior of sand-liquid slurries in a horizontal linear flow model. The quantitative data were restricted primarily to the thickness of sand deposits formed at the bottom of the fracture. An earlier paper provided basic data on the flow of sand in horizontal fractures. This study was designed to yield specific quantitative information on rate of advance of sand particles and pressure behavior under various flow conditions. A comprehensive photographic study was undertaken in a 10-ft windowed flow cell to provide the necessary qualitative and quantitative data. Since the number of potential variables far exceeded the capacity of the initial study, emphasis was placed on the effects of sand concentration, oil viscosity and oil flow rate. A detailed description of these experiments and the results are described in Ref. 9. However. the implications of this work on the fracture design calculations were not discussed. An analysis of these data as well as new data is provided in the following sections. EXPERIMENTAL RESULTS The primary objective of the experimental investigation was to provide information on the rate of advance of the solids in sand-liquid sturries. A 10%-ft long transparent

Jan 1, 1966

By W. A. Gibson, A. D. Trujillo

The Santa Rita copper deposit first served as a source of native copper for Indian implements and weapons. In 1801 Santa Rita copper, trans- ported by mule train to Chihuahua, began to be used commercially for making coins. It was a small- scale operation made difficult by Indian raids, transportation costs and depletion of the higher grade ores. There was, moreover, no known method for economically extracting copper from the under- lying low grade sulfide ores. Originally, copper and copper oxide ores were mined by underground methods, but there was virtually no development of the deposit-which became United States territory after the Treaty of Guadalupe Hidalgo in 1848-for almost a hundred years. A number of small underground mines had been started during that time, but it was not until 1873 that all the principal claims were consolidated under one control.

Jan 1, 1966

By J. C. Melrose

The contact angle is one of the boundary conditions for the differential equation specilying the configuration of fluid-fluid interfaces. Hence, applying knowledge concerning the wettability of a solid surface to problems of fluid distribution in porous solids, it is important to consider the complexity of the geometrical shapes of the individual, interconnected pores. As an approach to this problem, the ideal soil model introduced by soil physicists is discussed in detail. This model predicts that the pore structure of typical porous solids will lead to hysteresis effects in capillary pressure, even if a zero value of the contact angle is maintained. The model is generalized to situations in which the contact angle takes on values between zero and 40 degrees. For the imbibition branch of the capillary- pressure function, the model predicts a considerable departme ftom the usually assumed cos 0 relationship. In fact, according to the model, it is possible that a displaced wetting phase will not be able to reimbibe, even when the contact angle does not exhibit hysteresis. INTRODUCTION The subject of capillarity in porous media has long been of interest in many branches of engineering and applied science. The earlier investigators were those concerned with the physics of soils. 1-5 More recently, petroleum engineers and others dealing with the problems of petroleum production from reservoir rock have given much attention to the subject.6-10 Also, important additions to the literature of capillarity have been contributed from the field of chemical engineering.11-12 These attest to the wide range of industrial applications in which capillary phenomena play a role. The present paper is concerned with the role which wettability plays in capillary action in porous media. As is well known, capillary-rise (or capillary pressure) phenomena have frequently been interpreted with the aid of a model employing the concept of a cylindrical capillary tube. This approach has enjoyed a certain degree of success in correlating experimental results.13 The generalization of this model, however, to situations which involve varying wettability, has not been established and, in fact, is likely to be unsuccessful. In this paper another approach to this problem will be discussed. A considerable literature relating to this approach exists in the field of soil science, where it is referred to as the ideal soil model. Certain features of this model have also been discussed by Purcel114 in relation to variable wettability. The application of this model, however, to studying the role of wettability in capillary phenomena has not previously been attempted in detail. In the present paper, additional features of the model are introduced. These features are critical in determining the quantitative behavior of the model. GENERAL FEATURES OF CAPILLARY HYDROSTATICS BASIC PRINCIPLES When the interstices of a typical porous solid are occupied by two or more immiscible fluid phases, the fluids are microscopically commingled. Hence, fluid-fluid interfaces are found within a certain fraction of the pore openings. The fundamental equation of capillarity specifies the configuration of these fluid-fluid interfaces. This is known as the Laplace equation, when derived from mechanics, and as the Gibbs-Kelvin equation when derived thermodynamically.15 Given two fluid phases, a and ß, in hydrostatic equilibrium, the Laplace equation states that the respective fluid pressures in regions close to the interface are related by Here oaß is the surface or interfacial tension and rl and r2 are the principal radii of curvature of the surface or interface. The pressure difference, Pa-Pß, is the capillary pressure, PC. As Buff15 has shown, Eq. 1 states the condition for hydrostatic equilibrium within the two-phase confluent region, which is referred to as the interface. It thus can

Jan 1, 1966

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

By R. L. Loofbourow

The objectives of this report are (1) to indicate what is known of natural rock temperatures as deep as 5000 ft below the surface, (2) to indicate methods for estimating such temperatures in localities where directly measured temperatures are un- available, and (3) to outline as far as possible, areas in the 48 states within which rock temperatures are likely to be low and high. The usual object of temperature measurements in deep mines and wells is to learn the temperature of air or another circulating medium. This usually differs significantly from the undisturbed rock temperature. The subject is complicated further by the variety of the thermal properties and structures of rock masses, the migration of fluids within them, the diversity of the history of rock masses and of the processes which have been at work on them during that history.

Jan 1, 1966

By Allen L. Hatch

Perhaps the greatest challenge to the non-ferrous extractive industry, aside from discovery and cost reduction, is the continuous effort toward the total recovery of metals. This is-and must continue to be-a joint endeavor on the part of the miner and the smelter, the first to trend constantly toward that form of concentration most amenable to the smelter and the latter to upgrade technical efficiency on the product delivered. What follows will be a smelter's-eye view of the relative value of gold, silver, copper, lead and zinc as they each are contained in copper, lead and zinc concentrates. The variety of combinations is obviously of such a magnitude that it will be necessary to restrict the observations to fundamental concepts.

Jan 1, 1966

By R. Simon, D. J. Graue

This paper presents correlations for predicting the solu-lility, swelling and viscosity behavior of CO2-crude oil sys8i.m~. The correlations were developed from experimental data obtained by the authors. These data are also presented. The data were determined by measuring the properties of mixtures of CO, and nine different oils. Experiinental conditions covered a range of 100 to 250°F and pressures up to 2,300 psia. Properties predicted by the correlations have average deviations, expressed as per cent of experimental value, of 2 per cent for solubility, 0.5 per cent for swelling and 12 per cent for viscosity. INTRODUCTION Interest in CO2 injection as an oil recovery process has led to the development of performance prediction methods which can be applied to specific reservoirs.1 iS To use these performance prediction methods, it is necessary to know the solubility, swelling and viscosity properties of CO2-crude oil mixtures at reservoir conditions. Some information on these properties has appeared in the literature; however, this information did not cover the range of different oils and conditions needed to prepare generalized correlations for reservoir engineering purposes. Consequently, an experimental program was undertaken to collect the data needed. The data obtained and the correlations developed from the data are described in the following sections of this paper. SOLUBILITY OF CO2 IN CRUDE OILS CO2 solubility data in the literature come from six principal sources. The solubility prediction method of Welker and Dunlop3 is limited to 80F. The information in Ref. 4 is of two types: the first includes binary and ternary mixtures of CO, and light hydrocarbons (C1 to C6), and the second gives data for CO1 and heavy hydrocarbons for a temperature range of 40 to 90F. Ref. 5 contains a KCO2 chart for systems whose convergence pressure is 4,000 psia. The KCO2's are based mainly on CO2-natural gas mixtures. Poettmann's work covered CO2 solubility in one condensate and one crude oil6,7 Ja-coby and Rzasa measured CO? solubilities as a function of pressure and temperature for two natural gas-absorber oil mixtures and two natural gas-crude oil mixtures.'8 CO, concentration in these four systems was fixed at 5 mol per cent. The work reported in this paper extends CO, solubility data to a variety of different crude oil types in a temperature range from 110 to 250F and pressures up to 2,300 psia. The experimental procedure used by the authors to obtain the solubility data consisted of combining known amounts of pure CO, and crude oil in a visual cell at a fixed temperature and measuring the bubble point of the mixture. Measurements were made for a total of 40 different CO2-oil mixtures and the results are shown in Table 2. The mixtures included nine different oils (seven crude oils and two refined oils) whose properties are listed in Table 1. All nine oils had vapor pressures less than 1 atm at the experimental temperatures. Consequently, analysis of the bubble-point vapor showed a CO, concentration over 99 mol per cent. At no time during these experiments was a second, more dense, liquid phase observed. The solubility correlation which was developed from the data in Table 2 is presented in Figs. 1, 2 and 3. In these figures, solubility is expressed as xco2 the mol fraction of CO, in the CO2-oil mixture. Fig. 1 shows solubility as a function of CO2 fugacity6 and temperature. Fig. 2 shows the same solubility data expressed as a function of saturation pressure and temperature. The solubility shown in Figs. 1 and 2 is for an oil whose UOP characterization factor is 11.7. UOP characterization factors of crude oils can be determined from Ref. 10 if the viscosity and APT gravity of the oil are known. Fig. 3 gives the solubility correction factor for oils whose UOP characterization factors differ from 11.7. The solubility correlation in Figs. 1, 2 and 3 predicted

Jan 1, 1966

By C. D. Stahl, S. M. Farouq Ali

This investigation attempts to describe and simulate the alcohol displacement process by means of a cell model, as employed in chemical engineering practice. The proposed model is more simple than previously proposed models. and utilizes parameters chosen on a theoretical basis. The model successfully reproduced the formation of the stabilized bank and the breakthrough of alcohol, the latter depending on one of the model parameters, which may be correlated with the length of the porous medium. Moreover, the effects of the phase behavior of the liquid system involved, as observed in experimental studies, were reproduced. Several variations of the basic model were devised and tested on a digital computer. These included the cases in which: (1) the actual value of fractional flow was used in cell-to-cell computations; (2) the number of cells was varied within the same run; and (3) incomplete rather than complete phase equilibrium was assumed within each cell. The proposed cell model clarifies the basic mechanism of the process. Detailed concentration profiles obtained for each cell, for instance, showed the mechanism of bank formation in relation to the phase behavior characteristics. The results obtained indicated a varying degree of phase equilibrium concommitant with changes in the velocities of the phases in an actual alcohol displacement. This condition was approximated by changing the number of cells during the simulation. Interesting information was obtained on the influence of path length on the efficiency of alcohol displacement, which has been the subject of some controversy. Certain limitations preclude the use of the proposed model as a substitute for experimental studies. The results obtained were, nevertheless, of value in interpreting the experimental results. INTRODUCTION During recent years considerable effort has been directed toward an understanding of alcohol displacement, the process whereby oil and water are recovered from a porous medium by the continuous injection of a solvent. The complex nature of the physical process involved has so far defied a complete mathematical treatment. Other methods of approach, amounting to an overall material balance, have been proposed, yielding useful information on certain aspects of the process.1-3 Taber et al, in particular, defined the displacement mechanism in terms of the phase behavior of the alcohol-oil-brine system invoIved.2 Wachmann reported a mathematical treatment of alcohol displacement subject to certain simplifying assumptions? Donohue proposed the use of a "cell model" for simulating alcohol disPlacement.5 The nature of the assumptions involved limited the utility of the model. The present work attempts to examine the variables involved in the simulation of alcohol displacement and discusses several possible versions of the basic cell model. Under certain conditions the model results are similar to the experimental results. In particular, the spontaneous formation of of the stabilized bank and the effects of the system phase behavior were successfully reproduced. PREVIOUS WORK ON CELL MODELS Cell models and the theoretical plate concept are often used in solving chemical engineering problems in which an explicit mathematical solution may be difficult or impossible to obtain. Examples of such applications occur in distillation, gas-liquid chroma tography,6 reactor technology, absorption, etc. In petroleum engineering, such a model was used by Attra7 to describe non-equilibrium gas drive, and by Higgins and Leigh ton8 to calculate sweep efficiency in water flooding. Aris and Amundsen pointed out the equivalence between the diffusion model and perfectly mixed cells connected in series.9 Deans proposed a three-parameter cell model to simulate two-component

Jan 1, 1966

By R. F. Nielsen, R. G. Bentsen

Jan 1, 1966

By R. D. Vaughn

The analysis of laminar flow of power-law non-Newtonian fluids in narrow, eccentric annuli is employed in this paper to discuss the problems of lubricant flow in journa! bearings and of errors introduced by eccentricity in experimental studies with concentric annuli on extruders and wellbore annuli. The velocity profile and pressure loss-flow rate equations are developed for the laminar flow region. In addition, the expected error in flow rate and pressure-loss measurements for concentric annuli as a result of eccentricity is determined. For example, a 10 per cent displacement of the core of an almost concentric annul us would cause a 1.8 per cent decrease in the observed pressure loss for a fluid with a power-law exponent n of 0.25. The corresponding increase in the observed volumetric flow rate would be 7. 5 per cent. INTRODUCTION Non-Newtonianism and eccentricity occur simultaneously in two engineering problems: (1) flow of lubricants in journal-bearings and pressure-reducing bushings, and (2) flow of non-Newtonian fluids in plastic extruders and wellbore annuli. The lubricants used for moving parts are often non-Newtonian in character — often they are plastic in behavior. A solution to the problem of flow of non-Newtonian fluids in narrow eccentric annuli is particularly pertinent to this problem. In all experimental studies of laminar flow of fluids in concentric annuli, such as in extruders and well casings, the error due to eccentricity must be estimated or studied. A number of publicationsl-4,6,7 have dealt with this problem for Newtonian fluids; however, I am not aware of work for non-Newtonian fluids. This work is directed to the non-Newtonian problem. Before the solution to the problem is given, the pertinent conclusions from the work on Newtonian fluids will be reviewed. Heyda 1 and Redberger and Charles 2 have published general solutions to the problem of the laminar flow of Newtonian fluids in eccentric annuli, apparently without knowing of the earlier work of Ca ld well 3 and Bairstow and Berry,4 which is reported by Dryden, et al. 5 Although several mathematical routes are encompassed by the work of these authors, the results appear to be equivalent. Redberger and Charles2 show that the error caused by eccentricity in concentric annuli is negligible for small diameter ratios (k < 0.5); however, for large diameter ratios (k + l), the error in the predicted flow rate can be as great as 100 per cent or more. Partial solutions to the problem are available from the work of Dryden,5 Tao and Donovan6 and Piercy, et a1.7 Tao and Donovan examined the case of flow in narrow, eccentric annuli (k + 1) with and without rotation of the annular core. These authors also reviewed previous work on this subject and verified their approach with experimental data. Dryden gives the solution for the limiting case of complete eccentricity or tangency. Piercy, et a l. published an early solution to the problem of narrow eccentric annular flow. The conclusions of Redberger and Charles2 and the experimental proof of Tao and Donovan6 both suggest that the region of large diameter ratios (k + 1) is of main interest and that the parallel planes approximation to the solution in this region is satisfactory. This method will now be extended to the laminar flow of non-Newtonian fluids in narrow eccentric annuli. THEORETICAL SOLUTION The geometrical aspects of the problem are illustrated in Fig. 1. To represent the non-Newtonian fluid the power-law model was selected. This model has many disadvantages which have been pointed out; nevertheless, its simplicity, its frequent and wide applicability justify its use in this work. Fredrickson and Bird8 and Savins9 have used it as a basis for a theoretical study of laminar flow of non-Newtonian fluids in concentric annuli.

Jan 1, 1966