Search Documents

By B. H. Caudle, A. R. Khan

The oil production performance of thin-oil-column type reservoirs producing entirely by natural water drive are discussed in this paper. The dimensionless production histories were obtained in a series of scaled porous-medium laboratory-model tests under varying effects of rate of production, mobility ratio across moving boundary, ratio of oil- and water-sand thickness and mobility ratio across fixed boundary including the effects of different residual oil saturations in the water-invaded regions. This scaled-model study of thin oil columns indicated that greater recovery of oil results for lower values of the following: rate of production, mobility ratio across the moving boundary, water-sand thickness relative to the oil-sand thickness, and mobility ratio across the fixed boundary or residual oil saturation in the water-invaded region of oil sand The results of the study are given as dimensionless production history plots which may be used to predict the performance of the field they represent. INTRODUCTION Oil fields producing under natural water drive can be classified into two major groups: (1) an "edge-water drive" reservoir group in which the water is encroaching from the reservoir periphery, and (2) the group of reservoirs in which the oil-sand is entirely underlain by an aquifer. Reservoirs in the second group are further subclassified as "bottom-water drive" reservoirs and "thin-oil-column" reservoirs. The latter reservoirs are characterized by their relatively equal oil- and water-sand thicknesses and by having both the oil and water zone relatively thin. When oil is produced from a water-drive reservoir by a partially penetrating well, water breaks through sooner or later before all of the displace-able oil is recovered. Many analyses of this general situation have been made with varying degrees of realism in their approaches, both theoretical and experimental5.6 in nature. Pertaining to the scaled laboratory model studies of natural water drives using porous media, detailed investigations have been made by Henley, et al.7 on bottom-water drive and by Caudle8 on edge-water drive reservoirs. No laboratory model studies, however, have been made to investigate the performance characteristics of thin-oil-column reservoirs. This paper presents the results of a study on scaled models of thin-oil-column reservoirs, operated under varying conditions corresponding to the actual field practices. SCALING PARAMETERS The object of scaling is to produce a one-to-one relationship between physical and dynamic parameters of the model -and prototype. The most comprehensive treatment of scaling techniques for oil reservoirs is given by Geertsma, et a1.,9 which covers two methods for arriving at dimensionless groups for scaling. Table 1 lists the scaling parameters suitable for creating kinematic and dynamic similarities between the model and the prototype in the variables considered for the thin-oil-column reservoirs. In order to achieve similarity it is essential that the dimensionless groups listed in Table 1 have equal values in the model and the protorype simulated. Since it has been established that the role of interfacial forces is neglible in gross fluid movement in a reservoir no attempt was made to scale capillary forces- Although there are reports describing Scaling criteria for these forces,10,ll consideration of the necessary but restrictive conditions for

Jan 1, 1970

By T. K. Perkins, W. W. Krech

Jan 1, 1969

By A. G. Weber, R. L. Ledbetter, A. L. Hicks

A procedure which makes use of both analog and digital computers has been developed for predicting the pressure-production behavior. of water-drive reservoirs. The electric analyzer is used for matching the past pressure and production history of the reservoir and inferring the physical characteristies of the aquifer. When the history has been matched, a .set of base curves called influence functionce" are derived and the future behavior of the reservoir is calculated from these curves. INTRODUCTION Electric analog methods are particularly well-suited for complex water-drive problems. Like all other methods of analyzing water-drive reservoirs, however, analog methods are subject to certain limitations. One of the more obvious and serious of these results from inaccurate or insufficient knowledge about physical characteristics of the reservoir and aquifer. However, it is not insurmountable because the unknown characteristics can usually be inferred from the reservoir's previous behavior. Simons and Spain have shown that even with scant reservoir data and limited production history, electric analyzer studies can provide useful and reliable results. Aside from this limitation on analyzer predictions, the main disadvantage has been lack of a technique for keeping predictions abreast of latest production forecasts. Pressure behavior of a reservoir depends on the rate at which fluid is withdrawn; it is therefore necessary to specify the future withdrawal schedule before the future pressure behavior can be calculated. Conventionally, the electric analyzer predicts pressure behavior for a specific production schedule. If this schedule is abandoned in favor of some new one, as is frequently the case, either the study must be rerun or results adjusted in some manner. The usual procedure is to predict pressures for several different forecasts covering a wide range of production rates. Pressures for specific rates other than those studied can then be approximated by interpolating between the analyzer predicted curves. However, reliability of this type of prediction is doubtful. This paper describes a technique for predicting future behavior of water-drive reservoirs that retains all advantages of electric analyzer analysis and eliminates the main disadvantage. That is, the analyzer is used conventionally to infer unknown reservoir parameters from the pressure history, but not to predict the future behavior directly. Instead, a derived set of base curves, called influence l'unctions, can be used to calculate future behavior for any production schedule.

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 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 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 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 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 W. T. Bell, M. P. Lebourg

The perforating of multiple tubingless completions, in which two or more strings of 27/8-in. OD casing are installed in the same borehole, presents two basic problems. First, good completion practices require an eficient perforator that leaves no debris to interfere with subsequent completion operations. Second, and of more complicated nature, a system is needed for controlling the direction of fire so that adjacent strings are not damaged. Development of of a 2-in. OD, steel retrievable, shnped-charge gun has solved the first problem. The second was resolved through development of three different devices to provide directional perforating for the two types of completion methods being employed today. These include a mechanical orienting device and two self-orienting radiation devices. Conceivably, the latter methods could be adapted to the directional perforating (of .upper zones) of conventional dual and triple completions without removing production tubing and packers. INTRODUCTION The recent trend to multiple tubingless completions,' wherein two or more strings of 27/8-in. OD casing are installed in the same borehole, has presented two basic perforating problems. The first problem to be overcome was the design of a gun that would result in efficient perforations without causing debris to bridge in these small-diameter casings. The second problem was, of course, to devise a system for controlling the direction of fire such that adjacent strings would not be damaged. The problem of gun debris in slim-casing completions has been well defined in the numerous 27/8-in. single tubingless completions which have been effected during recent years with the expendable-type perforators. Failure of expendable-gun debris to settle out properly, even with good gun break-up, has resulted in bridging inside these small-diameter casings. This, in turn, has interfered with completion and re-conditioning operations. The past two years have clearly shown, in fact, that debris problems are not confined to this small range of casing with most designs of expendable-type guns. A high incidence of bridging has even been reported in 4%-in. OD casing. The perforating of multiple tubingless completions has entailed, therefore, the development of a completely debris-free gun of adequate performance, as well as the necessary devices to provide positive directional perforating in the two completion methods presently being employed. THE GUN —DESIGN CONSIDERATIONS The design of a shaped-charge gun for slim-casing perforating has presented a maze of interrelated prob-lems. Since the primary specification has been leveled at efficient charge performance with complete elimination of gun debris, it has been necessary to develop a retrievable steel-carrier gun of tubular construction to house the charges. The carrier must be retrievable through 2.29-in. API pump-seating nipples which are usually installed in 2 7/8-in. strings. To provide for recovery of the gun from the well, swelling of the carrier on firing must be carefully controlled. Control of swelling automatically means control of explosive weight in the shaped charge which, in turn, influences charge performance. Carrier swelling also depends upon the confinement pressure at the level being perforated, which imposes a hydrostatic pressure-depth limitation. Burr interference on the wall of the carrier, resulting from passage of the jet through the carrier wall, must be eliminated by suitable design. Clearly then, a satisfactory shaped-charge gun design for slim-casing operations represents a fine balance between carrier swelling, hydrostatic pressure, explosive weight and charge performance. Additionally, to lend itself to directional perforating, the gun design must provide a positive locking mechanism for individual charges within the gun carrier. This fixes charge orientation within the carrier and prevents

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

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

By G. R. Countryman, G. H. Thomas, I. Fatt

Procedures are described for scaling up turbulent friction pressure drops of inelastic fluids and laminar friction pressure drops of inelastic/elastic fluids in circular ducts. The laminar scale-up calculation method does not require the use of arbitrarily defined flow models, and can be performed directly from laminar capillary viscometry data, or by means of a numerical solution of the integral equation for Poiseuille flow, inputing true shear rate-shearing stress data from other simple shearing experiments. e.g. Couette and plane Poiseuille flows. The turbulent scale up procedure, based on the Dodge-Metzner correlation for inelastic fluids, requires evaluating the characteristic rheological parameters at the existing wall stress under laminar flow conditions. Since the wall stress is the quantity sought in a pipe flow scale-up problem, a trial-and-err or solution is indicated, and therefore a calculation method has been developed and programmed for a medium-size magnetic-drum-memory digital computer. Laminar flow data from other somces may also be adapted to this turbulent scale-upmetbod for inelastic fluids. Comparisons between predicted and experimentally observed friction losses for generalized Reynolds numbers ranging from 3 — 32.000 indicate an error band width of less than 10 per cent in circular tubes ranging in size from 0.5 to 1.89 in. in diameter. INTRODUCTION Considerable effort has been devoted to the problem of establishing an accurate procedure for predicting the laminar and turbulent flow behavior of non-Newtonian materials in pipes from data obtained with small scale laboratory or pilot plant equipment, as evidenced by the variety of reviews on the subject.1-3 Although there does not appear to be a universally accepted scale-up method for laminar flow, the various methods currently employed can be divided into two principal categories depending upon whether a functional relationship between shear rate and shearing stress is assumed, or whether a generalized correlation is adopted. The methods of Hedstrom4 and of Weltmann are typical of the first category, while the generalized correlation of Metzner and Reed6 illustrates the second category. Under appropriate conditions either approach yields a reasonably satisfactory description of laminar flow behavior. In contrast, the prediction of non-Newtonian transitional and turbulent friction losses is much more complex. For example, certain types of materials have been observed to sustain laminar flow at much higher flow rates than Newtonian liquids. When turbulence does develop, the level of agitation is damped to such an extent that the friction factor is lower than that of the original low viscosity Newtonian medium. Extremely broad transition regions and "diameter effects" have also been noted in the turbulent flow of certain chemically complex materials.' Aside from the latter complications, which are attributed in part to viscoelasticity, the existence of a variety of scale-up procedures for predicting transitional and turbulent friction losses is due principally to the chioce of viscosity parameter used in the definition of the Reynolds number. In contrast, Dodge and Metzner have demonstrated that the friction factorReynolds number relationship for inelastic fluids can be satisfactorily generalized provided the fluid property parameters are evaluated at the existing wall shearing stress.* However, the quantity usually sought in a pipe flow scale-up problem is the wall stress, and hence a trial-and-error procedure is required when dealing with non-power law materials. In this report we discuss two calculation methods: (1) for determining the wall stress, and hence the friction pressure drop, for inelastic fluids under conditions of turbulent flow by means of an iterative procedure, and (2) for predicting laminar friction pressure drops without recourse to arbitrarily defined flow models.

By O. D. Gaither

This paper is an analytical study of the flow of fluids through small vertical conduits. Small conduits are defined as 11/4-in. nominal diameter tubing size and smaller, and approximately twice this area for annular conduits (i.e., 1- X 21/2-in. annulus and smaller). Experimental data are presented for the 1-X2-in. and 11/4- X 2%-in. annuli, and the I-in. and 11/4-in. tubing, since these represent the small conduit sizes and configurations generally encountered in oilfield applications. Data have been gathered for these conduits for single-phase water, single-phase gas and two-phase water-gas mixtures, with particular emphasis on high gas-liquid ratios. Water rates in excess of 2,000 BID and gas rates in excess of 2.5 MMcf/D, and two-phase flow ratios in between these two, represent the scope of the data gathered. Existing equations have been applied to predict flowing pressures and compared with experimental data. New correlations have been developed. INTRODUCTION The increased economic pressure on the domestic oil industry in the United States has constantly required the use of new techniques and equipment designed to reduce the cost of finding and producing oil and gas. Since tangible items are most readily apparent in economic analysis, the advent of lower-cost well completions was inevitable. One of the methods used to reduce costs which has received widespread attention is the slim-hole completion technique where tubing is used as the well casing and in which small conduits are used for tubing if necessary. Small conduits, defined by Kirkpatrick1 as "11/4-in. diameter nominal tubing and smaller for tubing flow and less than twice the 11/4-in. diameter nominal tubing internal flow area for annulus flow", have also found widespread usage as siphon strings for de-watering gas wells and as "kill" strings in deep high-pressure oil and gas wells. The growing use of small-diameter tubing has resulted in an increased need for development of improved methods to measure or predict flowing bottom-hole pressures since the physical dimensions generally preclude the use of subsurface-recording pressure gauges. Even in the cases where small bombs are available, the relatively high velocities encountered at nominal flow rates make it necessary to use excessive weight bars or special hold-down devices. Attempts to use recognized correlations to accurately predict flowing or gas-lift performance in wells equipped with small conduits have been generally unsuccessful. Insufficient field data were available to allow the development of a correlation on this basis, and an experimental approach was applied in an attempt to obtain a workable relation. The experimental approach used to obtain the data presented in this paper was actually a compromise between a field installation and a laboratory study. A test well 1,000 ft in length was used to obtain flow data on single-phase liquid, single-phase gas and two-phase water-gas flowing mixtures. Liquid rates up to 2,200 B/D and gas rates up to 3 MMcf/D were used in the single-phase flow studies. Two-phase flow rates from 100 to 600 B/D with gas-liquid ratios from 500 to 8,000 cu ft/bbl were recorded. Experimental data were obtained for single- and two-phase flow through 1-in and 11/4-in. nominal tubing, and through the annuli between 1- and 2-in. and 11/4- and 2%-in. nominal tubing strings. Experimental results for the two-phase flow are compared to the Poettmann-Carpenter correlation2 which is widely used as a comparative standard for development of multiphase flow predictions in flowing and gas-lift wells. Correlations developed by Tek,3 Baxendell and Thomas" were also investigated. The experimental data recorded herein fell in between the two flow regimes as defined by Ros," and this correlation also failed to yield satisfactory results. The fact that existing correlations failed to confirm the experimental data led to the need for development of a new correlation. Although a two-phase flow study was the primary objective of this investigation, data were also recorded for single-phase flow of water and gas, and constants were developed relating to pipe roughness and equivalent diameters for annular flow. These single-phase studies assisted materially in the development of certain of the two-phase flow results. Considerable previous work has been published which presented relationship of surface measurements to bottom-hole condition. The works of Buthod and Whiteley,6 Jones,' Poettmannb and the Texas Railroad Commission" are classic examples of the successful use of mathematical relationships which allow acceptable predictions of subsurface pressures, when gas is the flowing fluid. Darcy and others have derived relationships which may be used with minor modifications to predict subsurface flowing conditions in injection and water-supply wells. As previously stated, the application of the single-phase flow relationships