Search Documents

By H. L. Hartman

The drop tester has proved an invaluable tool for the investigation of percussion drilling in the laboratory in "slow motion". It has allowed the process of rock penetration by impact to be studied a single blow at a time. In the present work, the relation of cratering to indexing has been determined under conditions simulating those in the bottom of the hole at atmospheric pressure. Indexed blows form craters which are influenced by others adjacent to them. Two unique effects are involved: (1) the provision of additional free faces in proximity to the point of impact and (2) the creation of subsurface damage by the previous blow(s). Both have a pronounced influence on the volume of rock removed per blow, which governs the rate of penetration in actual drilling. The emphasis in this study was placed on determining effects of index distance and impact energy on crater volume when blows were struck on a previorusly "drilled" or damaged surface, a situation most representative of down-hole conditions. The results indicated that (I) the optimum index distance on a damaged surface is greater than on a fresh, undamaged rock surface; (2) the volume of rock broken at optimum index distance is about the same for both surface conditions with die-shaped chisels but tends to be greater on a cratered surface with wedges; (3) the optimum index distance and maximum crater volume are proportional to the energy level; (4) at a given energy level, the volume of rock broken by dies generally exceeds that by wedges; and (5) maximum crater volume varies, approximately, inversely with width of die but is nearly independent of included angle of wedge. INTRODUCTION This is the third of a series l,2 of papers reporting on one phase of a continuing program of drilling research being conducted in the Rock Mechanics Laboratory at The Pennsylvania State U. The phase under consideration is crater geometry — blow energy relations in percussion drilling. While the previous papers have discussed mainly craters produced by single blows, indexed-blow studies are the subject of this paper. The basic mechanism of drilling under study here — that of penetration of a chisel-shaped bit subjected to impact — is the one employed in probably 90 per cent of all rock drilling in mining, percussion drilling. It also forms the fundamental action in nearly all oil well drilling with the so-called rotary (roller-bit) method. Here, the action is more complex, however, and may in actuality comprise rotary-percussion drilling. Certainly, the down-hole air and mud hammer drills, which superimpose percussion on rotation, belong in this latter category. In any event, impact blows constitute one of the principal mechanisms in the vast majority of drilling, mining or petroleum, and it is essential in drilling research to study percussion in all its basic aspects. Percussion drilling can be simulated in the laboratory by means of an impact drop tester. Striking one blow at a time with chisel-shaped bits, the drop tester permits study of most of the parameters in the drilling process except blow frequency and the interrelationship of certain variables. Single-blow craters, the simplest to produce, allow the determination of basic relationships, such as volume vs blow energy; while indexed craters most closely resemble the action in real drilling and provide data for the establishment of optimum operating criteria. Indexed-blow craters are those which are influenced by others in proximity. The simplest case to visualize, and the one studied in detail here, is that in which a given crater is formed adjacent and parallel to a previous crater, such that the second crater breaks into the first. A more complex situation usually exists downhole, where indexing is achieved by rotation of a radial-winged or rolling-cutter bit. But, admittedly, the phenomenon of indexing is more likely to be understood eventually if it is first studied in its simplest form, that related to parallel craters. Limited research has been conducted on indexed blows. The most complete study was that by Simon, 3,4 who proposed universal indexing curves to relate the important variables. His work, however, was conducted entirely with rows of craters formed on a

By Richard J. Fruehan

ThE thermodynamic activities of nickel and chromium in liquid Ni-Cr alloys have not been previously reported. There have been several investigations on the solid alloys.1-3 In the present investigation the activity of chromium in liquid Ni-Cr alloys at 1600° C was determined by the use of the following galvanic cell, The use of galvanic cells employing solid oxide electrolytes for determining the activity of metals in solid alloys was first demonstrated by Rapp and Maak4 on Cu-Ni alloys between 700° and 1000°C. Recently Schwerdtfeger5 tested the performance of the ZrO2(CaO) electrolyte at 1600°C and found it satisfactory for oxygen pressures as low as 10"12 atm. He also demonstrated the use of solid oxide electrolytes for determining the activities in liquid metal systems in measuring the activity of manganese in Fe-Mn alloys at 1500°C. Extensive investigations have been carried out in this laboratory, exploring the oxygen partial pressure lower limit at which ZrO2(CaO) electrolyte can be used satisfactorily in oxygen galvanic cells at 1600°C. In a paper to be published at a later date,9 it is demonstrated conclusively that ZrO2(CaO) exhibits insignificant electronic conductivity at oxygen pressures down to about 3 x 10"13 atm which is lower than the oxygen partial pressure in equilibrium with Cr-Cr2O3 at 1600°C. Therefore, the activity of chromium can be related to the measured reversible electromotive force of the cell in the present investigation since no electronic short circulating will occur. Fig. 1 is a schematic diagram of the cell used. The Cr-Cr2O3 electrode was made of equal weights of chromium and Cr2O3 powders which were mixed and packed into a recrystallized alumina crucible. A Pt-wire lead was embedded into the Cr-Cr2O3 compact to insure good contact. A ZrO2(CaO) crucible, supplied by the Zirconium Corp. of America, contained the Ni-Cr melt equilibrated with solid Cr2O3. The ZrO2(CaO) crucible was placed on the Cr-Cr2O3 electrode and additional Cr-Cr2O3 mixture was packed around it. Contact with the Ni-Cr alloy was made with a Cr2O3 tube, which is a weak electronic conductor. A platinum wire lead was attached to the Cr2O3 tube. The galvanic cell was lowered slowly into the hot zone of the furnace tube in order to avoid too large a thermal shock which could break the ZrO2 (CaO) crucible. The furnace was con- tinuously flushed with purified argon during the experiment. The temperature was measured with a Pt-Pt/10 pet Rh thermocouple and maintained at 1600° ± 5°C. The electromotive force was measured with a recording potentiometer. The furnace was turned off when measurements were being made in order to avoid induction effects and the reversibility of the cell was demonstrated in the usual manner. The electromotive force was measured several times over a period of approximately 1 hr. During that period the electromotive force remained constant within ± 3 pet and did not continuously drift in any one direction. After the experiment the alloys were analyzed for chromium in the usual manner. The measured electromotive forces of the cell at 1600°C, (2912°F), are listed in Table I. The activity of chromium, with solid pure chromium as the reference state, (a), was calculated using Eq. [l], , n F E I, n lna RT- U] where E is the measured electromotive force, F is the Faraday constant, and n is equal to 3, the number of equivalents per mole, (nF/RT = 18.59 volts'1 at 1600°C). The resulting activities are listed in Table I. The activity of chromium relative to pure liquid chromium was also calculated. The free energy of fusion of chromium at 1600° C was calculated using

Jan 1, 1969

By Reynold Q. Shotts

CHEMICAL methods, based on the relative rates of oxidation of fusain, bright coal, and dull coal by nitric acid, have been devised to determine these coal components.1-4 Results obtained by oxidation methods for fusain have been checked against results obtained from microscopic methods5, 6on duplicate samples of the same coals, but to the author's knowledge this has not been done for bright and dull coal components. For this reason it is not certain that the two methods of analysis identify essentially the same chemical or physical units. It would be highly desirable to see results of the application of both methods to duplicate samples, but in the absence of any such data the author has attacked the problem indirectly. From the U. S. Bureau of Mines samples were obtained of three Alabama coals which the USBM had analyzed optically and reported on over the past 30 years. These samples were subjected to analyses by oxidation rate methods. Results of this work, and comparisons with the USBM analyses, have been published." This was the first indirect approach to the problem. The present report attempts a second indirect approach by way of internal validation. By nitric oxidation samples of two whole coals were carefully analyzed for fusain, bright coal, and dull coal. One coal was analyzed in duplicate. Duplicate portions of each of the coals were divided into three density fractions by means of heavy liquids. A duplicate portion of one of the coals was divided by sieving into three size fractions. Each fraction was analyzed by oxidation and its percentage composition calculated in terms of fusain, bright coal, and dull coal. Because the weight percent of each fraction was known, a material balance calculation for the whole coal was also made. The resulting reconstituted analysis of the whole coal could be compared to that determined by direct analysis. In addition, specific reaction rate constants were determined for each component and for each whole coal or fraction. Arbitrary reactivity indexes were calculated by dividing by 100 the sum of the products of the percent of each component in the coal and its specific reaction rate constant. The resulting figure was an average reactivity index for each coal or fraction. Weighting the reactivity indexes for each fraction by the percent of the fraction in the coal gave a reactivity index for each whole coal which could be compared to that calculated directly from the whole coal analysis. If oxidation analyses really delineate definite physical entities within the coal, or even definite groups of similar entities, reconstituted analyses calculated from fraction analyses should check closely those made of the whole coal. It probably is true that optical methods identify and describe a greater variety of components than do chemical methods and that variations are wide in the appearance and quantity of those components identified optically. Chemical methods based upon differences in oxidation reaction rates would of necessity be less discriminating, as between similar components, than would optical methods. The procedures followed for oxidizing the samples, analyzing the residues, plotting and calculating the percentage of each component, and calculating its specific reaction rate constant have been fully described.'-' In brief the method originally proposed by Fuchs et al.1 for the determination of fusain consists of the oxidation of small samples of coal in boiling 8N nitric acid in a condenser-fitted flask. After boiling for periods of 1/2 to 4 hr, the unoxidized residue is filtered and washed. The washed residue is treated with normal sodium hydroxide, diluted, and allowed to stand for several hours. The resulting brown liquor is removed, and the filtered residue dried, weighed, ignited, and weighed again. The ash-free residue is expressed as a percent of original dry, ash-free coal. The percent residue is plotted against time, and the extrapolation of this line to zero time gives the percentage of dry ash-free fusain present in the original sample. The shape of the resulting time plots has been explained' by the assumption that they are the result of two different types of reaction, the first part representing a first order reaction with rate a func-

Jan 1, 1956

By L. M. Pidgeon, K. T. Aust

There has been considerable interest in the possible use of titanium in magnesium alloys.' Zirconium has shown some promise in this connection2 and its general similarity with titanium suggests that the latter might act in a similar manner. A literature survey revealed that quantitative data on the Mg-Ti system was unavailable. Several patents3 have claimed that titanium additions from 0.2 to 4 pct to magnesium alloys were possible, but no mention was made as to the form in which the titanium existed in the alloy. Kro114 succeeded in introducing only traces of titanium into magnesium by bubbling TiCl4 through the metal under argon or by reacting it with sodium titanium fluoride. The application of theoretical data given by Carapella5 based on Hume-Rothery's principles, involving atomic size factor, crystal structure, valency and the electro-chemical factor, suggests that a Mg-Ti alloy is a favorable case, and the system appeared to warrant experimental examination. Experimental Procedure and Results THERMAL ANALYSIS If titanium is appreciably soluble in magnesium, a change in the melting point of the magnesium might be detectable using standard cooling curve methods. Magnesium was melted in graphite crucibles under an argon atmosphere, the assembly being enclosed in a silica tube. Graphite thermocouple protection tubes served also to stir the melts. The apparatus was very similar to Fig 1, with the addition of a refractory and baffle system to prevent undue heat losses from the top of the crucible. Chromel-alumel thermocouples were calibrated using Al of 99.97 pct purity. Dominion Magnesium Limited sup- plied redistilled high purity magnesium of the analysis given above. Titanium was added in three different forms: 1. Titanium powder —100 mesh, from the Titanium Alloy Manufacturing Co., Niagara Falls, N. Y. 2. Sheet titanium from the U.S. Bureau of Mines, produced by Mg reduction of TiCl4. 3. Magnesium —50 pct titanium master alloy from Metal Hydrides Inc., Beverly, Mass. The melting point of the high purity magnesium used was measured experimentally as 651.0°C. More than a dozen tests were conducted using titanium from the three sources referred to above, in calculated additions up to 20 pct titanium, at temperatures between the melting point and 1000°C and holding periods up to 6 hr. In no case was evidence obtained of solubility of titanium in magnesium, using inverse-rate and time-temperature curves. The melting point of the magnesium was unchanged within the accuracy of measurement, namely -+0.5°C; and no other thermal arrests were detected. Metallographic investigation of the thermal analysis billets indicated that the titanium additions were apparently mechanically entrapped in the magnesium in segregated areas. Consequently, these samples were not analyzed for titanium. The master alloy proved to be a mechanical mixture of titanium particles in a magne- sium matrix. These results indicated that the titanium solubility, if such existed, could not be obtained by the usual thermal methods. X RAY DIFFRACTION INVESTIGATION In an effort to detect solubility of titanium in magnesium, samples were investigated using both the Debye-Scherrer and the Focusing Back-Reflection methods. Filings from samples of the thermal analysis billets and from pure magnesium were annealed in argon one hour at 350°C to relieve mechanical strain. Measurements made of the interplanar spacings showed no difference between the Mg-Ti samples and pure magnesium. The interplanar spacings could be measured to within 0.0002A, and the greatest variation found was 0.0004A, in the back-reflection method. The diffraction lines for magnesium were not shifted by the titanium additions indicating that the solid solubility of titanium in magnesium is of a very low order—less than 0.5 pct. From both diffraction methods, a d or interplanar spacing of 0.817A was obtained for the redistilled high purity magnesium. This latter value is not given in the standard X ray diffraction cards for magnesium metal or vacuum distilled magnesium. Theoretical calculations for a close-packed hexagonal space lattice for magnesium indicate that the planes {2134) should give a line which was found. The relative intensity for this reflection at 0.817A is slightly less than that at 0.870k for magnesium. SOLUBILITY OF TITANIUM IN LIQUID MAGNESIUM The Mg-Mn system was examined by Grogan and Haughton6 who were

Jan 1, 1950

By J. E. Drury, W. M. Byers, E. J. Lynch

Measurements made with starch ~~~~s indicate that the Presence of starch serves to depress the streaming potential of a clay mud system. Measurements with pure starch slurries show that the addition of NaCl or CaO reduces the streaming potential significantly. The addition of NaOH, however, causes only a minor reduction in potential even when there are large changes in mud resistivity. INTRODUCTION In their excellent study of the streaming potential, Hill and Anderson' proposed that the presence of conductive solids within the filter medium will reduce the generated potential. This theory was demonstrated with measurements on shaly sandstone cores having various "b" factors; the cores with the higher negative b factors showed markedly lower potentials. (The factor b was originally defined by Hill and Milburn2 in an expression relating the resistivity of a shaly sandstone to the resistivity of the interstitial water: log F, = log F,, + b (log 100R,)'. 'he factor b, therefore, is a measure of the contribution of the conductive solids to the total conductivity of the sample. It has been found experimentally to be directly related to the amount of chemically exchangeable ion per unit pore volume.) In the logical extension of this hypothesis to the streaming potential created by filtration through mud cakes, they suggested that the presence of nonconductive starch in a drilling mud would have the opposite effect; i.e., it would cause the streaming POtential to be increased. The high streaming potentials that they observed with lime muds were attributed to this cause. A cross-plot of some of their results which are relevant to this discussion is given in Fig. 1. The lime-mud curves and untreated-mud curve are taken from curves given in their paper. The gyp-mud line was estimated from their statement that the measured streaming potential for gyp muds in the range of 1.5 to 2.0 ohm-m was approximately 18 mv/1,000 psi. These gyp-mud potentials were the lowest that they reported for any of the bentonitic-clay muds, and were the only ones which were less than the untreated-mud potentials. This very low streaming potential of gyp muds has also been observed in field studies.' By way of comparison, a U. of Texas publication4 gives the approxi-ate compositions for lime and gyp muds shown in Table 1. These materials are added to a bentonitic-clay mud base. It can be seen that both of these muds contain significant amounts of calcium salts and that the gyp mud has the higher concentration of starch. It would appear from this and from the data shown in Fig. 1 that starch may well be the cause of a streaming-potential reduction, rather than a potential increase. EXPERIMENTAL A series of streaming-potential experiments has been performed using a cell of a design similar to that reported by Wyllie.5 Calomel electrodes were used in both the mud slurry and the filtrate. Pressure in the cell was generated with high-pressure nitrogen. The primary objective of the investigation was to determine the effect of starch on the potential. Five different slurries were made which ranged in composition from untreated ben-tonitic mud (Baroid Aquagel) to a pure starch suspension (Baroid Im-permex). The intermediate slurries were prepared by mixing together a 9-lb/bbl clay mud with a 12-lb/bbl starch suspension in the ratios 1:3, 1:1 and 3:1. When starch was used alone, it was necessary to use an 18-Ib/bbl mixture, since leaner mixtures would not form a suitable filter cake. The clay mud used in the experiments was aged for three months before use. In addition to this series of measurements on the starch-clay mud mixtures, other experiments were performed to determine the effect on the potential of the absolute starch and Aquagel concentration, and the effect of adding various chemicals. RESULTS The effect of starch on the streaming potential is shown clcarly in Fig 2. Increasing percentages of starch resulted in n continually decreasing

By Joseph L. Gillson

IT is a somewhat curious circumstance that the newest of the Institute's several awards should be conferred in the oldest of our several professional fields—for there is little question that geology and mining have antedated metallurgy and petroleum. Perhaps so did geophysics. After all, one geophysical instrument, the magnetic compass, has been successfully employed for many centuries. Young in years the Jackling Award may be, but in distinction—whether in terms of the title it bears or in terms of the first three medalists (Reno Sales, E. D. Gardner, and the late Father MaceIwane)—the Jackling Award is already one of the ranking awards of the Institute. The name of the 1957 medalist will bring added luster to an already illustrious roll. J. L. Gillson was born and brought up in Evan-ston, Ill., in an area characterized by singularly flat and uninspiring topography, and devoid of ore deposits and even of bedrock—a most unlikely environment from which to expect a future Jackling medalist to emerge! On graduating from high school he enrolled as a chemistry major at Northwestern University. At the end of his frkshman year his scholastic prowess had won him a prize of two chemistry books. But these went unread, for about this time he came under the spell of U. S. Grant, that distinguished and enthusiastic geologist who for many years was head of Northwestern's geology department. From thence forward Joe was a dedicated geologist. After a Navy interlude in World War I, he went to MIT for graduate work, substituted for the late Professor Palache for one year at Harvard in mineralogy, received his doctorate at MIT and stayed on as a meniber of the MIT faculty, where the geology department was headed by that great dean of all economic geologists, Waldemar Lindgren. In 1926 Joe spent a summer looking into barite deposits for the du Pont Co., and since 1928 he has been du Pont's chief geologist. Do men who lead such active lives in our exploration organizations generally have permission to publish their observations? Do they generally find time to publish, even if they have permission? Alas, no. But in Joe Gillson's case, yes. His bibliography runs to well over two score papers, including one of the very first papers in this country on the applications of petrography to problems of Portland cement. He has written on such topics as the origin of talc and the genesis of ilmenite, fluorspar, and alkaline rocks. Each year for a number of years his annual reviews of industrial minerals in the Mining Congress Journal have been eagerly awaited. - Do men who do so much to advance their science have anything left to give to the advancement of their profession? In Joe Gillson's case, yes. He has been a vice president of the Mineralogical Society of America and has been active in the councils of the Society of Economic Geology, particularly in their research programs. This year he is president of the American Geological Institute, currently perhaps the most challenging post to which a geologist may be elected. His contributions to AIME scarcely need documentation, so I shall mention only three—his chairmanship of the Committee on Democratization of the Institute, a committee whose work and recommendations brought a new and hopeful leaven to AIME at a time when this was sorely needed; his chairmanship, a few years ago, of the Industrial Minerals Division; and his current place on the Board of Directors as Vice President of the Institute. Such' a catalogue of contributions and achievement, even as condensed as I have had to make it, certainly merits the Jackling Award. Yet, implicit although not often expressed in the conferring of most awards is the hope that besides recognizing a distinguished career, it will provide stimulus for further effort. In the present instance, such an expectation is already well on the way to fulfillment. Only recently Joe Gillson has undertaken the editorship of a new and greatly revised edition of the Institute's Seeley Mudd volume, Industrial Minerals and Rocks. And so, Mr. Chairman, it is with pleasure and enthusiasm that I present Dr. J. L. Gillson for the 1957 Jackling Award.

Jan 1, 1958

By Rudolph G. Wuerker

THE scarcity of values of tensile strength of rocks has been explained by the lack of successful testing procedures. In the case of mine rock a description is given' of the difficulties encountered in testing a cylindrical specimen, such as a core, by conventional methods. Over a period of years the following method has given definite and reproducible results with the weakest as well as with the strongest rocks. It does not completely supersede the use of cores with special fixtures but is a supplement in all cases where cores cannot be obtained, as from soft rocks, or in cases where it is less expensive to prepare test specimens by cutting them out of the rock instead of drilling cores. Principle and equipment are the same as for the test for tensile strength of hydraulic-cement mortar.' The test specimen, Fig. 1, has the shape of a briquet. While in the original cement mortar test the briquet is cast in a special mold, it is prepared from rocks in different ways, depending on how easily they can be cut and shaped. Soft rocks, which cannot be core-drilled with a carboloy or diamond bit, are simply hand-cut. Only two dimensions need be watched. The first is the 1-in. diam at the narrowest cross-section of the briquet. The other critical measure is the radius of curvature of the waistline, as the roller supports in the grips have a fixed distance. This radius is ground out of the solid by means of a carborundum grinding wheel having a 3/4-in. radius. Medium hard rocks can be core-drilled with a carboloy bit. The resulting core can be used for nondestructive sonic testing first, and after that for any destructive test. By using an EX-bit and by carefully placing the coreholes, preferably by using a tenplate such as shown in Fig. 1, it is possible to obtain from the rock a punched sample from which numerous tensile briquets can be made. The outside radius of the EX-bit differs from the radius of curvature of the briquet by 1/8 in., but this still permits placing and aligning the specimen in the grips. In the case of bedded rocks the core might have bedding planes normal to the plane of the briquets, and rocks can be tested in any arrangement of the bedding planes desired. Hard rocks, limestones, igneous, and metamor-whic rocks can only be diamond-drilled or diamond-cut. Here the method of getting the tension briquets by accurate placing of EX-drill holes is especially economical. The tops of the briquets made from hard rocks cannot be rounded; they are straight cuts made with a diamond cut-off saw and rounded off on a polishing wheel. Results: As long as specimens broke over the waistline the results were considered acceptable. Further statistical treatment of the tests' showed a satisfactory percentage of standard deviation. The tensile strength values obtained by this method do not represent true values because of the stress concentration caused by the curvature of the side of the piece and because of the closeness of the grips. The ratio of maximum to average stress at the plane of failure has been determined to be about 1.75." All tensile strength values listed in Table I are corrected accordingly. To avoid this stress-concentration, if there are a sufficient number of cores, tensile strength can be measured by imbedding the cores in mortar in the two outer briquets in the gangmold.4 Strain-Measurements: The applicability of the briquet specimens for strain observations was tested in the case of sandstone and shale. Two element Rosette SR-4 strain gages were used. Young's modulus and Poisson's ratio, both in tension, were computed and found to be different from those in compression, determined during the same test series and from the same rock, see Table I. References 1 L. Obert, S. L. Windes. and W. I. Duvall: Standardized Tests for Determining the Physical Properties of Mine Rock. U. S. Bur. Mines R.I. 3891 (1946). 2 Test for Tensile Strength of Hydraulic-Cement Mortar, ASTM Standard C 190-44. S F. O. Auderegg, R. Weller, and B. Fried: Tension Specimen Shape and Apparent Strength. Proc. ASTM 11939) 39, pp. 1261-1269. 4 API Code 32.

Jan 1, 1956

By R. F. Jones, N. E. Truitt, M. Prats

In most published discussions and theories of in situ combustion, the combustion fronts are assumed to be vertical. However, evidence from field tests leaves no doubt that combustion fronts often advance more rapidly along the top than near the bottom of a formation as a result of difference in density between injected air and formation liquids. The approximation proposed in this paper to determine the movement of the resultant tilted combustion surfaces states that the vertical rate of movement of combustion surfaces is proportional to the horizontal oxygen flux. Where self-ignition is possible, the proposed method demands that a secondary combustion surface exist around production wells which produce some oxygen. These secondary combustion surfaces may be formed long before the primary combustion surface can advance from injection to production wells. Heat liberated near production wells at these secondary combustion surfaces can contribute to an early increase in production rate. Results indicate that significant oil recoveries cannot be obtained from the usual flood patterns (five-spots, seven-spots, etc.) without producing large volumes of unused oxygen. Ideally, to increase oxygen- consumption efficiency, well patterns should allow oil production from a first line of production wells and gas production from more distant lines of producers. However, it may be desirable to produce some gas at all wells to support (and benefit from) active secondary combustion surfaces. Results indicate that the well spacing through which combustion can be advanced is larger than that predicted by other methods. A large number of production wells may still be desirable to take quick advantage of gravity drainage. From a comparison with results at South Belridge field, California, it appears that this method adequately describes oxygen concentration and temperature histories and combustion-front shapes. However, this method does not accurately locate the most advanced point of the combustion surface. There is some field evidence to substantiate the actual presence of secondary combustion surfaces at South Belridge. Use of the proposed method appears warranted at this time when lay-over of the combustion surface can be anticipated. INTRODUCTION The assumption of vertical combustion fronts has been embodied in all previous publicationsl-6 which use the movement of combustion fronts away from injection wells to determine the temperature and fluid distributions in the reservoir. The only paper concerned with a mathematical model of the combustion process in which a nonvertical combustion front is used was written by Gottfried.11 Actually, nonvertical combustion fronts have been observed in most in situ combustion field tests for which adequate data are available.7, 8, 6 In practice, the typical vertical extent of the burned zone decreases with distance from the injection well, and this burned zone is at or near the top of the sand body. In some cases, such as at South Belridge field near Taft, Calif., the combustion surface is almost horizontal over a very large area.8 Thus, for some years an obvious and serious gap has existed between theory (vertical fronts) and practice (tilted fronts). This is indicated in Fig. 1. Tilted combustion fronts such as observed at South Belridge8 sometimes result from the natural tendency of injected gases to rise to the top of an oil sand. An initially present gas cap, even if thin, would make apparent the effect of gravity segregation from the first, particularly since formations considered for in situ combustion field operations are generally selected for their

Jan 1, 1969

By C. D. Stahl, S. M. Farouq Ali

This investigation attempts to describe and simulate the alcohol displacement process by means of a cell model, as employed in chemical engineering practice. The proposed model is more simple than previously proposed models. and utilizes parameters chosen on a theoretical basis. The model successfully reproduced the formation of the stabilized bank and the breakthrough of alcohol, the latter depending on one of the model parameters, which may be correlated with the length of the porous medium. Moreover, the effects of the phase behavior of the liquid system involved, as observed in experimental studies, were reproduced. Several variations of the basic model were devised and tested on a digital computer. These included the cases in which: (1) the actual value of fractional flow was used in cell-to-cell computations; (2) the number of cells was varied within the same run; and (3) incomplete rather than complete phase equilibrium was assumed within each cell. The proposed cell model clarifies the basic mechanism of the process. Detailed concentration profiles obtained for each cell, for instance, showed the mechanism of bank formation in relation to the phase behavior characteristics. The results obtained indicated a varying degree of phase equilibrium concommitant with changes in the velocities of the phases in an actual alcohol displacement. This condition was approximated by changing the number of cells during the simulation. Interesting information was obtained on the influence of path length on the efficiency of alcohol displacement, which has been the subject of some controversy. Certain limitations preclude the use of the proposed model as a substitute for experimental studies. The results obtained were, nevertheless, of value in interpreting the experimental results. INTRODUCTION During recent years considerable effort has been directed toward an understanding of alcohol displacement, the process whereby oil and water are recovered from a porous medium by the continuous injection of a solvent. The complex nature of the physical process involved has so far defied a complete mathematical treatment. Other methods of approach, amounting to an overall material balance, have been proposed, yielding useful information on certain aspects of the process.1-3 Taber et al, in particular, defined the displacement mechanism in terms of the phase behavior of the alcohol-oil-brine system invoIved.2 Wachmann reported a mathematical treatment of alcohol displacement subject to certain simplifying assumptions? Donohue proposed the use of a "cell model" for simulating alcohol disPlacement.5 The nature of the assumptions involved limited the utility of the model. The present work attempts to examine the variables involved in the simulation of alcohol displacement and discusses several possible versions of the basic cell model. Under certain conditions the model results are similar to the experimental results. In particular, the spontaneous formation of of the stabilized bank and the effects of the system phase behavior were successfully reproduced. PREVIOUS WORK ON CELL MODELS Cell models and the theoretical plate concept are often used in solving chemical engineering problems in which an explicit mathematical solution may be difficult or impossible to obtain. Examples of such applications occur in distillation, gas-liquid chroma tography,6 reactor technology, absorption, etc. In petroleum engineering, such a model was used by Attra7 to describe non-equilibrium gas drive, and by Higgins and Leigh ton8 to calculate sweep efficiency in water flooding. Aris and Amundsen pointed out the equivalence between the diffusion model and perfectly mixed cells connected in series.9 Deans proposed a three-parameter cell model to simulate two-component

Jan 1, 1966

By R. G. Gifkins

Mean boundary migration depended linearly upon creep extension under various conditions. Prolonged annealing or air-casting increased the proportion of immobile boundaries at a given stage of tests. Grain growth during annealing followed a power law. The linking of migration and sliding in creep and models for migration are examined; activation energies are calculated from these models. Segregation of oxygen to boundaries appears to be a controlling factor. THE idea that boundary migration and boundary sliding during creep are linked has often been suggested,'-5 especially when evidence has been obtained for an alternation between the two processes. In simple bicrystals sliding can occur without any migration6 but in more complex bicrystals7 it is observed. McLean' and Couling and Roberts 5 made measurements of both sliding and migration; McLean concluded that there was no linking of the two processes whereas Couling and Roberts supported the idea. McLean' also developed a model for migration during creep and annealing. The present work was initiated as part of a study of the recrystallization of lead during creep9 but has been extended to examine some aspects of the mechanism of migration. MATERIALS AND METHODS Specially refined high-purity lead* was used; the total of impurities, excluding oxygen, was 0.0008 wt pct. Oxygen was removed as far as possible by melting in vacuo, treating at 550" to 600°C with hydrogen and then casting at 450°C under a vacuum of 10-2 mm Hg; the oxygen was then believed to be 0.00005 wt pct. Specimens were made from extruded strip and then annealed at the required temperature in an air oven, controlled to 5 lac. They were then polished in a mixture of acetic acid and hydrogen peroxide. At this stage some specimens were etched to give grain contrast and their grain sizes determined by the intercept method. Migration of boundaries during creep was measured on the surface of polished (not etched) specimens, directly loaded at room temperature (22°C) or 50°C (52"C) to initial stresses of 300 to 700 1b per sq in. Twenty six vacuum-cast hydrogen-reduced specimens were used and, in addition, two specimens were from material air cast after hydrogen reduction ("partial" air cast) and two from air-cast material not treated with hydrogen ("full" air cast). Duplicates were made in the majority of tests. At intervals through the creep tests, specimens were unloaded and the migration measured, normal to the direction of migration at the point of measurement, for each boundary along standard traverses. Usually between 110 to 180 boundaries were measured at each of some twelve stages of the tests. Migration was measured on the projected image (X100) using prism illumination. Fig. 1 illustrates the method of measurement and shows the kind of variation of migration found from boundary to boundary and along particular boundaries. Fig. 2 gives an example of more extensive migration. The mean value of distance migrated was generally reproducible to *10 pct but occasional points gave a scatter of *20 pct. The correctibn applied by McLean8 to allow for the contribution sliding makes to the apparent distance of migration could not be used, for the sliding at most boundaries was too small to measure with the equipment available. There were two reasons for supposing the correction unnecessary in the present work;

Jan 1, 1960

By John H. Schemel

Large Zircaloy-2 hammer or press forged bars did not exhibit the uniform excellent corrosion resistance to steam normally expected of the alloy in wrought form. Weight gains of coupons cut from forged bars were 40 to 60 mg per sq dm compared to the 28 mg per sq dm obtained on hot rolled sheet and strip. The mechanism of this corrosion and its relation to microstructure is discussed. After these observations, a solution heat treatment followed by rapid cooling was tried on coupons from eight heats of forged bars. The basket-weave structure resulting from the quench did not show agglomeration of intermetallic compounds. These coupons showed uniformly good corrosion resistance with weight gains very close to the 28 mg per sq dm expected of Zircaloy-2 after the 14-day steam test. A large forging was solution heat-treated and tested for structure, mechanical properties, and corrosion resistance. Corrosion resistance was excellent and uniform throughout the section. The normally anisotropic mechanical properties were changed to completely iso-tropic. Strength levels were raised while there was a small loss in ductility at the service temperature. Tempering of the quenched structure below the ß transus did not improve the ductility at the test temperature of 600°F. ZIRCALOY-2 end caps for nuclear fuel elements Zircaloy-2 is an alloy of 1.5 pct Sn, 0.1 pct Fe, 0.1 pct Cr, and 0.05 pct Ni, balance hafnium-free zirconium patented by Westing-house Electric Co. used in pressurized water reactors are machined from large hammer-finished forged bars. The bars forged for this application range from 3 1/2 in. to 7 1/2 in. squares. One of the tests used to evaluate Zircaloy-2 for this service is a 14-day exposure to very high-purity steam at 75oF and 1500 psi pressure. Normally wrought products will exhibit a black lustrous film of zirconium oxide after the test and show a weight gain of approximately 28 mg per sq dm. Coupons representing the forged bars exhibited a wide variety of corrosion results ranging from acceptable black coupons to some covered with white crystals of zirconium oxide that had weight gains of nearly 100 mg per sq dm. One of the most common effects was a white corrosion product that looked like a stain to some observers and the outline of massive metal grains to others. Very careful specimen preparation prior to the corrosion test did not affect the result and metallographic examination did not reveal a structural feature of a size and shape that would correlate with the corrosion product. In addition to the relatively poor corrosion resistance, the forged product was not as ductile as desired. The problem of obtaining uniform and more acceptable properties in these heavy-forged sections seemed to be a function of the microstructure of the metal and, in particular, the distribution of the intermetallic compounds. None of the four alloying elements however are appreciably soluble in a zirconium and are present as compounds which usually appear scattered randomly throughout the structure. Moudryl showed that the distribution of these intermetallic compounds in Zircaloy could be related to a ostringero corrosion failure. Grozier et al. discussed another structural defect that causes a similar corrosion effect. In this case, an elongated gas-void was determined to be the cause. M. L. Picklesimer held in his discussion of the Grozier paper that at least a part of the observed =stringersn were caused by the distribution of the intermetallic compounds. He proposed a

Jan 1, 1962

By R. D. Young, E. A. Harre

Nitrogen exists in two broad categories commonly designated as elemental nitrogen and fixed nitrogen. Elemental nitrogen is found in nature as a diatomic molecule and constitutes about 78%, by volume, of the earth's atmosphere. Fixed nitrogen is a common term used for nitrogen that is chemically bound to other elements. Although elemental nitrogen is very abundant, nonleguminous plants cannot utilize nitrogen unless it is in the chemically combined form. The content of this chapter is confined to fixed nitrogen since the uses and properties of elemental nitrogen and the fixed nitrogen compounds are different. Certain bacteria associated with legumes are able to fix atmospheric nitrogen, but not to the extent required by succeeding nonleguminous crops. The need to use fertilizers containing fixed nitrogen compounds to grow food crops has been recognized for hundreds of years, and from antiquity farmers have used organic wastes such as animal manures, vegetable wastes, sewage sludge, fish scraps, etc., containing these compounds for fertilizers. These organic materials have a low nitrogen content, however, and cannot support large-scale food production. Known commercial-size deposits of natural nitrates containing much higher nitrogen content exist only in Chile, and the possibility of the discovery of other significant deposits is remote. The Chilean deposits were the major source of nitrogenous fertilizers until the Haber-Bosch process for the production of synthetic ammonia was developed. In the Haber-Bosch process, nitrogen reacts with hydrogen at high temperatures and pressures and in the presence of a catalyst to form ammonia. Modern plants use air, water, and natural gas or other hydrogen sources to produce the ammonia that is now the key building block of the fixed nitrogen industry. Several major problems that developed in the 1970s continue with the nitrogen industry into the 1980s. These are (1) the continuing sharp escalation in the cost of feedstocks for ammonia production as well as all energy sources, and (2) the very steep increase in the investment cost of production facilities. More than 75% of the fixed nitrogen produced domestically is consumed as fertilizers. The types of materials used as fertilizers have changed significantly during the past years. There has been a trend from the use of low-analysis fertilizers to types which contain high percentages of plant nutrients. The average plant nutrient content (N + P2O5 + K2O) of fertilizers used in the United States has increased from 18% in 1930 to 32% in 1960 and today approaches an average of 44%. In the case of nitrogen fertilizers, the low analysis natural nitrates (16% N) gave way to ammonium sulfate (20% N) after about 1920. After World War II, ammonium nitrate (33.5% N) previously used in munitions, became established as the leading nitrogen fertilizer. Urea use as a fertilizer (46% N) began in the 1960s, and in the past 10 years has become the world's leading solid nitrogen fertilizer. Nitrogen Materials Anhydrous Ammonia-Anhydrous ammonia is a colorless gas at ambient temperatures and pressures and has a sharp pungent odor. Anhydrous ammonia is stored and transported as a liquid in high-pressure cylinders and tanks at ambient temperature, and in low-pressure tanks, barges, and ships at -28°F. Anhydrous ammonia is applied directly at 4 to 7-in. depths in the soil as a high-pressure liquid by special applicators, after which it expands to a gas that dissolves in the soil moisture or is absorbed on soil colloids. In other systems ammonia is metered into irrigation water and is thereby introduced into the soil.

Jan 1, 1983

By C. Scala, M. Gondouin

Published laboratory data have established that very significant streaming potentials can exist across mud cakes subjected to pressure differentials such as exist between a mud column and formation fluids. Since shale is similar in kind and in properties to a well-formed mud cake, it appeared possible that streaming potentials might also arise across shales. A practical consequence would be that the streaming potential contribution to the SP would be smaller than might otherwise be assumed since the SP is measured with respect to the "shale bare line". The resultant contribution would be the difference between the streaming potential across the shale and that across the mud cake. In order to investigate this possibility, a laboratory apparatus was devised for measuring streaming potentials under pressure conditions normally encountered in wells. Reported data on streaming potentids across mud cakes were verified. More important, appreciable streaming potentials were found across shale samples. As in the case of mud cakes, the shale potentials varied with the concentration and nature of the impregnating electrolyte. Dependence on pressure was In a different rela-tion than that for mud cakes, possibly because of the lower compressibility of shales. The differences in streaming potentials between the shales md mud cakes were generally small, although they might vary considerably in practice depending on the relative electro-kinetic properties of the particular shale-mud cuke combination. A "poorly eflective" shale with a "very eflec-tive mud cake" would produce the largest negative addition to the electrochemical component of the SP. Results of field experiments with pressurized wells were considered in the light of the laboratory results. INTRODUCTION Field experiments wherein the pressure of the mud column in a borehole was varied have often indicated the presence of a streaming potential, but its magnitude was not easily predictable and its sign sometimes changed from one formation to another in the same well.' In the interpretation of SP logs from standard surveys, nevertheless, the correct value of R,, is very frequently obtained, within the limits of field accuracy, using only the simple electrochemical expression of the SP provided that proper corrections are made for bed thickness, ion activity, etc.' Accordingly, it was apparent that the streaming potential component of the SP in a borehole was considerably more complex than commonly assumed and clarification was indeed necessary The close similarity in nature and constitution of a shale to a well-compacted mud cake suggested that a streaming potential might also exist across shales as pointed out by M. P. Tixier in 1953 and as recently emphasized by Wyllie." e practical consequence would be that the streaming potential across the mud cake would be canceled, at least in part, on electrical logs since the SP is recorded with respect to the shale potential as a base line. This appeared to be a plausible explanation for many of the effects observed but no supporting data were available. In 1955 preliminary experiments at low pressures were conducted in our laboratories which supported the hypothesis of the shale streaming potential and formed the basis for later studies at borehole pressures.* EARLY LOW PRESSURE EXPERIMENTS The early experiments were conducted primarily to determine the possibility of a streaming potential across shales. The shale used was a compact type designated as No. 7 from Reagan County, Tex. (Permian). Cut samples were saturated for over one month in .2 to 2.0-N NaCl solutions. In order to eliminate all electrode polarization effects, a dynamic method was used. A transient pressure of about 1 atm was applied across

By C. S. Samis, J. C. Mitchell

Sodion activity data have been obtained for the Na-A1 system for sodium concentrations in aluminum between NN, = 25 . 10-and 300 x 10'. This concentration intert-a2 encompasses those sodium concentrations nor1)zally encountered in aluminum reduction cell operation. The sodium activity data were measured by a two-phase equilibration of sodium between a lead and an aluminam rich phase. The Na-Pb system was emplojed as a reference for the determination of the actiz7itj 0.f sodium. The activity of sodiin in aluminuwz exhibits large positi2.e deviations frorn ideality and a strong dependence, on concentration below NN, = 100 x 1K6. The inteated data indicate that the activity of aluminium is equal to its mole fraction over the range of miscibility of the Na-A1 sgstem. The data on the Na-A1 binary have been combined with experimental data on the cryolite-A1 system at 1025°C to obtain the activities of Na F and AlF, in liquid cryolite oz2er the cornposition 1-ange encountered in aluminum cell operation. The results are in agreement with activities calculated from the NaF-AlF, phase diagram. ALUMINUM is produced commercially by the electrolytic reduction of aluminum oxide dissolved in molten cryolite. There is considerable controversy over the actual mechanism of the discharge processes, but sodium potentials in the systems appear to play an important role in the reduction mechanism. The metal-salt reaction involving NaF and aluminum to produce A1F3 and sodium is responsible for the presence of sodium in the aluminum metal phase. Additional sodium may be produced from the simultaneous deposition of sodium and aluminum at the cathode as a result of a cathodic overvoltage. runert' has suggested that sodium gas at atmospheric pressure would be the result of reaction between cryolite and aluminum at 1000°C. Subsequent investigations have shown that production of sodium at atmospheric pressure does not occur. These later investigations can be conveniently separated into two categories. The first involves sodium, aluminum, and cryolite, in which the sodium activity was measured. The second includes a lead phase for sodium activity measurement. Jander and errmann,' studied the reaction: at 1090°C in alumina crucibles under an atmosphere of dry, oxygen-free nitrogen. The equilibrium constant for the reaction was expressed in terms of concentration. Their results in Fig. 1 show an increase in sodium content of the aluminum with increasing NaF/A1F3 weight ratio. Pearson and waddington3 performed similar experiments at 1000°C. Their results are included in Fig. 1. Dewing and ollinshead have measured the sodium content of aluminum as a function of the NF/AF, ratio at 1025"C, Fig. 1, and ewin' measured the sodium content of the aluminum in a system involving the solid reactants NaF, cryolite, and liquid aluminum. This reaction is different from those discussed above. Two solid phases, NaF and Na3A1F6, and a liquid aluminum phase were present. The sodium concentration in the aluminum was determined by the temperature. Dewing's measurements covered the temperature range 679" to 896°C. Omo, Matsushima, and 1t02" investigated the equilibrium 6NaF ential thermal analysis. Stokes and Frank6 determined sodium activities by

Jan 1, 1970

By C. G. Dunn, W. W. Martin

Various X ray diffraction techniques have been developed for determining orientations of crystals. Transmission and back-reflection Laue methods1,2,3 in particular have been found to be very useful. In many applications, however, a need still exists for a simplification of the details of analysis so that orientations can be determined more rapidly, and operators with little or no training in crystallography can be easily instructed in the necessary procedures. Filling this need in part is a system of analysis that has been in use for several years in the Metals Section of the Laboratory, General Electric Co., Pittsfield, Mass. This method is an extension of the one of Majima and Togino,4 which uses fifty-five indexed standard transmission Lauemams to cover the entire range of crystal orientations. Greater rapidity in determining the orientation of a crystal is achieved by using a slightly larger set of standard Lauegrams and by systematizing the steps followed. This system of analysis in its present form constitutes the subject matter of the present paper. Although the method uses transmission X ray Laue patterns, many of the operations are applicable to the analysis of back-reflection Lauegrams. Very briefly the method is based on the fact that Bragg and azimuthal angles for two or more identified Laue diffraction spots completely determine the orientation of a crystal, and standard Laue photographs of the Majima and Togino type provide a ready means for identifying Laue spots. After a preliminary analysis, which consists of identifying the pattern in terms of one of the standard patterns, a correctly marked Lauegram is selected to carry out the final analysis. From measured Bragg and azimuthal angles, the normals (poles) of the reflecting planes for the selected Laue spots are located on a stereographic projection and the orient.at.ion of the crystal determined (in terms of the positions of cube poles for example) by a few additional operations on a \I1ulff net. The actual procedure followed will be illustrated in detail after a brief treatment or X ray problems connected with orientation determinations and a discussion of useful principles and techniques. X Ray and Plotting Problems RELATIVE ADVANTAGES OF BACK-REFLECTION AND TRANSMISSION LAUE DIFFRACTION METHODS In the back-reflection Laue method it is easy to make a direct measurement of the tilt of a zone of Laue spots. Specimen thickness also is generally no problem. In the case of thin specimens, however, the time of exposure is less for the transmission Laue method. For example, with silicon ferrite this advantage holds up to a thickness of 0.6 mm if tungsten radiation at 40 kv peak voltage is used. Furthermore, the transmission photograph represents the entire thickness of the crystal whereas the back-reflection picture represents only a surface layer (a thickness of less than 0.05 mm for silicon ferrite), the deeper layers supplying a negligible amount to the diffraction pattern. If a crystal is distorted, the transmission Laue pattern gives more information particularly with regard to the orientation of the crystal and the overall amount of spread in orientation. Finally a 3 1/4 X 4 in. film suffices for transmission photograms; a 4 X 5 in. or 5 X 7 in. film with a central hole is needed for back-reflection patterns Not only is the manipulation of film easier because of the size of film, but

Jan 1, 1950

By M. Metzger, T. A. Read

The strengthening effects of hydroxide and synthetic plastic films in the creep of cadmium crystals were studied. The results were broadly consistent with the naive mechanical model. The dislocation-film reactions and the significance of the macroscopic stresses associated with certain films are discussed in the paper. ALTHOUGH a number of studies have been made of the strengthening of metal single crystals produced by surface films usually less than Ir in thickness, the mechanisms underlying this phenomenon are still not well established. Roscoe' first noted the critical shear stress and the stress-strain curve of cadmium-crystal wires to be raised markedly by the presence of oxide films. Subsequent investigators2-' have confirmed the effect for several metals and films including a number of plated poly-crystalline metal films. Harper and Cottrel1" with a thin film on zinc presumed to have formed as a result of etching and washing treatments) did not detect a strengthening at a rate of extension of 10 ' sec-' although at lo-' sec-' the film raised the flow stress by 7 pct at 0.0005 extension and by 22 pct at 0.01. The marked strengthening at moderate strains thus appeared to be the result of a dynamical effect dependent on the rate and amount of plastic flow. Other observations'. " are consistent with this view. The large effects of films on the tensile creep rate of single crystals were demonstrated when the observation of Andrade and Randall," that the creep rate of cadmium crystals increased on immersion in aaueous solutions of certain cadmium salts. was shown' * to hold only when a film was present initially and was explained by the attack of the solution on the film and the removal of its strengthening effect. Changes in the creep rate interpretable in terms of the removal or formation of various films have been observed by other investigators.3,8,10,31 Certain other experiments have been reported in which surface films or liquid environments have affected the plastic behavior of metals, but these are not within the scope of the present work which deals with the effect of solid films on single-crystal slip. Possible Mechanisms—The film does not significantly reduce the load on the crystal in the elastic range (the elastic moduli of film and crystal are not too greatly different and the area of the film is so small that its share of the load is negligible) so that the initial operation of internal Frank-Read sources is not affected. However, as Fisher" has pointed out, a source terminating on the crystal surface requires only half the critical stress for activation as an internal source of the same length and a relatively strong film (even a very thin one) would pin the free ends of the surface sources so that, if surface sources are responsible for plastic flow in a clean crystal. such a film would increase by roughly a factor of 2 the stress necessary for generation of dislocation loops at a given rate. The other general possibility is that the film resists the passage to the crystal surface of internally generated dislocations. If the film is adherent and if the specimen is to exhibit slip of the same general character as a film-free specimen, the film must deform or fail in shear at the active slip planes. The naive mechanical model, where shear-stress relaxation across the active slip plane in the metal transfers load to the film until it deforms and slip can

Jan 1, 1959

By M. King Hubbert, David G. Willis

A theoretical examination of the fracturing of rocks by means of pressure applied in boreholes leads to the conclusion that, regardless of whether the fracturing fluid be of the penetrating or non-penetrating type, the fractures produced should be approximately perpendicular to the axis of least stress. The general state of stress underground is that in which the three principal stresses are unequal. For tectonically relaxed areas characterized by normal faulting, the least stress should be horizontal; the fractures produced should be vertical with the injection pressure less than that of the overburden. In areas of active tectonic compression, the least stress should be vertical and equal to the pressure of the overburden; the fractures should be horizontal with injection pressures equal to or greater than the pressure of the overburden. Horizontal fractures cannot be produced by hydraulic pressures less than the total pressure of the overburden. These conclusions are compatible with field experience in fracturing and with the results of laboratory experimentation. INTRODUCTION The hydraulic-fracturing technique of well stimulation is one of the major developments in petroleum engineering of the last decade. The technique was introduced to the petroleum industry in a paper by J. B. Clark,' of the Stanolind Oil and Gas Co. in 1948, and since then its use has progressively expanded so that by the end of 1955 more than 100,000 individual treatments had been performed. The technique itself is mechanically related to three other phenomena concerning which ah extensive literature had previously developed. These are: (1) pressure parting in water injection wells in secondary-recovery operations, (2) lost circulation during drilling, and (3) the breakdown of formations during squeeze-cementing operations, all of which appear to involve the formation of open fractures by pressure applied in a wellbore. The most popular interpretation of this mechanism has been that the pressure had parted the formation along a bedding plane and lifted the overburden, notwithstanding the fact that in the great majority of cases where pressures were known they were significantly less than those due to the total weight of the overburden as determined from its density. Prior to 1948, this prevalent opinion had already been queried by Dickey and Andresen,2 in a study of pressure parting; and by Walker,"' who, in studies of squeeze cementing, pointed out that the pressures required were mostly less than those of the overburden, and inferred that the fractures should be vertical. J. B. Clark,' in his paper introducing hydraulic fracturing, and later Howard and Fast," and Scott, Bearden, and Howard," also of Stanolind, postulated that the entire weight of the overburden need not be lifted in producing horizontal fractures, but that it was only necessary to- lift an "effective overburden," requiring a correspondingly lower pressure. Hubbert,7 in discussing the paper by Scott and associates, pointed out that the normal state of stress underground is one of unequal principal stresses; and in tectonically relaxed areas, characterized by normal faults, the least stress should be horizontal. Therefore, in most cases, fracturing should be possible with pressure less than that of the overburden and, moreover, such fractures should be vertical. Harrison, Kieschnick, and McGuire,8 also on the expectation that the least principal stress should be horizontal, argued strongly in favor of vertical fracturing. Scott, Bearden, and Howard6 observed that, when

Jan 1, 1958

By M. Sheffield

This paper presents a technique lor predicting the flow of oil, gas and water through a petroleum reservoir. Gravitational, viscous arid capillary lorces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. INTRODUCTION Since the introduction of high speed computing equipment, one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance, 193-5,7,8 The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. DISCUSSION The method of predicting three-phase compressible fluid flow. in a reservoir can be summarized briefly by the following steps. 1. The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. 2. Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations.1,4,8 The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary forces are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and differenceequations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method of solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. The section

Jan 1, 1970

By R. G. Davies, R. H. Richman

Several techniques, X-ray diffraction, transmission electron microscopy, conventional, metallography, and mechanical property measurements are combined to investigate decomposition in the suPersaturated solid solution Fe-27 at. pct Be. The following sequence of structures for aging temperatures below 400°C is established: massive short-range order, pseudo-FeBe (L2o) order, Fe3Be (DU3) order, sinusoidal modulation of composition, general or continuous precipitation, and discontinuous or cellular precipitation. In addition to this sequence itself, several features of the individual stages are noteworthy; for example, when coherent modulations occur, their development is parallel to (100) everywhere in the specimen apparently without incubation, and both the wavelength and the composition amplitude of the modulations increase with aging time. Two types of discontinuous reaction are observed, depending upon the aging temperature. At low temperatures (<550°C) the boundary-nucleated variety of discontinuous precipitate predominates; but at higher temperatures (>600°C) the crystallographic variety occurs instead. Hardness and flow-stress measurements clearly reflect the quoted structure sequence. Some qualitative aspects of the structure changes and of the relationship between mechanical properties and structure are discussed. BERYLLIUM in solid solution markedly strengthens iron.&apos;" A recent study demonstrated that not only were supersaturated single crystals of Fe-27 at. pct Be very strong, but they deformed primarily by mechanical twinning.3 Such behavior, coupled with the decrease in solid solubility of beryllium in a Fe from about 33 at. pct at 1165°C to about 10 at. pct at 500°C4 as shown in Fig. 1, leads naturally to questions about property changes with aging in this system. In 1929 roll,1 and later Seljesater and Rogers,5 reported that useful increases in strength could be obtained by aging quenched Fe-Be alloys, but it was not until quite recently that any of the details of precipitation in these alloys were studied. Richman6 observed that both an Fe3Be superlattice (with a DO3 crystal structure the same as Fe3A1, illustrated in Fig. 2) and a modulated structure preceded the appearance of a discrete second phase in the 20 at. pct Be composition. Here, a modulated structure refers to a periodic composition fluctuation within one crystal, while its lattice remains coherent. Subsequent work on a Fe-Be alloys by Heubner7 and von Heimendahl and Heubner8 confirmed the existence of the modulated structure, although only an FeBe (L20) ordered structure was detected (also shown in Fig. 2). Both Rich-man and Heubner corroborated the large increases in hardness during aging reported earlier, but other mechanical properties were not examined. Since modulated systems are encountered relatively infrequently, the influence of coherent modulation upon mechanical properties is not well-known. With the accumulation of information about age-hardening systems, it has become clear that many precipitation reactions cannot always be described by the concept of heterogeneous nucleation at some lattice discontinuity followed by diffusional growth. Rather, complex sequences involving prior "preparation" of the matrix are quite common; that is, a transition structure occurs in the parent phase, instead of in some "nucleus" form. Delineation of such subtle sequences appears to be a worthwhile step toward the understanding of solid-state precipitation. Since not many structural progressions of comparable nature have been reported, the objectives of the present survey are to identify and elucidate the various

Jan 1, 1967

By R. Ehrlich, F. E. Crane

A consolidated porous medium is mathematically modeled by networks of irregularly shaped, interconnected pore channels. Mechanisms are described that form residual saturations during immiscible displacement both by entire pore channels being bypassed and by fluids being isolated by the movement of an interface within individual pore channels. This latter mechanism is shown to depend strongly on pore channel irregularity. Together, these mechanisms provide an explanation for the drainage-imbibition-hysteresis effect. The calculation of steady-state relative permeabilities, based on a pore-size distribution obtained from a Berea sandstone, is described. These relative permeability curves agree qualitatively with curves that are generally accepted to be typical for highly consolidated materials. In situations where interfacial effects predominate over viscous and gravitational effects, the following conclusions are reached. 1. Relative permeability at a given saturation is everywhere independent of flow rate. 2. Relative permeability is independent of viscosity ratio everywhere except at very low values of wetting-phase relative permeability. 3. Irreducible wetting-phase saturation following steady-state drainage decreases with increasing ratio of nonwetting- to wetting-phase viscosity. 4. lrreducible wetting-phase saturation following unsteady-state drainage is lower than fir steady-state drainage. 5. Irreducible nonwetting-phase saturation 101lowing imbibition is independent of viscosity ratio, whether or not the imbibition is carried out under steady- or unsteady-state conditions. Experiments qualithtiveIy verify the conclusions regarding unsteady-state residual wetting-phase saturation. Implications of these conclusions are discussed. INTRODUCTION Natural and artificial porous materials are generally composed of matrix substance brought together in a more or less random manner. This leads to the creation of a network of interconnected pore spaces of highly irregular shape. Since the geometry of such a network is impossible to describe, we can never obtain a complete description of its flow behavior. We can, however, abstract those properties of the porous medium pertinent to the type of flow under consideration, and thus obtain an adequate description of that flow. Thus, the Kozeny-Carmen equation, by considering a porous medium as a bundle of noninterconnecting capillary tubes, provides an adequate description of single-phase flow. With the addition of a saturation-dependent tortuosity parameter in two-phase flow to account for flow path elongation, the Kozeny-Carmen equation has been used to predict relative permeabilities for the displacement of a wetting liquid by a gas. 1-5 It has long been recognized that relative permeability depends not only on saturation but also on saturation history as we11.6,7 Naar and Henderson8 described a mathematical model in which differences between drainage and imbibition behavior are explained in terms of a bypassing mechanism by which oil is trapped during imbibition. Fatt9 proposed a model for a porous medium that consisted of regular networks of cylindrical tubes of randomly distributed radii. From this he calculated the drainage relative permeability curves. Moore and slobod,10 Rose and Witherspoon,11 and Rose and Cleary12 each considered flow in a pore doublet (a parallel arrangement of a small and large diameter cylindrical capillary tube). They concluded that, because of the different rates of flow in each tube, trapping would occur in one of the tubes; the extent of which would depend upon viscosity ratio and capillary pressure. They limited

Jan 1, 1970