Search Documents

By E. A. Gulbransen, K. F. Andrew

THIS paper. will present the results of our studies on the kinetics of the gas phase reactions of co-lumbium and tantalum with O2, N2 and H2. Studies on zirconium and titanium have been previously reported.9,11 Due to their high melting points and their many excellent physical and chemical properties, colurn-bium and tantalum have been finding their places as useful metals in science and industry. Chemically the metals are very inert to gas and liquid phase corrosion at room temperature. However, at moderate temperatures of 250°C and higher the metals react readily with oxygen and hydrogen. It is of interest to study the gas phase reactions of these metals for the following reasons: (1) to compare the rates of the several reactions with those of other metals and to correlate the rate data with theory and structure predictions, (2) to add basic information on the practical difficulties of the reduction, refining and working of columbium and tantalum, (3) to anticipate new uses for the metals and (4) to understand the role that the metals play as components in high temperature alloys. Literature Survey A number of papers exist in the literature on the gas phase reactions of columbium and tantalum with O2, N2 and H2; however, only fragmentary data have been reported on the kinetics of these reactions. Although the free energy of formation of the oxide Ta2O520at 25°C has been determined and some thermodynamic work has been reported on columbium and its oxides, the data are not sufficient to make equilibria calculations at elevated temperatures. Oxidation: 1. Kinetics: The rate of oxidation of columbium and tantalum has been studied by McAdam and Geil23 using the interference color method. Tantalum is reported to oxidize more slowly than columbium, zirconium and iron. The amounts of oxidation of columbium after 20 hr in air for several temperatures are reported by the Fansteel Metallurgical Corporation.36 Columbium starts to oxidize in air at about 200°C. The oxide is said to be adherent and prevents further oxidation until the temperature is raised. Columbium is reported not to become brittle on heating in air for short periods like tantalum because the oxide film prevents further reaction. The columbium oxides dissolve into the metal when heated in a vacuum at temperatures between red heat and 1200°C. Above 1200°C the oxides are reported to evaporate. Columbium, in the oxide free state, is an active getter in vacuum tubes. 2. Preparation and Structure: G. Grube, O. Kuba-schewski and K. Zwiauer15 have studied the decomposition of Cb2O5 in argon. They find that CbO2 is formed slowly at 1150° and rapidly at 1300°C. Cb2O5 can be reduced to CbO2 by moist H2. The reduction of this dioxide to the lower oxides is rapid. By varying conditions CbO2, Cb2O3, CbO and Cb2O can be formed. CbO2, CbO and Cb2O have independent lattices while Cb2O3 is a stoichiometric mixture of CbO and CbO2. In a later paper O. Kuba-schewski" has shown that CbO is cubic and that CbO2 has a rutile type of lattice. The most exhaustive work on the oxides of columbium is that of G. Brauer.8 The number and constitution of the solid phases in the binary system Cb-O were investigated by the preparative, analytical and X ray methods. He indicates that only the oxides Cb2O5, CbO2 and CbO exist with limited regions of homogeneity. The solid oxide Cb,O, exists

Jan 1, 1951

By M. A. Dayananda, R. E. Grace

Vapor-solid diffusion couples were employed in a study of ternary diffusion in the single-phase copper-rich corner of the Cu-Zn-Mn system Interdiffusion coefficients were measured at three different compositions corresponding to the intersections of diffusion paths determined at 850°C. Intrinsic diffusion coefficients were determined at two compositions by the Philibert-Guy method. New relations between the various intrinsic diffusion coefficients are also presented. MOST of the investigations of ternary diffusion have been made with "sandwich" couples formed by welding together two homogeneous solid solutions with or without markers placed at the interface. With such couples Kirkaldy and his coworkers1,2 have determined interdiffusion coefficients which describe ternary diffusion in the substitutional alloy systems, Cu-Mn-A1 and Zn-Cu-Al. Guy and Leroy3 have reported intrinsic diffusion coefficients for Co-Ni-Cr alloys as a function of composition using the analysis of Philibert and Guy.4 The use of semi-infinite diffusion couples, for example, vapor-solid couples, where the diffusing species are supplied from a vapor source in contact with a diffusion disk, has so far been confined to studies of binary diffusion. Balluffi et a1.5,6 have successfully determined the interdiffusion as well as the intrinsic diffusion coefficients in several binary metallic systems employing vapor-solid couples with markers initially placed at the vapor-solid interface. They have also discussed the general advantages and limitations of such couples. Vapor-solid diffusion couples can successfully be employed in ternary diffusion studies. The importance of such couples in ternary diffusion lies in the fact that inert markers lie close to the vapor-solid interface, and the composition of the marker plane is virtually known before experimentation begins. Another minor advantage is that the method of Philibert and Guy4 may be used with three alloys instead of four. Lastly, the diffusion paths need not be S-shaped, as pointed out by Meijering7 for "sandwich" couples. The Cu-Zn-Mn system is convenient for ternary diffusion studies because of the large single-phase field in the copper-rich corner. The vapor pressure of zinc above such alloys is known, and X-ray measurements indicate that the molar volume is nearly constant over a wide range of compositions.8,9 It is the purpose of this paper a) to extend the analysis of Balluffi and seigle5 to permit an evaluation of interdiffusion coefficients in ternary vapor-solid couples, b) to use the method of Philibert and Guy to determine intrinsic diffusion coefficients with such couples, and c) to present experimental data on both interdiffusion and intrinsic diffusion coefficients for the Cu-Zn-Mn system at 850°C. INTERDIFFUSION IN TERNARY VAPOR-SOLID COUPLES The analysis of Balluffi and seigle5 is extended in this section to a vapor-solid couple where two solutes diffuse from a vapor source into a ternary solid solution. To take into account the expansion of the diffusion disk caused by the absorption of mass from the vapor phase, the total fluxes are expressed on the basis of a coordinate system that moves with the velocity of the vapor-solid interface; thus. where D3ij are the intrinsic diffusion coefficients,4 v(x, t) is the local Kirkendall velocity of mass flow, and F(t) is the negative velocity of the vapor-solid interface (x = 0). Both these velocities are based on a fixed coordinate system far removed from the diffusion zone. Using the continuity relations and Darken's method1' for a ternary system in which components 1 and 2 diffuse from the vapor phase into a solid of negligible vapor pressure, one obtains

Jan 1, 1965

By W. A. Backofen, A. T. English

The kinetics of re crystallization were determined metallographically for a 3-1/4 pcl Si-Fe rapidly compressed at temperatures of 710° to 911°C, and held for various times at the working temperature. As in re crystallization following cold wwk, the time for re crystallization was reduced by increased temperature and strain. Nucleation occurred chiefly at grain edges, although other pre-existing interfaces, notably inclusions, also served as nucleation sites. The importance of intragranular nucleation to grain refinement in hot working is pointed out, together with some possible techniques by which nucleation might be induced within the grains. The observed isothermal growth rates are quite accurately described by the equation G = 4 x 10-'/t (cm per sec), independent of temperature and strain at all but the shortest times. Changes in G during recrystalliza-tion can exceed two orders of magnitude. It is argued that release of stored energy associated with lattice defects cannot adequately explain the time dependence. A suggested alternative involves processes of solute segregation and desegregation in the unrecrystallized matrix and at the moving boundary, respectively. EFFORTS to improve hot-working practice are hampered by an incomplete understanding of the fundamental metallurgy involved. Research directed at establishing the required basic foundation falls into two categories, the first concerned with measurement of mechanical properties under the conditions of high temperature and strain rate which characterize such processing, the second focused on structure and its relation to the deformation conditions. Rossard,' Rossard and Blain,2 and Hardwick and Tegart have conducted many high-speed, high-temperature torsion tests on various metals. One of their objectives was interpretation of stress-strain curves in terms of the structures of specimens rapidly quenched following interruption of twisting. They showed that appreciable hardening occurred in the early part of the test, and that this hardening reflected the concurrent development of an increasingly dense substructure within the grains. The detailed nature of this substructure was observed to depend on the ease of polygonization which, at least in fcc metals, varies with the stacking-fault energy.= By holding specimens after interruption but before quenching, the nucleation and growth of new grains could also be studied. Some success has been realized in relating the microstructure of hot-rolled steel to its temperature:strain rate:strain history with such testing.' The same general observations have been made by others investigating rolling4f5 and drop-hammer forging.~-~ Incident to these experiments, some efforts were made to determine the kinetics of re-crystallization as well as the velocity of the migrating boundaries. The results have been largely qualitative, owing to difficulties of interpreting hardness measurements5 and microstructures. The work described here was undertaken to obtain more quantitative information from tests on a material having suitable metallographic properties, and, in this way, to provide further insight into the fundamentals of this practically most important kind of deformation processing. I) MATERIALS AND EXPERIMENTAL TECHNIQUES Materials. Electrolytic etching of a Si-Fe alloy reveals dislocations and substructure in unrecrystallized grains.8'9 This unique characteristic makes possible a precise determination of the extent of recrystallization. For that reason, it was the material selected for all experiments in the present work. The two particular alloys that were used are identified in Table I. A major difference between "A" and "B" is the number of inclusions, as noted in the table. Most of the experiments were made with material "A". It was received as hot-rolled plate 1 in. thick, sawed to 2-in. widths, hot-forged at 1200°C to 5/8-in. rounds, and finally annealed in a mechanical vacuum at 1000°C for 24 hr. The result was substructure-free grains of mean diameter 0.05 cm.*

Jan 1, 1964

By Seymour G. Epstein

Using a modified capillary-reservoir technique, electromigration of dilute solutions of silver in liquid bismuth has been measured with relatively good precision. The effects of experiment duration, temperature. current density, and silver concentration have been investigated. Chemical diffusion coefficients of silver in liquid bismuth were also measured at temperatures ranging from 300°to 600°C. The effective valence of silver in bismuth was calculated as a function of temperature from a comparison of electric mobilities and diffusion coefficients. In an alloy containing 1 wt pct Ag, the electric mobility of silver increased from 1.65 x X sq cm per v-sec at 300°C to 2.38 x 10-4at 600°C; the diffusion coefficient in-creased from 3.0 to 17.5 x 10-5 sq cm per sec over the same temperature range. Calculated values of effective valence decreased from 0.28 at 300°C to 0.10 at 600°c. ALTHOUGH there has been considerable progress in the theoretical treatment of electron-transport properties in liquid metals, the mechanisms of mass-transport phenomena in liquid metals are still not clearly understood. This is especially true for the case of electro migration, alternatively termed elec-trotransport or electrodiffusion, which is mass transport induced by an applied electric field. A thorough knowledge of the mechanisms governing electromigration will considerably aid in understanding the electronic state of an impurity atom in a liquid metal. Electromigration has been long observed and studied in liquid metals and alloys. However, no satisfactory theory to explain the observed effects has been developed. The disparity in the present theoretical treatments, summarized by Verhoeven,1 has been attributed to the difficulty in describing the momentum transfer between electrons and ions and to the lack of precise measurements in liquid metals, especially at elevated temperatures. Prior to the last decade most of the electromigration studies were qualitative, in that only the direction of migration of alloy constituents was determined with any degree of reliability. These results have been summarized by schwarz,2 who was able to measure the electric mobility of several elements in mercury. More recently electromigration in liquid metals has been measured by a variety of techniques, with a varying degree of success. Mangelsdorf3,4 made precise measurements on dilute amalgams at room temperature, using a technique relating changes in composition to changes in alloy resistivity, but was unable to make measurements at elevated temperatures. Using a modification of this technique, Verhoeven and Hucke5 measured electromigration and resistivities in liquid Bi-Sn alloys. However, the resistivity technique is not sensitive to small composition changes and is impractical for very dilute alloys. Verhoeven and Hucke6 also measured electromigration of several solutes in liquid bismuth by adapting the capillary-reservoir method for determining diffusion coefficients in liquids. The reported precision of these measurements was poor, however, primarily due to the large uncertainties in the chemical analyses to determine concentration changes. Belashchenko7 has summarized the recent Russian investigations of electromigration, many of which, however, are of questionable accuracy. The present study was initiated to develop an experimental technique which would be both versatile and sufficiently sensitive to yield precise measurements of electromigration of solutes in dilute metal alloys at elevated temperatures. Using this technique, a comprehensive study of silver in liquid bismuth was made to determine the effects of experimental variables on the apparent electric mobilities. Chemical diffusion coefficients of silver in liquid bismuth were also measured over the same temperature range for comparison with the electromigration data. EXPERIMENTAL TECHNIQUE Electromigration Measurements. Electromigration of silver in liquid bismuth was measured by the capillary-reservoir technique, similar to that used by Verhoeven and Hucke,6 but modified to yield more

Jan 1, 1967

By V. J. Berry, J. D. Pendergrass

Well pressure transient tests provide a means for directly obtaining information about formation pressure and reservoir flow capacity. Such tests have also been proposed for determining presence and location of faults or other reservoir closures and for measuring oil in place. For mathematical convenience, most theoretical studies have considered the reservoirs to be homogeneous. Definitive information is not yet available to show whether the actual presence of nonuniformities will make pressure transient behavior different from that of a uniform reservoir. The conclusions reached from actual transient tests are questionable, therefore, insofar as they rely on the original assumption of homogeneity. One type of nonuniformity commonly assumed to exist is that of stratification. In most reservoirs the strata are thought to be in vertical communication. Equations for the transient flow of a single-phase, compressible fluid in a one-well, bounded, circular reservoir have been solved for several situations involving crossflow between multiple strata of various thicknesses and permeabilities. The results show that except for the very early flow period, which usually is too short to be analyzed, the transient performance observed at the well is substantially identical with that of a homogeneous reservoir having the same dimensions and having the same steady-state flow capacity. Thus, stratification does not adversely affect interpretation of well transient tests. This conclusion holds for all commonly encountered combinations of reservoir thickness and external radius. Deviations are observed for unusually thick reservoirs whose outer radii are relatively small. The results of these studies also show that the presence and the amount of stratification cannot be simply diagnosed from reservoir pressure transient data when there is crossflow between strata. INTRODUCTION The last decade has brought wide acceptance of the transient well pressure test for determining reservoir parameters. Following the original work of Hurst and van Everdingen,1 the mathematical theory was thoroughly explored. Numerous authors have suggested how to determine static reservoir pressure, permeability-thickness product, original oil in place and reservoir limits for different reservoir geometry. The same mathematical techniques have been used to predict the transient performance of a reservoir over a long period of time. Most of the theoretical work has been for homogeneous, isotropic systems. Some results have also been presented for a homogeneous, anisotropic reservoir.2 Petroleum reservoirs are not homogeneous. The deposition process seems to favor creation of a stratified formation. This concept is sufficiently well accepted so that the most natural extension of the transient flow theory beyond the homogeneous case is to a stratified formation. Results for a stratified reservoir with no vertical communication between layers can be obtained from the results for a homogeneous reservoir. Lefkovitz, et a1,3 have given a thorough treatment of the two-layer case. Two recent papers4,5 have treated the case of a two-layered reservoir with vertical communication or crossflow, between layers. Russell and Prats4 find that, after a relatively short time, the two-layered reservoir with crossflow exhibits a simple exponential pressure decline. From this time forward, the behavior is not distinguishable from the behavior of a homogeneous reservoir having the same steady-state flow capacity. The results of Katz and Tek 5 are equivalent. Russell and Prats4 also speculate that a multi-layered reservoir with crossflow will behave as a homogeneous system after long enough production time "... providing the contrast in kh between layers is not too great". They also suggest that, at intermediate times, "... the relative positions of the layers with respect to each other will have a great influence on the production behavior and on the time at which the previously mentioned large-time approximation might be valid". Katz and Tek5 remark upon the mathematical difficulty of treating a reservoir having many layers or strata. If relative position of the layers is also

By Wallace W. Wilson

This DaoerAA calls to attention for the first time, to the writer's knowledge, a purported recovery of oil by substantially horizontal gas cycling that is considerably in excess of 50 per cent of the oil originally in place. The average angle of dip along the path of travel of the gas is less than one degree of arc. The best recovery by gas cycling previously brought to the writer's attention1 occurred in the case of the Mile Six Field in Peru, where about 68 per cent of the original oil in place has so far been recovered. In this case, however. the recovery was greatly aided 1)y the high dip angle of about 17.5". An approximate estimate of recovery was made by the method developed by the writer' for the case of the Cook Ranch Field, using an average relative permeability ratio function dictated by past experience with a number of good quality sandstones. A result of about 40 per cent recovery of the original oil in place wa* obtained after the 8.5 pore volumes of cumulative gas cycling undergone by the Cook Ranch Field. Tile writer desires to emphasize what appears to him as the sound engineering practiced at Cook Ranch, and also the care attending the taking of data and the high quality of the calculations made. It is impossible to take issue on factual grounds with any of the conclusions drawn in the paper. It is, however, desired to point out the difficulty of explaining the unexpectedly favorable oil recovery on ordinary theoretical grounds. It would be of interest to know whether relative permeability ratios have been measured on the Cook sand. The operator is inclined to attribute some of the excess benefit to settling of the oil across the bedding planes into the lower part of the sand. It is difficult for the writer to appreciate the importance of such an effect, inasmuch as the concomitant opening of excessive vacated pore space in the upper part of the sand would he expected to cause excessive bypassing of the gas. REFERENCE 1. Welge. H. J.: "A Simplified Method for Computing Oil Recovery by (Gas or Water Drive." Trans. AIME, (1952), 195, 91. AUTHOR'S REPLY TO MR. WELGE Welge has poined out clearly the difficulty of applying material-balance procedures in calculating the Performance of low-pressure, dispersed gas injection operations. There can be no argument that even with the most favorable assumptions of K,/K,, and with cumulative gas cycling in excess of normal expectations, the calculated oil recovery will be considerably less than the 72.5 per cent recovery realized from this project to Jan. 1, 1951. There probably are at least two reasons for the unusually high recovery from this reservoir, the effects of which would be difficult or impossible to calculate by known methods: 1. The operator did not use the same injection wells for the entire operation, and exercised considerable control over producing GOR by changing the locations of injection wells whenever well performance indicated excessive bypassing. Fifty-six wells were used for injection, not over 29 of which were in use simultaneously. This procedure resulted in frequent shifting of the fluid movement patterns in the reservoir and undoubtedly increased the ultimate recovery. The writer knows of no satisfactory procedure for calculating the effect of this practice. 2. Vertical movement of oil across the bedding planes of the reservoir rock hy gravity drainage tended to maintain a relatively high oil saturation, and hence a favorable K,/K, ratio, in the lower sections of the reservoir. As pointed out by Welge, this condition has aggravated the tendency for high producing GOR, but at the same time it has resulted in increased oil recovery. As may be seen from the performance charts for the various battery areas, operating (GOR'S have been quite high for the last several years of reported history.

Jan 1, 1952

By Nicholas J. Grant, Robert H. Kane, Bill C. Geissen

The platinum-rich section of the Ta-Pt Phase diagram was worked out by metallographic and X-ray techniques, using twenty alloys; solidus temperatures and solid solubilities were determined. Besides s phase, four intermediate stoi-chiometric phases were found: TaPt (high-temperatzrre phase, structure unknown); TaPt2 (ortho-rhomhic, new type); a TaPt3 (orthorhombic, TiCu3 type); ß TaPf3 (monoclinic, new type). TaPt2 and TaPt3 melt congmentlv; two eutectic, two beritec-tic, and one eutectoid reactions occur in the inuestigated range. THE Ta-Pt system was examined as part of a study of the phase diagrams and the intermediate phases of niobium and tantalum with Group VIII transition metals. Since the main scope of this work was the relation among the ordered close-packed phases occurring in the concentration range from AB to AB3 (A = Nb, Ta; B = Pd, Pt, and so forth), only the platinum-rich side of the diagram was worked out thoroughly, while spot checks were made to determine the approximate outline of the tantalum-rich portion. The system has not been thoroughly treated in the literature; Greenfield and Beck1 in their survey of refractory metal intermediate phases found a o phase between -18 and 37 at. pct Pt which coexists at higher platinum concentrations with an unidentified intermediate phase. Browning2 worked on the range between 55 and 100 at. pct Pt; however, this preliminary study did not cover temperatures >1500°C, and included no conclusive crystallo-graphic data. More recently, a study of the platinum-rich section of the Ta-Pt system has been carried out by Ray,3 while Hartley, Parsons, and steedly4 have studied the tantalum-rich section. The crystal structures of TaPt2, a TaPt3, ß TaPt3, and others have been reported in abstracted form by Giessen and Grant.5 EXPERIMENTAL METHODS The methods of alloy preparation and processing and of the metallographic and X-ray procedures follow very closely those described by Ritter, Giessen, and rant' for other refractory metal-noble metal systems, and will not be repeated except where significant differences exist. Starting materials were: tantalum rod of 99.95 pct purity, supplied by National Research Corp., and platinum sponge 99.9 pct, supplied by The International Nickel Co. An oxygen analysis of the tantalum rod showed 30 ppm O; other impurities were C, N, Si, Nb, Mo, none of these exceeding 50 ppm, and 100 ppm W. Homogeneous buttons were readily prepared as follows: ingots of 5 g were weighed and arc-melted under welding grade argon, followed by two remelts after inverting the button. Between 40 and 100 at. pct Pt, weight losses during the first melt rarely exceeded 0.1 pct, in which case the alloy was rejected, and were nil on subsequent remelts. Since the vapor pressure of platinum far exceeds that of tantalum, the small evaporation loss could be assumed to consist of platinum; therefore, a corresponding amount of platinum was added on the second melt. Since it was confirmed on pure tantalum samples that no weight gain through oxygen or impurity pickup occurs (see Ref. 71, it could be assumed that the uncertainty in the composition did not exceed ±0.1 at. pct for the range 40 to 100 at. pct Pt. For the four control melts from 0 to 40 at. pct Pt, an error of ±0.5 at. pct was calculated. Microscopic examination of the whole buttons after melting showed no inhomogeneity other than local coring. This was eliminated by a homogenization heat treatment. A total of twenty-three alloys were prepared. Heat treatments were carried out in a tantalum

Jan 1, 1965

By Robert D. Pehlke, Robert G. Blossey

The solubility of nitrogen in liquid cobalt and a number of dilute cobalt alloys Was measured in the temperature range 1550" to 1700°C at nitrogen pressures to 1 atm. At 1600°C and 1 atm nitrogen pressure, the solubility in pure cobalt is 0.0047 wt pct [3.8 cu cm (stp) per 100 g Co]. The standard free- energy change for the solution reaction 1/2 N,(g) = gout pct in Co) is iF - 10,000 + 5.24 T cal per g atom N. The alloying elements Cv, Cu, Mo, Nb, Tu, W, and V increase the nitrogen solubility, while Al, Fe, Ni, and Si decrease the nitrogen solubility. THE solubility of nitrogen in cobalt was studied previously by Sieverts, Busch and Dodd,2 and Schenck, Frohberg, and Graf.3 Sieverts investigated the solution below the melting point of cobalt, and Busch and Dodd worked with the liquid metal; however, neither study found any solubility for nitrogen in cobalt. Schenck et al, reported a solubility of 0.005 wt pct N in cobalt at 1600°C. Seagle et al. 4 reported residual nitrogen in cobalt after electron-beam melting, although no reasons were ventured for the presence of this interstitial impurity. This investigation was undertaken to clarify these inconsistencies, and to provide additional data on atomic interactions in dilute metallic solution. EXPERIMENTAL PROCEDURE This study was well-suited to the Sieverts' method, which has not changed significantly since its introduction by Sieverts. 5 The apparatus was of the same design as that described by Humbert and Elliott, 6 with a reaction bulb similar to that reported by Weinstein and Elliott. 7 The buret was not water-jacketed, but the measured gas volumes were corrected for small changes in room temperature. The melt temperature was measured by a disappear ing-filament optical pyrometer sighted through a prism onto the melt surface. Temperature calibration was made based on a measurement of the melting point of pure cobalt taken as 1495°C. The difference between the reciprocal of the observed and the true temperatures was assumed to be a constant, neglecting the changes in emissivity with small additions of alloying elements. The total error in gas-volume measurement caused by accuracy limitations in reading the mercury level in the buret, variation in ambient temperature, and variation in melt temperature was about 0.5 cu cm. This represents approximately 10 pct of the total solubility, and is a larger relative error than in previous studies6-' because of the small nitrogen solubilities encountered in the present investigation. The source and purity of the materials used are listed in Table IA. Recrystallized alumina crucibles were used throughout the investigation. Examination of the crucibles after solidification of the melts showed no reaction or penetration. Alloying of the cobalt melt was accomplished with the dropping technique described by Humbert. 6 The investigation was started using cobalt from African Metals Corp. in shot form, but this material was not of sufficient purity and did not yield consistent solubility results. Sherritt Gordon Mines Ltd. cobalt was then used throughout the investigation. In experiments involving alloying elements which could react with oxygen dissolved in the melt, the cobalt was further refined by annealing at 1400°C in a dry-hydrogen atmosphere. This treatment, described by Woolley and Pehlke, 8 markedly reduced the oxygen content. The oxygen and nitrogen analyses of the cobalt used in this study are presented in Table IB. RESULTS Hot Volume. The hot volume was measured prior to each run and was 50 to 55 cu cm (stp) at 1600°C. It increased linearly with increasing pressure and decreasing melt temperature. Pure Cobalt. Measurements of the solubility of nitrogen in pure cobalt and the temperature coefficient of solubility are summarized in Table 11. At 1600°C and 1 atm nitrogen pressure, the solubility was found

Jan 1, 1967

By J. B. Coher, M. B. Bever

COLD work destroys long-range order, as was first observed by Dehlinger and Graf.1 Dahl2 showed that the mechanical disordering caused by cold work produces changes in those properties that are affected by long-range order. The existence of order has a marked effect on tile manner in which an alloy deforms and therefore on its mechanical properties.L' The relation between order and deformation behavior however, is only beginning to be understood. Seemann and G1ander4 suggested that cold work reduces an ordered alloy to a kind of gruel of antiphase domains. According to Cottrell,5 superdis-locations and pairs of dislocations in a slip plane connected by a strip of anti-phase domain boundary provide deformation mechanisms which do not change the state of order. In general, however, slip may be expected to produce new domain boundaries and destroy order if dislocations move on intersecting planes or cut through existing domain boundaries. The yield strength of an ordered alloy may pass through a maximum with decreasing domain size.5"7 The reduction of domain size resulting from cold work may, therefore, lead to a maximum in strength properties. Flinn8 suggests that dislocations can climb by short-range diffusion during deformation and thus can generate domain boundaries not only on slip planes, but also on other planes which will be those with the lowest domain-boundary energy. It has long been known that the rate of strain hardening of an alloy in the ordered state is larger than in the disordered state.3'9 Also the slip bands in ordered CU3, AU10 and Ni3Mn" are finer than in the disordered alloys. The stress which causes initial yielding in Cu,Au passes through a maximum with decreasing domain size.6'12 A maximum has also been observed in the hardness of ordered Ni3Mn as a function of strain.' These observations suggest, but do not prove, that domain size reduction occurs during deformation. According to Fisher13 slip destroys short-range order and thus the strength of an alloy is related to the degree of short-range order present. Rudman and Averbach14 and Averbach et a1.,15 demonstrated that cold work reduces short-range order in copper-gold and gold-silver alloys. Destruction of long-range and short-range order by cold work implies slip on an appreciable fraction of planes. Yielding of an alloy in which short-range order exists or antiphase domain boundaries are present should be discontinuous because the first dislocations moving across the slip plane meet with the largest resistance.5 Such discontinuous yielding has been observed in CU3AU6,16 and Cuzn.l7 This phenomenon, however, can also be explained by mechanisms not involving order.b The published investigations of the effects of cold work on alloys with long-range or short-range order have not been covered thoroughly in reviews of order-disorder phenomena. These investigations, therefore, will be listed here. They have been concerned with strength properties> 2,3,6,8,9,12,16-21 electrical resistivity,2,4,18,19,21-26: thermoelectric power, magnetic properties,2,18,21,26 and X-ray strain'1,2,14,15,27,28 The slip systems and the appearance of the microstructure after deformation in ordered and disordered alloys have also been investigated. 6.10.ll,l7,23,29,30 The alloys used in these investigations were copper-gold. 1-4,6,9,10,14,16,18 23-26,28,29 copper-zinc. L7. 19. 20,30 nickel-manganese,2 11,21 iron-aluminum. iron-nickel,' copper palladium,4,18,22 gold-silver.4,18,22 and palladium-silver-copper.4 In the work reported here. the changes in properties of polycrystalline Cu3Au were investigated as functions of strain in wire drawing and rolling. Specimens were heat treated before deformation

Jan 1, 1961

By Douglas J. Harvey

The surface tension of lead-tellurium alloys (in the range 0 to 6.70 at. pct Te) ad lead-arsenic alloys (in the range 0 to 10.53 at. pct As) has been examined by the maximum bubble pressure method. The surface tension of high purity lead can be expressed over a range of temperatures (M.P. to 750°C) by the relationship = 512 — 0.133t. Arsenic and lellurium were both found to be surface active in lead. The surface excess concentration at 700°C was .found to be constant at 0.74 x 10-10 moles per sq cm for tellurium. The surface excess of arsenic at 700°C varied between 0.26 x 10-10 moles per sq cm at 0.16 at. pct to 1.5 x 10-10 moles per sq cm at 8.3 at. pet. The surface tension of lead has been investigated sporadically over the last 90 years. The first recorded study was made by Quincke1 who used a rather crude application of the drop-weight method. Hogness2 obtained the first reliable value for the surface tension of lead by measuring the pressure necessary to force the liquid metal from the end of a small capillary of known radius. Hogness gives the empirical formula: ?(surface tension)= 444 - 0.077 (l-327) which describes the relationship between surface tension and temperature ("C); This work was followed by that of Bircumshaw3 who applied the two-capillar maximum bubble pressure method of Sugden.4 The values published by Bircumshaw are in good agreement with those found by Hogness. Matuyama5 examined the surface tension of lead by the drop weight method. However, the values found are somewhat higher than those previously reported. More recent investigations by Greenaway6 and by Melford and oar' have resulted in data which are not in close agreement, but which are well within the expected range at most temperature levels. The study of the surface tension of lead alloys began with Drath and sauerwald8 who examined the lead-bismuth system. They found that the surface tension values of the alloys in this system deviated only slightly from the ideal mixture rule. Matuyarna5 examined the lead-antimony system. Bircumshaw the lead-tin, and the lead-indium system was studied by Hoar and Me1ford.10 The results of these investi- gations are summarized in Fig. 1. It may be seen that the data vary only slightly from the ideal mixture rule. As none of the metals studied and reported in the literature has a very interesting effect on the surface tension of lead,*the next logical step might be to examine the effect of metals of low solubility or nonmetals having some liquid solubility. Baes, and Kellogg11 have pointed out that the elements most likely to be surface active in liquid metals are ones of limited solubility in the liquid state and possessing weak intermolecular bonding forces in the solid state. Therefore, Group V, VI and VII nonmetals are most likely to be surface active. It was decided to examine the effect of a borderline metal, arsenic, and a nonmetal tellurium. Arsenic forms a eutectic with lead having a melting point of 290°C and a composition of 6.9 at. pct As.13 Pure arsenic sublimes when heated above 610°C. Even though alloying with lead lowers the vapor pressure, fumes containing arsenic are evolved at a very noticeable rate when lead-arsenic alloys are heated above 650°C. Tellurium raises the liquidus temperature of lead sharply to 904°C,14 the melting point of PbTe (lead-telluride). The solid solubility of tellurium in lead was determined by Greenwood and worner15 to be 0.0007 at. pct at room temperature. EXPERIMENTAL METHOD Surface Tension Measurement—The two capillary maximum bubble pressure method for measuring surface tension developed by Sugden,4 and used later

Jan 1, 1962

By N. H. Harrison, J. K. Rodgers, S. Regier

This paper presents comparisons of data obtained from a laboratory reservoir study and from a calculated behavior prediction with the actual production history of a condensate reservoir. A small non-commercial discovery was depleted under closely controlled conditions and the well fluids were sampled at frequent intervals. Data on the reservoir and production variables were accumulated on a fixed schedule. A laboratory reservoir study war made using the initial well fluid samples as charging stock. The production procedures and operating conditions were held constant throughout the study wherever possible and in general paralleled the field work. The well fluid compositions and the cumulative recoveries ar a function of the reservoir pressure were also calculated using conventional flash vaporization procedures and equilibrium constants. Comparisons based on the composition of the well fluid show good agreement, the laboratory study agreeing within experimental accuracy with the field work and the calculated data comparing equally well. The gas-oil ratios are also in good agreement, but with somewhat greater deviations at the higher pressures. In the overall picture, it is believed that a model st,~tiy can predict within experimental accuracy the production history of a condensate reservoir. Better equilibrium constants for the heavier hydrocarbons are needed in order to attain improved composition accuracy by calculation. INTRODUCTION In Aug., 1955, a gas condensate well was completed in San Juan County, Utah, that was initially thought capable of good commercial production. These conclusions were derived principally from core data and electric logs, which indicated good permeability, porosity and gas content. However, after the usual series of potential tests it was found that the reservoir pressure had declined some 22 per cent, and it was obvious that the zone tapped was but a small pocket or trap. It became apparent that, with a controlled depletion of a small reservoir, a unique opportunity was available to compare laboratory and calculated studies with an actual field depletion and to further the present knowledge of condensate reservoirs. FIELD WORK The Coalbed Canyon Well No. 1 was conventionally completed in the Paradox limestone formation to a total depth of 5,912 ft. The producing zone from 5,762 to 5,806 ft was perforated with four jet shots per ft. The wellhead and field equipment were also conventional, the major items consisting of a two-pass indirect fired line heater, a high- and a low-pressure separator with the necessary controls and accessories, gas meters, back-pressure regulators, flare stacks and condensate stock tanks. The initial testing of this well con-sisted of a series of flow potential and pressure build-up tests during which some 30 MMcf of gas was produced. The reservoir pressure declined from an estimated 2,300 to 1,782 psig during this period, from which it was concluded that the reservoir was very small. In order to approach steady-state conditions in the reservoir and so provide optimum conditions for making comparisons, the field depletion was programmed to approach, if possible, constant production conditions. Bi-hourly readings were taken of the tubing pressure, the pressure and temperature of the separators, oil and gas rates, and other pertinent operating data. The gas rate, as indicated by the orifice meter, was held constant by the adjustment of the choke in the line heater. The temperature of the first stage separator was held constant by adjustment of the line heater jacket temperature. Practical considerations of production made the maintenance oi a constant gas rate impossible. The test started with a gas rate of 4 MMcf/D and a separator pressure of 250 psig. This rate was maintained until the choke was fully opened. The gas rate then declined with the falling tubing pressure and production was continued until the rate was about 2

By L. Adler

Any design of a mine roof in bedded deposits which ignores differential deflections at the supports can quickly lead to dangerous overstressing. As illustrated by the typical case presented on page 1029, if one support deflects only 0.1 in. more than its neighbor, stresses 50 pct higher than anticipated will result. Previous studies1-3 have assumed that the mine roof acts like a beam clamped at both ends, with no deflection or slope' at the supports (Fig. 1). Essentially this assumes:5* ?A = ?, = ?b and dA = dB = 0. However, a host of conditions can cause differential sag or deflection at one support, and even an apparently trivial deviation can have major consequences in the beam: variations in spacing, size, and material of supports; local settlement due to floor heaving or lower mine openings; non-uniformity in wedging or bolt tensioning.8 In one particular instance concrete pillars, poured in place, shrank away from the roof upon setting.' When understood and controlled, differential sag becomes a definite advantage, permitting the transfer of loads from one support to another. Where one support is a bolt or prop and the other a pillar or rib, the load on the former can be reduced while the load on the pillar or rib is increased. This means two things: 1) savings can be effected by using smaller artificial supports such as props and bolts and 2) loads can be increased in the rib to induce eventual failure for longwall caving. In this last respect, another factor underlying coal bumps becomes understandable. The design of a mine roof by beam theory assumes that when a bedded formation is homogeneous, iso-tropic, and elastic, it will act like a uniformly loaded beam, fixed or clamped at both ends,6 s in Fig. 1.'" These are important assumptions. The uniform load is the weight of the beam itself or its body load. For a beam of unit width, this load becomes per foot: w = yhb. w = uniform load y = specific weight, lb per ft3 h = depth or thickness of beam, ft b = width of beam, ft Since b = 1, w - yh [ 1 ] A beam fixed at both supports is statically indeterminate to the second degree, that is, two equations in addition to those for static equilibrium are required to determine the reactions at the supports. These equations result from the geometry. From statistics: ^F,= 0;RA + RB — wL = 0 [2] wU S,M0 = 0 ; MB H------------MA — RBL = 0 [3] it Using superposition in which the left support is kept intact and the right one resolved into an active force and moment, we obtain from Fig. 2:" ?b = ?m + ? b2 + ?B3 = 0 dB = dB1 + dB2 + dB3 / wU Ma RBL \ V d°=(-r+-2—r)-zr [5] Where: Subscripts refer to position on axis of beam R = force, lb M = moment, ft-lb w = uniformly distributed load, lb per ft L = span, ft 0 = slope, radians d = deflection, ft E = Young's Modulus, lb per ft" 1 = Moment of inertia about the neutral axis, ft4 Solving Eqs. 2 through 5 for the four unknown reactions, we obtain: 2 L3 Fig. 3 is obtained by plotting MA and M, vs d, and then plotting RA and R, vs d, and adjusting their scales to coincide with MA and M, vs d,. It should be noted in Eqs. 6 and 7 that as dB increases, MA increases the same amount that Mn decreases. Eqs. 8 and 9 show a similar relationship between RA and R,. Proof that MA is the maximum absolute moment occurring in the beam is as follows: where M, = moment at x wx-When where M,., = maximum positive moment Eqs. 6 and 8 into above: For small dn: MA> M,,,

Jan 1, 1960

By V. M. Dobrynin

Experimental data demonstrate that physical properties of porous rocks change under pressure. In this paper an assumption is made and proved that under pressure the changes of physical properties such as porosity, density, permeability, resistivity and velocity of elastic waves are controlled to a large extent by the pore compressibility of rocks. It is also shown that the pore compressibility of rocks can be determined, within the range of pressures from 0 to 20,000 psi, by knowing the maximum pore compressibility and the magnitude of the pressure. Mathematical equations were developed which permit one to define changes in physical properties of porous rocks under pressure. These equations were verified by experimental data obtained from the study of sandstones. INTRODUCTION In studying the behavior of porous rocks under pressure in the field of petroleum technology, the most interesting aspect is the observation of those properties which characterize the rocks as possible reservoirs — for example, porosity, permeability, resistivity, density and the velocity of elastic waves. The literature dealing with this problem mainly contains data concerning the study of only one or at most two of these parameters, but not of the group as a whole. An attempt is made in this paper to find general equations involving each of these parameters, which will permit the study of the behavior of rocks under pressure. All experimental data used here were obtained from the investigation of consolidated sandstones. EXPERIMENTAL In addition to the use of published experimental data, an experiment was carried out which studied the main physical properties of sandstones under pressure. Two homogeneous quartz sandstones were chosen for this purpose: (1) the Torpedo sandstone from Kansas, and (2) the Medina sandstone from Ohio. The porosity of the Torpedo sandstone was 20.2 per cent, and that of the Medina 8.7 per cent. Permeabilities were 45 md and less than 1 md, respectively. Each sandstone contained about 5 per cent clay minerals, consisting mostly of kaolinite and chlorite, which were distributed quite evenly throughout the samples. One cylindrical sample 2 in. in diameter and 5 in. in length was cut from each sandstone and then saturated in a vacuum with a 3N solution of NaC1. This high concentration was used in order to obtain true formation factors and to decrease the swelling of the clay minerals. The methods of mounting the samples and measuring the changes in porosity and resistivity were practically the same as those described by Fatt and Mann.1,2 Changes of resistivity under pressure were studied for sandstones with 100 per cent water saturation, and for sandstones with the irreducible water saturation. The irreducible saturation was obtained by enclosing the saturated rock samples in relatively fine silicate powder so as to remove the water by capillary action. This procedure is described by Orkin and Kuchinski.3 Changes of permeability with pressure were determined at room temperature using nitrogen as the flowing medium. In studying the effects of pressure, one series of measurements was made using an internal pore pressure Pi equal to the atmospheric pressure, while the overburden pressure Pe ranged from 0 to 20,000 psi. A second series of measurements was used over the same range of overburden pressure, but with an internal pore pressure of 1,800 psi. When the results were compared on the basis of net overburden pressure (Pe — 0.85 Pi ), there was practically no difference for these two sandstones. The origin of the factor 0.85 in the expression for net overburden pressure is given by Brandt,4 Fatt5 and Geertsma.6 The relatively large bulk volumes of the sample (about 15 cu in.) made it possible to observe a

By L. A. Roe, E. A. Elevatorski

Wollastonite, named after William H. Wollaston, an English chemist, is a calcium metasilicate, CaSiO3; CaO: 48.30%, SiO2: 51.70%. It has a short history as an industrial mineral. The earliest production of wollastonite is reported to be from a deposit near Code Siding, located north of Randsburg, CA. At this locality small tonnages of wollastonite were quarried during 1933-34 and 1938-41, and processed into mineral wool. This operation was largely experimental and virtually no United States production was again reported until the 1950s when a large deposit near Willsboro, NY, was developed by the Cabot Corp. A processing plant was placed onstream in 1953, with nearly continuous production to date. It is currently operated by NYCO, a division of Processed Minerals, Inc. Since 1958, wollastonite deposits in the Little and Big Maria Mountains of Riverside County, and in the Panamint Range of Inyo County, both in California, have operated intermittently for production of both ornamental and commercial wollastonite. During 1980, the United States was the major producing country, furnishing about 75% of the world's output. Current production comes from Finland, Mexico, India, and Kenya. Small amounts have been shipped intermittently from the USSR, New Zealand, Republic of the Sudan, Republic of South Africa, and Namibia (South-West Africa). The principal use of wollastonite is in the manufacture of plastics. Other uses are for paints, ceramics, adhesives, fluxes, glazes, thermal insulation board, and refractory products. Mineralogy Pure wollastonite, CaSiO3, has the composition of 48.3% CaO and 51.7% SiO2. However, it is seldom found in the pure state due to the ease with which it takes into solution the metasilicates of manganese, magnesium, iron, and strontium. Predominantly, wollastonite occurs as a contact metamorphic deposit forming between limestones and igneous rocks. Commonly associated minerals are garnet, diopside, epidote, calcite, and quartz. It has a specific gravity of 2.8 to 3.0, and hardness of 4.5 to 5 on Mohs' scale. When pure, it has a brilliant white color, but with impurities it may be grayish or brownish. Luster is vitreous to pearly. Melting point of wollastonite is about 1540ºC. Wollastonite occurs in coarse-bladed masses, rarely showing good crystal form. It is usually acicular or fibrous, even in the smallest of particles. The most unique property of crushed and ground wollastonite is its cleavage. Fragments of crushed wollastonite tend to be needle-shaped, imparting a high strength, and this property is the basis for many of its uses. The fiber lengths are commonly in the ratio of 7 or 8 to 1, length to diameter. The average diameter of wollastonite is 3.5 µm. Some crystals of wollastonite fluoresce under short-wave or longwave ultraviolet light, or both, with colors ranging from yellow-orange to pink-orange. Specimens may also show phosphorescence. The acicularity of wollastonite is a property of considerable importance to the marketplace. The plastics industry makes utilization of high aspect ratio grades of wollastonite (20:1) for reinforcing thermoplastics and thermoset polymer compounds. The naturally high pH of 9.9 (10% water slurry) is a prime property to the coatings industry. The coatings industry uses milled grades of wollastonite as a pH stabilizer in interior and exterior PVA and acrylic latex systems. Processed wollastonite can have a G.E. brightness of 90 to 93. Chemically, wollastonite is inert and this property makes it useful as a filler and reinforcing agent. There are two polymorphs of calcium silicates: wollastonite, a low temperature form, and pseudowollastonite, a high temperature

Jan 1, 1983

By Herman Dykstra, Irving Fatt

Relative wetting phase permeabilities calculated from capillary pressure-saturation data are compared with measured relative permeability data. The equation relating relative permeability to capillary pressure-saturation data is derived by assuming that a porous medium is analogous to a bundle of capillary tubes. The equation includes a term to correct for the difference between fluid path length and length of core. Relative oil permeability data in the water-oil-gas system are presented, and it is shown that if the water is considered part of he rock matrix, the relative oil permeability curve is typical of a wetting phase relative permeability curve. Apparatus for measuring two and three-phase relative permeabilities are described. INTRODUCTION The difficulty in measuring relative permeability of cores has made it desirable to have a correlation between relative permeability and some more easily measured property of porous media. Such an easily measured property is the capillary preisure-saturation relation. In the past there has been a tendency to separate the capillary pressure concept from the very complex pore geometry and to consider that the capillary pressure-saturation curve give? only some characteristic distribution of interfacial curvature between two fluids.' Recently, however, it has been realized that if capillary pressure data are to yield information concerning pore size distribution2 and fluid flow in porous media. the pore must be assumed to have a simple shape such as a cylinder or a sphere. Childs and George' showed that relative water permeability could be calculated from the capillary pressure-water saturation curve by assuming that the pores are cylinders in which the fluid flowing obeys Poiseuille's Law, and that the capillary preisure curve indicated the size and number of pores. Purcell5 derived an equation relating the permeability of a porous medium to the capillary pressure curve by assuming the porous medium to he analogous to a bundle of capillary tubes. Gates and Lietz6 calculated relative permeability by an equa-tion derived from Purcell's equation. In this paper an equation is derived giving the relative water permeability as a function of water saturation by as-sulning that the core sample can be represented by a bundle of capillary tubes in which the fluid path length is not the same as the bulk length, and in which the fluid path length varies with saturation. THEORY For a bundle of N capillary tubes the flow through dN tubes will be dQ = qav dN ..........(1) where Q is the total flow rate through all the tubes and q., is the aberage flow rate through the tubes in the interval dN. If the interval dN is made small the average flow rate through a tube in the interval dN is given by Poiseuille's Law, 8 nl where ?P is the pressure drop across the tube of length 1 and radius r, and µ is the viscosity of the fluid. Substituting for qav in Equation (1), dQ =pr4 ? P/8µl dN.........(3) Darcy's Law for linear flow of an incompressible fluid in a porous medium is Q = KA?P/µL...........(4) where K is the permeability and ?P is the pressure drop across a bundle of tubes which has replaced the porous medium of length L and cross-sectional area A. Differentiating Darcy's Law with A, ?P, µ, and L constant, dQ = A?P/µL dK.........(5) Equating (3) and (5) gives dK------------dN.........(6) 8Al If the pores are assumed to be cylinders dV =pr2ldN..........(7) where V is the volume of flowing fluid in the pores. Substituting for dN in Equation (6) gives dK = —dV..........(8) 8AI' By definition the saturation, S, of the core is v v S =V/Vp..........(9) Vv -PAL

Jan 1, 1951

By Irving Fatt, Herman Dykstra

Relative wetting phase permeabilities calculated from capillary pressure-saturation data are compared with measured relative permeability data. The equation relating relative permeability to capillary pressure-saturation data is derived by assuming that a porous medium is analogous to a bundle of capillary tubes. The equation includes a term to correct for the difference between fluid path length and length of core. Relative oil permeability data in the water-oil-gas system are presented, and it is shown that if the water is considered part of he rock matrix, the relative oil permeability curve is typical of a wetting phase relative permeability curve. Apparatus for measuring two and three-phase relative permeabilities are described. INTRODUCTION The difficulty in measuring relative permeability of cores has made it desirable to have a correlation between relative permeability and some more easily measured property of porous media. Such an easily measured property is the capillary preisure-saturation relation. In the past there has been a tendency to separate the capillary pressure concept from the very complex pore geometry and to consider that the capillary pressure-saturation curve give? only some characteristic distribution of interfacial curvature between two fluids.' Recently, however, it has been realized that if capillary pressure data are to yield information concerning pore size distribution2 and fluid flow in porous media. the pore must be assumed to have a simple shape such as a cylinder or a sphere. Childs and George' showed that relative water permeability could be calculated from the capillary pressure-water saturation curve by assuming that the pores are cylinders in which the fluid flowing obeys Poiseuille's Law, and that the capillary preisure curve indicated the size and number of pores. Purcell5 derived an equation relating the permeability of a porous medium to the capillary pressure curve by assuming the porous medium to he analogous to a bundle of capillary tubes. Gates and Lietz6 calculated relative permeability by an equa-tion derived from Purcell's equation. In this paper an equation is derived giving the relative water permeability as a function of water saturation by as-sulning that the core sample can be represented by a bundle of capillary tubes in which the fluid path length is not the same as the bulk length, and in which the fluid path length varies with saturation. THEORY For a bundle of N capillary tubes the flow through dN tubes will be dQ = qav dN ..........(1) where Q is the total flow rate through all the tubes and q., is the aberage flow rate through the tubes in the interval dN. If the interval dN is made small the average flow rate through a tube in the interval dN is given by Poiseuille's Law, 8 nl where ?P is the pressure drop across the tube of length 1 and radius r, and µ is the viscosity of the fluid. Substituting for qav in Equation (1), dQ =pr4 ? P/8µl dN.........(3) Darcy's Law for linear flow of an incompressible fluid in a porous medium is Q = KA?P/µL...........(4) where K is the permeability and ?P is the pressure drop across a bundle of tubes which has replaced the porous medium of length L and cross-sectional area A. Differentiating Darcy's Law with A, ?P, µ, and L constant, dQ = A?P/µL dK.........(5) Equating (3) and (5) gives dK------------dN.........(6) 8Al If the pores are assumed to be cylinders dV =pr2ldN..........(7) where V is the volume of flowing fluid in the pores. Substituting for dN in Equation (6) gives dK = —dV..........(8) 8AI' By definition the saturation, S, of the core is v v S =V/Vp..........(9) Vv -PAL

Jan 1, 1951

By J. Lohrenz, R. J. Francis, G. C. Clark

Compositional material balance methods have been extended to solution gas-gas cap-water drive reservoirs. Partial or complete pressure maintenance as well as reservoir repressuring by gas and/or water injection may be considered. Increrments of reservoir withdrawals and injections equivalent to a given time increment are used as the basic iterative step to allow the extension. An example field study shows how the reservoir volumes. extent of the water drive, and relative permeability ratio function may be estimated. Knowing these, the methods can be used to predict reservoir behavior. INTRODUCTION The compositional material balance is similar to conventional material balances in many aspects. The assumption of a pressure-homogenous reservoir with no transients is common to all conventional material balances1-1 although numerical methods are different. The compositional material balance also contains this assumption. Inclusion of phase equilibria directly in reservoir calculations is the key difference and improvement in the compositional material balance. A continuous account of the compositions of the reservoir fluids and effluents is determined. Calculations of fluid phases and properties can be based on these compositions. Excepting fluid phases and properties determined in this manner, the compositional material balance uses procedures exactly analogous to conventional methods. Brinkman and Weinaug5 and Jacoby and Berry6 have given compositional material balances for simple solution gas drive reservoirs; subsequently Jacoby and Berry' extended their procedure to include pressure maintenance by gas injection at constant reservoir pressure. This paper extends compositional material balance techniques to a more general combination drive reservoir with various types of pressure maintenance operations. The procedure for this extended compositional material balance and an example of its field application are presented. PROCEDURE The numerical calculations of this compositional mate- rial balance are performed on a high-speed digital computer. The most comprehensive case the present compositional material balance program can handle is a solution gas-gas cap-water drive reservoir with (1) water production, (2) water influx, (3) separate or commingled production from the oil and gas zones, (4) any surface processing arrangement including LPG and gasoline production, and (5) gas injection into the gas and/or oil zones." This description of the procedure applies to this most comprehensive case. Appropriate deletions in the procedure can be made for more simple mechanisms. ASSUMPTIONS In common with conventional material balances, the compositional material balance makes the following assumptions: (1) the total reservoir is at a uniform pressure at a given time; (2) the volumetric ratio of gas to liquid withdrawal at reservoir conditions is directly proportional to their mobilities; and (3) instead of assuming a pressure-dependent total gas-oil equilibrium, the cornpbsitional material balance makes a sianificantly different and more rigorous assumption as follows. Thermodynamic equilibrium exists for each component in the liquid-gas mixtures in the reservoir and in surface processing. These three assumptions are fundamental in the compositional material balance procedure. Other assumptions involve arbitrary simplifications of the reservoir mechanisms. This compositional material balance uses the following simplifying assumptions: (1) no bulk mixing occurs between the gas and oil zone; (2) total mixing occurs between any injected gas and the contents of the injection zone; (3) no in situ hydrocarbon liquid is produced from the gas zone; (4) no reservoir gas is retained behind any net bottom-water influx; (5) the residual saturation of hydrocarbon liquid at reservoir conditions behind any net bottom-water influx is constant; and (6) fluids retained behind any net bottom-water influx exert no further influence on the reservoir. These six assumptions are not fundamental. Modifications of the present procedure can be made for different assumptions of the reservoir mechanisms. THE TIME INCREMENT Previous compositional material balances have used a pressure increment for simple solution gas drive reservoirs"" and an increment of injected gas for pressure maintenance.' These increments are suitable for these simple cases. They are, however, inadequate for reservoirs with time dependent drive mechanisms including water drive. A time increment removes this inadequacy and is

By R. D. Begley, A. W. Khair

Introduction The technique being investigated here uses a physical mine model which consists of a plexiglass room-and-pillar model 11.4 x 11.4 cm (4.5 x 4.5 in.). It is placed at the bottom of a larger plexiglass box 25.4 x 30.5cm (10 x 12 in.). The box is then filled with a model material representing geologic overburden. Sand has been selected because of its deformability under gravitational forces. The mine model is capable of creating a cavity and induces movement on the surface of the model material. This surface disturbance can be measured using a photographic technique called holographic interferometry. This paper discusses results of an investigation of the possibility of generating subsidence criteria using physical mine models and laser holographic interferometry. Problem Statement Subsidence appears to be widespread in two states - Pennsylvania and Illinois - although serious subsidence is found in localized areas of numerous other states (HRB Singer Inc., 1977). The US Bureau of Mines estimates that more than 1618 km2 (400,000 acres) of urban areas are threatened by subsidence. Pennsylvania and West Virginia contain 57% of these threatened areas, with the remainder being spread out over 16 states (Johnson and Miller, 1979). Events of subsidence have been reported in urban areas, including Pennsylvania (Gray et al., 1977) and Illinois (Mahar and Marino, 1981). One location has observed three episodes of subsidence beginning in 1950 (Bruhn et al., 1981). These subsidence events all occurred over abandoned room-and-pillar mines. Some mines have been abandoned more than 50 years (Bauer and Hunt,1981) . Available research does not offer much on the problem related to time dependent subsidence, especially those associated with the failure of remnant coal pillars. The physical modeling technique developed at West Virginia University is capable of measuring both instantaneous and time delayed subsidence. Laser Holographic Interferometry Laser holographic interferometry, or holometry, uses the interference properties of laser light to measure small movements of objects. The object requires only optical access and no surface instrumentation. This technique is closely related to photoelasticity. Both techniques rely on the uniqueness of the light reflected off or through an object in any given state of stress. The entire measurement process is done remotely without any surface instrumentation and without using a camera. The interferogram is constructed and visual observation will yield fringes on all areas that yielded or moved. The displacement represented by one fringe was determined by Wilson et al. (1977) and is given by the relationship: [d ='1(1) 2 cos 8 where d = displacement of one fringe 6 = incident angle of laser light A = wavelength of light (632.8 x 10-9 meters for Helium-Neon lasers)] In summary, holometry is capable of detecting very small movements of objects under study. It is a very sensitive process that requires laser light and vibration free experiments. Experimental Program Lab Facilities The primary components of the facilities were the laser and the optical table. The light source was 50 milli-watt, Helium-Neon laser. The optical table was a granite slab weighing nearly one ton. All experiments were performed on the optical table. It consists of a large table sitting on air cells similar to car innertubes. Six air cells were placed on top of two layers of wood and Styrofoam. Figure 1 is a schematic diagram of the vibration free table. Even with these precautions, the experiments still had to be performed at night to reduce potential traffic related vibrations. Room-and-pillar mine models consisting of three different pillar sizes were constructed out of plexiglass. The models were designed to represent three different extraction ratios. Figure 2 shows the models. The extraction ratios were determined by the tributary loading concept (Peng, 1978). These models were the primary component of a seam extraction simulator. A schematic of this simulation device is shown in Fig. 3.

Jan 1, 1986

By R. W. Fountain, J. F. Libsch

ONSIDERABLE time and effort have been ex- pended recently in research designed to provide a better understanding of the solid state transformations leading to the permanent magnet qualities of many commercial alloys. The knowledge derived from these investigations has resulted in the development of new magnetic materials as well as improvements in some of the older materials. The permanent magnet qualities of the newer alloys having high residual inductions and high coercive force values, have been attributed to the precipitation of a second phase or an order-disorder transformation in a solid solution. The ternary Co-Fe-V alloys provide an interesting group of magnetic materials having a wide range of magnetic properties. The low vanadium alloys' possess soft magnetic properties similar to the binary Co-Fe alloys, while at higher vanadium contents (about 10 pct) interesting permanent magnet properties are obtained.2,3 Previous work by Nesbitt' on the ternary Co-Fe-V permanent magnet materials, resulted in the conclusion that the precipitation of a second phase, ?, was responsible for the permanent magnet qualities of these alloys. However, unlike the conventional precipitation hardening alloys, the ? phase, which is stable at high temperatures, is precipitated in the low temperature a phase upon aging. Recent work by Geisler and Martin"." has established that the order-disorder reaction inherent in the binary Co-Fe alloys is also present in the ternary Co-Fe-V alloys. On the basis of their investigation, it was proposed that ordering of the a phase is also a possible explanation of the mechanical and magnetic hardening in these alloys. Thus, it appears that the Co-Fe-V alloys are rather unique in that two solid state transformations are taking place which may account for the properties of these alloys. In view of these facts, the present investigation was undertaken to study the mechanism of mechanical and magnetic hardening in the 10 pct V-Co-Fe alloy with particular reference to the role of the ? phase; and to correlate the results of this study with other magnetic properties. In addition, the influence of increasing vanadium contents on the transformation of Co-Fe alloys was also undertaken to supplement the main objective of the investigation. Experimental Procedure The Co-Fe-V alloys used for study consisted of: 1—a commercial lot of Vicalloy I, a 10 pct V alloy, and 2—alloys prepared by powder metallurgy and containing 52 pct Co and 0 to 14 pct V, balance iron. The commercial alloy was received in the cast and forged condition and was used primarily for the detailed study of hardness, electrical resistivity, and magnetic properties of the 10 pct V alloy. The powder metallurgy alloys were prepared from electrolytic iron, commercial cobalt, and ferro-vanadium powders. The processing cycle was as follows: 1—Elemental powders were weighed and mechanically mixed for 24 hr. 2—Mixed powders were pressed at 60,000 psi into the form of bars 5/8 x 5/8 x 5 in. 3—Bars were sintered in an atmosphere of purified dry hydrogen for 24 hr at temperatures ranging from 1380" to 1340°C, depending upon the vanadium content. The as-sintered bars were cut in half; and one half was machined to 7/16 in. diam for thermal analysis. The remainder of the bar was hot forged to approximately 0.050 in. for X-ray diffraction analysis. Additional details of sample preparation are presented subsequently according to the type of test performed. Thermal Analysis Measurements One of the simplest and most widely used tools for determining phase relations as they change with temperature and composition is to study the rate of change of temperature of a material as heat is supplied or extracted at a constant rate. The thermal analysis method used in this investigation is the inverse-rate method described by Smith.'

Jan 1, 1954

By F. Brons, V. E. Marting

In the past, other authors.1, 2 have studied the influence of a skin effect on the productivity of a well. This skin effect, expressed by the skin factor S, is considered to be caused by a thin layer of impaired permeability immediately around the wellbore and extending vertically over the whole productive interval penetrated by the well. The skin factor S is defined as follows. Based on this, the impairment in productivity caused by a skin can be expressed by the fractional loss in productivity I, which is the loss in productivity divided by the total unimpaired productivity. For compressible flow in a stabilized well which drains a circular area of radius re, 1 is given by The present paper deals with a different kind of productivity impairment. Consider a well in which part of the productive formation is blocked off completely, either by incomplete penetration or by exclusion of parts of the productive zone by blank casing. In Fig. 1 (A, B and C), three examples are shown. Fig. 1(A) shows the situation where a well only partially penetrates the formation. This often is done to combat the actual or imagined danger of bottom-water coning. Fig. 1(B) shows a well producing from only the central portion of a productive interval. This type of completion is sometimes used where both water and gas coning are a problem. Although the case of a well producing through perforated casing cannot be treated in a manner similar to the previous two cases (where radial flow in the horizontal plane is assumed), Fig. 1(C) shows several intervals open to production and qualitatively describes this case (as will be discussed later). To study the loss in productivity in all these cases, two parameters are introduced which fully determine the impairment. The first is the penetration ratio "b," i.e., the total interval open to fluid entry divided by the total thickness of the productive zone. The second is the ratio hlrw. In this ratio, rw is the wellbore radius. The definition of h is more cumbersome. In Fig. 1(A), h is the thickness of the total productive interval. The streamline configuration for this case of partial penetration is basic to the other two cases considered (Views B and C of Fig. 1). It will be obvious from Fig. i (A) that the flow lines in the uppermost portion of the formation will be essentially horizontal, while those in the lower portion will curve upward toward the well. In Fig. 1(B), with only the middle portion of the zone open to production, the streamline configuration. of the upper half will be an exact mirror image of that in the lower half of the zone. Hence, for the case illustrated in Fig. 1(B), h is defined as one-half the total sand thickness. It follows, then, that in Fig. 1(C) h is one-half the distance between corresponding points in adjacent intervals. (In gun-perforated casing, h would be one-half the distance between perforations.) In a paper by Nisle4 and a paper given by the present authors5, the mathematical theory was developed for the cases under consideration. In the present publication, the emphasis is put on the results of these studies; consequently, the equations derived in Refs. 4 and 5 will be omitted, for the most part. The pressure drop p in a well producing from only a portion of the total formation thickness can, in analogy with Eq. 1, be expressed by Numerical solution of Eq. 5 by use of the IBM 650 leads to the following important conclusions. First, during a short period after starting production (usually on the order of a few minutes), the function p(b) is given by