Search Documents

By Oren C. Baptist, Eliot J. White

Samples of four reservoir sands having different degrees of water sensitivity were subjected to several laboratory tests to determine differences in capillary behavior attributable to clay-mineral effects. The degree of water sensitivity was estimated by the difference between gas and water permeability. Unsaturated pore volumes determined with mercury injection were comparable to irreducible water saturations determined with the semipermeable bar-rier method only for the non-sensitive sands that contained small amounts of non-swelling clay minerals. Irreducible saturations were independent of salinity of the water, which suggests that most residual water is adsorbed 0.7 a film over solid surfaces and that internal retention by clay minerals is relatively insignificant. The water-sensitive sands imbibed water to higher saturations than did non-sensitive sands, and this result is believed to be a combination of the effects of clay minerals and preferential wettability; however, the relative importance of each of these effects is riot determinable. Thus it is concluded that the semipermeable barrier method is the best of the laborafory methods used for estimating formntion water .saturations in the sands studied. INTRODUCTION This investigation is a continuation of studies of the role of clay minerals in petroleum production in the Rocky Mountain area. Previous work has shown large differences in gas and water permeability in numerous reservoir sands; these water-sensitivity effects observed in the laboratory have been correlated with well behavior in the field, and the degree of water sensitivity imparted to reservoirs by the various clay minerals has been determined.'.' The effect of interstitial clay minerals on formation water3 saturations is not so well known. It has been reported that formation water saturations probably are higher in clayey than clean sands,' but this effect has not been evaluated in the laboratory, and the relative effect of the various types of clay minerals on formation water saturation has not been determined. The clay-mineral groups of most importance to this study and to petroleum production in general are kaolinite, illite, and montmorillonite. The minute layers of montmorillonite are loosely bound to one another so that water enters between them, causing great expansion of the particles. Kaolinite and illite do not exhibit interlayer hydration, but they do immobilize considerable amounts of water owing to their large surface areas, which are of the order of 15 to 97 m2/gm of clay.5 The four Wyoming sands used in this study were selected to represent a wide range of water sensitivities. These sands, with the number of wells and fields from which the samples were obtained, are: Second Frontier, five wells in two fields; Newcastle, six wells in one field; Tensleep, two wells in two fields; and Lower Muddy, five wells in one field. EXPERIMENTAL METHODS AND PROCEDURES Water-displacement data were obtained by the semipermeable barrier method," using both fresh water and brines having either 16,500 ppm sodium chloride or calcium chloride. A displacement pressure of 35 psi was found to be sufficient to obtain the practical irreducible water or brine saturation (IWS, IBS) in most samples. Mercury-injection data were ob-tained with equipment and pro-cedures similar to those described by Purcell.7 Exceedingly high pressures are required to attain complete mercury saturation by this method; however, since saturations in most samples change very little with large pressure increases above 1,000 psia, the maximum pressure used was 1,400 psia. To facilitate comparison of mercury-injection and water-displacement results, it is convenient to express mercury-injection data in terms of unsaturated pore volume. The abbreviation, (100 - SHg)min, is used to denote the pore volume unsaturated by mercury at 1,400 psia injection pressure. Gas permeabilities, ki, were deter-mined by the method proposed by Klinkenberg.8 Methods of determining types and amounts of clay minerals by X-ray diffraction analyses, as well as methods of determining water permeability, have been described previously.' The imbibition characteristics of

Jan 1, 1958

By R. S. Shoemaker, D. L. Harris

PST investigations into the possibility of concentrating the low-grade iron ores of the Gogebic Range in Wisconsin have been hampered by the complex association of the constituent minerals. In part the problem arises from the fact that the iron occurs essentially as two minerals, hematite and goethite, having nearly the same chemical composition but exhibiting different physical properties. Individual grains of these minerals are usually so small that the ore must be ground to a size finer than 325 mesh if substantial liberation of the minerals is to be secured. Appraisal of beneficiation tests on this iron ore requires, therefore, some method of determining quantitatively the contents of constituent minerals in the various products of the beneficiation. Because of the chemical similarity of the two iron minerals, chemical analyses do not provide a means of differentiation. Although the minerals are different physically, they are both opaque and friable, and as the ore must be ground to a very fine size for liberation to be secured, quantitative microscopic methods of analysis are extremely difficult at best. In an effort to overcome these analytical problems, a quantitative X-ray diffraction analysis method has been developed at the University of Wisconsin for the research program on iron ore beneficiation. The X-ray diffraction analysis of chemical compounds of crystalline form has proved of inestimable value for structure studies and rapid qualitative analytical purposes in a number of fields.' In some cases X-ray diffraction analysis has been found suitable for quantitative determinations. H. P. Klug2 found that quantitative procedures for quartz powders when mixed with calcium fluoride were reproducible within ±1.0 pct. Talvenheimo and White3 reported quantitative analyses of clay minerals to be accurate within 10 pct and also detected 1 pct of bentonite in 99 pct illite. In other investigations,'" powders of quartz and mica have been analyzed quantitatively with accuracies from 1 to 5 pct. A powder of tungsten carbide was analyzed similarly by Redmond7 and Rooksby8 reported the determination of small amounts (0.11 pct) of calcium oxide in magnesium oxide. In the most generally practiced use of X-ray diffraction, traces of the diffracted beam are recorded as lines on a film much as lines of elements are recorded on an optical spectrograph film. As with the spectrograph, the darkness of the lines denotes the amount of the constituent in the sample. If the film is replaced by a Geiger tube assembly and a graphic recorder or a scaling unit, the intensity of the diffracted beam can be registered as a peak on a chart or energy impulses per unit of time on the scaler. Several methods of quantitative measurement have been investigated and reported in literature. Klug2 used manual operation of the goniometer and a counting technique to determine the difference in counting rates at the diffracted angle and for the background. In another case3 the area under the recorded diffracted peak was the basis for the determination. In other instances the height of the recorded diffraction peak above background, measured while the goniometer was moving, has been the basis for the measurement. The sizes of particles most suitable for diffraction techniques are reported in the literature as follows: for tungsten carbide a maximum of 40 microns7 and for quartz a maximum of 5 microns.' Very fine sizes of crystals are known to cause a broadening of the diffused beams; hence some workers have set a lower limit of 0.1 micron for quantitative results. Guinier10 has indicated that broadening of the diffracted X-ray beam is not appreciable for crystallites larger than 0.02 to 0.03

Jan 1, 1956

By Herbert H. Kellogg

Graphs representing the standard free-energy of formation as a function of temperature for 21 fluorides are presented, along with estimated values for the standard free-energy of formation of 20 atedadditional fluorides. A few of the many possible uses of these data in metallurgical calculations are discussed, including the fluorination of oxides, sulphides, metallurgicaland chlorides, and the reduction of metal fluorides. RECENT papers which have presented free-energy data for oxides, sulphides, chlorides, carbonates, sulphates, and silicates in the form of free energy vs. temperature plots,1-5 have found important uses in teaching and metallurgical development work. This method of presentation of thermo-dynamic data serves to make available in readily usable form thermodynamic constants which are scattered throughout the literature, and to show by graphical means the stability relations which exist between the compounds of a given type. It is certainly desirable to have these free-energy-temperature diagrams for as many types of compounds as are of present or potential use in metallurgical processes. Thus, even though there are but few uses of metallic fluorides in extractive metallurgy today, the potential uses of fluorine processes are of great importance and serve as the principal justification of this paper. Recent developments in the production of fluoro-carbons have resulted in vastly improved methods for the production and handling of fluorine and corrosive fluorides, with consequent reduction in costs.".' The extractive metallurgist should no longer overlook possible fluorine processes when development of a new process for a relatively high unit-price metal is being considered. In a previous paperv he author outlined in some detail the importance and use of free-energy data in metallurgical calculations and also discussed the properties of a standard free-energy vs. temperature diagram. The same method of presentation is adopted in this paper, and the reader is referred to the earlier work for explanations of the use of free-energy data. Thermodynamic Data This paper presents standard free-energy vs. temperature diagrams, Figs. 1 and 2, for those fluorides for which the major thermodynamic constants are known. In the calculations of the standard free-energy of formation (AGO), values of the enthalpy of formation at 298°K (?Ho,) were taken from "Selected Values of Thermodynamic Constants"8 most of the values of standard entropy at 298°K (S°298) were taken from "Selected Values"," but in a few cases the values given by Kelley9 ere preferred; data for the Phase changes came from "selected Values" and several of Kellley's publics- tions;0-12 and most of the specific-heat data were taken from Kelley.l0 The appendix gives a detailed list of the sources of the data. As explained in an earlier paper,6 he accuracy of some of the standard free-energy relationships is not high. An error of +1 kcal at 1000°C is expected in most of the curves of Figs. 1 and 2. The error may be as high as +3 kcal in some of the curves. Rather than attempt to estimate the reliability of each curve, the source of the data used is reported in the appendix, and the reader may compare these data with newer values as they become available. Figs. 1 and 2 fail to show values for the fluorides of many important metals. In most cases this is because the entropy of the metallic fluoride is unknown. For preliminary calculations it is often useful to have available the best estimate of an unknown datum, rather than no information at all. With this in mind, Table I lists estimated standard free-energy values for a number of important metallic fluorides not included in Figs. 1 and 2. The values reported in Table I were, calculated by Brewer et al.13 and are based on estimated entropy values. The estimated uncertainty of the values for each fluoride were calculated from the uncertainty values given by Brewer1' for ?H°, and the entropy function.

Jan 1, 1952

By W. M. Rutherford

A knowledge of the limits of miscibility between reservoir oil and possible injection fluids is required for selection of the optimum miscible-injection fluid. Limits of miscibility can be estimated from the results of equilibrium phase-behavior experiments. They can also be determined by means of displacement experiments conducted in a high-pressure sandpack. This paper describes the equipment and procedure which have been developed for determining miscibility conditions by stable displacement. A systematic series of displacements of a West Texas reservoir oil was carried out. The results indicate that, at constant pressure, miscibility is a function only of the pseudocritical temperature of the injection gas. This fact, together with improved experimental methods, makes the displacement technique a rapid, reliable means for determining miscibility conditions. In conjunction with the displacement experiments, phase diagrams were constructed for the oil with dry gas and propane and with dry gas and ethane. Phase behavior of the methane-ethane-propane system was determined at 110°F. The experimental work demonstrates the feasibility of using ethane-rich gases to reduce cost and pressure requirements for miscible displacement. INTRODUCTION In recent years, interest in the miscible displacement of oil by light hydrocarbon mixtures has been high. Many pilot and a few field scale projects have been started.1 These projects have made use of various methods for achieving miscibility: (1) the LPG-slug process, (2) the enriched-gas-drive process and (3) the high-pressure gas-drive process. Some field projects have been successful; the results of others are debatable. In general, projects which have performed best have involved the injection of an appreciable fraction of a pore volume of miscible material. Economical application of miscible displacement depends strongly on the cost of the miscible-injection fluid. If an appreciable fraction of a pore volume of material is required for successful application of these methods, a precise knowledge of the minimum requirements for miscibility in terms of composition of injection fluid is essential. Therefore, reliable experimental methods for determining miscibility conditions are important, and a procedure for estimating these conditions from the composition of the reservoir fluid is highly desirable. The subject of this paper is the problem of determining conditions which result in miscible displacement of oil by light hydrocarbon mixtures. Miscibility conditions can be estimated by means of equilibrium experiments conducted in a PVT cell, or they can be determined by means of high-pressure displacement experiments. This paper describes the equipment and procedure which have been developed for the determination of miscibility by high-pressure displacement experiments. These methods have been applied to the displacement of a West Texas reservoir oil with mixtures of dry gas, ethane and propane. In conjunction with the displacement experiments, triangular phase diagrams have been constructed for mixtures of the oil with dry gas and propane and with dry gas and ethane. The effect of injection-gas composition on conditions for miscible displacement in high-pressure sandpacks and cores has been the subject of several published papers.2-5 The experimental methods used in these investigations were such that displacements were unstable, and the effects of fingering and/or gravity layover are clearly evident in the results. Miscibility conditions were probably correct in spite of the instability phenomena, but the experiments evidently were time-consuming, and limited data were reported. Systematic high-pressure flow studies which would support a correlation of miscibility conditions have not been reported; however, Wilson5 has proposed the use of pseudocritical temperature of the injection gas as a parameter and Benham, et al,6 have based a miscibility correlation on observed and calculated equilibrium data.

By L. R. Mabson, F. M. Russell

INRODUCTION Roan Consolidated Mines Limited (RCM) operates five underground mines - Mufulira, Luanshya, Baluba, Chambeshi and Chibuluma (figure 1) . Associated metallurgical plants provide for the operations of milling, smelting two refineries for copper and a cobalt plant. Primarily a producer of high quality copper (242,000 tonnes during 1980) the Company is also the world's second largest producer of cobalt metal. RCM's remaining ore reserves at 31st March 1980 were 275 million tonnes at 2.84% copper. Monthly ore production is in excess of 1.2 million tonnes. Sublevel open stoping was the choice of mining method for the first underground operations on the Copperbelt for ore with dips steeper than 30 degrees. This was for use at the Roan Antelope mine, subsequently to be called Luanshya mine. The sublevel open stoping system was planned based on the methods used in many of the iron mines of Michigan and New York and in the copper mines of Tennessee. Peterson (Ref. 1) gave the following reasons for employing sublevel stoping, reasons which are still valid today: (1) the working place could always be barred and made safe, and the sublevel exit was always at hand; (2) the working conditions were similar day by day and drilling operations were repetitive; (3) pillar spacing can be varied as conditions demand; and (4) stoping operations could be shut down completely or opened up to full capacity at short notice, without affecting grade or causing wastage of reserves. Such was the success of the method that it was adopted as the basic mining method at all the RCM mines in their initial stages. Although cut and fill, block caving, and sub- level caving methods have been used, sublevel open stoping in its many variants still accounts for nearly 90% of total ore production in the Company. The following descriptions outline some of the applications of sublevel open stoping within RCM. LUANSHYA Geology There are two main ore deposits at Luanshya: the Roan (Luanshya Mining Area), in which mining began in 1930, and Baluba, in which mining production started in 1973 (figure 2). Both orebodies occur in synclinal and variably folded basins. These basins are underlain and completely surrounded by much older schists and granites. The ore formation, an argillite, lies near the base of the sediments and is underlain by conglomerates, quartzites and arkoses. Copper mineralisation is mainly found in two well defined zones, separated by a varying thickness of pyrite mineralisation. In the eastern part of the Roan deposit only the upper orebody is present, but west of 18 shaft and at Baluba the lower orebody predominates. Locally the subjacent quartzites of the footwall formation are also mineralised with copper. The copper mineralisation is found princi- pally as disseminated sulphides, mainly chalcocite, bornite and chalcopyrite. Near the surface the sulphides have been changed to oxides. At Baluba the cobalt mineral carrollite is found with the sulphides. Mining methods Sublevel open stoping is the basic mining method. In its standard form this method is best suited to steeply inclined orebodies, with simple structures and with a strong foot- wall and hangingwall. In recent years, as more complex folding and flatter dipping orebodies have been encountered, considerable variations have been introduced. Stope drilling is effected with drifters and extension steel.

Jan 1, 1981

By A. M. Adair

The effect of cathodically charging zone-refined iron with hydrogen and the resultant yield stress-grain size relationship has been investigated by tensile testing at -31°C. The resulting Petch plots (using a least-squares analysis) reveal a significant reduction in the slope and an increase in the intercept for specimens cathodically charged with hydrogen, regardless of whether the hydrogen is present during testing or is removed by evacuation. Other than surface blistering, the major effect of cathodically charging this material is to create occasional grain boundary cracks and thus presumably a large number of free dislocations at grain boundaries. These results are consistent with the hypothesis that the slope of the Petch plot is related to the stress required to activate dislocation sources at or near pain boundaries. While there is no significant effect of the presence of hydrogen on the lower yield stress of this material, the ductility is considerably reduced. After hydrogen removal by evacuation, the ductility recovers substantially and almost complete recovery of ductility is experienced by the smaller grain size specimens. THE influence of hydrogen on the mechanical properties of ferritic materials has received much attention throughout the history of the metallurgy profession. This interest has stemmed from the obvious importance of preventing the catastrophic failure of parts in service, and from the desire to determine what mechanism is responsible for the observed behavior. The various mechanisms proposed for hydrogen em-brittlement have been discussed in detail by Tetel-man,'who concludes that a combination of the applied stress and hydrogen pressure within voids best accounts for the many experimental observations reported. From a quite reasonable correlation of the strain aging and delayed-failure behavior of a variety of materials,2-7 Troiano and coworkers are of the opinion that the delayed-failure mechanism (slow strain rate embrittlement as well) is better explained as a lattice effect; solute elements diffuse to regions of high triaxial stress, form uncondensed atmospheres there, and effect a reduction in the cohesive strength of the lattice. Since hydrogen yield points8,9 have been observed in ferritic materials, and since the hydrogen embrittlement process may be related to strain aging, one further contribution to embrittlement could be hydrogen-dislocation interactions which make crack blunting by plastic deformation a more difficult process. The etch" relationship offers a means of comparing the relative ease of initiating and propagating plastic deformation in metals and alloys. In the present study, this relationship was used to compare the behavior of zone-refined iron in three conditions: uncharged, cathodically charged with hydrogen, and cathodically charged with hydrogen and evacuated. These Petch plots and the observed ductility parameters allowed a separation of the effects of the very presence of hydrogen and the act of cathodically charging hydrogen into this material. MATERIALS AND EXPERIMENTAL PROCEDURE An ingot of zone-refined iron was obtained from Battelle Memorial Institute which contained the following (typical) interstitial impurities in ppm: carbon, 10 * 5; nitrogen, 1; oxygen, 22; hydrogen, 0.9. The ingot was cold-reduced to a 1/2-in.-sq cross section, recrystallized 1 hr at 575°C, and then straight-rolled to a thickness of about 0.105 in. Tensile specimens were milled from the strip having a 1-in. gage length

Jan 1, 1967

By S. A. Spachner

An analysis of the powder patterns of the omega transition phase in three binary titanium alloys containing 8 pct Cr, 13 pet Mo, and 15 pet V, respectively, showed that all could be indexed on an orthorhombic unit cell. The unit-cell dimensions were all quite similar. There is evidence for the occurrence of more than one transition phase in Ti-Cr and Ti-V alloys. FOLLOWING the discovery of a new brittle phase in certain titanium-chromium alloys by Frost, Parris, Hirsch, Doig, and Schwartz,' named ? by these investigators, much interest has developed concerning the identification of alloy systems in which this phase exists, kinetics of the transformation in such systems, and the crystal structure of the ? phases appearing. Omega has been found in alloys of titanium-manganese, titanium-molybdenum, titanium-iron,' and titanium-vanadium.' In the present work, a study was made of the powder patterns of the ? phase occurring in alloys of Ti-8Cr, Ti-13Mo and Ti-15V in an effort to obtain a basis of comparison for ? phases appearing in these different alloys. Comparison of Alloys for ?-Phases Two hundred mesh powders of Ti-8Cr, Ti-13Mo and Ti-15V were prepared. A portion of powder sufficient for three x-ray powder specimens was placed in a Vycor tube, evacuated to less than 0.1 µ pressure, filled with argon, and sealed. The tubes were then heated to 900°C, held at this temperature for 30 min, and water-quenched. Following this transformation to the ß phase, the specimens were quenched in lead or solder pots, held at the desired temperature for varying lengths of time and water-quenched. The heat-treatments used and the phases identified are listed in Table I. The intermediate water quench was used in lieu of direct downquenching from 900°C to the bath temperature desired to circumvent the possibility of transformation of the powder specimen from ß to (ß + a) while it was being moved from the homog-enization furnace to the lead bath. To check the possibility that upquenching to the bath temperature might bring into play a different mechanism than downquenching to this temperature, resistivity measurements were compared on Ti-8Cr rod specimens which were transformed isothermally at 400°C by downquenching from the solution-treatment temperature and those which were transformed at 400" C by upquenching ß specimens from room temperature. The results are shown in Fig. 1. No difference in transformation mechanism was noted. It soon became apparent that sufficiently detailed Debye-Scherrer photograms of the ? phase in the alloys which were to be used could not be obtained by the use of copper radiation. Although satisfactory photograms of ß and ß + a phases were obtained with this radiation, the ? phase appeared to exhibit mild fluorescence, making interplanar spacing computation virtually impossible. Molybdenum radiation appeared to offer promise in cutting down this fluorescent effect, the mass absorption of titanium for molybdenum radiation being 23.7 in contrast to 204 for copper radiation. The reduction in the accuracy in determining the location of the photo-gram lines caused by the shorter wave length of molybdenum characteristic radiation did not appear to be a significant factor. Utilizing a camera of 57.3 mm radius, best results were obtained for exposure times ranging between 30 and 40 hr. Tube potential

Jan 1, 1959

By Hsun Hu

Direct observations were made on the recrystalli-zntion of a cold-rolled (110) [001] crystal of 3 pct Si-Fe by annealing thin sections of the specimen inside the electron microscope during examination. Subgrain formation, presumably by polygonization, occurs rapidly, followed by a slow increase in the subgrain size before recrystallization can be detected. No migration of the subgrain boundaries was observed, and the formation of recrystallization nuclei by coalescence of subgrains is suggested. Re-c~ystallization characteristics of various texture components were individually examined. AS reported in a previous paper: orientations of the various deformation texture components in a 70 pct cold-rolled (110) [0011 crystal of 3 pct Si-Fe were completely confirmed by transmission electron diffraction patterns taken at various areas of the specimen. These orientations can be summarized as follows: Over most areas of the deformed crystal (111) planes are parallel to the surface. These areas correspond to either $- the two main defor: mation textures, the (111) [112] or the (111) [112].2-5 Within a few areas, (100) planes parallel to the surface are frequently found. These areas represent the minor deformation texture component, the (100) [Oll] type, which is derived from rotated Neumann bands.2'4'5 Only in very rare instances are (110) planes found parallel to the surface. This orientation corresponds to the very weak texture component of the (110) [001] type. In a pole figure, this (110) [0011 orientation is located within the orientation spread between the two main (111) [1121 type textures. Usually the electron diffraction patterns taken at various areas of the specimen are very simple ones, i.e., for any one area the pattern consists only of sharply defined spots of a single orientation. In some areas, however, more or less continuous rings of weaker intensities are superimposed on a simple pattern, indicating the presence of a wide range of additional orientations. In order to understand in detail the recrystallization of this cold-worked crystal, the recrywtallization characteristics of each of the various deformation texture components must be examined. This can best be done by using a high-temperature stage and transmission electron microscopy techniques, whereby small areas corresponding to the various texture components can be selected individually for examination. The present paper presents results obtained by this method. MATERIAL AND METHOD Specimens were prepared from a previously studied1 single crystal of a commercial grade 3 pct Si-Fe alloy with an initial orientation (110) [00l], cold-rolled 70 pct from 0.014 in. to 0.004 in. in thickness. Thin foils were prepared by electrolytic polishing as described in the earlier paper.' Specimens of this thin foil were then mounted in a high temperature stage similar to the one designed by Whelan.6 A direct current was used as the heating source. The temperature of the heating grid was calibrated by means of a Pt vs Pt-10 pct Rh thermocouple in a vacuum chamber (0.01 to 0.02µHg). Maximum variation within the temperature range of 250o to 000°C was about ± 25°C. After the desired areas were selected for examination, the specimen was heated to a low temperature for some time while the selected area was constantly watched. When no further change in structure and orientation could be observed, a higher temperature was applied. Photomicrographs and diffraction patterns were taken from time to time during examination, with the heater temporarily turned off* while the photographs were

Jan 1, 1962

By L. G. Jones

An approximate method of calculation is developed in this paper which allows duplication of radial unsteady-state gas flow computer results where Darcy's law applies, such as those reported by Aronofsky and Jenkins and Bruce, Peaceman, Rachford, and Rice. Moreover, the new calculation method can be used to obtain results for radial unsteady-state gas flow obeying the quadratic flow law proposed by Duwez and Green. Means are discussed for predicting well behavior at single or superimposed flow rates in finite or infinite reservoirs, determining reservoir rock properties from well-test data, reproducing and interpreting back-pressure test data, and determining the radial extent and reserves of gas reservoirs from well-test data. Example calculations are presented for gas flow following both Darcy's law and the quadratic flow law. INTRODUCTION Since the publishing of U. S. Bureau of Mines Monograph 7,1 most gas-well testing methods have been based on the equation q = y (pl2 - pw2)a, where q = production rate, pl and pw, are formation pressure and sandface pressure, respectively, and y and a are constants to be obtained from test data. These methods, used for predicting both deliverability and "open-flow" capacity of gas wells, have been useful and accurate in many cases but unsatisfactory in others. Even at best, however, they do not supply information about the formation or lead to an understanding of nonsteady-state gas flow in porous media. Many theoretically based studies2-4 of gas flow obeying Darcy's law have been made.' Since the partial differential equations which result from combining Darcy's flow law with the continuity equation are nonlinear, most of the published research consists of either numerical solutions or analytical solutions for linear approximating equations. Such solutions have been of limited value in field work due to their unhandy form and their failure to correlate most field data. There is evidence5,6 which indicates that Darcy's law is inadequate to describe gas flow at some flow rates of practical interest. A quadratic flow law: which reduces to Darcy's law at low rates, is more successful in accounting for experimentally observed behavior. This flow law has been applied successfully to a few hypothetical reservoir cases in work which has recently been published.10 However, these numerical solutions of the equations involved have been successful only on a special analog computer. Routine applications to field cases would be awkward and have not been attempted. The present paper describes an approximate method for computing nonsteady-state gas flow solutions which has been completely successful in predicting the results for both Darcy flow and quadratic flow obtained by elaborate numerical methods. The new calculation method allows determination of the observable variables in gas-well testing at constant rates. It is similar to the scheme of using a succession of steady states suggested by Muskat7 in that it makes use of steady-state and material-balance equations. It also is similar to the work of Aronofsky and Jenkins2 in that the new method includes Aronofsky and Jenkins' stabilized flow equation as a special case. It improves upon both of these calculation schemes in that it accurately describes all portions of reservoir history and suggests means of determining reservoir rock properties from well-test data. This paper deals only with production from a reservoir, in which case the rate is defined as being negative. The reservoir studied here is a homogeneous, disk-shaped porous body of uniform thickness, with all boundaries sealed except the inner radial boundary of the well.

By John Johnston

WHEN I was honored by being invited to give the Howe Memorial Lecture, I decided to read Howe's book, "The Metallography of Steel and Cast Iron," published in 1916-that is, about 25 years ago-in search of a text. I found the book written in fine, clear, attractive, yet precise style, as I expected, knowing that his father ¬Dr. Samuel Gridley Howe, his mother-Julia Ward Howe, and all of his four brothers and sisters had lived in a New England literary atmosphere and had written books. This leads me to suggest that many authors of metallurgical papers and books could well study Howe's writings with a view to the improvement of their own. I failed to find a specific text, but did find a subject, for in reading his book I was struck by the few instances in which time is considered explicitly as a factor in the phenomena examined, and by the fact that many of the hottest arguments would have been settled had time been taken properly into account. Indeed, Howe hardly mentions time except in a rather incidental way, as in a few references to lag; for instance: The lowering of temperature at which the transformation occurs by rapid cooling . . . is a phenomenon of undercooling or lag. It is in part explained by the fact that the transformation itself is a time-consuming reaction which, setting in during a rapid fall of temperature. is naturally protracted till the temperature has sunk far below that at which it is due. He does discuss, quite fully, the iron-carbon equilibrium diagram, as known at that time, introducing it by stating: As he is a reasonable man who is not deterred from using a hand truck for moving his trunk by the consideration that he has thus to move a trunk plus a truck, so am I reasonable in putting into your hands this useful, indeed indispensable, tool for mastering the ABC of iron metallurgy. Somewhat later, in introducing a discussion of the phase rule, which "it would be hardly proper to pass . . . by without an attempt to outline its meaning," for "a clear exposition of its applicability or jurisdiction is needed to restrain the half-initiated from misusing it and the uninitiated from being misled by these misuses," he writes-and I quote his statement because we cannot properly discuss the influence of time on a process until we know the state of equilibrium which would ultimately be reached by the system under consideration It is a most remarkable and valuable generalization, its conceptions help greatly toward getting a broad outlook on metallography, but its misconception has brought out a flood of obscuring writings. It tells us about the constitution towards which alloys tend, that which they reach when equilibrium is complete, when all tendencies have corrected themselves and have been complied with completely. Its application, in absence of true equilibrium, is like determining density with a telescope or even with a tuning fork; it is like the mad tea party at which the butter put into the works "f the watch was the best butter. If the reader learns no more than to beware of attempts to decipher the conditions of inequilibrium by a law which touches only the conditions of exact

Jan 1, 1942

By Arthur C. Green

A group of 13 cities situated in Los Angeles and Orange counties in Southern California is engaged in constructing an aqueduct to carry water from the Colorado River at a point near Parker, Arizona, to a distribution point near Riverside, California. The aqueduct will have a total length of approximately 241 miles, and is designed to deliver 1605 cu. ft. of water per second, or a little more than one billion gallons per day. The aqueduct crosses a barren, desert waste and the project includes approximately 91 miles of 16-ft. tunnel. The large amount of tunnel driving that has been done in this country in the past few years has resulted in some remarkable developments in equipment and methods. The building of the Colorado River Aqueduct,, with its 91 miles of tunnels, has called on the best engineering knowledge in the country. New and better equipment has been designed for all phases of the work, and this equipment has been coordinated and adapted to the various conditions with a high degree of efficiency. The result has been that records for tunnel driving have been set under adverse natural conditions that exceed the estimates for progress under good conditions. All of the tunnels in the main aqueduct, except the Valverde and Bernasconi at the western end, have a finished diameter, inside the concrete lining, of 16 ft. Fig. 4 shows a section of the tunnel with the concrete lining in place. Fig. 5 illustrates the finished tunnel in three typical sections. The standard tunnel grade is 3.432 ft. per mile in the 16-ft. tunnels, but the Valverde and Bernasconi tunnels have a steeper grade and a finished diameter of 15 ft. 3 in. The rough bore of the large main aqueduct tunnels averages 12.95 solid cubic yards of material per foot of advance in timbered sections and 10.76 solid cubic yards per foot in unsupported sections. Tunnel Driving While the Conway mucker has been used on tunnel work for years, a larger and more powerful Conway was designed to meet the require-

Jan 1, 1937

By J. Marin, H. A. B. Wiseman

This paper presents results of an experimental study dealing with the plastic stress-strain relations of aluminum alloy 14s-T6 subjected to combined biaxial tension and compression stresses. Plastic stress-strain relations for both constant and various types of variable stress ratios were determined and a comparison made with the simple flow theory of plasticity. THIS investigation was undertaken for two main purposes: 1—to obtain plastic stress-strain relations for Alcoa 14s-T6 when subjected to various combinations of biaxial stresses, and 2—to determine the validity of the flow theory in predicting plastic stress-strain relations for combined stresses. The biaxial stresses used were tension combined with compression. These stresses were produced by subjecting a thin-walled tube to axial tension and torsion. Various types of combined stress conditions were investigated. To provide control test data and information on the influence of biaxial stresses on the strength, the usual constant stress ratio type tests were made. These tests showed that the biaxial yield strength agrees approximately with the distortion energy theory. However, the Prager semi-empirical theory agrees best with the test results. For these constant stress ratio tests the plastic stress-strain relations were found to agree approximately with the flow theory. (For constant stress ratio the flow and deformation type theories are identical.) In view of the fact that the constant stress ratio tests cannot distinguish between the flow and deformation theories, a number of variable and special biaxial stress tests were made. These tests showed that the flow theory is inadequate, since large differences exist between the experimental and theoretical results. Special tests were also conducted to check the validity of the isotropic yielding assumption. This assumption is made in the linear-type flow theory. In these tests specimens were loaded in tension to predetermined values beyond the proportional limit stress. The tension load was removed from the specimen and the specimen was subsequently loaded in torsion to fracture. Other tests were applied with the order of torque and tension loading reversed. The nominal stress-strain curves for these two tests approximately coincide indicating that the assumption of isotropic yielding is valid. In addition, these tests verified the requirement by the slip theory of plasticity that prestraining below 140 pct of the proportioned limit stress in either torsion or tension did not influence the subsequent plastic stress-strain relations in tension or torsion. Another type of test was conducted to check the assumption that the axes of principal stress and strain coincide during plastic flow and variable stress conditions. These tests were conducted by applying combinations of axial tension, internal pressure, and torsional loads such that the axes of maximum principal stress could be rotated through 90" during the loading path. The results of those tests showed that the angle between the axes of principal stress and strain varied greatly and that the assumption of coincidence between axes of principal stress and strain is not valid. Introduction In various metal processes including forming of sheets, rolling of bars, and extrusion of rods, metals are subjected to stresses beyond the yield strength of the material. Often these stresses are not simple stresses acting in one direction, but are combined stresses acting in more than one direction. Structural and machine members are often subjected to combined stresses. To adequately determine the factors of safety for these members, it is necessary to know the plastic stress-strain relations of the material used for various combined stress conditions. Various theories have been proposed for predicting the plastic stress-strain relations for combined stress in terms of the simple tension plastic stress-strain relation for the material. To determine which theory might be adequate, test results must be obtained for various combined stress conditions, in order to compare the actual with the assumed theoretical behavior. Although there is considerable test data available for combined stress conditions in which the principal stress ratios are constant, relatively little test information is available for vari-

Jan 1, 1954

By R. R. Angel

Drilling rate is a parameter that should be considered in determining the volume requirements for air and gas drilling. The use of past methods which ignore the effects of the solids content upon the pressure and velocity of the annulus flow stream can result in undercalculation of the required volurne by as much as 50 per cent. A vertical-flow equation is presented for determining volume requirements. This equation includes the effect of the solids that are transported up the annulus in the flow stream by incorporating the drilling rate as one of the parameters. The effect of down-hole temperature on required circulation rates is also ana-lvzed. A simple approximate method of determining volume requirements is presented. This method is more accurate than the methods used in the past. Hole cleaning difficulties are analyzed for a recent air drilling job where past methods indicated that - excess air was being used. Sample curves of calculated boitom-hole pressures are presented for air and gas drilling in several hole sizes. INTRODUCTION In certain areas the use of air or natural gas as a circulating medium for drilling oil and gas wells is becoming a common practice. Large increases in penetration rate and bit life are achieved through the use of these media in preference to mud or water. Drilling rates as high as 90 ft/hr have been obtained in shales. The importance of maintaining adequate circulation rate is generally recognized; however, much disagreement exists among drilling operators as to what constitutes "adequate" circulation rate. In quarry drilling, where annular velocities can be accurately determined, an annular velocity of 3,000 ft/min of standard air is required for best results in rocks having approximately the same density as those commonly penetrated in drilling oil or gas wells'. Although this standard air velocity has proven satisfactory for quarry drilling, some oil and gas well drilling operators believe that an equivalent annular velocity of more than 4,000 ft/min is required; others believe that as little as 2,000 ft/min . is sufficient. Much of this disagreement results from determining the required circulation rates with methods which fail to incorporate the drilled solids in an equation which is applicable to vertical flow. Hughes Tool Co. Bulletins No. 23' and 23-A3 present data for determining circulation rates based on the Weymouth formula. These data do not include the drilling rate as a parameter and, therefore, neglect the effect of the solids being transported up the annulus. In spite of this apparent defect and the fact that the Weymouth formula is not applicable to vertical flow, the Hughes data have well served the drilling industry in many areas and are important and timely contributions to the science of air and gas drilling. The Hughes data purposely omit a correction for increasing down-hole temperature. At slow drilling rates this effectively compensates for the use of a formula which is not valid for vertical flow; however, volumes determined by the Hughes method are not sufficient to support rapid drilling rates at moderate and great depths. For example, Phillips Petroleum Co.'s Cauthorn "D" No. l in the Vinegarone field for Val Verde County, Tex., was air drilled from 1,500 to 9,300 ft using a compressor delivering 1,400 cu ft/min. The 8 3/4-in hole was drilled with 5-in. drill pipe, and drilling rates as high as 90 ft/hr were obtained between 7,000 and 9,300 ft. No water or caving hole was encountered. At 7,728 ft it was necessary to wash-out 60 ft of cuttings to reach bottom after a trip to change bits. At 8,130 ft a twist-off occurred and the drill col-lars were stuck in drill cuttings. These difficulties indicate that the 1,400 cu ft/min of air was not sufficient to keep the hole clean. Hughes data indicates that 1,180 cu ft/min is sufficient to produce an annular ve-

Jan 1, 1958

By P. E. Busby, C. Wells, M. E. Warga

Fundamental data on the rate of diffusion of boron in austenite and solubility of boron in the a and y phases of iron and steel have been obtained from deboronizing experiments and provide partial explanations for some of the phenomena observed in boron steels. The rate of diffusion of boron is about the same as carbon in austenite. The solubility of boron in austenite at normal heat-treating temperatures is less than 0.001 pct. A partial tentative Fe-B phase diagram in the important low boron concentration ranges and an equation representing the diffusion of boron in austenite are presented. A LTHOUGH a search of the literature revealed ALTHOUGHthat boron diffusion had not been studied quantitatively and systematically prior to 1948, a few qualitative observations which were made by following the diffusion of boron into iron packed in ferro-boron1-3 suggest that boron obeys normal diffusion laws and that the rate of diffusion increases with increasing temperature. The claim of Cornelius and Bollenrath1 that boron does not appear to influence the rate of carbon diffusion has recently been substantiated by Wells, Batz, and Mehl.4 Calculations based on data from the paper by Campbell and Fayv indicated that the diffusion coefficient (D) for boron in iron at 900°C is approximately 3x10 sq cm per sec; a value of 2x10-7 sq cm per sec at 1038°C has been reported by Digges et al in connection with decarburizing experiments on a commercial boron steel containing 0.43 pct C. During the period 1948 to 1950 several diffusivity constant (D) values for the diffusion of boron in the y phase of iron and steel and a number of solubility values for boron in both a and y phases were reported. It was concluded tentatively as a result of Metals Research Laboratory studies that D values for boron in y and a phases are about the same as for carbon in the comparable phasesa and that carbon up to 0.4 pct did not affect the rate of diffusion of boron in y iron above 1000°C within the limits of experimental error. Whether saturation of the y or a phases with carbon would significantly affect rates of diffusion of boron through them was then and still remains in doubt. The best estimate of Q (activation energy) was reported to be about 25,000 cal, somewhat lower than that for carbon in y iron. The solubility of boron in y iron at 1000°C was thought to be higher than 0.004 pct B but not as high as a more recent analysis of available data shows it to be. The solubility of boron in a iron at 700°C was reported to be 0.0004 pct B, but it would not have surprised the authors if the true solubility value is actually lower than this. Breaks in diffusion curves which were not understood when first obtained prior to 1950 have now been recognized as indicating solubility limits, and these are included in the present paper. Experimental Procedure Two methods were employed in an effort to determine diffusion coefficients (D) for boron in austenite. Early experiments involved the use of welded couples prepared in accordance with the procedure of Wells and Mehl,' but when the results of these experiments proved to be unsatisfactory for calculating D values, studies were continued by means of deboronizing experiments. In both cases, boron analyses following the diffusion anneal were carried out spectrographically and all D values were computed using the Grube solution of Fick's law. The precision of the spectrographic method used proved to be about equal to that reported by Corliss and Scribner" —average deviation from the mean about 5 pct of the amount present when the boron content is 0.003 pct. Concentration-penetration curves obtained from spectrographic analyses of welded couples generally indicated that practically no boron was transported across the weld interface during the diffusion anneal. Furthermore the concentration-penetration

Jan 1, 1954

By G. M. Schwartz, G. A. Thiel

THE first record of the quarrying of dimension stone in Minnesota dates back to 1820 when limestone was quarried locally for part of old Fort Snelling. Limestone quarries were operated at Stillwater, Mankato, and Winona as early as 1854. Granite was quarried first at St. Cloud in 1868, and within a few years thousands of tons were shipped to widespread points. Rough dimension stone for large buildings furnished the first important market, but beginning in 1886 paving blocks were in demand. The largest shipment was in 1888, when 1925 cars were shipped from the St. Cloud area. Quartzite was quarried first at New Ulm in 1859 and somewhat later at Pipestone and elsewhere in southwestern Minnesota: The productive dolomite quarries at Kasota were opened first in 1868 and have continued as large producers of a variety of stone to the present time. At present, the industry is controlled by relatively few operators, and for that reason detailed figures on dimension stone are not released for publication. A general idea may be obtained from the data in the Minerals Yearbook for 1948. The figures for total stone produced in Minnesota are 1,804,000 tons valued at $5,090,652. Probably the largest item in the latter figure is received from dimension stone. A better idea of the situation, in relation to the country as a whole may be gained by using the data for 1930 when more companies were operating in Minnesota, and complete figures were published. In that year Minnesota produced granite valued at $2,668,119 and ranked third among the states in value. Minnesota's production of granite was almost exclusively for dimension stone. In the same year Minnesota produced 300,000 tons of limestone (dolomite) valued at $840,860, and this likewise was mainly dimension stone. In finished limestone Minnesota ranked second among the states in 1930. Sandstone and minor amounts of quartzite are the only other dimension stones that have been produced in Minnesota, but the quarries are now inactive. The commercial stones of Minnesota have been described in two reports by Bowled and by Thiel and Dutton.2 The early history of quarrying in Minnesota and extensive notes on the various rocks are given by N. H. Winchell.3 Small limestone and dolomite quarries were numerous throughout the area of Paleozoic rocks in southeastern Minnesota. Early production was largely dimension stone. With the increased use of Portland cement, most of these ceased production, and today only those at Kasota and Winona remain in operation. In recent years many quarries have reopened and new ones started, but these are devoted to the production of crushed rock and agricultural lime. As the application of modern quarrying and finishing methods increased, small companies in the granite business have dropped out, and the remaining companies have modernized their plants, purchased old quarries, and opened up new ones, thus furnishing a wide variety of granites suitable for most of the customary uses. It is the purpose of this review to present notes on' the geology and operations of each of the quarries now operating within the state. Granites and Related Igneous Rocks The term granite as used in this report includes granites, gneisses, diorites, gabbros, and other igneous rocks. The granites of greatest economic importance are found in three widely separated regions, see Fig. 1. 1-Central Minnesota in the region of the city of St. Cloud, 2-the upper Minnesota River valley region, 3-the northeastern portion of the state, commonly referred to as the Arrowhead region. The St. Cloud Region: The rocks of the St. Cloud region are mainly granites and related rock types such as monzonites and quartz diorites. The stones may be grouped into three major types, namely, pink granite, red granite and gray granite. Most of the pink granite occurs in the area to the southwest of St. Cloud. The rock is best described as stone with large pink crystals set in a finer grained black and white background. The minerals of the matrix occur in remarkably uniform sizes, and the pink crystals are sufficiently uniform in. their dis-

Jan 1, 1952

By Harlan W. Nelson, James B. Purdy

A study of data obtained during laboratory tests to determine the suitability of bituminous coals for use in residential underfeed stokers of the clinkering type has led to the following general conclusions: 1. Performance data obtained during the test furnish useful and practical information upon which to judge the performance of stoker coals. 2. Conclusions derived from the results of evaluation tests are in general agreement with use experience of the coals in the field. The method is of practical use in this regard. 3. Few direct relationships were found between data furnished by the evaluation test and determinations by standard laboratory tests. This is particularly true with regard to coke formation and clinker formation, two of the most important factors in determining the degree of satisfaction to be obtained in applications of the small underfeed stoker. It has not been the intention in this paper to prescribe this evaluation test as the best method for the evaluation of coals for use in residential stokers. That there are other test methods having similar objectives serves to point out the necessity for full-scale laboratory tests. At the present time, a standard test procedure for the evaluation of coals for use in residential stokers is not available. For the past several years a joint committee, appointed from the membership of the Residential Stoker Committee of Bitulninous Coal Research, Inc., and from the Engineer- ing and Research Committee of The Stoker Manufacturers Association, has been working to develop a standard procedure for testing and evaluating bituminous stoker coals. This procedure has been completed to the extent that it has been drawn up in tentative form. It is now planned to install the equipment in several laboratories and to run check tests, using portions of identical lots of coal, to determine the degree of reproducibility offered by the method. Completion of this standard test procedure should greatly facilitate the evaluation of coals for use in residential stokers. Acknowledgments The cooperation of coal operators sponsoring the coal evaluation tests, in granting permission to use data and results from their tests, is gratefully acknowledged. Acknowledgment is also made to Carroll F. Hardy, Chief Engineer, Appalachian Coals, Inc., for his interest and timely suggestions on the preparation of this paper. Special thanks are due to Ralph Sherman, Assistant Director, Battelle Memorial Institute, who was largely responsible for development of the evaluation procedure, for helpful suggestions and advice. References I. H. A. Sherman: The Evaluation of Coal for Use in Domestic Stokers. Univ. of Ill., Eng. Experiment Station, Circular No. 39 (1939). 2. . Q. Shotts: The Relation of Free- swelling Indexes to Other Characteristics of Some Alabama Domestic Stoker Coals. TP 2314, Coal Tech., (Feb. 1948); Trans. AIME (1948) 177, 502. 3. . J. Helfinstine and C. C. Boley: Correlation of Domestic Stoker Combustion with Laboratory Tests and Types of Fuels 11. Combustion Tests and Preparation Studies of Representative Illinois Coals. Ill. State Geol. Survey, Re1. of Zntlestigations No. 120 (19467 ) p. 35. DISCUSSION E. R. Kaiser*—AS a former fuel engineer assigned to residential stokers at Battelle, the writer has read with considerable interest the paper by Messrs. I'urdy and Nelson. It is gratifying to note that earlier test techniques have been developed and that criteria found important in the years 1935 to 1938 are still important in judging the performance of stoker coals for residential heating. The authors have plotted a number of the primary coal data singly against performance data of the stoker in an effort to establish relationships. Unfortunately, wide scattering of points usually resulted, which did not illustrate what experience generally indicates to be the case. For example, Fig 3 indicates almost no trend in the percentage of released ash converted to clinker with change in ash-softening temperature. Other factors, such as fuel bed temperatures and zones of heat release must have influenced the results. In the present state of our knowledge, we cannot explain all of the reasons why suitable stoker coals perform satisfactorily, nor why one satisfactory coal may be better than another. It would therefore

Jan 1, 1950

By G. B. Thomas

Several methods of analyzing pressure build-up data in wells have been presented by various authors. This paper reviews the theory and method of D. R. Horner and presents example calculations performed on data obtained by testing several different types of wells. These calculations include, (1) graphical estimation of final static pressure, (2) determination of the productive capacity of the pay away from the well bore and. (3) the degree to which the formation adjacent to the well bore has been damaged by completion or other causes. The methods of testing and precautions which should be taken to assure the best data possible are discussed. Limitations and reliability of calculated results are also treated. INTRODUCTION Pressure testing of wells is generally limited to the determination of producing and static mean formation pressures. The so-called "static" pressure determination.. along with PVT, electric log and production data, enable the reservoir engineer to determine, within reasonable limits, the drive mechanism of the reservoir and in some cases. the amount of edge water encroachment. Producing pressure tests enable calculation of productivity indices and allow the engineer to plan the systematic production of a pool for optimum conservation of suhsurface energy. The radial flow formula advanced by Muskat' has been based on the assumption of incompressible radial fluid flow. It has been known that reservoir fluids do not behave in an ideally incompressible manner. For example, incompressible flow theory indicates a simple logarithmic relationship between the difference of the instantaneous and static well pressures when plotted against time. 'The latter stages of this type of plot of pressure build-up data generally show a marked devia-, tion from the earlier straight line trend. which deviation may be shown to be due to the compressible flow of fluids toward the well bore. Shut-in times of 24, 48. 72. or at most 96 hours are currently in wide use for determining so-called "static" reservoir pressure. Due to the continuation of compressible flow of fluids into the well bore long after this arbitrarily taken shut-in time, the determination of static pressures has almost invariably resulted in lower than equilibrium values. Materials balance calculations made early in the life of a reservoir often result in a calculated reserve which later observations prove to be too low. Failure to obtain reliable "static" reservoir pressures within the prescribed 24 or 48-hour build-up period has undoubtedly been a major factor in obtaining these low estimates. Comparison of theoretical and actual productivity performance has indicated that formation damage and not being able to attain true static pressures have been partially responsible for the observed discrepancies. Research into the theory of compressible flow behavior has resulted in methods of applying these theories to the testing of wells. Application of the developed theories to pressure build-up performance indicates that in many cases the following can be estimated. 1. The static reservoir pressure. 2. The in-place effective formation permeability away from the well bore. 3. The degree or extent by which the formation has been damaged adjacent to the well bore, either through completion methods or subsequent damage due to fluid entry. REVIEW OF THEORY Horner2 has shown that the pressure build-up within a "point source" well is approximated by the following formula: Equation (1) is the approximate "point source" solution to the radial compressible flow equation advanced by Muskat. The solution assumes the following: 1. A point source well is producing at constant rate from the center of an infinite reservoir with a constant pressure at its external boundary. 2. The fluid flowing is present in one phase only. 3. The compressibility and absolute viscosity of the fluid remain essentially constant over the range of temperature and pressure variation encountered. 4. The well i.; shut-in at the sand face arid there is 110 after production into the well bore. 5. The formation permeability is homogeneous in the direction of flow. From Equation (1) it can be seen that if an ideal well were shut-in while producing from a reservoir under the conditions assumed, the pressure build-up would be a logarithmic func-

Jan 1, 1953

By A. Romero-Juarez

This paper is concerned with the effect of propping-agent concentration on flow capacity of a fracture in the case in which there is embedment of the propping agent. Previous published studies have shown definitely that there is a relationship between fracture flow capacity and prop ping-agent concentration, and it has been shown that the theoretical results are confirmed by laboratory experiments. The problem of directly finding sand concentration for maximum flow capacity of a sand-propped fracture is solved, and formulas and charts are given to obtain this concentration under various conditions of effective overburden pressure, medium sand-grain diameter and rock properties. It is shown that, for the same effective overburden pressure and the same rock characteristics, optimum sand concentration in pounds per gallon does not depend on medium sand diameter. For conditions met in bydraulic fracturing operations, it is found that sand concentration in grains per square inch for maximum flow capacity varies within a wide range of values; this indicates the convenience of using data of fracturing pressure and rock characteristics for calculating sand concentration in order to achieve the best results in fracture treatments. INTRODUCTION It has been pointed out in the published literature1 that one important factor controls the success of a hydraulic fracturing operation — the propping of the fracture. The main effect of the propping agent is to hold the fracture open by means of the reaction forces that oppose the pressure due to the overburden. It is assumed in this paper that the propping agent is sand, that the grains are spherical and uniform in size, and that this sand is distributed in a monolayer in the fracture. The existence of a sand concentration value for maximum fracture flow capacity has been recognized since the publication of one of the first studies on the subject.2 In the reference cited, flow capacity was assumed to be proportional to the cube of the fracture volume per square inch not filled by sand. The obtained experimental curves relating to flow capacity and sand concentration had the same appearance as the curves obtained with the afore-mentioned assumption. Recent work3 has shown that fracture permeability, in the case of a monolayer, can be calculated accurately by means of a modified Kozeny-Cacman equation; thus, it can be shown that fracture pqmeability is a function of free volume in the fracture and of the rock and sand wet surfaces. The flow capacity of a fracture is dependent on the permeability of the fracture and on the fracture width. If all excessive quantity of sand is present in the fracture, its width may be large but its permeability i:; low due to the greater flow resistance; too few grains of sand per square inch are present in the fracture, sand embedment in the rock is greater and the width becomes small enough to appreciably decrease the rate of flow. To determine optimum conditions of flow, therefore, it is necessary to know the dependence of fracture width on density or concentration of propping agent. This relation is given in a study4 in which the approach to embedment is similar to the one made in a penetration hardness test of metal. With the equacions presented by Darin and Huitt,3 the number of sand grains per square inch that results in maxirnum fracture capacity can be determined graphically by plotting fracture capacity vs sand concentration. Since the relationship between fracture capacity and sand concentration may be considered known, one is led to formulate the analytical problem of directly determining sand concentration for maximum fracture flow capacity, in order to study the effect of the various quantities that come into the process. The solution of this problem appears to be interesting both from the viewpoint of fluid mechanics and from the practical viewpoint of oil production in hydraulic fracturing treatments. EQUATION FOR FRACTURE FLOW CAPACITY It can be shown3 that the fracture flow capacity of a monolayer is given by the following expression,

By L. D. Jaffe

From the limited experimental data in the literature, preliminary values were derived for the thermal diffusivity of titanium alloys and for the quenching severity of various mediums used in heat treating them. The thermal diffusivity, under quenching conditions, was found to be approximately 0.006 sq in. per sec. The quenching severity H for brine is approximately 8, for water 4, for oil 0.9, for air 0.08, and for argon 0.02 in. -' The quenching severity of the Jominy-Boegehold end quench water jet against titanium alloys appeared effectively infinite. Using these results, graphs were prepared showing ideal round sizes and Jominy positions having cooled conditions equivalent to those at the center and surface of titanium alloy rounds, plates, and sheets quenched in various media. THE importance of proper quenching of steel parts has long been recognized, and much quantitative work has been done on the subject. Recently, it has become evident that for titanium alloys, too, the rate of cooling from high temperature is important in obtaining high strength and hardness. Maximum hardness and strength may not be attained in the as-quenched condition but rather after a subsequent tempering or aging treatment. Nevertheless, if cooling from temperatures near the ß transus is too slow, for the hardenability of the alloy used, a soft a-ß structure will be produced which will not harden on subsequent aging.'-" The quenching of parts of practical interest can usually be approximated by simple shapes, such as cylinders (rounds) and plates (or sheets) of various thickness and combinations of these. A standard that has been used extensively for comparison purposes is the center of an ideal round: a cylinder quenched in a medium that instantly cools its surface to room temperature." For steel, considerable practical use has been made of the Jominy-Boegehold end quench hardenability specimen,".' and its value for titanium alloys is becoming evident. It therefore appeared worthwhile to attempt comparison of cooling conditions in titanium alloy rounds and plates of various sizes quenched in various media with cooling conditions in ideal round and Jominy quenches. These comparisons are based upon measurements reported in the literature. Quenching Severity of Jominy Water Jet and of Brine As one starting point, there were available the data of Phillips and Tobin8 for positions on a Jominy specimen having the same hardness as the centers of round bars of various size quenched in 5 pct NaOH brine, shown in Table I. The Jominy speci- men is a round bar, normally 1 in. diam and 4 in. long, quenched by a water jet at one end and air cooled on the other surfaces.' After quenching, hardness is measured vs distance from the quenched end on longitudinal flats ground 0.015 to 0.10 in. deep. The data of Table I were based on a subsized specimen of 3/4 in. diam. Since the positions tested are close to the jet quenched end, cooling there may be taken as due solely to the jet, and air cooling can be neglected. These positions are then equivalent to various locations near one surface of an 8 in. plate jet quenched on both sides. Assuming that the thermal diffusivity is constant and that Newton's law of cooling holds, temperature vs time curves at corresponding positions in Jominy specimen and round bar then may be calculated for various assumed quenching severities of the water jet and of the brine. Cooling curves are generally taken to be metallurgically equivalent when they show the same time to reach a temperature halfway between the initial (heat treating) temperature and the final (quenching medium) temperature.5 On this basis, using Russell's tables,9 upplemented by heat flow tables derived in similar fashion by the author for positions near the surface of severely quenched plates, it was found that the experimentally observed equivalence, Table I, was best matched if the quenching severity H of the Jominy water jet was taken as infinite (H = co). This means that the jet quenched surface should be considered as instantly cooled to the temperature of the water

Jan 1, 1956

By J. B. Cheatham, J. G. Yarbrough

Although adequate removal of cuttings from beneath a drill bit is important for efficient drilling operations, very little basic data are available relative to the fundamentals of chip removal by hydraulic jets. A discussion is presented in this paper of an experimental investigation of the jetting action of hydraulic jets in removing loose particles from the bottom of a cylindrical hole. Conditions for which the jet is no longer capable of removing chips from the bottom of the hole are determined. This situation represents equilibrium between the chip removal force and chip Bolddown forces such as gravity and pressure. In most of the tests loose particles were jetted with water or a water-glycerine mixture to determine the dependence of chip removal on hole size, jet size, height of jet off bottom of hole, flow rate, particle density and fluid viscosity. ,A test with a pressurized mud system indicated that relatively low pressures can completely overcome the removal action of a hydraulic jet. Although the system studied herein is not directly applicable to a rotary drill bit, the work with such simplified systems can provide a better understanding of the chip removal action of jets, and with logical extensions it may provide a reasonable basis for the best use of fluid jets in drilling. INTRODUCTION The primary deterrent to maximum drilling rates is the inability of the drilling system to remove rock cuttings efficiently enough to prevent interference with the drilling action.1-2 The objective of chip removal studies is to permit predicting and controlling removal forces under downhole drilling conditions. Conditions at the bottom of a hole during rotary drilling are exceedingly complex and are not likely to be described in a quantitative way by investigations in terms of the total drilling action until a better understanding is developed of the simplified components of the problem.3 The present study is concerned with the elementary condition of removal of chips by a single central jet. Even this relatively simple model provides mathematical difficulties because of the turbulent nature of the flow from the jet and because of the shape of the bottom of the hole beneath the jet. Theoretical and experimental studies have been made of turbulent jets impinging normally on an infinite body and deductions based on analytical solutions to simplified problems can give some insight into the problem of cutting removal by a jet. However, because of the present lack of understanding of the behavior of the interaction between the fluid jet and the chips being removed, an experimental approach was chosen for the present study. Methods have been developed for maximizing hydraulic horsepower, impact force and jet velocity; but whether maximizing these parameters maximizes chip removal with present drilling bits has not been demonstrated. Simplifying the problem of chip removal may make it possible to develop some understanding of the manner in which the jet velocity is dissipated. Better understanding of a simple case should materially assist in extending analysis to more complicated cases. Thus, we are not concerned in the present study with the rock fracturing process itself but only with the removal of the debris from the bottom of the hole. A problem which is quite similar to the chip removal problem is the suspension of solids in stirred vessels. This problem has been studied by the chemical industry and correlations have been obtained by dimensional analysis which permit the design of mixing vats7 An approach similar to that used in the mixing vat problem is used in the analysis of the jetting data in the present paper. EXPERIMENTAL PROCEDURE The test equipment arrangement shown sche-matically in Fig. 1 allows the jetting action to remove particles until an equilibrium height is attained for each combination of hole size, jet size and flow rate.*** Equilibrium conditions require that the removal force is unable to remove additional particles. This balance between holddown and removal forces implies a relationship between the two forces which is constant for the particular system. When the holddown forces are constant,

Jan 1, 1965