Search Documents

By E. B. Brauer, W. Hurst, J. D. Clark

Because of drilling, completion, and workover practices, the permeability around a wellbore generally is different from the permeability of the formation. The zone with the altered permeability is called "skin" and its effect on the pressure or flow behavior of the well is called "skin effect". Hawkinsl has shown that the radius and the permeability of this zone are related to the skin by: If the permeability in the skin zone is less than that of the formation, the skin is positive; if it is more than that of the formation, the skin is negative. If the two permeabilities are equal, s is zero; that is, there is no skin. Van Everdingen2 and Hurst3 have given mathematical solutions to the case of a zone of reduced permeability around the wellbore. This skin effect is illustrated in Fig. 1. They treat the positive skin as a zone of reduced permeability of infinitesimal thickness around the wellbore. When applied to a well with negative skin, however, their solutions lead to the calculated flowing well pressure which is smaller than the formation pressure, i.e., an injection situation. This is a physical contradiction. In this paper, it is shown that this mathematical difficulty can be overcome by assuming an effective wellbore radius larger than the actual wellbore. Existing solutions2"4 are modified to include the effect of an enlarged wellbore radius, and thus enable the engineer to deal with negative as well as positive skins. Theoretical Development Van Everdingen and Hurst,4 through the use of Laplace transformation, obtained the following expression for transient fluid influx into a wellbore: The values for QtD, reported extensively by van Everdingen and Hurst,4 are obtained from the following integral:

Jan 1, 1970

By W. T. Malone, J. R. Williams, R. L. Tiner, J. M. Tinsley

Hydraulic fracturing methods for production stimulation have become a common procedure in the oil and gas industry. Fracturing treatments are performed on wells of various potentials to help increase the production rate. In addition to the desire for increased production is the need to predict what increase might be expected. This knowledge is useful in economical treatment planning to accomplish the desired production goals for the well. Theoretical production increases due to horizontal radial fractures may be calculated using equations developed from mathematical analyses and electrical models.1,2,6 It also is possible to predict theoretical production increases due to vertical fractures using curves developed from electrical3 and mathematical models.8 However, these studies assumed that the vertical fracture height and the formation height were equal, which may not be true." In addition, the final fracture height through which fluid is produced may not be equal to the created height. The development of vertical fracture proppant transport relationships has made unrealistic the assumption that propped fracture height is equal to created fracture height. Through the use of published proppant transport equations,4'12 it is possible to calculate the height of a bed of propping agent deposited in a fracture of a given height. Assuming that the fluid flowing in the fracture moves through only the propped portion of that fracture (i.e., that the un-propped part of the fracture heals) it becomes essential to know how production increase varies with the ratio of propped fracture height to formation height. It should be recognized that future investigation may determine that the fracture does not heal completely and that a highly permeable channel exists over the deposited bed of proppant. The effect of this channel was not considered for the present study. The following analysis of experimental data allowed development of curves to determine theoretical production increases considering fracture height as a variable. Discussion In designing a fracturing stimulation treatment a number of variables must be considered. The size of the treatment, type of proppant used, proppant concentration, treatment fluid properties and well conditions are all factors affecting the dimensions and permeability of the fracture. Production increase due to a vertical fracture is affected by propped length, height, width and permeability of the fracture for a given set of formation conditions. Production increase may be defined as the ratio of production rate after fracturing to the production rate before fracturing at the same pressure differential, pe — pw. The producing rate in an

Jan 1, 1970

By H. H. Keller, E. J. Couch

Damage to production wells subjected to the high temperatures associated with in-situ combustion projects has been a problem since the advent of this thermal recovery technique. Injection of water down the annulus of hot wells has been successfully utilized for cooling and prolonging the life of such wells. A mathematical model that is suitable for predicting the effectiveness of cooling water in a particular application is developed. Typical results for practical ranges of the various parameters involved are presented. INTRODUCTION In the normal course of in-situ combustion operations, highly damaging conditions develop in the production wells as the combustion zone approaches, and these wells are subjected to high temperatures (up to 1500F). The resulting well damage is recognized as one of the more serious operating problems unique to the in-situ combustion process. Even in operations where well completions have been specifically designed to cope with this problem, the question arises as to how this destructively high temperature can be kept within tolerable limits. A simple and economical approach to this problem is to circulate cooling water down-hole. This method has been used with some success in shallow reservoirs. At greater depths, however, the heat exchange with the overburden surrounding the well and between the fluids in the tubing and the annulus tends to equalize the temperatures of the two fluid streams so that the cooling capacity delivered to the bottom of the hole decreases with well depth. This method would provide no appreciable control over the bottom-hole temperature at great depths. In order to delimit the condition under which cooling by water circulation is effective, the wellbore temperature must be related to the water inlet rate and temperature; to the reservoir depth, production rate, and initial temperature; and to the combustion at at bot tom-hole depth. A model that meets these requirements and that will offer an approach for developing engineering solutions to practical design and operating problems is described in this paper. MATHEMATICAL MODEL Previous studies of temperature distributions include those for flowing wells1p2 and those for hot fluid injection. Adams and Haddenhorst 6 presented solutions for the down-hole circulation of hot oil in the treatment of paraffin-deposition problems. In this paper the thermal model presented by Ramey 4 is adapted to the flow system shown schematically in Fig. 1. Water at temperature Ti is injected at a rate wi and flows down the annulus outside the tubing. It mixes at bottom-hole depth D with the fluids produced at race w from the hot reservoir at temperature T. This mlxtureJ. of injected water and reservoir fluids flows up the tubing and is produced at the surface at TI,. In order to describe the temperature response in this system, we make use of the usual conservation of

Jan 1, 1969

By R. L. Robinson, P. DeFrank, R. F. Hatfield

Production equal to or greater than open-hole completions is possible through perforated completions if the flow paths throughout the perforations are free of obstructions.' Previous investigations have indicated that, in some cases, maximum productivity is not obtained due to plugging of the perforations with debris from the wellbore fluid and/or debris from the perforating process. Control of wellbore fluid and proper well completion techniques are very important in the elimination of perforation plugging.' Selection of properly designed perforating equipment also is quite important for maximum well productivity. Studies made with the use of the perforation flow laboratory, in which perforating tests are run under simulated well conditions, have shown that perforation plugging occasionally occurs. Plugs in perforations result from various causes. One of the contributors to plugging that can be controlled by the perforating device is the slug or carrot resulting from jet perforating. This condition has undergone extensive study, and previous developments have resulted in a decrease in carrot size—but not in complete elimination. It was felt that, if a method of completely eliminating the carrot could be found (with charge performance equal or better than previously obtained), a significant contribution would have been made. Such a charge has been developed and is reviewed in this paper. THEORETICAL CONSIDERATIONS OF CARROT ELIMINATION There are several popular theories of how a shaped charge jet is formed."," One of these is the hydrodynamic theory, wherein the liner is treated as a fluid under the impact of the explosive forces. Another is the spalling theory, where the interaction of shock waves in the liner causes small particles to be torn from the inner surface of the liner and projected at high velocity toward a common focal line to form the jet. It is not unreasonable to assume that both of these theories may be partially correct. For purposes of our discussion, it is simpler to use the second theory. Generally speaking, when a plate of metal is attacked by an explosive shock pulse, the pulse is transmitted through the metal at high velocity, as shown in Fig. 1(A). Pressures at the shock front may be 6 X 10' psi or higher. This shock wave (a compression wave) upon striking a free surface, Fig. 1(B), is partially reflected back as a tension wave. The interaction of the tension wave with the decaying compression wave causes a stress concentration in the metal; when this stress becomes great enough, a failure occurs and a small particle flies off and is projected away from the metal surface, as shown in Fig. 1(C). A new free surface is then presented to the shocks; the process repeats itself until there is no longer enough energy to cause failure, leaving a portion of the metal relativcly untouched. See Fig. 1(D). To apply this theory to a perforating shaped charge, consider the flat plate to be folded into a somewhat conical shape, with a column of explosive placed in intimate contact with the cone (Fig. 2). An explosive initiator detonates the charge, creating a shock wave

By J. T. Raleigh, D. L. Flock

This study investigated the nature of formation plugging with bacteria and attempted to relate its characteristics to physical rock properties. Fifteen core samples of four specific formation types were defined and plugged using a commonly occurring, uniform sized dead bacteria, Bacillus subtilis. Two of the formation types were fairly uniform-grained sandstones and two were heterogeneous carbonates. The injection rate and concentration of solids in the brine were held constant during the test runs. The pore geometry factor G was a significant petro-physical rock characteristic for correlation with depth of plugging, but lacked importance as a parameter for rate of plugging. Although the G factor may not completely describe porous rock geometry for plugging predictions, for reservoirs with multi-pore systems, there is evidence that the pore geometry factor might form a basis for establishing injection water specification in more llornogeneous reservoirs with only primary porosity development. INTRODUCTION GENERAL PURPOSE Injecting water into oil-bearing formations to increase recovery is common practice in the petroleum industry. One of the most serious operational problems involved with water injection wells is the tendency of the porous rock immediately adjacent to the injection wellbore to plug partially. This results in increased injection pressures and/or a decrease in the rate of injection. This plugging in injection wells can result from a number of causes, including bacteria contained in the injected water.1 The purpose of this work was to investigate the nature of injection plugging with bacteria and to attempt to relate its characteristics to physical rock properties. Many variables exist in research of this type, so, to reduce them, injection rate and concentration of bacteria in brine were held constant. The bacteria used for plugging was a uniform-sized, dead bacteria, Bacillus subtilis, a micro-organism commonly found in water, air and soils. Four specific porous rock types were plugged—two consolidated sandstones and two heterogeneous carbonates. The experimental work involved two distinct phases. The first concerned determining the physical rock properties and the second was the experimental work of actually plugging the selected core samples. PORE GEOMETRY FACTOR Capillary pressure is necessary to displace a wetting fluid from a capillary opening with a nonwetting fluid. The factors which govern capillary pressure are the surface or interfacial tension of the fluids involved, the system wettability and the equivalent pore radius. Thus: P,= 2ucososrd..........(1) Thomeer used the idea that the location and shape of the capillary pressure curve reflect characteristics of the pore structure of any porous media sample.' He presented the mathematical description of capillary pressure curves and of indicated differences in pore geometry of samples. where G = pore geometry factor PC = capillary pressure Pi = displacement pressure (V,)p, = fractional bulk volume occupied by the displacing fluid at any capillary pressure (V,), = fractional bulk volume occupied by the displacing fluid at infinite pressure or total interconnected pore volumes. The shape of the curve is defined by the pore geometry factor G. The curve shape is related directly to the sorting and interconnection of the pores of the sample. The location of an infinite number of curves can be defined by the same asymptotes. However, these curves differ in shape, each curve being described by a specific value of G. Thus, the description of a curve is unique when, in addition to the position of its asymptotes, its shape is defined. These curves can be used graphically to determine the characteristic parameters of any arbitrarily located curve. Those parameters that can be approximated from any capillary pressure curve are: (1) total interconnected pore volume, represented by (V,),- ; (2) extrapolated displacement pressure, the abscissa axis; and (3) pore geometry factor, by the curve shape. PROPERTIES OF BRINE AND BACTERIA The fluid used throughout this study was an aqueous solution of sodium chloride, 25,000 ppm by weight. The addition of salt to fresh water was intended to eliminate clay swelling in the sandstone cores and to retard leaching in the limestone cores. The viscosity of the brine solution at 100F, the temp-

Jan 1, 1966

By R. C. Davis, R. F. Berry, G. W. White, B. T. O’Connell, L. A. Stacha

P,, and v, may now be solved for by trial and error between Eqs. 2 and 3, by assuming time approaching zero and equating PI, to P12. Observe that the pressure under the slug, P12 of Eq. 3, does not remain equal to the P1, at the valve in Eq. 2 as the slug proceeds up the hole. To plot the profile of pressure vs depth, it is necessary to express P12 in terms of P11. To do so, the extent to which the original submergence length S changes with respect to time must be analyzed. At the end of time t after injection, the actual value of L in Fig. 3 is The fall-back v,(t) is now distributed in the gas and exerts pressure on the bottom. This, in addition to the weight of the gas, amounts to the difference in pressures from point to point. With the average pressure profile calculable, a true numerical value for the average pressure in the tubing may be obtained. The value of v, must be obtained experimentally, and since v,, S and D are known, the basic Eq. 1 may be solved for the amount produced. Knowing P, and the actual volume of tubing occupied by gas (vat), it is possible to calculate the volume of gas at standard conditions which is present in the

By J. S. Hegwer, P. M. Dunlap

An improved acid for the treatrrzent of sulfate-con-raining limestones and dolomites is described. The acid is designed to reduce he reprecipitation of dissolved calcium sulfate and the possibility of plugging permeable flow channels. In addition, this improved acid has a much lower reaction rate than that of regular acid; the advantages of a "retarded" acid are obtainable. Field testing of the acid has shown it suitable for use in sulfate-containing formations. Substantial improvements in productivity generally resulted. INTRODUCTION Acid treatments of limestones, dolomites and other formations bearing carbonate deposits are frequently unsuccessful when the calcareous formation contains sulfate, either as anhydrite (CaSO,) or gypsum (CaSO4. 2H2O). Preliminary dissolution in acid followed by redeposition of calcium sulfate appears to be a major factor contributing to poor well performance after acidizing. The precipitate is usually the gypsum form of calcium sulfate, but in higher temperature formations it may be anhydrite. The freshly precipitated crystals are nearly always very small and needle-like. They may occupy a gross volume many times that of the original anhydrite crystals and will obviously constitute an impediment to flow through newly enlarged flow channels. It is believed that the redeposition problem is most severe when anhydrite lines the fracture systems and large pores which supply the effective permeability of a formation. Microscopic inclusions of calcium sulfate also present large sulfate surface areas for dissolution in acid. In either case, great amounts of calcium sulfate may dissolve before the acid can be spent on formation carbonates. For regular spent acid (originally 15 per cent hydrochloric acid) the precipitate could be as much as 270 Ib gypsum/1000 gal acid. Two techniques have been applied by the industry for reducing sulfate plugging during acidizing. The method'.' commonly employed in the field is the attempted removal of a quantity of regular treating acid before it has reacted completely with the formation. This is practiced because the solubility of calcium sulfate is greater in a solution that is still acidic than in one which has been largely spent on the formation rock. The chance of precipitative plugging is therefore reduced if the withdrawal is successful. However, it is often impossible to get the acid out of the formation before precipitation occurs. A second possible method, which at first glance appears practical, involves addition to the acid of sequestering agents which form strong soluble complexes with calcium ions. These chemicals do increase the "solubility" of calcium sulfate in fresh acid, but to a lesser extent in spent acid. The sequestering agents have, therefore, proved unsatisfactory because the amount of sulfate eventually deposited from the spent acid may be greater than that from regular acid. Another logical approach to the problem of calcium sulfate reprecipitation is the prevention of the initial dissolution of calcium sulfate by the common ion effect. This may be accomplished by adding a soluble calcium salt to the fresh acid. The use of calcium salts in treating acids is not entirely new. An earlier suggested use of a soluble calcium salt in hydrochloric acid apparently failed to recognize the full extent to which the solubility of calcium sulfate could be suppressed. The present study extends this earlier work and adds certain improvements toward the development of a practical anti-anhydrite acid. LABORATORY DEVELOPMENT Calcium Sulfate Solubility Table 1 shows the results of a laboratory study performed to establish the effect of calcium chloride concentration on gypsum solubility. Because of the strong tendency of calcium sulfate to form supersaturated solutions, accurate solubilities are difficult to determine. These solubility data probably are reliable to within ± 15 per cent. The solubility of calcium sulfate in 15 weight per cent hydrochloric acid increases with increasing temperature. This trend is also followed in hydrochloric acid which cantains calcium chloride. This is contrary to the solubility behavior of calcium sulfate in water, wherein the solubility decreases with increasing temperature.

By W. G. Bearden, G. C. Howard, J. W. Spurlock

An investigation of the factors af-fec.ting the inter-zonal flow of fluids in the casing-wellbore annrc1rt.s of an oil or gas ivell is presented Laborntory tests revealed that failure of the caring-cement bond allows communi-cation of fluid at pressure differentials near 1.000 psi. Perforating with a 11011011: carrier gun does not shatter or crack the cement sufficiently to allow communication through these cracks and does not seriously affect the casing-cement bond failure pressure. The pressure to cause communication at the casing-cement bond can he increared to 3,000 psi by. limiting pipe deformation during cementing and subsequent operations. Communication can be prevented, however. up to pressures sufficient to cause the cement to fail or the casing to col-lapse, by attaching to the casing a seal ring of deformable ruhber molded between two steel flanges. When installed, one flange is attached rigidlv to the casing to serve as back-up for the gasket While tile other flange is movable. To date, these gasket seal rings have been field-tested in 27 Wells. Althotrgh communication has been noted in two wells, inter-zonal flow has been prevented in other wells in the same area by installing more seal rings between the communicaring zones. INTRODUCTION The inter-zonal flow of fluids behind casing cemented in oil wells during normal production or stimulation operations is a major operating problem. If it cccurs during production operations. it can seriously limit the pressure to which a reservoir may be depleted without harmful invasion of unwanted fluids. If it occurs during stimulation operations, it may preclude the possibility of stimulating the desired zone. This communication can be caused by failure of the bond between the casing and cement, failure of the bond between the cement and formation, or failure of the cement itself. This paper presents the results of an investigation of factors affecting com-munication at the casing-cement bond and methods which aid in preventing it. LABORATORY INVESTIGATIONS The laboratory tests performed in this study were divided into three phases. The first phase consisted of tests to determine if perforating would shatter and crack the cement sheath so severely that communication would occur through the cement itself. The second phase consisted of tests to determine more precisely the cause for the casing-cement bond failure. The third phase consisted of tests to evaluate the effectiveness of various de-formable seals attached to the casing in preventing communication. EFFECT OF PERFORATING, The effect of perforating was invest~gated in tests performed with apparatus shown on Fig. 1. Cement was poured into the annulur space hetween 5%-in. casing and 8%-in. casing and allowed to cure for 72 hours. The apparatus was then perforated with a hollow carrier gun containing eight duPont 26-A shaped charges. Bull plugs were placed over the perforations in the 8 5/8-in. casing and hydraulic pressure was applied inside the 5 1/2i-in. casing. The pressure to cause fluid communication past the cement sheath was observed and termed the "bond failure pressure". Several different cement compositions were included in these tests, the results of which are shown in Table 1. In Tests 1 through 4 of Table 1, bond failure occurred at- the inter- .

Jan 1, 1966

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

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

By 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 Wieland D. R., Hurst R. E., A. R. Hendrickson

Analysis of acidizing techniques, in correlation with reservoir data and a backlog of past treatments, has resulted in the development of a valuable engineering guide for planning acidizing treatments. Such treatments fall into three categories: (I) acid injection into the pores of the matrix; (2) acid injection into natural formation fractures at less than parting pressure; and (3) combination acidizing-fracturing treatments in which acid solutions (without propping agents) are injected at treating pressures sufficient to open and extend fractures through which the acid flows. Because the spending tirnze of acid during a specific well treatment does not change appreciably, maximunl penetration is attained when the first increment of injected acid is completely spent. Additional acid injection cannot be expected to further extend the benefits of the treatment. Depth of penetration will depend upon the reaction rate of the acid under treatment conditions, the injection rate of the acid into the matrix or fractures and the area-volume relationship existing in the flow channels. Based on Darcy's flow formula, extremely low injection rates must be used in order to keep bottom-hole injection pressures below formation fracturing pressure. As a result, only limited penetration of unspent acid will occur. Treatment records indicate that, in most acidizing treatments, formation parting pressures are exceeded, greatly extending acid penetration. Under these conditions, stimulation benefits are limited to the fracture area produced during the spending time of the first increment of acid injected into the formation. This area may be calculated from laboratory and well data to estimate depth of penetration. This, in turn, may be correlated with productivity data to assist The art of gas and oil well acidizing has been characterized by many changes in treating materials and techniques since its inception. These developments have been designed to provide greater production increases. prolong production declines and shorten payout time. Such improvements have been based primarily on data derived from laboratory research and field experience. As more of the variables influencing these treatments have been recognized and evaluated, acidizing has become less of an art and more of a science. Recent studies of fracturing treatments,' in light of individual well conditions and the results of thousands of fracturing treatments, made possible the formulation of an engineering guide that is now being used to select optimum treating techniques and to forecast probable results of such treatments. A similar analysis of the factors controlling acidizing treatments has been made and is the basis for this paper, The findings herein can be used as a guide in the selection of acidizing solutions and techniques, tailored to fit specific well conditions and to provide optimum stimulation per dollar cost. Acidizing treatments may be classified into three basic categories—(1) treatments in which the acid is injected uniformly into the pores and flow channels of the matrix, (2) treatments in which the acid enters natural fissures and fractures in the formation at less than fracturing pressures and (3) injection of acid into the formation at a pressure sufficient to open and extend fractures into the rock through which the acid penetrates (without the inclusion of a propping agent). TYPE 1 — MATRIX ACIDIZING This category consists of treatments in which acid solutions are injected into a homogeneous carbonate

By Brown, R. W., G. H. Neill, R. G. Loper

All facets of ball sealer behavior must be known and understood to design for their optimum use in well treatments. The down-hole factors including the inertial forces, drag forces and holding forces have been discussed. Equations have been developed and presented for these forces. The use of these equations in this process will contribute to the most efficient use of an economical and proven ball sealer process. The urocedure presented has been used to successfully analyze field results indicating problems involving the use of ball sealers. These problems have been isolated and corrected in future treatments by using the procedure presented to revise the perforating program, the injection rate, the design of ball sealer stages and/or the type and size of ball sealers. The procedure presented also indicates that dislodgement of ball sealers occurs primarily in ultra slim-hole completions. INTRODUCTION The perforation sealing process has been proven highly successful and economical since it was introduced to the oil and gas industry in early 1956.' Since it has proved so successful in multistage fracture treatments at a much lower cost than conventional packers, the ball sealer process has become a byword in well completions.' Not only has it changed the design of fracture treatments, it has changed the concepts of selectively perforating and has greatly aided the success of the single-point and limited-entry techniques." Yet with all the advantages, there are still facets of ball sealer behavior that are not widely recognized. These must be known and understood to better design and execute the optimum ball sealer treatment. The down-hole behavior of ball sealers is influenced by a number of factors including an "inertial force", a "drag force" and a "holding force". These are considered in two phases—the initial seating of the ball and the forces tending to unseat the ball after contact has been made (see Fig. 1). This paper presents a theoretical study of these factors and formulas for predicting their influence. The authors acknowledge the fact that the flow equations presented herein describe the behavior of Newtonian fluids and do not cover the so-called power-law fluids. The assumptions made, however, are reasonable and well within practical limitations of the conditions encountered during field application. For example, it is assumed that the perforations are round with a discharge coefficient of 0.82. While this may not be exactly true, it has been verified within practical field limitations.' The discharge coefficient would not vary significantly even if the perforation size or configuration varied pronouncedly. The approach of assuming uniform perforation size is not as strained as it might first appear. Normally, a perforator contains charges equally spaced with the same number in each vertical (or horizontal) plane. If the perforating gun is not centralized, the holes in one plane will be smaller than average, while the holes in opposing planes will tend to be larger than average. The net result is that the two effects tend to balance. FACTORS AFFECTING CONTACT WITH PERFORATlON Ball sealers are injected into well treating fluids for the purpose of contacting and sealing those perforations which are accepting the fluid flow. The efficiency of the sealers is primarily influenced by: (1) the velocity of the balls down the pipe, and (2) the fluid velocity through the perforations. To divert the sealer to the perforation, the inertial force of the ball must be overcome by the drag force created by the fluid velocity through the perforation. BALL VELOCITY The final, or stable, velocity of the ball is the sum of the actual fall velocity and the fluid velocity. This will be

By H. A. Wahl

This paper presents methods of designing hydraulic fracture treatments in formations saturated with slightly compressible liquids. Howard and Fast describe the fluid-loss control resulting from the viscosity and compressibility effects of the reservoir fluid in terms of a fracturing fluid coefficient. However, for horizontal fractures, this coefficient is applicable to infinitely thick reservoirs. These methods of evaluating fluid loss from horizontal fractures during treatment differ from previous methods in that: (I) formation thickness is incorporated in the velocity function describing the rate of fluid loss to the formation; and (2) for a given thickness the relative position of the fracture to the horizontal boundaries of the reservoir or any interrnediate permeability barrier is incorporated in the design procedure. This theory is applied to several actual and hypothetical field cases. The results of these calculations demonstrate that the eflect of formation thickness on fluid loss from horizontal fractures may result in fractures several times larger than predicted by other methods. INTRODUCTION Hydraulic fracturing has been the dominating well-stimulation technique during the last decade. From the time it was commercially available, fracturing has had a very favorable success ratio compared to other stimulation methods. However, numerous wells that were successfully stimulated remained economic failures. Only in the last few years has the oil industry recognized the importance of designing treatments for each well which will yield the maximum profit. This implies that for a given set of treatment conditions the size and general physical characteristics of the fracture can be determined. One of the outstanding contributions in the field of fracture design was presented by Howard and Fast.' In their paper, three flow mechanisms that control fluid loss from the fracture to the formation are defined: (1) viscosity effects of the fracturing fluid, (2) viscosity and compressibility effects of the reservoir fluid, and (3) wall-building effects of a treated fracturing fluid. The rate of fluid loss to the formation adjacent to a fracture is shown to be a function of time. For each flow mechanism the velocity function v(t) is expressed as a constant or fracturing fluid coefficient divided by the square root of time. This form of v(t) can be used in the Carter equation, as developed in the Howard and Fasc paper, to estimate the areal extent of fractures. This equation may be written as In the following sections the effect of formation thickness on fluid loss from horizontal fractures is evaluated when the controlling flow mechanism is the viscosity and compressibility of the reservoir fluids. The major assumptions inherent in the Carter equation are brought out in subsequent mathematical developments. THE EFFECT OF RESERVOIR FLUIDS ON FRACTURE EXTENSION When a fracture treatment is performed using a fracturing fluid identical to the reservoir fluid, the only flow mechanism controlling fluid loss is the viscosity and compressibility effects of the reservoir fluids. In developing a fracturing fluid coefficient for this mechanism, Howard and Fast considered a segment of the formation bounded by a unit area on the fracture face and extending to infinity in the direction of leak-off as shown in Fig. 1(a). The direction of leak-off z is normal to the fracture face. The equation governing the linear isothermal flow of fluid into this segment is where a = Eq. 2 is a form of the diffusivity equation. It is derived from the equation of continuity, Darcy's law and Eq. 3

By J. H. Henderson, T. E. Morrisson

An irvestigation has been made of the production of oil through horizontal fractures of high capacity and large radius placed at the base of producing formations. The specific aitns of the study were to (I) determine by model studies the reservoir performance when gravity drainage is the only producing mechanism, and (2) develop a method of predicting such perfortnance from commonly measured reservoir proper-ties. Although the theoretical approach tailed to develop a rigorous analytic sorution to the transient problem, it (lid lead to a practical method of predicting performance. Predictions (Ire in good agreement with rates and recoveries observed in scaled-model experinzents. Capillary effects and the influence of relative permzenbil-ity are taken into account. The experimental program was mainly concerned with the influenre of fracture radius, fracture capacily and capillary hold-up on the total recovery and the rate of recovery. Much of the experimental work was done on a triangular model which represented one-eighth of a square well pattern, so constritcted that horizontal fractures of varying radii and flow capacity could be studied. An electrolytic model was used to relate the initial production rates to fracture radius, formation thickness and the drainage radius of the well pattern. The conclusions reached as a resrtlt of the work are as follows. I. The production of oil under gravity force alone into a well-propped horizontal fracture, whose radius is of the order of hundreds of feet, is a recovery method which combines high recoveries with very acceptable production rates. 2. The mechanism is effective in "dead oil" reservoirs as well as in reservoirs which possess the energy of solution gas. 3. A method of predicting performance has been developed which gives good agreement with experimental behavior in cases of practicnl interest. Even for very thick beds and large fracture radii, where actual conditions do not conform to the assumptions made, the method gives order of magnitude predictions which would not otherwise be possible. INTRODUCTION Gravity drainage of oil from a formation which has continuous ver- tical permeability over a reasonable depth interval is a surprisingly effective method of recovery. However, for a conventional well the production rates which result from gravity drainage are so low that the economics are seldom attractive. The same formation may produce rapidly and profitably by gravity alone through a high capacity, extensive (i.e., > 100-ft radius) horizontal fracture placed near the bottom of the pay section. Fig. 1 shows a schematic diagram of the plan. Energy in the reservoir gas is not a necessary condition for this recovery process. With a sufficiently large fracture, attractive production rates will result from gravity flow alone, and the producing GOR need only be the solution GOR. It is, of course, necessary to keep the fluid level in

By J. L. Huitt, J. L. Pekarek, D. K. Lowe

In a theoretical study of hydraulic jetting, the velocity of the abrasive material relative to the velocity of the fluid in the jet stream is analyzed as the jet stream moves through the convergent and straight sections of the nozzle and the region between the nozzle exit and target. The results revealed that the abrasive material exits from the jet nozzle at a lower velocity than the fluid. The exit pmticle velocity can be increased by increasing either the density of the fluid or the length of the nozzle, and/ or decreasing either the particle density or particle diameter. In the divergent jet stream, there exists a point after which the particle velocity exceeds that of the fluid. The relative velocities were considered in the derivation of an equation to predict cutting rate of a circumferential notch and maximum notch depth. Data of a general nature and data which substantiate the theoretical results were obtained experimentally. INTRODUCTION The use of a fluid containing an abrasive material has been an established technique for cleaning and cutting for many years. In the petroleum industry, the early effort to use this technique1 to perforate and/or to overcome wellbore damage met with only limited acceptance because of the short life of the jet nozzle. With the introduction of improved perforating techniques, and later, hydraulic fracturing, the use of hydraulic jetting as a well completion technique became even less appreciated. It was only in recent years that interest in hydraulic jetting was revived. Once this interest was revived, the results of surface tests stimulated the interest of the industry even more than the state of the technology probably warranted because many of the tests were not appropriate for down-hole conditions. However, because of the stimulated interest, the development of the jet nozzle progressed very rapidly to the point where the nozzle life was no longer a prob- lem. With this accomplished, the use of hydraulic jetting in well completion became an accepted practice in a short time. The purpose of this paper is to present a theoretical analysis of the hydraulic jet stream as it passes through the nozzle and travels to its target. With a better understanding of the jet stream and the effects of various parameters, the performance of the process can be predicted more accurately. Equations are presented for cutting rate as applied to circumferential wellbore notching that relate the jet stream make-up, notch configuration and formation material. Also, experimental data are presented on some factors pertinent to hydraulic notching that are not theoretically analyzed. RELATED STUDIES Most of the studiesl-5 reported in the recent literature have pertained to the more practical aspects of hydraulic jetting; i.e., the effects of certain parameters as interpreted from experimental results, and the application of hydraulic jetting in well completion. In reviewing the effects of various parameters, it is interesting to note the reported depths of penetration obtained under various imposed conditions. In general, the depths vary from a few inches to several feet; however, a depth of penetration of less than 6 in., as reported by Thompson,4 seems more realistic for the usual field practice of hydraulic jetting with sand in water for a period of 20 to 30 minutes. In addition to the practical aspects, the study of Brown and Loper5 included a theoretical approach to hydraulic jetting. Their study resulted in the development of a theoretical expression for the maximum depth of penetration if jetting were continued for an infinite time. An analysis of the equations presented reveals that the initial cutting rate is infinite. The equation expressing centerline velocity is that of Forstall and Gaylord,6 which is applicable for a jet stream exiting in a large stationary medium. Since practically all of the fluid pumped into a perforation (or cut) must flow back through the perforation prior to re-entering the wellbore, a description of the medium as finite and non-stationary seems more reasonable. Thus, in this

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

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 T. C. Poulter, B. M. Caldwell, D. E. Brink

The flash X-ray has been used more than a decade to study the configuration of the jet from a shaped charge. The high-speed, rotating-mirror smear camera has provided time-distance graphs of detonations and shock fronts in transparent materials, and the high-speed framing camera has given pictorial representations of the progress of explosive phenomena. Adequate means for measuring the paths and velocities of all parts of the shaped-charge cavity Iiner during the collapse phase have not existed heretofore. A technique enabling the single-lens framing camera to make stereoscopic photographs of the cavity liner while it is collapsing has been developed. Analysis of this photographic record gives the directions and velocities of various parts of the Iiner surface, permitting direct quantitative measurement where this has previously been impossible. Significant improvements in shnped-clrarge design are expected to result. INTRODUCTION The familiar oilwell jet perforating charge, generally referred to as a "shaped charge", is related to the Mun-roe charge first described by C. E. Munroe in 1888 and later by Neumann in 1910. It differs in construction from the Munroe charge in that its cavity is lined with some inert material (usually metal)', and it differs in performance by projecting a fast jet of dense liner material against the target instead of a stream of expanding detonation products. Fig. 1 represents a typical shaped charge in axial cross section. The charge is circularly symmetrical about its longitudinal axis. The cavity at the right is lined with metal, usually copper; and, with the exception of air, the liner is otherwise completely empty. High-explosive material is intimately in contact with the convex surface of the liner and extends in the other direction to the booster. This is a pellet of high explosive having somewhat different characteristics and serves to couple the detonator to the main charge. When the detonator is fired, a small hemispherical reaction zone of very high temperature and pressure occurs at its output end. The booster responds, adding its own latent energy to the reaction, thereby initiating the main charge. The reaction progresses throughout the main charge at a velocity of 7 to 8 mm/microsecond until all the explosive material has reacted. With a reaction pressure of about 160,000 atm, the surrounding material is moved rapidly away from the high-pressure zone. As the detonation proceeds along the cavity liner, moving from apex to edge, the liner is collapsed toward the axis of the cavity, where liner material from all parts of the concave surface converges. A zone of extraordinarily high pressure results. with only one escape path—outward along the axis. The particles from the liner, thus, become a jet of very high velocity, with a density characteristic of the material from which the liner was made. The kinetic energy of the jet, concentrated in a small area of the target, has great penetrating power. It overcomes the strength of any target, and the target material crumbles or deforms plastically to escape from the path of the jet. An unlined cavity charge, on the other hand, projects a turbulent body of hot gaseous material against a target and erodes the exposed surface. With equal quantities of explosive, the penetration effectiveness of an unlined charge is a small fraction of that of a lined charge. Given an explosive charge with a detonator at one end and an axially symmetrical metal-lined cavity at the other, one has a perforator of considerable potency. To maximize the penetrating and other desirable effects,

By J. van Lookeren

In the case of a reservoir where the oil underlies a large gas cap and overlies bottom water, production can be inzproved considerably if wells are perforated below the water-oil contact rather than opposite the oil-bearing section. The theoretical arguments for this improvement, which in some cases may be tenfold, are supported by experimental evidence. INTRODUCTION Reservoirs frequently consist of one sand body containing an oil layer underlying a thick, gas-bearing section and overlying bottom water (Fig. 1). Wells in such reservoirs are customarily opened up opposite the oil-bearing interval or, rather, opposite the lower part of this interval (Fig. 2). More often than not, it is inefficient or prohibited to blow off large quantities of gas with the oil. This im-psses a severe limitation on the drawdown, which results in low oil-production rates. Because of the small difference in specific gravity between oil and water, the latter will easily cone upwards, so that the production of water together with the oil is not prevented by this completion method. In such a situation wells should be opened deeper, in the originally water-bearing zone. The gas can be shut off by not perforating over the upper part of the interval (Fig. 3) or by a tubing hung close to the bottom of the interval (Fig. 4). Initially much water will be produced, but when steady-state conditions have been established (as depicted in Figs. 3 and 4) larger drawdowns are permissible than in a conventionally completed well, and correspondingly higher oil-production rates will result. Water cuts are only slightly higher because of the small effect of gravity on the water-oil flow distribution. CALCULATIONS The ratio of improvement has been estimated by means of approximate calculations. The completion depicted in Fig. 4 is the most convenient basis for these calculations. By postulating that vertical flow gradients can be disregarded, the Appendix shows that the well-known Dupuit formula can be obtained except for some complication due to the presence of two fluids rather than one. The Appendix makes the reasonable assumption that, for suitable boundary conditions, the oil layer occupies a constant fraction of the total fluid height throughout the reservoir. The oil- and water-production rates can thus be related to the liquid levels in the borehole as follows:

Jan 1, 1966

By K. H. Coats

In tnany oil and gas fields, significant investments are required in piping networks for collection, injection and disposal of various fluids. This paper treats the problem of optimally locating the trunkline(s) for a given systet to nlinimize both the total length of piping connecting the individual wells or lease batteries to the trunkline and the length of the latter. Single and multiple trunklines, varying line sizes, and grouping of wells are considered. Some attention is paid to the case of a curved trunkline. A field case is treated and conclusions are that the problem of locating a single trunkline is relatively straightforward while the location of multiple trunklines requires further research. INTRODUCTION la many oil and gas fields, significant investments are required in piping networks for collection, injection and disposal of various fluids. A single field may concurrently require piping systems for gathering of oil, gathering of gas, injection of water or gas and disposal of salt water. Two problems arise in the consideration of any one of these networks. The first involves the sizing of the various segments of the network to minimize piping and compression costs. The second problem, treated herein, is that of optimally locating the trunkline(s) for a given piping system to minimize the total length of pipe connecting the individual wells or lease batteries to the trunkline. This minimization problem may arise in field automation if it involves replacement of individual lease batteries by a central field battery which will require trunklines and connecting pipelines to individual wells. The problem also arises in connecting a gas pipeline through a subsidiary trunkline to wells in a gas storage field. The first analysis given below treats the case of a single trunkline and of separate connecting lines for all field units.* The analysis is then extended to consider multiple trunklines and common connecting lines for groups of field units. Varying sizes, and therefore costs, of the connecting lines are taken into account. A simplified development is given to illustrate the manner in which line sizing calculations may be integrated with the line placement problem. This study is primarily directed toward the placement of straight trunklines, although for certain field shapes or well distributions curved trunk-lines might be preferable. Some mathematical consideration is therefore given to the placement of a curved trunk-line in a field. The determination of optimal locations for two or more trunklines involves some difficulties which are discussed below. The need for further research in this area is emphasized. The problem of a single trunkline is solved in relation to an actual oil field. The problem of optimally locating two trunklines in the same field is also treated and the results are compared to the single-trunkline case. All calculations reported were performed on an IBM 7072 digital compker. The importance of right-of-way, roads and terrain in locating trunklines and connecting lines is realized. The present analysis is advanced on the premise that in planning piping networks one should start with the optimal conditions whenever possible and proceed therefrom.