Search Documents

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.

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 P. J. Jacquard, C. Jain

Interpreting pressure measurements is one of the fundamental problems in the operation of oil reservoirs. The principal methods and means of study zuhich are currently used are reviewed. Emphasis is placed on the necessity of improving the quality of the information to turn it to the best possible account. One of the important aspects of the problem is the development of a system of automatic interpretation of the pressure measurements. The proposed method aims at this objective: starting from the relation which exists between a local perturbation of permeability and its effect on the velocity potentials in the wells, we determine by successive approximation the modifications to be made to the permeability map by fitting the least squares between the measured pressures and the pressures indicated by the mathematical model. The method has already been successfully employed in the case of heterogeneous reservoirs of radial- circular symmetry. Its application is extended here to the two-dimensional case, and is illustrated by two examples, for which the results are presented. The conditions required for the method to become fully operational are likewise analyzed. I — INTRODUCTION One of the basic problems of production and recovery from oil fields is interpreting pressure measurements made in wells in the single-fluid-phase regime. Simplifying assumptions, for example supposing the reservoir to be homogeneous, make it possible to use more or less simple methods which permit the estimation, from these pressure measurements, of a mean value of the permeability around the wells. Considering the actual heterogeneity of the majority of reservoirs, and the phenomena of interference between wells, these methods are only indicative. Much has been published on the studies of different cases of comparatively simple heterogeneity: some have discussed the radial heterogeneity around a single wel1,l or vertical heterogeneity.2 The results of these works have shown the influence of variations of permeability on the pressure, but unfortunately they are not very useful in practice because they all are concerned with the phenomena around a single well, without interference between wells, and they correspond to relatively simple and theoretical schemes. These studies of reservoirs or parts of reservoirs on a larger scale are made by the electric analyzer or digital model. We define a permeability map and a porosity map, then reproduce the production history of the wells on the wells of the model, and note the corresponding pressures. These pressures are then compared with actuality, and the discrepancies observed are analyzed. From this analysis we deduce the modifications to be made to the maps of permeability and porosity, and repeat until the coincidence between the pressures indicated by the model and the actually observed pressures is considered adequate. We have adjusted the model. Obviously, in this method of investigation, everything depends on the formula which permits the modification of the permeability to approach the adjustment. Considering the complicated nature of the problem, some oversimplified formulas are poorly adapted to reaching a good adjustment. For example, we have shown that the elementary operation of increasing the permeability not far from a producing well to compensate for too low pressures in the well can, in certain cases, express itself by an effect exactly opposite. 3 The project described here is for a system where we begin with the pressure measurements and the several known parameters of the reservoir, and develop the most probable permeability map, taking account of the measurements available and a subdivision of the reservoir fixed a priori. Ideally, this system is conceived as one in which no manual intervention would be necessary. The model simulating the reservoir is only one part of the whole system, and it is neither the longest nor the most complex part. The need to integrate the model into an aggregate of much more extensive processes has led to the adoption of the digital solution, even though a solution combining the digital and the analog is perfectly conceivable. The simulating model is in that case no more

Jan 1, 1966

By N. Mungan

Formation damage, i.e.. reduclion in permeability, has been generally attribuled to clay minerals which expand or disperse upon contact with water that is less saline than the connate water. Luboratory, studies show that penneahility reduction can also occur in formalions containing only nonexpandable clays such as illite or kaolinite, and can be caused also by changes in pH. Furthermore, pH changes can damage even formations that are essentially free of clays. It is suggested that permeability reduction is due to the small passages being blocked by particles, which may be dispersed clays, cemenlion material or other fine parricles. These particles are dislodged by dispersion of clays due to changes in salinity or by dissolution of calcareous cement by acids, or of silicaceous cement by alkaline solutions. In working with reservoir cores, it was found that extracted cores damaged more easily and extensively than nonextracted cores. The extent of damage depended also on tenlperatltre. INTRODUCTION Permeability is an important property of porous media and has been the subject of many studies by engineers and geologists. Many of these studies are conccrned with formation damage, i.e., reduction in permeability, resulting from exposure of oil-producing formations to water substantially less saline than the connate water. This effect causes understandable concern since during drilling, completion and production phases formations are often exposed to fresh water. The damage resulting from contact with relatively fresh water has been attributed to expansion and dispersion oF clay minerals. During laboratory investigation of the use of NaOH as a wettability reversal agent to increase oil recovery from oil-wet reservoirs, several cores used in the displacement studies suffered loss in permeability. Despite the traditional usage of NaOH for conditioning aqueous mud systems, the role of the caustic filtrate in wellbore damage seems to have been overlooked. Browning2 as recently reported on the effects of NaOH in dispersing clay minerals but he was concerned only with complications that may arise in drilling massive shale beds. The following study was made to examine the role of pH and salinity changes in core damage. Where cores from reservoirs were used, tests were performed with extracted and nonextracted cores both at room and reser- voir temperatures, since it was felt that the test environment and core condition may affect the results. Because of its limited coverage and exploratory nature, this study is not intended to provide answers to field formation damage problems. It is hoped that it will encourage research into new aspects of the permeability reduction problems, particularly those allied to new recovery and production processes. PROCEDURE In all permeability tests, fluids were pumped through the cores at a constant volumetric rate. Only deaerated fluids and reagent grade chemicals were used. The fluids were passed through two ultrafine filters before injection to remove any entrained particles. The cores, with the exception of the unconsolidated cores, were mounted in Hassler holders. Water was used to transmit pressure to the sleeve. The inlet endpiece had two entry ports which permitted scavenging one fluid with another to avoid any mixing in the small holdup volume. The cores were flushed with CO2 gas, evacuated for 5 to 6 hours and saturated with the first liquid at a pressure of 1,000 psi for 24 hours to eliminate any free gas from the cores. Pressure differences up to 20 psi were measured by transducers, calibrated in inches of water and continuously recorded. For greater pressure drops, gauges were used. All reservoir cores were cleaned with a light refined mineral oil, then with heptane, and finally dried with CO2. Compatibility tests showed that no precipitates formed when mineral oil and the crudes were mixed. Some cores were extracted in Dean-Stark-type solvent extractors using xylene and trichloroethane and dried in a vacuum oven at 450F. Each test consisted of a sequence of water, test solution, and again water flow. RESULTS AND DISCUSSION STUDIES IN BEREA CORES Salinity Contrast Berea cores 2-in. in diameter and 12-in. long were cut from sandstone quarried in Cleveland, Ohio. The clay minerals were identified by X-ray diffraction to be chlorite, kaolinite, illite and incerlayered illite. Flow of fresh water or 30,000 ppm brinc does not cause any permeability reduction (Fig. 1). However, after injection of brine the core is readily damaged by fresh water. Damage starts almost instantly as the fresh water injection is begun, and at a cumulative injection of 1.2 PV fresh water, the permeability has dropped from 190 to 0.9 md. Upon continued injection, the effluent contains clay minerals dislodged from the core. The final core per-

Jan 1, 1966

By S. G. Gibbs

A new method for predicting the behavior of sucker-rod pumping systems is presented. The pumping system is described by a flexible mathematical model which is solved by means of partial diflerence equations with the aid of computers. Polished rod and intermediate-depth dynamometer cards can be calculated for various bottom-hole pump conditions. The technique permits simulation of a wide variety of operating conditions, both normal and abnormal. The data generated with the new technique are useful in refining the criteria for design and operation of sucker-rod systms. INTRODUCTION Sucker-rod pumping systems are used in approximately 90 per cent of artificially lifted wells. In view of this wide application, it behooves the industry to have a fundamental understanding of the sucker-rod pumping process. Oddly enough, our understanding has been rather superficial. This is evidenced by the semi-empirical formulas which have been used as the basis for design and operation of sucker-rod installations. Thoueh- we have realized the limitations of our methods for many years, it has not been computationally feasible to use more refined techniques. With the advent and widespread use of digital computers, it is now possible to handle the mathematical problems associated with sucker-rod pumping. This paper summarizes a computer-oriented method which can provide greater insight into the sucker-rod pumping process. It is hoped that this technique, and techniques which may evolve from it, will prove to be the tool needed by industry to obtain the most efficient use of rod pumping equipment. THE MATHEMATICAL MODEL Prediction of sucker-rod system behavior involves the solution of a boundary value problem. Such a problem includes a differential equation and a set of boundary conditions. For the sucker-rod problem, the wave equation is used, together with boundary conditions which describe the initial stress and velocity of the sucker rods, the motion of the polished rod and the operation of the down-hole pump. Of these items, the wave equation, the polished rod motion condition and the down-hole pump conditions are of primary importance. Discussion of the mathematical model centers about these factors. ROD STRING SIMULATION WITH THE WAVE EQUATION The one-dimensional wave equation with viscous damp- is used in the sucker-rod boundary value problem to simulate the behavior of the rod string. This equation describes the longitudinal vibrations in a long slender rod and, hence, is ideal for the sucker-rod application. Its use incorporates into the mathematical model the phenomenon of force wave reflection, which is an important characteristic of real systems. The viscous damping effect postulated in Eq. 1 yields good solutions, even though nonviscous effects such as coulcomb friction and hysteresis loss in the rod material are present. Fortunately, the nonviscous effects are relatively small, so the viscous damping approximation used in the wave equation is adequate. The coefficient v is a dimensionless damping factor which is found in field measurements to vary over fairly narrow limits. For mathematical convenience the gravity term is omitted in Eq. 1. The effect of gravity on rod load and stretch can be treated separately, as will be noted later. Since Eq. 1 is linear, the legitimacy of this procedure is easy to demonstrate. POLlSHED ROD MOTION SIMULATION The motion of the polished rod is determined by the geometry of the surface pumping unit and the torque-speed characteristics of its prime mover. By determining the motion of the polished rod, we formulate an important boundary condition. From trigonometrical considerations it can be shown that the position of the polished rod vs crank angle 0 is given by (see Fig. 1) These equations are obtained from the general solution of the "four-bar" linkage problem and can be used to describe the kinematics of any modern beam pumping unit.' If prime mover speed variations are disregarded, the angular velocity of the crank is constant, and Eq. 2 can be used to predict the position of the polished rod vs time. However, the constant-speed condition leading to constant crank angular velocity is only approached in practice: hence, it is better to make provisions for prime mover speed variations in the mathematical model. The speed at which the prime mover runs is determined by its torque-speed characteristics and the torque imposed upon it. The torque that the prime mover "feels" is the net torque arising from the polished rod load and the opposing torque from the counterbalance effect. The

By K. H. Ribe

Water often enters an oil reservoir during completion or workover operations on a well and forms a partial "water block" to oil production. A mathematical study of radial two-phase flow, neglecting capillary effects, has been employed to study the formation of such a water block and subsequent re-moval by production from the well. The effects of reduced oil permeability about the well on the well productivity were studied. The fluid saturation distributions about the well during formation and removal of the water block have also been computed. Several relative permeability relations and viscosity ratios were employed. If water has invaded the formation, its influence through relative permeability effects alone can cause the following. 1. Oil productivity will be depressed for extended periods after production is resumed and will build up only gradually as the water is removed. 2. Oil injected for treatment of water blocking will delay rather than promote restoration of full well productivity by enlarging the region invaded by water. Thus, unless the specific action of chemicals contained in the oil is needed, oil injection appears undesirable. INTRODUCTION During oilwell workover operations, water may enter the oil-bearing formation from the wellbore. When production is resumed, oil must flow through the region invaded by this water. The presence of this region can cause both well productivity and oil production rate to be low and oil to be produced with high water-oil ratio for some time after production is initiated. This situation is sometimes described as a water block. The introduction of water into the formation may result in other actions which also lead to reduction in well productivity and which are also usually included in the connotation of the broad term, water block. Often considered, for example, are the possibilities of clay swelling by contact with fresh water and the formation of emulsions with the formation oil. If it is suspected that such specific actions have taken place, remedial treatments are undertaken which usually involve the injection of chemicals in oil. Since the introduction of water, even in the absence of specific interactions with the formation or oil, will cause a temporary water block (which might be misinterpreted as evidence of a more severe situation), it is of importance to evaluate the magnitude and duration of this blocking which results purely from the reduction in relative permeability to oil in the vicinity of the wellbore. It is also of interest to evaluate the effect of oil injection on the productivity of a well blocked by water in this manner. Inasmuch as this unfavorable condition may persist for some time, it may lead to premature condemnation of a workover or premature abandonment of a potentially productive pay zone. A quantitative evaluation of the influence of water entry on the oil productivity through changes in relative permeability was made by solving a radial form of the Buckley-Leverett equation. The distribution of water saturation around the wellbore during the entry of water was calculated and was followed by a similar calculation of the saturation distribution during the period of resuming production. At any stage in the removal of the invading water, knowledge of the distribution of its saturation permitted calculating the attendant loss in oil productivity. The influence of the shape of the relative permeability relationships was also evaluated by carrying out the calculations for two hypothetical cases. Further, the effect of the oil-water viscosity ratio was examined by repeating the calculations, for several ratios of unity and greater, with the same relative permeabilities. Fi-nally, results are presented to show how the length of time a well must be swabbed to resume production depends on the length of time it has been subjected to invasion by water. STATEMENT OF THEORY Differential Equations Water is assumed to enter a producing formation which is initially at the connate-water saturation and contains no gas. The water and oil are treated as in-

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 K. A. Slagle, D. K. Smith

weight obtained. Additives used in conjunction with salt in these slurries have included silica flour, calcium ligno-sulfonate and cellulose retarders, granular lost-circulation materials, bentonite and selected low-water-loss additives that are not significantly deteriorated by the presence of the chloride ion. In some other areas of South Texas, the salt-saturated slurries have been used quite extensively for improvement of the flow properties of slurries and attainment of better circulation characteristics at lower displacement rates. Concurrently with this property, the protection of shales and shaly sands is also realized, as well as useful retardation of the slurry. Resultant superior cementing jobs have been indicated by both communication tests and acoustic logs for bonding to pipe and formation.' In one section of Louisiana, a major oil company has been using salt-saturated API Class A cement with calcium lignosulfonate retarder for cementing through the Miocene at 9,000 to 10,000 ft. This is another of the situations where interbedding of sands and shales exists, creating difficulty in maintaining formation competency when a fresh-water slurry contacts the clay minerals of the formation. Further work has also been done in the shaly Miocene formation at 13,000 to 15,000 ft where fairly close water or gas contacts are encountered. Indications thus far are that better segregation of these various fluids is obtained by use of the saturated salt slurries because of their improved formation bonding characteristics. In addition to the properties of salt in this situation, attainment of turbulent flow at minimum displacement rates has also been accomplished by use of an additive to help provide exceptional dispersion and viscosity reduction of the slurry. Another oil company was encountering considerable expense in completing wells in Southwestern Louisiana due to extensive block squeeze requirements for effective separation of zones. A very effective mud program was being used to minimize washout in the shale sections and, apparently, a nearly gauge-size hole was being obtained. However, primary cementing results with fresh-water slurries were generally poor. On a few occasions when slurry was actually circulated to the surface, large pieces of shale formation were brought out of the hole with the slurry, indicating a severely water-sensitive, sloughing formation. Inhibition of shale heaving was being accomvlished in the drilling program, and immediately indone ipon circulation of the fresh-water slurry even though it contained a low-fluid-loss additive to reduce filtrate damage in the sands. The subsequent change to salt-saturated slurry yielded 11 successful primary cement jobs out of 12, compared to the previous success ratio of practically zero. Since these were deep, high-temperature wells in the range of 13,000 ft with high-pressure zones necessitating 17.5-lb/gal fluid densities, the slurry used was API Class E cement, silica flour, weight material, retarder, salt saturation (which also reduced the amount of weight material) and maintenance of low fluid loss by use of a salt-compatible additive. Other salt slurries have been used to a limited extent in this same general area for similar problems at depths ranging from 5,000 to 17,000 ft. In the shallower wells, the cement has usually been API Class A where the salt functions as a retarder, and in the deeper wells API Class E cement is used with the additional salt advantage being its increased slurry weight and inability to dissolve salt stringers. A considerable number of squeeze jobs have also been done on older wells using the salt-saturated slurries with very good results. MID-CONTINENT In North Texas, salt-water slurries have been used for cementing the Woodbine sand, Strawn sand, KMA sand and Pettit lime. Shales surrounding these formations have created the same difficulty in obtaining separation of producing zones that has been the problem in other areas. Depths in this area range from 3,400 ft for the Strawn to 7,400 ft for the Woodbine, and concentration of salt has varied accordingly. In the deeper wells, where retardation is desired, saturated salt-water cement is used; for the shallower wells, in order to provide shorter waiting-on-cement times, the amount of salt has been 18 per cent by weight of the mixing water. Results have been excellent with no reported failures on any of the salt cement jobs; where acoustic bond logs have indicated indifferent bonds previously, they are indicating very good bonding for the salt slurries. In Oklahoma various shales of Pennsylvanian age exhibit a high degree of sloughing in the presence of fresh water, causing severe washouts above and below sand formations which it is desired to isolate. This situation exists to some degree in practically all parts of Oklahoma and includes formations of other ages such as the Wood-ford shale. For the past few years, salt-saturated cements and displacement rates as high as practical have been used as a remedy for this problem, with very good results. The Layton and Bartlesville formations are two examples of shaly sands where saturated slurries have been helpful. In one area where five wells were drilled through this type of problem shale without obtaining a satisfactory primary cement job, a change was made to salt-saturated cement preceded by a suitable chemical wash for the drilling mud involved. Acoustic bond logging indicated excellent bonding, and final completion bore out this result by being trouble-free. This type of slurry has also been used extensively on squeeze jobs where shales have been heaving around the producing formations. Predominantly, the basic slurry has been either API Class A cement or a pozzolan cement—although, as deeper wells are being drilled, the use of salt in Class E cement is also increasing. Salt cement in West Texas has been used primarily to help prevent channeling through the shales around the Delaware sand, Queens dolomite and Hope lime. Many of the shaly and dirty sands of this area are sensitive to the filtrate from a fresh-water cement. Salt at 16 to 18 per cent by weight of mixing water has been added to centent, and has been effective in controlling formation damage and communication between zones in these formations. Use of these lower salt concentrations is dictated in this area by the relatively low formation temperatures where retardation of the slurry would create unduly long waiting-on-cement times. Also, quite a bit of cementing has been done in this area using salt concentrations in the accelerating range— that is, 2 to 5 per cent by weight of water. Specifically, these concentrations have been used in coiljunction with high percentages of gel to overcome the retarding effect of the calcium lignosulfonate dispersant, although there are probably several shales where these concentrations could provide some degree of formation stability. On several occasions, the salt has also been used to lower the critical velocity or rate for turbulence with the slurry, particularly in the pozzolan cements. On wells in the Hope lime, it has usually been necessary to squeeze the shale above and below the lime to get a water-free completion. Use of salt-saturated slurries has largely eliminated this

By H. A. Wahl, J. M. Campbell

This study extends our information on solid-liquid slurries to the flow of sand in horizontal fractures. Inasmuch as this is basically an unsteady-state process, a comprehensive photographic study was undertaken in a 10-ft windowed cell to determine if the basic flow regimes described for steady-state flow in pipes applied to the subject process. Since the number of potential variables far exceeds the capacity of a single study, emphasis has been placed on the effects of sand concentration, oil viscosity and oil flow rate. The extensive photographic evidence obtained has proven very valuable in gaining an insight into the basic flow mechanisms. Being able to follow visually the flow characteristics that accompany the quantitative data is valuable in the application of the results. Although the use of dimensionless parameters was carefully investigated. it was found that the data obtained could be more easily, and as accurately, correlated by judicious use of the dimensional variables investigated. However, a study into the feasibility of scaling slurry flow was made in the event this technique is justified in future investigations. The data presented show that the pressure behavior observed in solids transport in pipes basically applies to slurry flow in horizontal fractures. The roles of the parameters are altered but a basic equivalence exists. The most significant correlating parameter was the oil viscosity (µo) and the bulk velocity of the slurry (vn), expressed as ''µv" product. The most significant correlation expresses the rate of advance of the sand as a function of the variables investigated. There are many practical ramifications of this phase of the investigation that should aid in better treatment design. Evaluation of sand advance rates provides a means of estimating sand placement efficiencies during a treatment and the resulting sand distribution in the fracture. The results show that sand placement efficiencies are low under typical treatment conditions. A brief description of the effects of overflushing is also included. INTRODUCTION The flow of sand-oil slurries in fractures is an area in which little basic knowledge is available. This stems to some degree from the fact that it is impossible to duplicate fractures at the surface. They occur in various shapes and sizes with an infinite combination of irregularities. Unfortunately, we can never "see" these fractures except in cores and by indirect means of measurement. In spite of this inherent difficulty, it is desirable to develop some basic concepts that will provide a better understanding of the sand transport mechanism. An insight into the problem is provided by investigations of fluid flow in rectangular conduits. Several studies on the flow of liquids in non-circular conduits1,13 show that a Reynolds number-Fanning friction factor relationship can be written if the hydraulic diameter is substituted for the regular diameter in a circular pipe. This hydraulic, or equivalent, diameter is taken as four times the cross-sectional area occupied by the flowing fluid divided by the wetted perimeter. Eq. 1 expresses an extension of this same work when applied to infinite parallel planes b distance apart.' Eq. 1 is a theoretical equation expressing the friction factor as a function of the Reynolds number for laminar single-phase fluid flow. This expression has been verified experimentally. The equivalent expression for a smooth circular conduit differs only in that the value of the constant is 16 instead of 24. Numerous studies have related friction losses to Reynolds number in both circular and non-circular conduits. These results are widely used and are not reviewed here. Huitt' investigated the effect of surface roughness on fluid flow in simulated fractures. He concluded that fluid flow in fractures may be treated similarly to fluid flow in circular conduits. This work, together with that of Nikuradse,' shows that surface roughness has no appreciative effect upon the resistance to flow in the viscous flow region. In the region of turbulent flow, surface roughness is a prominent factor. Hydraulic conveyance literature is another important source of information. Durand3 has attempted to organize systematically the variables involved in hydraulic-solid transport in pipes. He has classified the modes of flow into three types according to the size of the particles in the mixture— homogeneous mixtures, intermediary mixtures and heterogeneous mixtures. With the usual concentrations and flow rates used in hydraulic transportation, particles with diameters of less than 20 or 30 microns form eszentially homogeneous mixtures with water. The data show, however, that even small materials will tend to settle out under laminar flow conditions. Mixtures containing solids over 50 microns in diameter do not achieve total homogeneity even under turbulent flow conditions. Particles from 50 microns to 0.2 mm in diameter may be transported in fully suspended flow at normal transport velocities although the concentration in the vertical plane is not uniform. Above 2 mm in diameter solid materials are transported along the bottom of the conduit at a velocity substantially less than that of the liquid itself. Between 0.2 and 2 mm in diameter, the particles tend to be in a transition zone between heterogeneous suspended flow and deposit flow at normal hydraulic transport velocities. The sand sizes used in fracturing usually fall in this size range. It is interesting to note that the grain size range designated by Durand for this transition zone corresponds closely to the transition zone between

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

By J. Handin

Triaxial compression tests with independently applied external confining pressures and internal pore pressures show that the ultimate compressive strengths of representative oil well cements are nearly linear functions of effective pressure — the difference between external and internal pressures on the jacketed cylindrical specimens (to 15,000 psi). The strengths are little affected by the test temperature to 350F (not to be confused with the curing temperature). At an effective pressure of 15,000 psi, strengths are in the rage of 30,000 to 50,000 psi, comparable to those of sedimentary rocks under similar conditions. The cements become very ductile even under low effective pressures; permanent shortenings of 30 per cent or more are attainable without rupture. INTRODUCTION Since the pioneering work of Richart, Brandtzaeg and Brown1 on the failure of cement under combined compressive stresses, it has been recognized that ultimate compressive strength is greatly enhanced by the application of confining pressure. More recently, Mc Henry2 showed that the strength of concrete was a linear function of the effective pressure (the difference between the external confining pressure on a jacketed specimen and the internal fluid pore pressure), at least for a range of 0 to 1,500 psi. The effect of temperature had not been investigated, and no previous systematic triaxial compression testing of naterials used for oilwell cementing seems to have been done. The present work was suggested by the late J. hi. Bugbee3 who stated that "consideration of the common application of high-pressure hydraulic fracturing to the initial completion or recompletion of wells, and the large pressure drawdowns in some producing wells, particularly those in abnormally high-pressure gas-condensate reservoirs, raises the question of what is a suitable cement strength for various completions. The intuitive belief exists that cement strength need be no greater than formation strength. Tests should, however, be conducted at downhole conditions." The ultimate compressive strengths of rocks penetrated by the borehole must rise several fold with increasing depth.4 This marked enhancement of strength is due to the influence of the effective pressure, the total weight per unit area of the overburden less the hydrostatic pore pressure. (The effect of temperature due to the geothermal gradients is relatively small for depths to 30,000 ft.) .A significant comparison of the strengths of rocks and cements at downhole conditions requires knowledge of the confined compressive strengths of cements as well. EXPERIMENTAL PROCEDURES The theory and technique of triaxial compression testing are fully discussed in earlier reports.4-6 Briefly, cylindrical specimens 1-in. long and 0.5-in. in diameter are jacketed in thin copper tubes of negligible strength, placed in the test chamber and subjected to an external confining pressure of kerosene and loaded axially by the piston at a strain rate of 1 per cent per minute. Pore pressures of water (or kerosene) are applied independently through the hollow piston and are maintained constant during the shortening of the specimen, Tests at sensibly 0 pore pressure are arranged so that any free water in the cement can escape to the atmosphere. (If egress of water were denied, pore pressure would rapidly attain the value of the external confining pressure because of reduction of pore space.) The test chamber can be heated for high-temperature experiments. Unless otherwise noted, the cement samples were air dried for about a week. Recorded during a test are pore and confining pressures, shortening and axial differential force (total force less the product of the confining pressure and the area of the piston). True stress-strain curves are constructed by dividing the measured differential force by the actual cross-sectional area of the specimen at a given

Jan 1, 1966

By J. C. Savins, G. C. Wallick, W. R. Foster

The dual differentiation-integration method of rbeological analysis is applied to Couette flow. Using machine processing techniques, a spectrum of characteristic derivative functions for a variety of ideal Generalized Newtonian and viscoplastic models has been developed; As for the case of Poiseuille flow the functions for Generalized Newtonian models as a class are strikingly different from the corresponding functions for the class of viscoplastic models. It is shown that this dual scheme of analysis is a highly sensitive analytic method for dotermining the applicable range if a rbeoloaical model. The unique characteristics of this function in the case of viscoplastics may lead to a more precise detection and evaluation of the yield point or yield stress in these materials. INTRODUCTION It has been shown 1,2 that the dual differentiation-integration method of analysis is a discriminating and flexible method for interpreting flow behavior and determining the applicable range of a particular model from Poiseuille flow experiments. In the preceding paper3 the basic equations governing the Couette flow experiment were reviewed and the significance of changes in boundary conditions on data analysis in the case of the viscoplaaic discussed. In this paper the dual method of analysis is applied to Couette flow using the suites of rheological models considered previously.