Search Documents

By T. D. Mueller, H. Dykstra

This paper describes a method for obtaining compositions of gas and oil phases in equilibrium with each other at a given reservoir temperature and pressure as lean or enriched gas is injected into a reservoir oil not in equilibrium with the gas. The method is based on calculation of equilibrium constants, K, from the equation log Kp = A + B F, where p is pressure and F is component characterization factor. A and B are constants that vary with a composition parameter. The paper explains how A and B could be determined with a minimum amount of data. The method consists of calculating equilibrium constants as a function of liquid composition and using the K's to calculate changes in compositiort of injection gas as it successively contacts original resetvoir oil at the gas-oil front and the change in composition of reservoir oil as it is successively contacted by the injection gas in the vicinity of an injection well. An iterative process is required. Two examples of the method are given which show how oil and gas compositions and properties are calculated with the aid of a computer. INTRODUCTION The calculation of recovery from lean-gas or enriched-gas injection into a saturated or under-saturated reservoir oil requires a knowledge of how phase properties change in a reservoir between injection well and producing well. To determine the phase compositions or properties experimentally would be difficult and time consuming, especially if one wanted to look at the effect of various injection gas compositions. For some fields an engineer would want to make preliminary calculations on the economics of gas injection before he makes the decision to spend a corrsiderable sum of money for laboratory work to obtain the necessary fluid properties. A method for calculatine compositions and properties based on a minimum of laboratory data is described. The minimum amount of data required is injection gas and reservoir oil composition and injection gas solubility as a function of pressure for an undersaturated oil, or the compositions of mixtures of reservoir fluid with two or more enriched gases, which mixtures have a bubble-point pressure that is the same as the expected average displacement pressure. On the other hand, if the economics appear to be attractive, more laboratory data can be obtained in the form of equilibrium constants for three or more mixtures of the reservoir fluid and injection gases. The method is based on the calculation of K's from the equation, log Kp = A + B F, where K is equilibrium constant, p is pressure and F is component characterization factor. This method of plotting K data by means of the above equation is described by Hoffman, Crump and Hocott.1 A and B are constants that vary with a composition parameter. It is shown how experimental data are used to determine how A and B vary with composition. BASIS OF METHOD EQUILIBRIUM CONSTANTS Equilibrium constant (K) data are required for this method. The K of each component is obtained from an equation that gives the product of K and pressure as a function of component characterization factor. The characterization factor F = b(l/Tb — l/T), where b is a constant characteristic of a given component, Tb is the boiling point in OR at atmospheric pressure of the given component, and T is the temperature of the system in OR. This method of plotting K data is described by Hoffman, Crump and Hocott.l The straight-line equation of the log Kp - F plot is log hp = log Kop +BFi, where Kop is at F =0, i designates the component, and B is the slope of the line. The factor F can be

Jan 1, 1966

By S. Pearson, L. Rotherham

Measurements have been made of the variation of internal friction with temperature for spectroscopically pure silver, and for o series of solid solutions of silver with cadmium, indium, and tin, using a Ke-type torsion pendulum apparatus. Some experiments have also been made to investigate the effect of nonmetallic impurity on grain boundary relaxation in silver. The effect of the alloying elements is to increase the grain boundary viscosity, and to raise the activation energy for grain boundary relaxation from 22,000 cal per mol for pure silver to 43,000 cal per mol for the solid solutions; the same value being obtained, within the limits of experimental error, for all the alloying elements and solute concentrations investigated. Results of the experiments show exactly the same trend as those obtained previously for a similar series of copper solid solutions. They are in agreement with the general theory of grain boundary relaxation developed by Zener and Ke, but do not seem to be in agreement with either of the mechanisms so far put forward to explain grain boundary slip. EXPERIMENTS described in this paper were made as part of an investigation into the effect of alloying elements in substitutional solid solution on grain boundary viscosity, measured by the internal friction method developed by Ke.&apos; The apparatus and experimental procedure have already been described.&apos; Preliminary work was done with spectro-scopically pure silver. Experimental Work Materials Used—All the materials were supplied in the form of % mm diam wire by Messrs. Johnson Matthey and Co. Ltd. The silver and all the alloying elements used were spectroscopically pure. Composition of the alloys is given in Table I. They were all melted and chill cast under the same conditions as the copper alloys. The ingots were rolled to 1/4 in. sq rod with intermediate annealing at 600°C in an atmosphere of cracked ammonia when necessary. These rods were then drawn to wire by standard procedure with intermediate bright anneals at suitable steps in the drawing operations. Specimens from the annealed wires were sectioned longitudinally. mechanically polished, and examined under the microscope. They were then etched and reexamined. Grain sizes were reasonably uniform for all wires, and no traces of a second constituent were observed. The wires were weighed before and after testing so as to measure any loss of the alloying element due to evaporation. The loss was less than the accuracy of measurement—< 0.1 pct of the total weight— for the Ag-In and Ag-Sn alloys, but appreciable loss of weight occurred in the case of the Ag-Cd alloys. The percentage loss is given in Table 11. X-ray diffraction photographs were taken of all thc wires after test to determine whether there was any marked preferred orientation. Results of the examination are given in Table 111. Experiments on Spectroscopically Pure Silver-— Experiments were made on specimens in which the grain diameter was small relative to that of the wire, and on one specimen in which the grains were grown by the usual strain-anneal method to extend across the full diameter of the wire. These large grained specimens will be referred to as single crystal specimens. The wires were annealed for 11/2 hr at 700°C, and experiments were made at two frequencies of vibration, 1.5 and 0.4 vibrations per sec. The variation of internal friction with temperature is shown in Fig. 1. Experiments on a second small grain specimen gave approximately the same results as for the first wire. It will be seen from Fig. 1 that the peak in the curve for the small grained wire, which is presum-

Jan 1, 1957

By E. A. Gulbransen, A. Pebler

Pressure-composition isotherms have been determinedfor the systems ZrMo2 -H2 between 0" and 900°C at hydrogen pressures between 10-4 and 760 Torr. Tkese studies plus X-ray diffraction analyses of selected compositions show that these intermetallic compounds form a complete series of hydrogen solid solutions within this temperature and pressure range. The partial molar enthalpies and entropies of solution of hydrogen in the three intermetallic compounds were calculated. The solubility of hydrogen in these Laves phases can be related to their free electron concentration. ThE reaction of hydrogen with inter metallic compounds of zirconium with Laves phase-type structure has both theoretical and practical interest. The in-termetallic compounds have structures different from any of the pure metals. In addition the electron to atom ratio of the intermetallics varies appreciably. In the nuclear materials field, these materials are of interest since intermetallic compounds may exist as precipitates in commercial zirconium alloys with transition metals. From studies on the hydrogen reaction, one can establish the stability of the intermetallic compounds relative to the hydrides of the component metals. The intermetallic compounds may also be used as getters for hydrogen and as materials for hydrogen storage. This paper will present the results of equilibrium and structural investigations of the ZrCr2-H2, the ZrV2-H2, and the ZrMo2-H2 systems for the temperature range of 0" to 900°C and the pressure range of Torr to 1 atm. A complete thermochemical analysis will be made of the measurements. Some of the results have been referred to in a survey paper&apos; covering the types of hydrogen behavior observed in zirconium alloys. LITERATURE ZrCr2 is the only intermetallic phase in the Zr-Cr system.&apos; It exists in two allotropic modifications. Between 900" and 1000°C the hexagonal MgZa (C14) structure transforms into the cubic MgCu2 (C15) structure. Beck4 has previously shown that the intermetallic ZrCrz takes up hydrogen into solid solution up to a composition H/ZrCr2 = 1.14 at room temperature and 1 atm hydrogen. A ZrCr2 hydride was stated to have a fcc structure. The binary Zr-V system has been investigated by williams5 although Rostoker and Yamamoto 6 had given earlier a partial phase diagram. ZrV2 was the only intermediate phase found in the system. wallbaum7 reported earlier that the compound ZrV2 hp the MgZm (CJ4) type structure with a = 5.288A and c = 8.664A. This was not comfirmed by Elliot8 nor by Matthias, Compton, and corenzwit9 who reported that ZrV2 has a cubic MgCua LC15) type structure. The latter authors gave a = 7.429A. Also, several authors do not agree on the temperature of the peri-tectic horizontal. It appears that ZrV2 may exist in several modifications. A high capacity for occluding hydrogen was reported by Beck4 who found H/ZrV2 = 1.38 at room temperature and 1 atm of hydrogen. ZrMoz is the only intermetallic phase in the Zr-Mo system.10 ZrMoa is formed by a peritectic reaction of the melt with a molybdenum-rich solid solution. ZrMo2 has the C15-type structure with a = 7.59A. MATERIALS High-purity chromium, vanadium, and molybdenum were used together with grade 1 crystal bar zirconium to prepare the alloys. Table I lists the chemical analysis and the sources for the metals. Ten-gram alloy compacts were levitated and melted in an inert argon

Jan 1, 1968

By F. D. Patchen

Changes in the properties of partland cement upon the addition of fine-ground silica are discussed. Data were collected from formulations cured for periods up to 60 days at temperatures varying from 180" to 350°F. Compressive strengths, permeabilities and thickening times are compared. It is demonstrated that addition of find-ground silica in proper proportion results in a hydrated cement with superior properties above 230°F; i.e., early compressive strengths are higher and increase with additional curing, permeabilities are significantly lnrver and thickening times are increased. X-ray analysis of hydrated cement samples have been used to determine more fully the complex chemistry of these systems. Dicalcium silicate alpha-hydrate was found in all neat cements cured above 260°F. A 10 weight per cent silica addition produced a maximum amount of this phase. Xonotlite and tobermorite have been identified as hydration products at 320°F in systems containing 35 or more weight per cent silica. Formation of these components appears to be primarily responsible for the gain in strength and reduction of permeability observed after the addition of silica. INTRODUCTION For several years it has been recognized that loss of compressive strength occurs when portland cement is cured at a temperature above 250°F. This phenomenon has been called "strength retrogression". More recently it has been recognized that increased permeability usually is associated with decreased compressive strength at these temperatures. At curing temperatures above 300°F, these factors reach serious proportions. A compressive strength loss of 50 per cent between curing times of 24 hours and 14 days has been reported,&apos; and permeabilities as high as 10 md have been measured in neat cements cured at 320°F for seven days.&apos; Earlier investigators,&apos; although unsuccessful in eliminating strength retrogression, stated ". . . .it appears possible that by the proper selection of additives for cement, strength loss with age at high temperatures can be either eliminated or considerably lessened". One solution suggested was the use of a pozzolanic material, hydrated lime, and a retarder. Menze15 in 1934 discovered that portland cement containing finely divided silica possessed high compressive strengths when the curing temperature was 350°F. The use of such a composition is now an established practice in the manufacture of steam-cured cement blocks. It is reasonable to assume that the same beneficial effect can be used advantageously in high-temperature oilwell cementing operations. Ludwig first mentioned the use of sand for this purpose. Later, a portland cement containing moderate additions of silica flour and retarded with CMHEC, a cellulose derivative. was recommended for high-temperature application." These formulations have been successfully employed in the field."&apos;" This paper presents the results of a research program to test the efficacy of finely ground silica in increasing the compressive strength and reducing the permeability of oilwell cement compositions containing portland cement. Included also are X-ray diffraction patterns of powdered samples which have been used to determine some of the important chemical reactions occurring in these systems. LABORATORY PROCEDURES AND MATERIALS All compressive-strength test specimens were prepared and cured from a slurry in accordance with API RP 10B, with two modifications. A final curing pressure of 5,000 psi was maintained for a one- or three-day period; but, for further curing, pressures were reduced to the vapor pressure of water at the curing temperature. Also, the water content was varied in samples containing silica to maintain the slurry viscosity within a range of 5 to 15 poise as measured on a consistometer. Test specimens at 180°T were cured in a water bath at atmospheric pressure. At higher temperatures the specimens were cured in a steam-jacketed high-pressure curing vessel at 5,000 psi for a one- or three-day period. cubes cured for longer periods were removed from the pressure-curing vessel after 24 hours, sealed in small bombs containing water and placed in an aging oven.

By C. J. Bechtoldt, H. C. Vacher

Quenched alloys, prepared by powder metallurgical techniques, were examined by microscopic and X-ray diffraction methods. The compositions and heat treatments were chosen so that the chromium and nickel solvuses could be determined and a range of composition and temperature, in which a eutectoid had been reported, could be explored in small intervals. The results showed no evidence of eutectoid reactions having taken place. ALLOYS containing iron, chromium, molybdenum, and nickel are extensively used in high-temperature applications. In order to provide basic information on the solid-state reactions that occur in such alloys, the National Bureau of Standards is engaged in an extensive study of the equilibrium phases in the binary, ternary, and quaternary systems involving these four metals. Attention was directed to the chromium-nickel system because a review of the literature showed conflicting data. A recent compendium1 of the constitution of binary alloys gives two diagrams for the chromium-nickel system, one a simple eutectic and the second a diagram showing a eutectoid and a eutectic. The eutectoid was the result of an allotropic transformation in a high-temperature solid-solution form of a chromium-rich phase called 0. The approximate temperature and composition of the eutectoid was 1200°C and 30 at. pct Ni. In the preliminary study, at the Bureau, no transformation was indicated in the chromium-rich phase of ternary and quaternary alloys in which chromium and nickel were major components. Similar results have been reported from time to time in the literature,2-8 In consequence of this situation, a systematic study was made of the constitution of heat-treated chromium-nickel alloys using metallographic and X-ray diffraction methods. The compositions of the alloys were selected so as to determine the boundaries of the chromium-rich and the nickel-rich phases and to explore in small intervals the ranges in composition and temperature in which the eutectoid had been reported. The results of this study are reported and discussed in the present paper. EXPERIMENTAL PROCEDURES Materials, Alloys, and Heat Treatments—The alloys were prepared by a powder metallurgy technique similar to that used in earlier work.9 Analyses of the materials used are given in Table I. Powders were mixed using a vibrator in 30-g batches of predetermined composition; 2-g pellets were made from the mixture by compacting at room temperature in a 1/2-in. diam mold under 60,000 lb per sq in. pressure. The pellets were sintered in dry hydrogen for 10 days at 1300°C. At the end of the sintering treatment the specimens were quenched in water, which caused a pale green film to form on the chromium-rich alloys and thin black scale to form on the nickel-rich alloys. The pellets after sintering were approximately 1/2 in. diam by 3/32 in. thick. The purification train9 was modified by inserting a tube of titanium turnings, heated at 700°C to increase the efficiency of the removal of oxgen and nitrogen. It was found in the earlier work that the sintering treatment lowered the carbon, oxygen, and nitrogen contents of chromium-rich alloys to 0.01, 0.01, and 0.004 wt pct, respectively, and that there was a loss of chromium. In order to ascertain the extent of this loss, the chromium contents of certain alloys were determined chemically. They included those alloys used to determine the parameter-composition curves and those of critical composition, used to determine the chromium-solvus microscopically; these analyses were made on samples after the final heat treatments. The mixed powder values for the chromium content of those alloys that were not analyzed were adjusted to conform to a smooth curve based on the chemically analyzed values for

Jan 1, 1962

By C. Laird, C. E. Feltner

An investigation into the effect of a material&apos;s slip character on its high strain fatigue properties has been carried out using copper and Cu-7.5 wt pct A1 as representative wavy and planar slip mode materials, respectively. In copper, the constant term and exponent in the Coffin-Manson equation, ?EpNf1/2 = const, are unaffected by the temperature of testing or the history of the metal. On the other hand, in Cu-7.5 wl pcl Al, while the exponent is invariant, the constant is increased with respect to copper and increased further with decrease in temperature. The life in Cu-7.5 wt pct A1 is double that of copper at room temperature and about four times at 78°K. Fractoflaphic studies indicate that these diffevences in life are associated with the early stages of cracking. In copper, severe surface rumpling provides crack nuclei at pain boundaries independent of&apos; testing temperature or the history of the metal. In Cu-7.5 wt pct Al, the planar slip mode of deformation prevents grain boundary folding; crack nucleation is therefore mainly transgranular and the increase in life with respect to copper is associated with a prolonged linking-up stage of such transgranulay cracks. At the lower temperature the increased difficulty of cross slip inhibits multiplication of trans-granular cracks and retards the linking-up process, giving a further enhaancement in life. It is expected that such behavior will be found general in materials of wavy and planar slip mode, respectively. STAGE II crack propagation in high strain fatigue takes place by a repetitive process of crack tip blunting in the tension part of a fatigue cycle followed by resharpening of the crack in the compression part.1-3 This mechanism has come to be called the "plastic relaxation process".3 At very high strains (lives to failure less than 103 cycles), most of the life of a specimen is spent in crack propagation by the plastic relaxation process.4,5 Since the differences in the cyclic strain hardening behavior of metals and alloys are comparatively small in this strain range,6,7 the amount of blunting at a crack tip depends primarily on the applied strain amplitude. The rate of crack propagation in this life range is therefore nearly constant with respect to material and temperature of testing. In the longer life portion of the high strain fatigue range (l03 to 103 cycles to failure), a considerable fraction of the lives is spent in crack initiation and Stage I growth. Grain boundaries have frequently4,5,8 been found to be the sites for this growth which is also associated with a surface folding process.4&apos;5 Such folding behavior may also be general in any material having a grain structure and capable of plastic deformation. It has therefore been suggested475 that the universality of the Coff in-Manson law,9&apos;10 constant,* where ?Ep is the plastic strain range and Nf the number of cycles to failure, is reflected in the approximate invariance of both the crack initiation and propagation behavior. It happens that, in most investigations of the phenomena associated with failure,4,5,8 materials with a wavy slip mode have been used as vehicles for the studies, or else "normally" planar slip mode materials have been tested under conditions giving rise to wavy glide deformation, e.g., stainless steels cycled at high temperature.&apos;&apos; The purpose of the present investigation was to determine whether or not different cracking phenomena would occur when the slip mode of a material is planar and if the associated life to failure would be different than that of a wavy slip

Jan 1, 1968

By Ulf S. Landergren

Diffusion coefficients and marker movements have been determined in brass using welded couples. Three different concentration ranges were employed at 750°C, while a fourth concentration range was measured at 500°, 600°, 700°, and 800°C. Marker displacements toward and porosity development in the high zinc side of the couple were observed in all cases. The results were interpreted as favoring a vacancy diffusion mechanism. IN the large amount of diffusion literature which has accumulated during the past half century, very few experimental data have been presented on body-centered-cubic metallic phases. The impressive combination of experimental and theoretical studies which has given the vacancy mechanism a fairly secure position as the preferred mechanism for most face-centered-cubic metals has no counterpart for body-centered-cubic phases, although the theoretical picture seems to be in a better state than the experimental work. The attempts to select diffusion mechanisms by the theoretical calculation of the activation energy for the various conceivable unit diffusion processes and the energy of formation of the type of lattice defect required have been based on models of two kinds. The first kind utilizes potential functions fitted for certain properties of pure crystalline copper; the second kind is fitted for sodium. The latter is an example of an open metal in which the metallic ion is small compared with interatomic distance,&apos; while the former is characteristic of the larger class of metals in which the ions have diameters about as large as the interatomic distance. brass should resemble copper in this fundamental property. It should be evident that results on the mechanism of diffusion obtained for one of these classes need not necessarily apply to the other. Huntington and Seitz, using a model of the copper type, showed that the activation energies to be expected for interstitial migration or a two-atom direct interchange were much higher than that for vacancy diffusion. Zener pointed out that a four-atom ring rotation would involve an activation energy considerably lower than the two unlikely mechanisms although still higher than the vacancy mechanism in the face-centered-cubic structure. Furthermore, he argued, the four-atom rotation would be relatively easier in a body-centered-cubic lattice and might be competitive with the vacancy mechanism. Paneth showed that a body-centered-cubic crystal-like sodium could form interstitial defects of special type, called a crowdion, by crowding an extra atom into < 111 > directions where the defects would affect a line of about eight atoms. The defect was constrained to move along the defining line. Unlike the face-centered-cubic case, therefore, the theoretical calculations for the body-centered-cubic mechanism do not provide a very clear-cut answer, and the absence of adequate experimental studies is the more serious. This investigation was undertaken for the purpose of obtaining useful information on diffusion in body-centered-cubic phases. brass was chosen for several reasons: 1—the phase field is sufficiently wide at high temperatures to permit the use of diffusion couples with adequate concentration ranges for satisfactory chemical analysis, 2-—thermodynamic activity data are available for the calculation of mobilities from diffusion coefficients, 3—the phase is as close to the theoretical model of the copper -type on which the calculations have been done as any body-centered-cubic metallic phase is likely to be. Measurements of the movements (if found) of inert markers during the diffusion process would indicate that a vacancy or Or mechanism played a major role in diffusion in brass. If no marker movement should be found, a ring mechanism would be indicated but not assured. Experimental Procedure The alloys used in these experiments, listed in Table I with analyses, were supplied by the Naval Research Laboratory. They were from ingots of 3x3 in. cross-section, heated to 760°C, reduced under a

Jan 1, 1957

By L. L. Handy

The most common method of identifying hydrocarbon-bearing strata in a well that penetrates many different formations involves measurement and interpretation of the electrical properties of the formations as determined by electrical logs. Even though this method is used extensively, and even though in a great many instances it is capable of indicating presence of oil or gas, situations arise for which it is extremely difficult, if not impossible, to deduce the presence of hydrocarbons. These situations may involve the following. 1. A thin formation, bounded by highly resistive formations, in which it is impossible to obtain the actual resistivity of the uninvaded zone with existing logging devices. 2. A formation in which invasion has been so extensive that a value for the uninvaded zone resistivity cannot be obtained. 3. A very shaly formation in which the resistivity index, I, is lower than that usually associated with productive formations. 4. Laminated formations comprised of thin productive sands separated by thin shale streaks in which the individual sand and shale streaks are too thin to permit measurement of uninvaded-zone resistivity with existing logging devices. 5. Productive formations in which the water saturation is high. To extend the utility of electric log interpretation to identification of hydrocarbons in all types of formations, there is strong incentive to find a method not subject to these limitations. Some time ago, in connection with research on the wettability of reservoir rock, an investigation was conducted in which the resistivities of cores were measured shortly after they were removed from a core barrel, and again after they had been extracted and restored to their original oil and brine saturation.&apos; The resistivities after extraction were generally lower. Other tests made on the cores indicated that they were more nearly water wet after they were extracted; thus, it was assumed that the observed changes in resistivities were due to a change in wettability of the cores. Other experiments&apos;," have shown that resistivities of rock samples are sharply dependent on wettability. These experiments have shown that oil-wet samples are more resistive than water-wet samples. To obtain an understanding of how the wetting properties of the surfaces of core material affect electrical resistivity, a series of experiments was conducted. Two groups of core samples were prepared for testing. One group contained brine, but no residual oil. The other group was saturated with brine, flooded with oil to a low water saturation, then flooded with brine to a final residual oil saturation. Resistivity measurements were made on each group. Both groups were then flooded with the original brine to which a chemical had been added that renders sand and clay surfaces preferentially oil wet, a so-called reverse-wetting agent. very little change in resistivity was observed in cores containing only water. The group containing residual oil, however, showed resistivity increases of 100 to 200 per cent. These experiments showed that the resistivity of a core containing oil could be altered by changing wet-tability of the core. Moreover, the possibility was introduced that reverse-wetting agents might be employed as the basis for a logging method for identification of oil-bearing strata. Since behavior of a porous rock containing gas and water might be expected to be similar to that of a rock containing oil and water, such a method should also be applicable to identification of gas-bearing zones. In principle the wettability of the invaded zone could be reversed without altering conductivity of the interstitial water or the hydrocarbon saturation therein. Those strata showing significantly increased invaded-zone resistivities would, therefore, contain hydrocarbons; those with no significant change would be filled only with water. Addition of a reverse-wetting agent to a hydro carbon -bearing zone which is, by nature, already preferentially oil wet would not result in an enhancement of its resistivity. It is generally believed, however, that most hydrocarbon-bearing strata are preferentially water wet.

By Keith R. Bock, Norman Parlee, Albert M. Barloga

Coils of pure iron and iron-carbon alloy wire (0.05 to 0.80 pct C) and sufficient sulfur to saturate the solid phase were equilibrated in evacuated or argon filled tubes. After rapid cooling, and removal of the outside nonmetallic layer, the wires were analyzed for carbon and sulfur and the data used to construct an Fe-S binary and isotherms of the Fe-C-S ternary in the range 1149° to 1427°C. THE solid solubility of sulfur in steel is of interest in connection with such phenomena as hot shortness, burning, and so forth. "Burning," the more or less permanent damage that some steels suffer when heated for forging or rolling, has been shown to be related closely to the behavior of sulfur and less closely to carbon and oxygen.&apos; Attempts to interpret burning phenomena in steels fail because of lack of data on the Fe-C-S and the more complex systems in this family. Rosenqvist and Dunicz 2 and Turkdogan, Ignatowicz, and pearson3 have largely elucidated the Fe-S diagram in the region of interest but no information on the Fe-C-S diagram in this region appears to be available in the literature. This paper deals with the elucidation of the Fe-C-S diagram in these interesting ranges. The method employed is different from those used by previous workers2&apos;3 on the Fe-S system. EXPERIMENTAL METHOD Pure iron wires (Ferrovac E) or iron carbon alloy wires of 1 mm in diam were cleaned with acid and acetone, coiled, and placed in silica tubes (7 mm OD and 5 mm ID) previously closed at one end. Enough sulfur was added to assure saturation of the solid iron phase. The filled tubes were either simply evacuated and sealed, or filled with argon at a reduced pressure and sealed. The argon was required at the higher temperatures to prevent collapse of the tubes. The filled and sealed tubes were placed in the middle uniform temperature zone of a Globar tube furnace and equilibrated at temperatures ranging from 2100°F (1149°C) to 2720°F (1493°C). After equilibration the tubes were removed from the furnace and quenched in air or water, the form of quenching being found to have no effect on the results. The tubes were broken open and the coils were placed in a 1 : 1 HC1 solution to remove the sulfide-rich layer. The coils were then cut into small pieces and analyzed for sulfur and carbon. In the early stages of the investigation different equilibration times ranging up to 17 hr were tried and the cores of the wires were analyzed to test for saturation. One hour appeared to be sufficient to reach maximum sulfur content at 1454°C and 2 hr sufficient at 1149°C. The practice adopted was to use at least 3 hr at the higher temperatures and at least 5 hr at the lower temperatures. The iron-carbon alloy wires used were made by carburizing pure iron wire with carbon monoxide gas in a one inch diameter ceramic tube at about 1204°C. Differing carbon contents were obtained by allowing fairly large coils to react with the gas for varying lengths of time at a flow rate of 800 cc per min. The reacting times ranged from 15 min to 4 hr depending upon the amount of carbon desired. The furnaces used were controlled by means of Leeds and Northrup Speedomax Type H instruments with D.A.T. Control attachment. Thermocouples were calibrated against the melting points of gold and copper. The temperatures recorded appear to be accurate within i2.7OC. Starting with run No. 85 and continuing with the

Jan 1, 1962

By W. M. Parris, P. D. Frost, H. A. Robinson

The effects of hydrogen content and strain rate on the static tensile and notch-rupture properties of the Ti-3Mn-complex alloy heat-treated to differed strength levels were investigated. The extent of embriulement was found to be proportional to the amount of hydrogen present and the strength to which the specimens were heat-treated. The research suggests the possible necessity of specifying lower hydrogen limits than are now allowed when titanium alloys are heat-treated to higher strengths than are presently used. The research also supports a hypothesis advanced earlier for a mechanism of hydrogen embrittlement. MUCH attention has been focused on the effects of hydrogen in in Most of the published in- formation has dealt with the effects of hydrogen on alloys in the annealed or lowest-strength condition. The work of Kotfila and Erbin3 is one of the few exceptions. In view of the fact that heat-treatments are now being used to produce high strength in titanium alloys,&apos; it is necessary that more attention be given to the effects of hydrogen on these alloys in the heat-treated condition. Research was initiated to establish the effects of hydrogen on the properties of the Ti-3Mn-complex alloy at various strength levels. This research, which supplemented that of Kotfila and Erbin, was started because it was found that high-strength tensile specimens of this alloy were strain-rate sensitive. The essential results of the research are summarized in this paper. MATERIALS AND EXPERIMENTAL PROCEDURES Material for this investigation was a 41/2-in.-diam bar forged from a single commercial heat of the Ti-3Mn-complex alloy having the composition given in Table I. The as-received forging was reduced further by forging at 1750°F to 3/4-in. square bars. About one fourth of this material was vacuum annealed at 1400°F for 24 hr to obtain a base hydrogen content of 25 ppm. The remaining bars were hydrogenated to levels of 90, 140, and 260 ppm in a large Sievert&apos;s apparatus in the following manner: Four 3/4-in. square bars, 5 in. long, were placed in the reaction chambzr. The chamber was evacuated, heated to 1400°F, and a measured amount of hydrogen was introduced. When absorbtion was complete, as indicated by the drop in pressure, the reaction tube was cooled and the bars removed. Each batch of bars was then sealed in individual evacuated Vycor containers and heated for 24 hr at 1400°F in order to distribute the hydrogen homogeneously. After vacuum-annealing or hydroge nation, the bars were rolled at a temperature (1400°F) below the ß-transus to 1/2 - in. - diam rounds. Groups of the rolled specimens representing each of the hydrogen levels were given the following heat-treatments: 1) One hr at 1300°F, water-quenched. 2) One hr at 1300°F, water-quenched, and aged 8 hr at 1100°F. 3) One hr at 1300°F, water-quenched, and aged 4 hr at 1000°F. 4) One hr at 1300°F, water-quenched, and aged 48hr at 800°F. Heat-treated bars were machined into standard 0.250-in.-diam tensile specimens having 1-in. gage lengths. Throughout the processing and testing of this material, hydrogen analyses were made by vacuum-fusion methods to check for accuracy of hydrogen additions, homogeneity, and possible losses during processing. A summary of these analyses is given in Table II It may be seen that the actual hydrogen contents were reasonably close to the desired levels. In addition, it was noted that the variation in hydrogen content within any one bar did not exceed

Jan 1, 1959

By H. F. Webster, C. G. Dunn

Rolling md annealing procedures are described for developing the (110) preferred orientation in tantalum for use as a thermionic-emission materinl where electrodes of a uniform high work function are needed. The development of a (110) [001] annealing texture is also reported. A preferential growth process due to low (110) gas-metal surface energy is considered in explanation of the observed growth of (110) grains. ThE growth of (110)-oriented grains in thin strips of rolled zone-refined tantalum, previously discussed briefly,&apos; is reported together with more recent information on the formation of a (110)[001.] annealing texture, i.e., one where the [001] direction of each grain lies close to the rolling direction. The development of tantalum sheet with a (110) preferred orientation (position of [001 ] directions unspecified) is of importance as a thermionic-emission material in thermionic-energy converters. Since therm ionic-emission properties of cesium-coated crystals vary over a wide range depending on (hkl) a2-7 a single (hkl) preferred orientation is more useful than either several preferred orientations or none. Furthermore, the (110) preferred orientation yields the highest thermionic-emission density for a cesium coverage of less than a monolayer.5-7 A consideration of the stability of the grain structure needed at high operating temperatures suggests that the preferred orientation itself should be produced by annealing processes. Information on annealing textures in tantalum, however, is rather limited. Pugh and Hibbard 8 obtained (111)[112] and (lll)[112] components in a 1400°C recrystallization texture in 95 pct cold-rolled foil; these components were enhanced on annealing at 2500°C. Muller9 rolled 99.97 Ta to 90 pct reduction in thickness in the final stage and noted only a slight sharpening of the rolling texture, with some intensification of the (111)[112] component upon annealing at 1200°C. The experimental procedure followed in the present work stemmed from the idea of using possible differences in (hkl) surface energies to provide a small driving force for growth of grains having the lowest (hkl) surface energy.l0* This would mean using high-purity tantalum processed into thin sheet form and annealed at high temperatures either in a noncontaminating atmosphere or in one producing additional sample purification. EXPERIMENTAL Zone-Refined Tantalum. National Research Corp. tantalum in the form of 1/8-in.-diam rod and of 99.98 purity (53 ppm 0, 15 ppm N, and 8 ppm C) was zone-melted four times. The traversing speed was 4 mm per min. Such treatment increased the resistance ratio, R293k/R20K , where R293K is the resistance of the specimen at 20°C and R20K the re-resistance at liquid-hydrogen temperature, from a value of 25 for the initial rod to 175 for the zone-refined rod (the increase in the resistance ratio provides a measure of purification11,12). Part of a zone-refined rod was rolled unidirec-tionally to a strip 0.003 in. thick and 0.28 in. wide. Samples were given anneals in the range 1300" to 2100°C using ac conduction heating in an oil-free vacuum system. Recrystallization to small grains appeared to be complete at 1300°C. Large grains were formed at higher temperatures, but these were not entirely strain-free, some grains showing asterism in Laue X-ray patterns and some displaying a range in orientation. Part of the difficulty was caused by thermal expansion of the tantalum strip as it was held between fixed supports. However, in a different arrangement with the tantalum freely suspended in a tantalum can, high-frequency induction heating in vacuo at about l0-6 Torr produced

Jan 1, 1964

By P. L. Raffo, T. E. Mitchell

Single crystals of Ta-Re alloys have been deformed in tension and compression as a function of solute concentration and temperature. The temperature dependence of the yield stress is decreased by alloying. At low temperatures there is a minimum in the yield stress (at about 2 at. pct) as a function of concentration. Three-stage hardening was not observed in the alloys. Linear hardening rates, where observed, were characteristic of stage II in pure tantalum. At low temperatures and high concentrations, work softening was often prominent at high strains. Electron microscopy revealed that dislocations are mostly of screw orientation, although heavily jogged and cusped, and often occurring in dipole pairs. The alloys become more prone to cleavage with increasing rhenium content, due mostly to an increasing yield stress rather than a change in cleavage strength. Alloying apparently decreases the Peierls stress of bcc metals at low temperatures by facilitating kink nucleation at solule atoms on screw dislocations. Higher solute concentrations further increase the yield stress. Quantitative calculations and comparisons with other bcc alloy systems are made to show that the predominant hardening mechanism is due to the interaction of dislocations with misfitting solute atoms close to their slip planes. The interaction is so localized that the core dilatation of dislocations could be important. CONSIDERABLE interest has arisen in recent years concerning the mechanical properties of unalloyed bcc metal single crystals.&apos;-7 These and similar studies on polycrystalline materials have resulted in a considerable understanding of the macroscopic yielding characteristics8, 9 and the process of work hardening.&apos;&apos; This latter advance is due in a large part to the increased use of transmission electron microscopy in P I RAFFO junior Membar AlMF is Metallyureist Iewis such studies and the general observance of three stages of work hardening in the shear stress-shear strain curves. The properties of substitutional alloy single crystals of bcc metals have not been determined in as much detail as in the pure metals. Insufficient systematic studies have been made on the effect of alloying on the magnitude and temperature dependence of the yield stress or the stress-strain curves. The alloy systems which have been studied to some extent have been Ta-W,11, 12 Ta-Mo," Ta-Nb,12&apos;14 Mo-Re,15 W-Re,16 and alloys of iron with nickel,17 platinum,17 vanadium,18 and aluminum.18 A large body of data exists on Fe-Si alloys, but these have mainly been a result of the experimental convenience of using this material in place of unalloyed iron, rather than a comprehensive study of the mechanical properties of the alloy system per se. The present authors have previously reported additional data on the properties of Ta-Re, Ta-W, and Ta-Mo alloy single crystals.1g Some of the data on Ta-Re alloys will be described later in this paper. Briefly, single crystals of Ta-Re (1.2 to 10.6 pct), Ta-W (6.2 and 10.2 pct), and Ta-Mo (2.9 and 9.0 pct) were grown by the floating zone method in vacuum. Compression tests were performed in the temperature range 77" to 623°K. Slip was found to occur in Ta-Re alloys at 77&apos; and 298°K on the (~ll)[lil] system and slip lines were coarser and straighter than in unalloyed tantalum. No three-stage hardening was observed in any of the alloy systems and the curves were nearly parabolic. At temperatures above 175°K the yield stress* increased monotonically with in- creasing composition while at 77" and 175°K the yield stress initially passed through a minimum after which it began to increase again for concentrations greater than about 2 pct Re. Similar minima were observed in Ta-Mo alloys and to a lesser extent in Ta-W alloys. The minimum was interpreted as being a result of a

Jan 1, 1969

By Richard H. Hopkins, R. W. Kraft

A new technique facilitating the rapid determination of the three-dimensional morphology of phases suspended in an opaque matrix is described. The vapidity of the technique is based upon continuous cinephotomicrographic recording of specimen microstructure as the specimen surface is removed under controlled conditions in an electrolytic cell. The advantages of the technique for determining spatial morphologies in opaque materials are illustrated by application of this method to the study of the configuration of faults in a specimen of directionally solidified CuAl2-Al eutectic alloy. A spatial model of the eutectic fault network deduced in this manner has substantiated the facts that faults in this alloy tend to lie parallel to one another in a direction close to that of the specimen growth axis and that the majority of faults tend to bridge lamellae rather than run parallel to them as growth proceeds. For the specimen studied fault density decreased as growth proceeded. ALTHOUGH metallographic observation is a powerful tool for investigating the correlation between structure and properties in metallic alloys, ordinary metallographic methods suffer the fundamental limitation that any conclusions concerning the three-dimensional structure of these opaque materials must rely upon inference from two-dimensional microsections. For this reason it is frequently impossible to make valid judgments regarding the true spatial distribution of the micro-constituents in specimens under investigation. Several techniques have been devised to circumvent the basic limitation of planar microsec-tioning and employed with varying success. These include construction of three-dimensional models of alloy microstructure from sequences of photomicrographs,1-3 quantitative deduction of spatial features based upon statistical studies of data from random and/or selected plane microsections,4-6 and selective dissolution of one phase in a two-phase alloy to reveal the topography of the phase remaining.&apos; However these techniques also have certain obvious disadvantages. Construction of three-dimensional models is a very time-con- suming procedure due to the tedious polish-etch-photograph cycle needed to obtain the necessary photomicrographs to portray the spatial microstructure. Statistical methods by their very nature give only average values of the quantities being studied and cannot therefore give a detailed physical picture of the microstructure to be observed. Finally, the technique of selective dissolution is only applicable in systems where the phases are of vastly different chemical reactivity; it is of little use in the study of single-phase alloys. Since the determination of the three-dimensional structure of opaque materials is hindered by two facts—the materials are opaque to radiation in the visible spectrum and construction of spatial models of microstructure is quite laborious—two solutions to the problem suggest themselves: 1) use radiation of energy great enough to render the alloy transparent, i.e., some microradio-graphic method, or, 2) devise a rapid means to section and record alloy microstructure. Stereomicroradiography has indeed been used to study internal arrangements of phases in several materials8,9 but the method, in general, is limited to very thin specimens (due to X-ray absorption) and simple microstructures. It would seem then that development of a rapid, simple, and generally applicable technique for studying the details of alloy microstructures in three dimensions can best be undertaken by following the second approach above. Using this idea a new technique for rapid evaluation of spatial microstructures was developed based upon cinephotomicrographic recording of the micro-structure of a specimen undergoing controlled electrolytic dissolution. It was further felt that application of this technique to the careful study of the spatial configuration of faults in a CuAl2-Al eutectic alloy would yield information that could be correlated with the extensive data3,6,10-12 already compiled concerning the solidification and crystallography of this eutectic. The practicality and usefulness of similar methods for studying three-dimensional features in alloys has recently been illustrated by the study of dislocation loops in zinc13 and spheroidized structure in steel.14

Jan 1, 1965

By H. J. McQueen, G. C. Kuczynski

ThE useful properties of k6-Alfenol and the thermenol alloys (summarized by J. E. Nachman and W. J. Buehler&apos;) have encouraged research into those alloys which, having compositions near 25 at. pct Al, 75 at. pct Fe, undergo an ordering transformation between 500° and 550°C. The equilibrium diagram of the iron-rich iron-aluminum alloys was first determined by Bradley and Jay.2,3 They distinguished three structures: the disordered, which is bcc; FeAl, which is simple cubic with one Fe atom and one Al atom per lattice site; and Fe3Al, which is fcc with three Fe atoms and one Al atom per lattice site. This last structure may be conceived as eight body-centered cubes with Al atoms at the centers of alternate ones, so that Al atoms have only Fe as first and second nearest neighbors. Later, Osawa and Murata4 determined the boundaries of the different superlattice regions in the iron-rich Fe-Al solid solution by dilatometric and X-ray methods. Two of the six phases, which they found between 0 and 50 at. pct Al, were never identified by other workers. A. Taylor and R. M. ones,&apos; concurrently with and independently of this research, redetermined this diagram by measuring lattice parameters at temperature. Further reference will be made to their results. EQUILIBRIUM DIAGRAM A) Experimental Procedure—The alloys were prepared from NRC gas-free iron (99.9 pct Fe) and 99.99 pct Al. They were melted in alumina&apos; crucubles under vacuum in an induction furnace. If the ingots were unsound or showed more than a 3 pct weight loss, they were rejected. The ingots were homogenized at 1100°C for 24 hr before hot-rolling at 700" to 900°C from 1 in. to 0.25 in. thick. It is interesting to note that alloys containing more than 28 at. pct Al were able to take an 80 pct reduction of area in the ordered FeAl region. The dilatometer specimens and specimens for chemical analysis were cut from these plates. A Leitz dilatometer was used to measure the coefficient of linear expansion as a function of temperature. The specimens were slowly cooled from 900" to 300°C where they were annealed for 10 hr before slowly cooling to room temperature. The expansion was recorded photographically as the alloys were heated and cooled at 1.5oC per min between 100" and 1000°C. The magnetic measurements determining the Curie points were performed on the dilatometer specimens in the dilatometer, by using direct current in the furnace windings both to heat the specimen at 1.5°C per min and to serve as the primary coil. The regular specimen holder was replaced by a thin alundum tube on which was wound, as a secondary, a closely spaced coil of fine copper wire. A reversal of current in the furnace windings produced, in the secondary, an emf which was proportional to the induction in the specimen, and which could be detected on a galvanometer. In this way, it was possible to determine the Curie points of alloys to within * 5°C. B) Results and Discussion—The coefficients of linear expansion, as obtained from the dilatometer curves, are plotted as a function of temperature and composition in Fig. 1. Since in many cases the heating and cooling curves were superimposed, it was decided to measure slopes of heating curves only. Definite high peaks, due to the transformations, were observed on each curve between 400" and 520°C, and, for alloys containing more than 23 at. pct Al, a second set of smaller peaks at temperatures above 660°C. The sharpness of the peaks indicates that ordering to Fe3Al is a first-order reaction. The compositions and the temperatures of these peaks are listed in Table I. When these temperatures are plotted as a function of the composition, the result is the constitution diagram presented in

Jan 1, 1960

By W. A. Tiller, S. O’Hara

Grain refinement in stirred melts has previously been shown to arise from dendrite segmentation. The present work discusses experiments capable of distinguishing between remelting and mechanical effects which result from fluid motion. Nucleation of still un-dercooled Sn-1 pct Pb melts with a single-crystal seed shows that the grain size of the resulting ingots decreases as the undercooling at nucleation increases up to about 1Z0C, then increases again for larger undercoolings. Pure melts exhibit little refinement as a result of nucleation at any undercooling up to 12°C. The presence of fluid motion gives rise to grain multiplication in pure melts, but the magnitude of this effect is partly obscured by solid-state pain growth. The results indicate that two grain-multiplication pvocesses can operate, one which involves dendrite remelt ing, and another associated with fluid motion. Small ingot experiments indicated that an equiaxed transition can be caused by fluid motion. THE presence of forced convection during recales-cence of undercooled Sn-1.3 pct Pb melts has been shown to give rise to extensive grain refinement during solidification of small ingots of this alloy.&apos;&apos;2 In Ref. 2, this refining was shown to result entirely from dendrite fragmentation and not from any nucleation mechanism. The purpose of the present work is to establish what processes lead to dendrite segmentation. Experiments were therefore carried out to determine the resulting structure of both pure tin and Sn-1 pct Pb melts, nucleated at various undercoolings by single-crystal seeds in the absence of forced or natural convection. By comparing the information obtained from these experiments with the results of Refs. 1 and 2, the effects of fluid motion on the ingot structure, and the details of the segmentation mechanisms, should become clear. The results of these experiments should also be helpful in understanding the solidification structures of larger ingots since natural convection, which is known to play an important role in determining ingot structure13 differs only in degree from the fluid-motion effects studied here. I) EXPERIMENTAL PROCEDURE The apparatus shown in Fig. 1 was constructed to permit controlled nucleation of stirred or unstirred supercooled melts using a seed oriented in the tin dendrite growth direction. Because the melts, in all cases, were capable of larger supercoolings than the one used in each experiment, the structure observed always resulted from a single nucleation event at the surface. The sequential procedure was as follows: Tin was melted under purified hydrogen and the oxide removed at a temperature of 450°C . Electromagnetic stirring was then used to clear the surface completely and during this time the guide tube was lowered underneath the melt surface. A pointed Czochralski-like crystal was then propagated from a specially prepared tin dendrite seed4 inside the guide tube while the temperature of the bulk melt was held just above the melting point, ensuring a perfectly clean seed which was epitaxial with the original dendrite. The dendrite and epitaxial seed, still in the guide tube, were withdrawn from the melt which, in turn, was undercooled to 2C° below the temperature at which the nucleation experiment was to be performed. The melt was held there for several minutes, both in stirred and in still-melt experiments. The temperature was next raised to the required value and the melt nucleated by the carrot seed inside the guide tube, the heater power being shut off just after nucleation. In experiments where quenching after nucleation was required, a high-density graphite crucible was used, but otherwise the procedure was as described above. The temperature gradient in the melt before nucleation was kept below 1C° per cm by carrying the experiment out just below the center of the heater shown in Fig. 1.

Jan 1, 1968

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By C. S. Matthews, P. Hazebroek, F. Brons

A method has been developed for calculating the average pressure in a bounded reservoir. The reservoir is first divided into the individual drainage volumes of each well, by using the criterion that at steady state each individual drainage volume is proportional to a well&apos;s production rate. The average pressure in each drainage volume is then calculated by a method developed in the report. By volumetrically averaging these individual drainage volume pressures, the average pressure in the entire reservoir is obtained. To calculate the average pressure in each drainage volume, a correction is applied to the ordinary extrapolated pressure, i.e., the pressure obtained by extrapolating to infinite time the linear portion of the graph of closed-in pressure versus log[?t/(t + ?t)], where At is the closed-in time and t the production time. The correction, which is a function of the production time, is presented in graphical form for different shapes of the drainage area (horizontal cross section of the drainage volume). INTRODUCTION It is important to be able to find the volumetric average pressure in a reservoir so that the size of the reservoir may be determined from material balance calculations. It is also desirable to be able to find the approximate distribution of pressure within a reservoir for detection of fluid movement. &apos;The purpose of this paper is to present a method for calculating both the average reservoir pressure and the approximate distribution of pressure within a bounded reservoir — that is, a reservoir with no water drive. In reservoirs where the pressure builds up rapidly after wells are shut in, the determination of average pressure generally poses little problem, for one often need only average the final buildup pressures. It is when pressure buildup is slow that difficulties arise. For practical and economical reasons, the time allowable for closing in wells is limited. If at the maximum allowable closed-in time the pressure has not reached a con- stant value (and this is more often the case than is generally realized), calculation of average pressure presents difficulties. One method which has been used in the past for obtaining an "average" pressure involves the extrapolation to infinite closed-in time of curves of pressure as a function of closed-in time. Contour maps of these extrapolated pressures are used in arriving at an average pressure. In some cases wells are shut in for arbitrary Lengths of time, and contours of these arbitrary shut-in pressures are used for calculation. All these methods are somewhat arbitrary and suffer from lack of fundamental theoretical basis. For example, no precise meaning can be given either to the extrapolated or to the arbitrary shut-in pressure. Furthermore, there seems to be no reason why either set of pressures should be contourable. It was because of the vagueness of these old methods and their lack of a sound theoretical basis that, the present method for determining average reservoir pressure was developed. SYMBOLS A = drainage area of a well, cm&apos; A, = drainage area of well j, cm2 A , = total drainage area of all wells in a reservoir, cm2 c = coefficient of compressibility of fluid at reservoir temperature, atm-1 f = hydrocarbin-filled porosity, fraction h = net thickness of formation, cm k = permeability of formation, darcy p = pressure at a well, atm p0 = initial reservoir pressure, atm p* = pressure obtained by straight-line extrapolation of linear portion of plot of p versus log [At/(t+ At)] p = volumetric average pressure inside drainage volume of a well q = volume flow rate at a well, cc/sec a: prevailing reservoir conditions q, = total volumetric flow rate from a reservoir, cc/sec at prevailing reservoir conditions rw = radius of well bore, cm T = kt/fµcA t = corrected time of production of a well, seconds, defined by

Jan 1, 1955

By R. W. Kraft

A comprehensive analysis of the crystallographic and metallographic structure of several Al-CuAl, lamellar eutectic grains solidified under equilibrium conditions has been performed. A consistent pattern was discovered which was interpreted in terms of a detailed atomistic model of an interface describable as: Interface II {111} AL 11 {211} <101>Al <120> ? Several arguments are presented to substantiate the theory. It is compared in detail with another interface model characteristic of the Wid-manstatten precipitate of CuAl2 in A1 at low Cu levels. In addition, a concept is proposed to help explain the structure of the compound CuAl2. It has been shown that the eutectic alloy between A1 and CuA1, can be forced to solidify as substantially parallel lamellae throughout a relatively large volume&apos; by appropriately controlling the solidification parameters of growth rate, thermal gradients, and impurity level and by causing solidification to proceed in a unidirectional fashion from one end of an ingot to the other. A typical ingot solidified in this manner shows three distinct zones along the length of the ingot, Fig. 1. At the head end is a region consisting of many different grains, some of which are elongated slightly in the direction of solidification. After a few centimeters of growth, this polygrained region gradually merges into the second zone which consists of a very few elongated coarse grains. The third zone at the tail end is a region of eutectic colonies caused by segregation of impurities. Microscopic examination of transverse and longitudinal sections of the ingots reveals that during the initial stages of growth (zone I) the lamellar structure is frequently highly distorted and twisted but that it straightens out in zone 11. This characteristic pattern suggests that nucleation of a eutectic grain begins randomly at the head end of the ingot, and that some particular orientation relationship between the phases is more favorable from an energy consideration than others for one reason or another, so that some grains grow and others which are less favorably oriented do not. The favorably oriented or "low energy" grains continue to grow until excess solute rejected at the advancing liquid-solid interface destroys the parallel lamellar structure by forming a cellular rather than a substantially planar interface. When this occurs, the parallel lamellae fan out and eutectic colonies are formed because the lamellae grow perpendicular to the local interface. Questions naturally arise as to the crystallographic nature of the "low energy" grains. What do the grains which persist into zone II have in common? Can this explain the stability of such grains during growth? To answer the first question it is necessary to determine the crystallographic orientation of the lamellae of each phase relative to one another. This has been done for Al-CuAl2 eutectic alloys by two different techniques on specimens prepared in different ways.2&apos;3 The results were different from each other and from a third orientation relationship between A1 and CuA1, reported for overaged A1-Cu alloys.4 The three orientation relationships are summarized below to illustrate their similarities and their differences. (001) ? 11 (001) A1; [l00] ? 11 [l00] Al Eutectic, Ref. 2 (001) 8 1) (001) A1; [1001 ?11 [1101 A1 Eutectic, Ref. 3 (001) ? 11 (001) A1; [l00] ? 11 [310] Al; Widmanstfitten, Ref. 4 It is not surprising that different orientation relationships were found by Ellwood and Bagley,&apos; and Takahashi3 since the solidification conditions were quite different. Nor is it surprising that Mehl, Barrett, and Rhines 4 found still another value since they were studying another type of alloy. Nevertheless, it would be expected that some logical and systematic pattern could be found if all the alloys were at equilibrium since the same two phases are

Jan 1, 1962

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system&apos;s behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell&apos; have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine&apos; have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.&apos; A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By J. Geertsma

The relationship between porosity and the speed of propagation of acoustic waves in fluid-saturated porous rocks as measured by the Sonic log and by ultrasonic tecbniques is analyzed. Biot&apos;s continuum theory 1 is used to explain the difference in acoustic wave propagation between a dry and a liquid-saturated porous material. The porosity is a variable in this theory. However, the acoustic wave propagation in the dry rock depends too on porosity, and this dependence is not predicted by the theory. Frequently in dry sandstones, a nearly linear relationship between reciprocal acoustic wave velocity and porosity is observed in the low-porosity range. The physics behind this behavior is outlined. An empirical relationship of the form, 1/v = A + F, applies accordingly for many porous dry rocks, provided the porosity is the only variable. The presence of a liquid in the pores changes the value of B, and this change is found to be in agreement with the Biot theory. The time-average relation introduced some years ago2 results in an equation of the same type — 1/V = F/Vf + (1 - F)/Vr — but is not based on a sound physical picture. Still, this relation sometimes predicts approximately correct A and B values. Carbonate rocks with their complicated pore structures sometimes show a different relationship between wave velocity and porosity, unfavorable for log interpretation. Examples are presented. The simultaneous presence of calcite, dolomite and anhydrite, with their different grain densities and matrix compressibilities, complicates acoustic-log interpretation in carbonate rocks still further. Other complicating effects of acoustic-log interpretation are discussed. They are related to the influence of shale streaks and natural fractures on the average wave velocity observed by the logging tool and to the effect of adsorption phenomena on wave propagation in unstressed rocks particularly in sandstones. INTRODUCTION The velocity of propagation of sound waves in porous sediments as a function of depth is a quantity frequently measured. Velocity loggers of various d-sign may be used, for which velocity is defined as the. distance between a wave generator and a wave detector (one-receiver system) or the distance between two wave detectors (two-receiver system), divided by the shortest time required for a vibrational pulse to cross this distance. Velocity loggers ale used to assist in seismic prospecting, to differentiate between the various types of sedimentary rock layers, and also to determine rock porosity and pore fluid content. This paper is especially concerned with the relation between velocity and porosity and the effect of the pore fluid. Most velocity-log interpretation in terms of porosity is based on a formula advocated by Wyllie, et al,2 suggested earlier by Hughes and Jones,3 and known as the "time-average relation". This formula app1ies for a model consisting of a layered system of parallel alternating slices of a solid and a liquid, crossed by the wave path perpendicular to the interfacrbs. This must be an unsuitable model for the derivation of wave propagation properties of liquid-saturated porous media. It suggests that only rock matrix and fluid properties influence the wave velocity. The surprising fact, nevertheless, is that applications of this formula to clean sandstones under specific conditions (sufficient effective stress) results in porosity predictions which are close to reality. Other authors, such as Gassmann,4 White and Sengbush,5 Brandt,6 and Hicks and Berry,7 make use of more realistic models to arrive at a wave velocity formula for porous sediments. Their models consist of various liquid-saturated packings of frictionless spherical grains and are, therefore, also highly idealized. A number of experimental results obtained in the laboratory can be better explained with these model theories than by the time-average relation. A study of a general character, not depending a priori on a special model of the porous structure, was still lacking.