Search Documents

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst's horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter's edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

By F. W. Schremp, J. W. Chittum, T. S. Arczynski

This paper describes a laboratory study of causes of external casing corrosion and the test work that led to the use of oxygen scavengers to prevent this attack. External casing failures are classified as water-line, casing-casing, collar and body failures. A corrosion mechanism based on principles of differential oxygen availability is developed that is consistent with facts known about each kind of failure. The field use of oxygen scavengers is depicted as a direct result of the laboratory study. A part of the paper is devoted to reporting on the field use of hydra-zine to control external casing corrosion. Results of field measurements made over a period of several years are presented as evidence of the efectiveness of the hydrazine treatment. The first conclusion reached is that the use of hydrazine materially reduces the cathodic protection requirements for treated wells. This result is interpreted to mean that a reduction is taking place in the amount of corrosion on the casing. Results indicate also that hydrazine shows its greatest usefulness within the first 12 to 18 months after a well is completed when pitting corrosion is likely to be most active. INTRODUCTION According to surveys sponsored by the National Association of Corrosion Engineers,' the cost of repairing casing leaks caused by external corrosion may exceed $4 million per year. In addition, well damage and lost production resulting from casing leaks probably costs the petroleum industry an additional $5 to $6 million per year. Concern about the cost of external casing corrosion led to an extensive laboratory study of factors causing this external corrosion and to the development of a new approach to its prevention. This paper presents a discussion of various causes of external casing corrosion, details of laboratory studies and the results of the field use of an oxygen scavenger in well cementing fluids to prevent the external corrosion of oil-string casing. Measurements on test wells over a period of several years show that cathodic-protection current requirements are greatly reduced when hydrazine is used in cementing mud. Reduction of current requirements can be interpreted to mean that removal of oxygen by hydrazine has greatly suppressed corrosion cells on the external surface of the casing and thereby, has reduced corrosion. To date, hydrazine has been used by the Standard Oil Co. of California in more than 200 well completions. KINDS OF CASING FAILURES A survey of a large number of casing leaks disclosed four types of external casing failures — water-line, casing-casing, collar and body failures. These types are identified largely by their location on the casing. Water-line failures are found just below the surface of water or mud in the casing annulus. Casing-casing failures occur on the oil string just below the shoe of the surface string. Collar failures are found in the threaded ends of casing joints where they are screwed into casing collars. Body failures may occur at any point on the body of a casing joint. Ex- amples of each kind of failure have some of the general characteristics that are shown in Fig. 1. Water-line failures usually result in the circumferential severance of an oil-string casing. The corrosive action causing a water-line failure usually is sharply defined and is limited to a short length of the casing. Casing-casing failures usually are accompanied by pitting corrosion distributed around the oil-string casing for distances up to 100-ft below the shoe of the surface string. Casing-casing failures may also sever the casing. Collar failures seem to start on the first thread at the bottom of recesses between collar and casing joint. Corrosion proceeds across the threads by what appears to be a normal pitting mechanism. Both casing and collar are severely attacked. Body failures are the result of highly localized pitting at any point on a casing wall. Besides the pit that perforates a casing, a large number of other pits usually are found along one side of the casing joint. The pits occasionally are filled with corrosion products consisting largely of oxides and sulfides.' Frequently, the mill scale is largely intact on the rest of the casing. Examination of a casing failure does not always reveal the cause of the failure. Frequently, the necessary details are destroyed when the failure occurs. For example, formation water flowing through a perforation at high velocity may enlarge the hole and destroy any remaining evidence of the cause of the failure. One way to obtain undistorted information about a failure is to study the nature of other pits on the casing in the vicinity of the failure. A study of such pits frequently suggests that they are characteristic of an attack resulting from the differential availability of molecular oxygen.

By F. G. Miller, R. C. McFarlane, T. D. Mueller

During the process of gas storage in pressure-depleted oil reservoirs, it has been observed that in some instances additional liquid oil is recovered and that the composition of the storage gas is materially altered. A mathematical study was made of the dynamic behavior of such a depleted oil reservoir undergoing gas injection. The important variable considered in this study, not included in previously published work, was that of compositional effects on the phase behavior of two-phase flow. Pressure, saturation and component composition profiles were developed for a linear, horizontal and homogeneous porous medium containing oil and gas but undergoing dry gas injection. Special new techniques were developed to overcome the problems of numerical smoothing which arise in the solution of the equations representing such systems. The method of solution includes the development of partial differential equations describing the behavior of the system, representing these equations by finite difference approximations, making certain simplifying assumptions and, finally, applying methods of numerical analysis with the aid of a high-speed digital computer. In an example calculation, results using the mathematical model are compared with field observations made on a gas storage project in Clay County, Tex. This field project involved a depleted oil reservoir used 'for gas storage and gas cycling purposes. As a result of these processes, the reservoir yielded substantial amounts of secondary oil, both in the form of stock tank oil and as vaporized products in the produced gas. The methods derived in this study may be applied to a variety of oil reservoir problems which are dependent on compositional effects. In recent years the number of oil reservoirs being used for gas storage purposes has increased greatly, and there has been at least one published account of additional oil recovery resulting from gas cycling a depleted oil reservoir after repressur-ing with dry gas for storage purposes. Additional oil recovery from oil reservoirs resulting from gas storage operations could become an important secondary recovery process. This is especially true since the use of natural gas in large metropolitan areas continues to increase and more gas storage volume near these areas is needed. These facts provided the motivation for the work reported here. This paper reports on a study of the inter-relations of composition, saturation and pressure changes which occur when hydrocarbon gas is injected into an oil reservoir system. From an understanding of the process, prediction methods may be developed for use in forecasting the secondary recovery products from gas storage operations in oil reservoirs and, consequently, .the economics of such projects can be developed. The process of storing hydrocarbon gas in an oil reservoir and the cycling process both involve the transient multi-phase flow of oil reservoir fluids with changing fluid compositions. There is a scarcity of work published on the subject. One of the first laboratory models developed to study a problem similar to the one posed here was published in 1948 by Standing, Lindblad and Parsons. l They developed a model to study gas condensate recovery by gas cycling, including the effect of heterogeneous reservoir permeability. Their results showed the compositional effects of injecting a gas of known composition into a two-phase condensate reservoir of known composition. A major finding of their study indicated that maximum oil recoveries from gas cycling projects result when low cycling pressures are used. Previous to their work it was generally considered

By Henry M. Howe

My attention has been drawn within a few days to a series of articles in Volume XVIII of the Engineering and Mining Journal, 1874, by Mr. J. A. Church, in which it is stated, among other

By C. F. Gates, H. J. Ramey

Recent literature shows that pronounced increases in oil recovery can result from the use of miscible systems in recovery operations. This literature also points out certain problems associated with maintaining miscibility, e.g., compositional changes resulting from alterations in pressure and temperature or from zone dilution. The purpose of the work described in this paper was to study miscible zones for possible use in waterflooding operations and examine the manner in which these zones break down. The fluid system selected for study was oil—terNary-butyl alcohol—water, which, because of its narrow range of miscibility, is ideally suited for investigation in short, linear systems. The laboratory-observed effects on oil recovery and zone stability of a number of variables are reported. Among the variables investigated were miscible zone size, viscosity ratio, length of travel, and interstitial water saturation. Results obtained with this system show that interstitial water adversely affects recoveries because of premature phase break. Salinity of the water further aggravates Ais situation. This effect, although pronounced for this system, would probably occur in any fluid system containing an alcohol. As expected, viscosity ratio has a decided influence on recovery efficiency and zone stability. The results also show that use of a suficieno zone size, even though the system has poor phase characteristics, yields higher oil recoveries than are obtained from straight water floods. Length-of-travel studies showed a square root relationship between zone growth and path length for favorable viscosity ratios. No such clear-cut dependence was observed under unfavorable viscosity conditions. INTRODUCTION As new oil reserves have become increasingly more difficult and expensive to find and develop, the oil industry has devoted more and more time and money to finding more efficient methods for exploiting known reserves. The increasing number of water floods is one manifestation of this attempt to improve oil recoveries from existing fields. Since considerable quantities of oil are by-passed by the waterflood process, it is only a partial answer to the problem. The search continues. Of late, several proposed processes have received the critical attention of investigators. Among them are gas drives, in situ thermal reactions and miscible displacements. The present paper is concerned with miscible displacements in conjunction with water floods. Recent literature1- reports that, under proper conditions, oil recoveries from LPG floods, propane sweeps, condensing gas drives, etc., approach 100 per cent of the oil in place in the region contacted by the flood. This literature also points out some of the problems associated with these processes, one of the main ones being the maintaining of miscibility. Compositional changes resulting from pressure and temperature variations or from zone dilution give rise to phase breaks and loss of miscibility, the conclusion being that, after loss of miscibility, much of the beneficial influence of the zone is lost.' The results of the subject work show that, if sufficient zone material exists ahead of a phase interface, considerable improvement in oil recovery can be obtained. In order to increase our understanding of the behavior of miscible systems as they degenerate into immiscible processes, a laboratory system having rapid deterioration characteristics was needed. The fluid system oil—ter tiary-butyl alcohol (TBA)—water possesses these properties and was the system investigated. Although these fluids are different from those generally employed in miscible studies, an understanding of their behavior can be applied to other miscible slug processes. However, under conditions where mass transfer maintains a miscible system through multiple contacts, e.g., enriched gas drive or high pressure gas drive, the observed behavior of alcohol zones would not be applicable. The variables of interest were zone size, viscosity ratio, length of travel, initial water saturation, rate and core characteristics. Of these, only the effects of the following on zone breakdown and oil recovery are reported in detail: (1) zone size, (2) path length, and (3) initial water saturation. The effects of viscosity ratio and core characteristics are discussed in conjunction with other variables. Although initial water saturation is not envisioned as important to the phase behavior of hydrocarbon systems, e.g., LPG, it is important in alcohol systems. The properties of the fluids and cores used are given in Table 1. As previously stated, this fluid system was selected because of the limited range of miscibility

The first session of the Institute was held at the Smithsonian Institution, on Tuesday evening, February 22d. The members were welcomed to Washington and to the Smithsonian by Prof. Joseph Henry. President Holley responded to the cordial welcome of Professor Henry, and then read an introductory address on the subject of Technical Education. Following the discussion on President Holley's address the following persons, duly proposed and indorsed by the Council, were elected members and associates of the Institute : MEMBERS. Emile R. Abadie,. .. .. St. Louis, Mo. R. J. Anderson,.....Pittsburgh, Pa. Commander L. A. Beardslee,. .. Washington, D. O, Charles Bender,. .. .. New York City. Charming M. Bolton,. .. . Richmond, Va. Thomas E. Bowman,. .. . Silverton, Colorado. Calvin E. Brodhead,. .. . Bethlehem, Pa. M. L. Brosius,. .. .. . Lewistown, Pa. Harvey H. Brown,. .. .. Cleveland, Ohio. James A. Burden,.. .. . Troy, N. Y. J. P. Clemes,. .. .. . Sonora, Mexico. N. H. Cone,. .. .. . Nederland, Colorado. W. F. Durfee,.. .. . New York City. Robert Frazer, Jr.,. .. . Ashland, Pa, Joseph E. Galigher,.. .. . Bingham, Utah. William Golding,.. .. . Steel Works P. O, Pa. Samuel E. Griscom,.. .. . Pottsville, Pa. Gen. E. Burd Grubb,.. .. Beverly, N. J, Arthur R. Guerard,. .. .. Charleston, S. C. Edward Hart,.. .. .. Easton, Pa. Prof. Lewis M. Haupt,. .. . Philadelphia. E. S. Hawley,.....Buffalo, N. Y. John W. Hoffman,.. .. .. Philadelphia. Nelson P. Hulst,.. .. . Milwaukee, Wis. H. C. Humphrey,.. .. . Philadelphia. Ed. P. Jennings,.. .. .Ithaca, N. Y. B. F. Jones,......Pittsburgh, Pa.

By J. D. Testerman

A statistical technique to identify and describe naturally occurring zones in a reservoir and to correlate these zones from well to well is described. The technique is particularly useful in describing a reservoir where cross flow between adjacent strata is important in determining reservoir behavior. Although it has been used primarily for permeability zonation, the technique is general and can be used to correlate any reservoir property or related data, such as the information contained in well logs. INTRODUCTION One of the first problems encountered by the reservoir engineer in predicting or interpreting fluid displacement behavior during secondary recovery processes is that of organizing and using the large amount of data available from core analysis. Permeabilities pose particular problems in organization because they usually vary by more than an order of magnitude between different strata. Due to the sheer volume, it is almost always necessary to group data and to use an average value to represent a number of measurements. Perhaps the most common method now used to group permeability data is the capacity-fraction technique, which ranks permeabilities in order of magnitude, regardless of the physical location of the permeabilities within the reservoir. The cumulative per cent capacity is plotted against cumulative per cent thickness. This plot is divided into an arbitrary number of zones, generally of equal thickness. Five zones (or averaged groups of data) usually are obtained, each of which is treated as homogeneous in subsequent calculations. The division so obtained has no physical meaning; strata in the same zone, calculation-wise, are usually not adjacent in the reservoir. Reservoir engineering techniques being developed will handle crossflow that occurs between adjacent communicating reservoir strata because of imbibition and gravity segregation. Since crossflow occurs between physically adjacent layers within the reservoir, a new zonation technique recognizing the actual location of strata within the reservoir is necessary. Similarly, the recognition of natural zones is important for predictions of oil recovery by processes involving diffusion. One such process is miscible displacement, where predictions of lateral diffusion within the reservoir must recognize the actual location of the invaded zones in relation to the rest of the formation. Natural zones must also be adequately recognized to account for heat transfer within the reservoir during thermal exploitation. Because of the complexity of the problem, statistics appear to offer the only practical hope of dividing a reservoir into physically-meaningful natural zones. This paper presents a statistical technique for identifying these natural zones and for ascertaining which ones are likely to be continuous between adjacent wells. The zones defined have minimum variation of permeability internally and a maximum variation between zones. The technique is general and can thus be applied to reservoir properties other than permeability. The method will guide the reservoir engineer in estimating which zones are likely to be continuous between wells. However, a statistical correlation based on permeabilities in two different wells is no guarantee that the zones so defined are, in fact, continuous. Rather, the assumption of continuity must be consistent with geological data concerning the depositional environment, as well as justified on the basis of engineering judgment in combination with statistics, just as judgment is required with conventional zonation methods. CALCULATION PROCEDURE The reservoir zonation technique is a two-step operation. The steps are individually described, and a sample calculation is presented in the Appendix. ZONATION OF INDIVIDUAL WELLS First, the set of permeability data at a single well is zone, into Zones. These zones are selected so that variation is minimized within the zones and maximized between the zones. The equations4,6 used to zone the data are where B = the variance between zones, , = the number of zones, i = the summation index for the number of zones, j = the Sumation index for the number of data within the zone, mi = the number of data in the ith zone, k,. = the mean of the permeability data in the ith zone, k . = the over-all mean of the data in the well, W = the pooled variance within zones, N = the total number of

By D. G. Russell

A simple method has been developed with which flowing bottom-hole pressure data from two-rate flow tests in oil or gas wells can be analyzed to estimate the formation permeability, skin factor and average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized producing rate of the well is changed to another, higher or lower, rate. The new method yields the same information as a conventional pressure buildup analysis, but eliminates the need for closing in the well. The analysis of a two-rate flow test is of the same degree of difficulty and requires about the same engineering time for application as a conventional pressure buildup analysrs. The extended closed-in periods experienced with conventional buildups because of long, low-rate afterproduc-tion periods are eliminated by fIow tests. Other anomalous pressure buildup effects, such as "humping" due to weli-bore phase segregation, can be successfully eliminated with the new method. Generally, flow tests of about 24 hours' duration run with conventional pressure measurement equipment are sufficient for interpretation purposes. Field testing of the two-rate flow test method has established that it is a reliable and economical method which can be used in many instances to complement or even replace conventional pressure buildup methods. INTRODUCTION The principal method for estimating formation permeability and well damage, or skin factor, in a producing oil or gas well is the analysis of shut-in bottom-hole pressure buildup data.' This familiar method has been used quite successfully by reservoir engineers for many years. It is based on the solution of the radial flow equation for constant rate conditions, and requires that the well be closed in for a sufficient period of time to obtain a clearly defined linear portion on the plot of observed t + ?t bottomhole pressure vs log t + ?/?t(where At is shut-in A? time, and t is producing time to the instant of shut-in). From the slope of the plot and other normally obtainable data, the permeability, skin factor, and reservoir pressure at infinite shut-in time (if the reservoir were infinite) can be estimated. Over the years several drawbacks have become apparent in the use of conventional shut-in pressure buildups for determining permeability and skin factor. The conventional pressure buildup interpretation theory assumes that a well is closed in at the sand face and that no production into the well occurs after shut-in. In practice, of course, the well is closed in at the surface, and inflow into the well continues until the well fills sufficiently to transmit the effect of closing in to the formation. This adjustment period is commonly referred to as the "after-production" portion of the pressure buildup. In the tight reservoirs, long, low-rate after-production periods frequently occur, and the well must be shut in for several days or, in some instances, even weeks to obtain an interpretable buildup curve.' Obviously, such long shut-in times can cause loss in current income, both from reduced oil production and from the fact that personnel and pressure measurement equipment are occupied with a single well for too long a time. In other cases, even long shut-in periods do not seem to be of much aid in obtaining an interpretable buildup." If there is considerable phase redistribution (liquid fallout or bubble rise) after a well is shut in, then curves with no interpretable portion are often obtained during the buildup. In addition to instances in which wellbore effects cause trouble, there are also cases in which the major objection to use of the closed-in pressure buildup is simply the fact that the well must be shut in. When there is no proration and when the well has limited producing capacity. closing in the well means loss of income. From the foregoing discussion it is apparent that it ib desirable to have an alternative method of obtaining the same information as that derived from a conventional buildup without the need of closing in the well. One possible solution which has been offered for this problem is the use of a bottom-hole shut-in tool' which isolates the major portion of the flow string from the formation face during the buildup. In this paper an alternative method which is frequently successful in avoiding wellbore effects and which does not require the use of special equipment is presented. A new, simple method has been developed with which the flowing bottom-hole pressure data from flow tests in oil or gas wells can be used to estimate permeability, skin factor and the average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized production rate of the well is changed to another, higher or lower rate. The need for closing in the well is eliminated, and pressure measurement periods of only 18 to 24 hours arc usually sufficient, even in tight reservoirs. Thus, the new

By C. S. Matthews, P. Hazebroek, H. Rainbow

It ha been suggested that lormation fractures created by well stimulation treatments will adversely affect sweep-out efficrency in injection operations. Fluid-flow model studies involving vertical fractures of various lengths and fluid systems of various mobility ratios have been carried out to study this subject. In addition, limited data have been obtained on one model containing a horizontal fracture. It was found that relatively long and highly conductive fractures (not generally obtained fracturing operations) were required to affect the sweep-out efficiency substantially. In a given case in the field an approximate distinction can be made between the presence of long conductive fractures and shorter or less conductive ones. This is done with pressure build-up analyses along with data on the relationship of fracture length and conductivity to well productivity. This type analysis shows that usually the fractures induced are either short or of limited conductivity and therefore do not damage sweep-out efficiency. INTRODUCTION Improved well productivity and in-jectivity can frequently be exploited in injection operations. Higher total throughput can yield improved economics. On occasion, the achievement of an increased productivity or injectivity in specific wells can bring about a more uniform sweep of the reservoir. Higher rates can be exploited particularly in water floods of "depleted" reservoirs where a rapid "fill-up" is desired and where low pressures contribute to low well productivity. The creation of fractures local to the wellbore is an excellent means for achieving these objectives. Even though fracturing has been employed in some floods with success,1"3 there still seems to be some reluctance to employ this tool for fear of undue damage to the flood pattern and ultimate recovery. We therefore need to examine the influence which fractures of varying length may have on flood performance and then determine the length of fractures which obtain in the field with conventional fracture treatments. A substantial influence of fractures on the recovery obtained at breakthrough of the injected fluid has been presented by Crawford and Collins for equal fluid mobilities for the line-drive pattern.4,5 In addition, the influence which a fracture has on the production performance after breakthrough and on the ultimate recovery warrants consideration in reaching a conclusion concerning the use of induced fractures in flooding operations. This report presents this type of data for the five-spot injection pattern for several fluid mobilities. In arriving at some conclusion concerning the fracture lengths obtained in field operations we must examine the performance characteristic most affected by the fracture. This is the change in the flow system as reflected in the change in productivity and pressure build-up behavior. A study of these changes is also presented in this report. RESERVOIR ANALOGS The X-ray shadowgraph technique, employing miscible displacement in porous models, has been used in the study of the influence of fractures on pattern sweep-out efficiency. The X-ray shadowgraph procedure is described in detail in an earlier report8. Fractures were represented by leaving the proper portion of the model surface exposed to either injection or production. This assumes the fracture resistance to be negligible compared to that of the formation. Actually the flow resistance in propped fractures obtained in the field is sometimes not negligible so that the results with this model indicate the maximum influence of the fracture. Two types of models were necessary to represent vertical fractures in a five-spot flood. These pattern elements are illustrated in Fig. 1. In studying the influence of the horizontal fracture, only one well spacing length to thickness ratio of 50 was used. The pool unit of a Carter Electric Analyzer was used in studying the influence of fractures on productivity and build-up behavior. The square drainage system of a well was represented by a network of 576 elements of equal volume and resistance. One-fourth of this drainage area was studied as a model unit using a network of 144 condensers and resistors. Vertical fractures were represented by a shunt directed from the well perpendicular to the drainage boundary. Horizontal fractures were represented by a circular shunt. Resistance of a shunt was varied in representing different conductivities for the fractures. FRACTURE DIRECTION AND LENGTH The direction at which a vertical fracture extends into the formation and its length and conductivity influence the sweep behavior. The two

By R. L. Templin

The increasing use of aluminum and its alloys in commercial fields tias demanded a better understanding of their machining properties. This fact is exemplified by problems that have arisen in the automotive and airplane industries, but many in other fields might be cited. As pure alunlinum and its alloys in their various commercial conditions show appreciable differences in their machining properties, it is not surprising that quite divergent solutions have been offered for the machining problems encountered. However, if the fundamental requirements of the most suitable cutting tools for these metals are understood, these machining problems lend themselves more readily to satisfactory solutions. Since the machining of free cutting brass and mild steel is understood by most persons accustomed to working these metals, it may serve our purpose better to first make a general comparison of the tools more commonly used in machining these metals with the tools most suitable for machining aluminum, then proceed to a more specific discussion of the individual tools. Comparison oF Cutting Tools Cutting tools commonly used for machining free cutting brass usually have little, if any, top and side rake; they are ground on a medium to coarse abrasive wheel and used without any cutting compound or with a cutting compound that has a paraffin base. Those ordinarily used for steel have some top and side rake, are usually ground on a medium to fine abrasive wheel, and are often used with soluble-oil cutting compounds. The proper tools for aluminum and its alloys should have appreciably more side and top rake than the tools for cutting steel; should have very keen edges obtained by grinding with fine or very fine abrasive wheels supplemented in many cases by hand stoning with an oil-stone; and should be used with suitable cutting compounds whenever possible. In many cases, tools suitable for machining aluminum and its alloys are not appreciably different from tools commonly used for cutting hardwoods. The front clearance of a tool most suitable for machining aluminum and its alloys should be about 6°, the top rake from 30" to 50°, making the total angle of the cutting edge of the tool from 35' to 55". A side

By J. L. Thompson

There are two basic forms of the LPG-slug process: the gas driven and the water driven. The pressure required for miscibility between gas and LPG prohibits the use of the gas-driven LPG process in shallow reservoirs. The water-driven LPG slug process normally exhibits good sweep efficiency. However, displacement of the LPG by water is poor. An improvement in this process appears possible by injecting a slug of carbon dioxide between the LPG slug and the water drive. Laboratory experiments were conducted in homogeneous linear systems to determine the effect of pressure on the various displacement zones. In the carbon dioxide-LPG displacement zone, the effect of pressure on the displacement efficiency is negligible above the miscibility pressure. However, the displacement efficiency decreases linearly with decreasing pressure below the miscibility pressure. With water displacing carbon dioxide, it was found that the recovery of carbon dioxide after water breakthrough can be predicted from its solubility in water. A displacement test was conducted with LPG and carbon dioxide slugs large enough to avoid interference between the oil-LPG, LPG - carbon dioxide and carbon dioxide - water displacement zones. Under these conditions, essentially complete oil and LPG recovery was obtained. However, a substantial amount of carbon dioxide was left in the core at water breakthrough. INTRODUCTION It has been shown in homogeneous linear systems that essentially complete oil recovery can be obtained by displacement with a miscible fluid such as LPG. Economic considerations dictate that the valuable LPG also be recovered. Methods proposed for recovering the LPG include (1) displacement with natural gas at a pressure high enough to make the two fluids miscible1 (the disadvantage of this method is that the required pressure of about 1,500 psi is too high for shallow reservoirs) and (2) displacement of the LPG with water2 (principal disadvantage is that a large amount of LPG remains in the formation because of poor displacement efficiency). This paper treats an improvement in the water-driven LPG slug process3,4 where carbon dioxide is injected between the LPG and water. Carbon dioxide is completely or at least partially miscible with propane and soluble in water; thus, displacement efficiency is improved at both ends of the carbon dioxide slug. To minimize effects of other variables, all experiments were conducted in linear core systems. GENERAL PROCESS DESCRIPTION The displacements of propane by carbon dioxide and carbon dioxide by water are considered separately, assuming enough of each is injected to prevent mutual interference. Displacement descriptions are confined only to the part of the pore volume containing mobile fluids unless otherwise specified. The schematic representation is given in Fig. 1. Starting at the producing well, a zone containing only oil is followed by a zone containing a single-phase mixture of propane and oil. The concentration of propane increases in the direction of the injection well and ultimately grades to pure propane. Then follows a zone containing a mixture of propane and carbon dioxide in continuously varying proportions until 100 percent carbon dioxide is reached. In this zone, the concentration of carbon dioxide in the immobile water phase also varies from zero to full saturation. Whether the displacement of propane by carbon dioxide is of complete or partial miscibility depends on pressure and temperature (Fig. 2). If the horizontal line representing the prevailing pressure intercepts the boundary curves of the phase envelopes for the prevailing temperature (i.e., if the pressure is lower than the critical

By F. C. Pittman, D. W. Harriman, J. C. St. John

This paper mentions briefly the history of hydraulic jetting as applied to perforating and fracture initiation. It points out the advantages of hydraulic perforation and undercutting as an aid for creating a fracture at the point desired rather than depending upon the weakest point in the formation for breakdown. It describes early experimental work with jet nozzles in search of better nozzle materials. The effect of splash-back damage and its subsequent influence on jet body and tool design during this work is discussed. A series of cutting-rate curves for jets cutting steel and steel-cement-formation combinations is presented to show the effect of hydrostatic head and the point of diminishing returns with respect to cutting time. There is a series of photographs showing various types of rock formation in which perforation and undercutting tests were made. These stones were drilled and the casing cemented in place as in an actual well. The casing was perforated and circularly cut as if preparing for a fracturing job. The conclusion reached in the paper is that hydraulic jetting with sand-laden fluid can be used for perforating and undercutting casing, cement and formation rock for the intended purpose of inducing the formation to fracture at a desired point. INTRODUCTION Hydraulic jetting as a means of cleaning a formation during acidizing has been used for many years. The acid-jet gun was used to clean formation and for drilling cement from casing with acid. This tool was dropped to a seating ring installed on the tubing and could be retrieved by reverse circulation or wire line. In this process, the fluids were directed against the formation to clean and penetrate beyond the mud-damaged area. After the development of the hydraulic fracturing process, interest in jetting the formation with abrasive-laden fluid was renewed. The cutting action of abrasive fluid on pumping equipment, which is a well known problem, has been utilized to perforate and undercut casing, tubing and formation. In this process, surface areas are provided in the formation on which the frac- turing pressures may work to produce a fracture at the location desired. EQUIPMENT The tools used consist of bodies containing two, three or four jets in a single horizontal plane equally spaced around the circumference. The body is threaded to the end of a string of tubing and run into the well. In addition, an expendable tool is available that may be left in the hole after perforating and produced through until necessary to remove the tubing. The bottom of each tool is equipped with a ball-type back-pressure valve. The ball can be reversed to the surface and recovered. This allows reverse circulation so that loose sand which accumulated below the tool while jetting may be removed before the fracturing job. The present jetting tools are of such length that, when screwed together, the jets fall on 12-in. spacing. ANALYSIS OF THE PROBLEM Early work with nozzles indicated that hardened steel was not a suitable material for jets and that more durable material would be required. Ceramics were tried and found to be superior to steel but still too shortlived for most applications. Tungsten carbide, which was tried and found to be the best material available, is now being used for jet nozzles. Analysis of the problems presented by this process indicated that the following information was needed. 1. The effect of submergence on a jet stream and the effect of submergence on the cutting rate of a jet stream. 2. The cutting rate and depth of a jet stream in steel. 3. The cutting rate in various formations and point of diminishing returns. 4. Size of hole or slot cut in a pipe by a given jet. 5. Fluids and sands used. 6. The flow of various weight fluids through the jet nozzle and the discharge coefficients. 7. The type cut formed in the formation by the jet stream. 8. The stretch and contraction of tubing due to temperature and pressure. The cutting effect of a jet stream on a target in the atmosphere is rapid and spectacular as most of the

By J. M. Smith, D. Kunii, P. Adivarahan

Effective thermal conductivities were measured for seven samples of porous rocks through which gases or aqueous salt solution were flowing, parallel and countercurrent to the flow of heat. The results are analyzed in terms of the contributions of the solid and fluid phases to the total energy transfer. The solid contribution for the sandstone samples studied was found to be a Junction of the porosity and an effective thermal conductivity for the solid phase. The fluid contribution was a function of the Peclet number. Based upon the equations presented, the effective conductivity of the system can be estimated from the properties of the fluid, the flow conditions, and the porosity and pore-size distribution of the rock. INTRODUCTION Heat-transfer characteristics of porous media containing fluids are needed for the design and evaluation of thermal methods of petroleum production. The study reported here considers the problem of countercurrent flow of fluid and heat in sandstone rocks, and is a continuation of previous studies concerned with parallel and radial flow of energy and fluid in both rocks517 and beds of fine particles.6,10 In addition to these earlier works, Asaad,1 Kimura,4 Preston,8 Zierfuss and vlietll and waddams9 have reported effective conductivities for systems filled with stagnant fluids. There appears to be no published data on rocks through which fluids are flowing. The purpose of this study was to present such information and, particularly, to investigate the effects of properties of the rock and the flow rate on the results. The rock samples were naturally occurring sandstones from different locations and were similar to petroleum-bearing formations. The properties of the seven samples showing the range in porosities and permeabilities are given in Table 1. The mean pore sizes were calculated from mercury penetration data by averaging the radius with respect to pore volume. The details of this method and comparison with other procedures for estimating 7 are given elsewhere. Four fluids, nitrogen, carbon dioxide, helium and 10 weight per cent sodium-chloride solution were used with each sample. For Bartlesville sandstone, BA-5, runs were made at mean pressures from 15 to 250 psig using nitrogen. The other measurements were carried out at a mean pressure of 100 psig. The maximum temperature range eacountered was 70° to 300°F. The flow rate of fluid ranged from about 1.0 to 130 Ib/(hr)(sq ft), and this corresponded to a modified Reynolds-number range of 5 x 10-4 to 1.0. The variation in Prandtl number was 0.70 to 4.1. EXPERIMENTAL WORK The equipment consisted of apparatus to supply a constant flow of gas, or salt solution, at constant pressure to the bottom of a vertical test section, and metering and flow-regulation apparatus after the test section. The flow of salt solution was

By P. H. Dudley

In giving a brief account of the experiments in progress to inquire into some of the facts in regard to "railway resistances," recently commenced upon the Lake Shore and Michigan Southern Railway, with the dynagraph (which is essentially an instrument recording upon paper the force required to draw a train, at all speeds, marking the time in short intervals, so that the force absorbed due to any

By A. R. Gregory

A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy. S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King, 1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. 5- and P-wave velocity measurements on Bandera sandstone and Solenholen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2 and by Gregory 3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fate sand by Gregory, 3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King 1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8 A new shear wave velocity laboratory apparatus has been developed for testing rock samples under elevated pressures. This apparatus employs liquid-solid interfaces as mode converters. The technique for generating pure shear waves at liquid-solid interfaces was developed from studies of critical angle methods.6 Two identical transducers are placed in contact with opposite ends of the rock sample, one transducer acting as an ultrasonic energy emitter or generator and the other acting as an ultrasonic energy receiver or motor. The P-wave energy emitted from the piezoelectric crystal in the generator transducer is converted into pure S-wave energy at the Drakeol-aluminum

By A. G. Weber, K. H. Coats, M. H. Terhune, R. L. Nielsen

Two computer-oriented techniques for simulating the three-dimensional flow behavior of two fluid phases in petroleum reservoirs were developed. Under the first technique the flow equations are solved to model three-dimensional flow in a reservoir. The second technique was developed for modeling flow in three-dimensional media that have sufficiently high permeability in the vertical direction so that vertical flow is not seriously restricted. Since this latter technique is a modified two-dimensional meal analysis, suitably structured three-dimensional reservoirs can be simulated at considerably lower computational expenses than is required using the three-dimensional analysis. A quantitative criterion is provided for determining when vertical communication is good enough to permit use of the modified two-dimensional areal analysis. The equations solved by both techniques treat both fluids as compressible, and, for gas-oil applications, provide for the evolution of dissolved gas. Accounted for are the effects of relative permeability, capillary pressure and gravity in addition to reservoir geomety and rock heterogeneity. Calculations are compared with laboratory waterflood data to indicate the validity of the analyses. Other results were calculated with both techniques which show the equivalence of the two solutions for reservoirs satisfying the vertical communication criterion. INTRODUCTION Obtaining the maximum profits from oil and gas reservoirs during all stages of depletion is the fundamental charge to the reservoir engineering profession. In recent years much quantitative assistance in evaluating field development programs has been provided by computerized techniques for predicting reservoir flow behavior. Because of the spatially distributed and dynamic nature of producing operations, automatic optimization procedures, such as those now in use for planning refining operations, are not now available for planning reservoir development. However, present mathematical simulation techniques do furnish powerful means for making case studies to help in planning primary recovery operations and in selecting and timing supplemental recovery operations. A number of methods have been reported which simulate the flow of one, two or three fluid phases within porous media of one or two effective spatial dimensions>-4 However, applying computer analyses to actual reservoirs have been limited mostly to two-dimensional areal or cross-sectional flow studies for two immiscible reservoir fluids. To obtain a three-dimensional picture of reservoir performance using such two-dimensional techniques, it has been necessary to interpret the calculations by combining somehow the results from essentially independent areal and cross-sectional studies. To the authors' knowledge, the only other three-dimensional computational procedure, in addition to those presented here, was developed by Peaceman and Rachford8 to simulate the behavior of a laboratory waterflood. Two computational techniques which may be used to simulate three-dimensional flow of two fluid phases are described in this paper. The first method, called the "three-dimensional analysis", employs a fully three-dimensional mathematical model that treats simultaneously both the areal and cross-sectional aspects of reservoir flow. The second method, called the vertical equilibrium (VE) analysis", is applicable

By H. A. Kendall, W. C. Goins

Several investigations in recent years have shown that drilling rates are increased significantly with increased hydraulic horsepower. But, there has been no over-all method of designing jet-bit programs that efficiently uses the surface power. A study of present practices indicates that frequently as little as 50 per cent of the possible effects at the bit are used. Some observers have indicated that the best utilization of hydraulic horsepower (maximum effect on drilling rate) occurs when the bit hydraulic horsepower is maximum; others have stated that jet impact force is more important, and others have believed that maximum jet velocity is required. Limited efforts to date have shown some optimum conditions for bit hydraulic horsepower and impact, but these conditions cannot exist during drilling of a large part of the hole and do not provide a basis for designing a complete jet-bit program. This paper shows the maximum obtainable bit horsepower, impact force and jet velocity at all depths, taking into account the limitations of the pump, piping, hole and minimum circulating rate for adequate cuttings removal. Ranges of operation are developed; and flow rates, surface pressure and bit pressures are specified for each range to provide a maximum of any one of the desired effects. It also is shown that, by proper selection of nozzle sizes and by following the rules presented, the maximum obtainable quantities can be effectively utilized from surface to total depth. Finally, a simple graphical method of selecting nozzle sizes and flow rates is presented which can be used with familiar bit-company hydraulic tables and calculators to design jet-bit programs for maximum bit hydraulic horsepower, impact or jet velocity, as desired. These programs make most effective use of the pumps. Heretofore, there was no method available for designing field tests which adequately separated the effects of bit horsepower, impact and jet velocity. The programs and procedures developed in the paper are dissimilar and, when used in future field testing, should demonstrate which program is the most important in obtaining the fastest drilling rate. INTRODUCTION During the past decade, rig hydraulics has come into increasing prominence. There has been a definite trend toward providing higher horsepower pumps, jet-type bits have had increased use, numerous investigators"" have reported increased drilling rates as a result of increased hydraulics, and bit manufacturers have provided tablesa-" and calculators that are now commonly used 10 design jet-bit programs. Opinion has varied as to the hydraulic quantity which has the great- est effect on drilling rate. Papers and data have been presented that show pump horsepower,'9 it hydraulic horsepower and jet impact force,' each to be the most significant factor affecting drilling rate. Examinatibn of jet-bit programs of the bit companies indicates emphasis on jet velocity. Only pump horsepower can be eliminated because it can be used to produce any one of the bit effects which, a priori, must be more relevant factors. This contradictory state of opinion and practice regarding the bit effects is unfortunate, but several published references have been concerned with making one or another of the factors maximum; and, because these in each case have given results applicable to only intervals of the hole drilled, there seems to be ample reason to complete the previous efforts. It also is believed that the differences in programs for each effect, where they exist, should be delineated so that future use may determine which hydraulic effect is the more relevant. It is the purpose of this paper to: (1) show the theoretical maximum bit hydraulic horsepower, jet impact force and jet velocity available at all depths, taking into consideration all necessary restrictions on operating conditions; (2) illustrate procedures by which the maximum available horsepower, impact force or veIocity may be obtained; and (3) present a graphical method for rapid selection of jet-nozzle sizes and flow rates to be used with conventional procedures to design jet-bit programs for maximum bit horsepower, impact force or velocity as desired.