Search Documents

By L. V. Pirela, S. M. Farouq Ali

Some interest has been expressed recently in the application of solvents in conjunction with a thermal drive, such as a steamflood. At least one field project of this type has been reported. This paper makes the first attempt to provide some of the basic information necessary for choosing a solvent for such an application. Results of displacements conducted at elevated temperatures are also discussed. This paper shows that, in the case of alcohol-hydrocarbon-water or brine ternary systems, the system miscibility may increase, decrease, or both increase and decrease within the same system when the temperature increases. The distribution coefficient was consistently found to increase in favor of the oleic phase with an increase in temperature, which is advantageous from the standpoint of oil recovery. Results of solvent displacements in a sandstone core showed that even a relatively small temperature effect in the favorable direction (increase in miscibility and distribution coefficient) can lead to a 50 percent increase in oil recovery. On the other hand, if the system miscibility decreases with the temperature, even though the distribution coefficient increases, the oil recovery at the higher temperature may not increase. INTRODUCTION Miscible displacement has been the subject of many investigations. Miscible-phase solvent flooding, employing solvents such as alcohols that are miscible with both oil and water, has received considerable attention. Investigations conducted by Gatlin and Slobod,l Taber, Kamath and Reed,2 Holm and Csaszar,3 and Farouq Ali and Stah14 have helped in understanding the mechanism of solvent flooding. While the use of alcohols and similar solvents as oil recovery agents is economically unattractive, these solvents nevertheless provide a means of understanding the mechanistic aspects of important miscible-phase oil recovery processes such as the Maraflood process (discussed by Gogarty5). As shown by Taber et a1.,2 the phase behavior of a solvent-oil-water system is of .utmost importance in determining the efficiency of the oil displacement. This paper makes the first attempt to obtain the ternary phase behavior data for nine alcohol-hydrocarbon-water/brine systems at elevated temperatures. Results of core tests at elevated temperatures are also presented which corroborate the effect of temperature on oil recovery, as judged on the basis of phase-behavior data. The use of solvents at elevated temperatures has been suggested earlier by Farouq Ali.6 At least one field test of this type has been reportedO7 With the widespread use of thermal recovery techniques, it is possible in some situations (especially in steamflooding) that the use of a suitable solvent in conjunction with the heat carrier may be economically feasible. This paper attempts to provide some of the basic information needed for judging such feasibility. LOCATION OF THE PLAIT POINT The phase behavior of a typical ternary system consisting of a solvent, hydrocarbon and brine (or water) can be represented by a triangular diagram such as the one shown in Fig. 1 (for temperature TI). Point Pon the binodal curve represents the plait point, being the limit of tie lines such as Y1 Y2 The compositions of the coexisting oleic and aqueous phases are given by points Y2 and Y1, respectively, for any mixture composition along Y1Y2. The characteristics of ternary diagrams have been discussed by Findlay8 and others. Taber et aL3 have shown that the position of the plait point in a particular ternary system plays the decisive role in determining whether or not theoleic phase would be displaced in a continuous manner during an alcohol displacement. The area above the binodal curve also is important, since it represents

Jan 1, 1969

By D. H. Polonis, Gary A. Dreyer

This investigation involved a study of the reactions occurring during aging of supersaturated a phase in a CIL-Si alloy. The aging processes at temperatures below 552°C were studied by means of metallography, microhardness, and electrical resistivity measurements. During the initial stage of aging in the temperature range 275o to 525oC an anisotropic phase, which is visible in the form of striations in the micro-structure formed prior to the appearance of the -phase. This transition phase has been designated k' due to its structural similarity to the k phase. The rate of K' formation during aging increased with aging temperature. At 450o and 525°C the reaction appeared to reach completion before significant amounts of precipitation were observed. The overall precipitation reaction in this system involved four stages: formation of formation of at a grain boundaries, growth of particles in 7-egions, and coagulation of into large particles. The relative role of each of these processes is dependent on the temperature of aging. The formation of phase from supersaturated a during aging increases the electrical resistivity and slightly increases the hardness of a 4.9 pct silicon alloy. The formation of during aging decreases the electrical resistivity and increases the hardness until coalescence of particles begins during the final stages of aging. The coalescence of results in a decrease of both electvical resistivity and hardness. THE available literature does not record a detailed analysis of the aging processes in Cu-Si alloys. Such an analysis has been difficult due to both the limited data available and misleading microstructural observations which have been reported. Previous work suggests that a rather complicated series of processes may be involved in the aging of supersaturated a. The present investigation is concerned with an analysis of the sequence of these aging processes in a 4.9 pct Si alloy. The currently accepted version of the Cu-Si equilibrium diagram, Fig. 1, was published by C. S. Smith in 1940.' This diagram shows a eutectoid at 5.2 pct Si with the high temperature phase, having its composition range entirely to the right of this value. By quenching alloys in the composition range 4.6 to 5.2 pct Si from the region, Smith was able to retain the phase. However, this phase exhibited a striated microstructure which was anisotropic to polarized light. These striations were suggested by Smith to be duplex twins. Using an oscillating crystal X-ray method Barrett2 studied alloys containing 4.0 to 5.4 pct Si quenched from the and regions. The striations observed in quenched specimens were attributed by Barrett to quenching stresses. He found that striations could also be formed by deformation after quenching and identified them as stacking faults in alloys containing 4.0 to 5.0 pct Si. Barrett related the striations in alloys containing 5.0 to 5.4 pct Si both to stacking faults and to a metastable phase which he called The striated structure reported by smith' was also obtained by Hoffmann, et al. 3 By deeply etching their specimens the anisotropic phase on the surface was removed, leaving a structure which was microscopically similar to pure copper. By means of X-ray analysis before and after deformation at room temperature it was concluded that the anisotropic phase was metastable This conclusion has been questioned by Barrett who attributes the discrepancy to superposition of k and lines on the powder patterns of Hoffmann, et alm3 The decomposition of supersaturated according to smith' and Anderson,4 is very sluggish. However, Nony and Dreyer,5 and Hoffmann, et al.' found that deformation increases the rate of decomposition. X-ray studies by Hoffmann, et al on

Jan 1, 1962

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 R. E. Collins

The existence of fluid migration across fixed boundaries in oil and gas reservoirs has been known for many years. Several techniques have been developed in the past for estimating The rate of migration across fixed boundaries as an aid in planninng field development and in the valuation of oilfied properties. The principal de/errc7nt to the use of these techniques lies in the rather extensive reservoir and field data required for ca1cula1ion.s of fluid migration. For this reason, a new, simplified procedure has been developed which makes possible the calculation of fluid migration with a minimum of field and resrvoir data. This new technique is based on certain. solutions of the differential equations describing flour in the reservoir which assume, for the portion of the reservoir of interest, that the formation can be approximated by a homogeneous rock of uniforrn thickness, that only a Single mobile fluid phase exists, and {hut fluid production at the well is solely a result of expun.rioti of rile resersoir fluds. The results of the, present work are compiled in a set of curves. These curves can he used to calculate both rule and cummulative fluid migration when the abovc assumptions ore ju.r/ifiril. The only data required for such calcu1ations are the production hisrories of all wells in the field, the permeability and porosity of the reservoir., the, compressibility and viscosity of the fluid. and the shape of the reservoir, INTRODUCTION The migration of fluids across fixed boundaries in oil and gas reservoirs has long been recognized as an engineering and economic problem, but no methods for estimating the extent of migration during the development stage of a reservoir have been published. The investigation reported here has been directed toward providing a simple technique for estimating the extent of such migration and has therefore been restricted to the simple case of a homogeneous reservoir of uniform thickncss containing compressible liquids. It is assumed that only a single mobile fluid phase exists and that fluid production at the wells is solely by expansion of the reservoir fluids. Obviously, few, if any. reservoirs conform to these assumptions during their entire productive life. On the other hand, most reservoirs approximate fluid expansion reservoirs during their initial stages of primary production. Consequently, the present work should be applicable during the field development of most reservoirs, regardless of whether or not they ultimately are gas cap, dissolved-gas, or water-drive reservoirs. The mathematical analysis of this problem has been outlined in the Appendix, while the results of the investigation are presented in the following discussion. DISCUSSION OF CALCULATION PROCEDURE SLNGLE Well, Ini-iimite Reservoir The procedure for calculating fluid migrations is based on solutions of the partial differentia1 equations describing the flow of a homogeneous compressible liquid in a homogeneous reservoir of uniform thickness. This differential equation has the form of the well-known diffusion equation. The equation, as applicable to our problem, is where p is fluid density, is porosity, EL is viscosity, c is compressibility,* and k is permeability; .r, y and t are the space coordinates and time. Derivation of this equation is given in the Appendix. The analytical solution of this equation for a single well of negligible radius producing at a constant rate. q, from a reservoir of infinite areal extent for a time, 1. can he readily obtained. From this solution the rate of flow and also the cumulative flow across a straight line segment of length, 1, at a perpendicular distance. 6. from the well can be computed. One end of the line must be on the perpendicular joining the well to the

By R. G. Garlick, H. B. Probst

Tungsten single-crystal specimens of various orientations were deformed in tension at room temperature. Slip traces indicated both (112)(111) and (110) (111) slip; however, about 10 pct plastic dejormation was required bejore these traces could he seen. Resolzled shear-stress considerations indicated that at least some of the particular systems observed were not those on which slip initiated hut were secondary systems. This is highly probable, considering the high plastic strains involved. Analysis of the stress data with respect to possible slip systems indicated that slip must have been initiated at stresses below the proportional-limit stresses measured. ALTHOUGH there is general agreement that slip deformation in fcc and hep metals takes place on the plane of close packing and in the direction of close packing, there is no such agreement for slip in bcc metals. It is generally agreed that the close-packed direction (111) is the slip direction, but investigations have led to conflicting descriptions of the slip plane. As Hoke and Maddin1 point out, there are essentially four points of view. 1) Deformation in bcc metals is composite shear on (112) and (110) planes or on nonparallel (110) planes, and any apparent slip on other planes can be accounted for completely by slip on these planes. Chen and Maddin have reported this to be the possible case for molybdenum." 2) Slip occurs in the (111) direction and is limited to {110), {112), and (123) planes (the three closest packed planes of the lattice). 3) Slip occurs in a (111) direction, but not necessarily on planes of close packing (banal slip). 4) Slip occurs in those directions and on those planes for which ß, the ratio of Burgers vector to interplanar spacing, is a minimum. The slip systems of the bcc metal tungsten have never been fully determined, although several investigators have reported pertinent results.3"7 The present work was initiated to determine the active slip systems in tungsten at room temperature, the critical resolved shear stresses associ- ated with these systems, and the orientation dependency of operative systems. These features of tungsten deformation at room temperature have never been reported in the literature. The authors find this somewhat surprising, particularly in view of the availability of tungsten single crystals since the advent of electron-beam zone-melting techniques and the widespread interest in refractory metals. It would seem apparent that such work has been attempted but the results have not been reported, possibly due to the conflict with deformation theory. The results of the present work conflict with any straightforward description of deformation based on the proposed explanations of slip in bcc metals. Although these conflicts have not been resolved, a presentation of the problem seems in order and should be of value to others working in this field. MATERIALS AND PROCEDURE Single-crystal tungsten rods of 1/8-in. diameter were grown from undoped commercial rod by one-pass electron-beam zone melting in apparatus previously described by witzke.8 Representative analyses of the rod before and after melting are given in Table I. There was no detectable zoning of impurities during melting. Buttonhead tensile specimens, as shown in Fig. 1, were machined from single crystals. Two flats 90 deg apart were ground in the gage length of each specimen. These flat surfaces facilitated the analysis of slip traces. No attempt was made to position these flats with respect to the crystal orientation. The worked surface layer of the machined specimens was removed by electrolytic polishing in a NaOH solution. Tensile-axis orientations were determined by standard Laue back-reflection techniques. The sharpness of the Laue spots indicated that the crystals were free of strain.

Jan 1, 1964

By K. Takahashi, I. Iwasaki

The methyl stretching vibrations for a series of normal xanthates, isoxanthates and normal amines were determined with the aid of an infrared spectro-photometer, and the results were compared with published information on the corresponding alcohols. The frequency shifts of methyl stretching vibrations may be attributed to the inductive effect of polar groups. The hydrogens in the shorter chain alkyl groups are somewhat polar and can be expected to interact with water molecules by hydrogen bonding. Contact angle, surface tension lowering, solubility, and critical micelle concentration appear to be related to the induced polarity of the alkyl groups. Traube's rule states that in homologous series of compounds solubility decreases and surface activity increases regularly with the size of the hydrocarbon group. Analyzing the work of transfer of a solute from solution to surface, Langmuir' interpreted this rule to mean that each CH2 group added to a molecule produced its effect independently of the others. It is also well known2 that the maximum equilibrium contact angle with sulfhydryl collectors is independent of the polar groups and of the sulfide minerals, and that it increases to a certain limit with the number of carbon atoms in the alkyl group. This observation is thought to demonstrate the orientation of the collector molecules at the sulfide surfaces. The dependence on the hydrocarbon chain length is usually explained as due to the screening of the force field of polar mineral surfaces.3-5 Here again, Langmuir's principle of independent surface action appears to be implied. Heteropolar compounds consist of nonpolar and polar groups, the former being the hydrocarbon group and the latter containing such electronegative atoms as O, N, S, P, or halogens. According to Palit6 the hydrogens in the alkyl groups are more electroposi- tive than in the paraffins owing to the negative inductive effect of the polar groups and a possibility that a hydrogen bond with water molecules ensues. Goto and Takenaka7 provided direct evidence of the inductive effect by measuring the shift of the C-H stretching vibrations of alcohols, acids, and alkyl bromides with an infrared spectrophotometer. Their approach appeared to be directly applicable to the studies of surface activity both at air-water and water-mineral interfaces encountered in flotation systems. The present article attempts to relate the C-H stretching vibrations of alkyl xanthates, alcohols, and amines with such flotation chemistry parameters as contact angle, surface tension lowering, solubility, and critical micelle concentration. MATERIALS AND METHODS Normal xanthates from methyl to octyl inclusive and the isopropyl, isobutyl, and isoamyl xanthates used in this investigation were prepared in the usual manner.8 Normal amine hydrochlorides with carbon chain lengths of 1, 2, 4, and 6 were from Distillation Products Industries, of 8 from Archer-Daniels-Midland Co., and of 10, 12, 14 and 16 from Armour and Co. Normal propylamines, pentylamines and heptylamines from Distillation Products Industries were used both as received and as hydrochloride formed by converting the amines according to the procedure described by Li9 Perkin-Elmer Model 521 and 621 infrared spectro-photometers were used to investigate the spectra in the frequency region of 3050 to 2800 cm-' with frequency scale expansion. The accuracy of the frequency readings was within + 1 cm-1. The KBr micropellet technique was used to prepare the xan-thate samples for measurements, and the halocarbon mull technique was used for amines and amine hydrochlorides. EXPERIMENTAL RESULTS Fig. 1 shows the CH3 asymmetric stretching vibrations as a function of the number of carbons in the alkyl group for normal xanthates, isoxanthates and amines. Also included in the figure are data for normal alcohols and isoalcohols taken from the article

Jan 1, 1970

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 H. J. Ramey

A general solution is presented for the transient temperature distribution caused by radial movement of a cylindrical heat source through a homogeneous medium of infinite extent. This problem represents a highly simplified model of the movement of a combustion front during the thermal recovery of oil. Numerical solutions are presented for a heat source moving at a constant velocity, and a velocity inversely proportional to the radial location of the heat source. Numerical solutions are presented for both finite and infinite vertical thickness of the heat source, i.e., with and without vertical heat losses. The numerical solutions are used to estimate the fuel concentration needed to maintain a combustion front during the thermal recovery of oil. The sensible heat carried to the combustion front by the injected gas stream is discussed in regard to the fuel concentration. Numerical solutions for a heat source (combustion front) of finite vertical thickness indicate the field conditions that may be necessary to sustain the temperature of a combustion front above the ignition temperature of the fuel over considerable distances from an injection well. The results also have implications in regard to the quantity of heat required during ignition. The results of this computation apply to a highly idealized model of a thermal recovery process. But the results may be used as a guide in engineering consideration of the thermal recovery process. INTRODUCTIO N The thermal recovery of crude oil has received considerable attention since the publication by Kuhn and Kochl of laboratory and field tests of this method of oil recovery. Other publications2,3,4,5,6, since that time have dealt with various features of the thermal recovery of oil and have illustrated the extremely complex nature of this process. Vogel and Krueger7 described an electric analog computer designed to solve a moving heat source problem which was analogous to movement of a combustion front during the thermal recovery of oil. The heat source was maintained at constant temperature and was a cylindrical source of infinite height (no vertical heat loss) moving radially at either constant velocity, or at a velocity inversely proportional to the radial location of the heat source. Jenkins and Ramey2 presented analytical solutions to a similar heat conduction problem and pointed out the possible importance of vertical heat losses. The following presents a general solution to the transient heat conduction problem introduced by Jenkins and Ramey4 as well as numerical solutions for a wide range of conditions possibly similar to those that might exist during field operation of the analogous thermal recovery process. DESCRIPTION OF PROBLEM The problem considered in this paper is determination of transient temperature distribution caused by a cylindrical heat source of infinite or finite vertical height, moving radially through an isotropic medium of infinite extent. It is further assumed that the heat flux generated at the surface of the moving heat source may be a function of time such that the fuel concentration required may be constant, or permitted to vary, and that heat is generated at the surface of the heat source only. In regard to the thermal recovery of oil the assumption that heat is generated at the surface of the front is equivalent to assuming that the thickness of the combustion zone is infinitely small, or that the reaction rate between fuel and oxidant gas is infinite. The last assumption does not appear stringent. See laboratory information published by Martin, et al3 Under these assumptions the temperature distribution throughout an infinite medium caused by a heat source moving radially with an arbitrary velocity may be described by The heat generation function or heat flux at the surface of the source, +(t), is Note that the fuel concentration may be defined as a function of time. Eq. 1 is derived in the Appendix.

By Charles Stein

The effects of several different prior treatments an the creep behavior of 99.9995 pct aluminum at 260°C and 1000 psi canstant stress are compared with annealed specimens. Quenching from 538oC, irradiation with 2 mev electrons, and tension-compression room temperature fatigue damage were used to change the substructure of the specimens prior to their creep testing. The quenched and the irradiated specimens showed a larger primary and transient stage contribution to the creep curve than those specimens which were only annealed prim to creep testing. The specimen receiving prior fatigue damage at room temperature showed no first stage or transient creep when tested under identical conditions as the above specimens and had an average creep rate seven orders of magnitude lower than that of either the annealed, the quenched or the irradiated creep specimens. In a previous electron transmission microscopy investigationl of the substructure developed during the creep of pure aluminum at 0.57 Tm, it was noted that a large concentration of vacancy loops and dislocation loops were present in and near the subboundaries while the interior of the subgrains contained few of these defects and had a low dislocation density, see Fig. 1. Yim and Grant2 and Hazlett3 have shown that the presence in nickel of a substructure developed by prior cold work effectively reduces the primary and transient contributions to the creep curve. Hultgren,4 McLean,5 Gervais, Norton, and Grant,6 Chang and Grant,7 and others have shown the same effect with the higher stacking fault energy material, aluminum. However, the specific defect responsible for the change in the creep behavior developed by the prior cold work was not established. Specifically, what effectively interfered with dislocation motion in these materials at elevated temperatures? Was it jogs on dislocations produced by their interaction with other dislocations or with vacancy or dipole loops, see Fig. 1, or did the subgrain walls determine the mean free path for glissile dislocations? This would be a realistic possibility only if subgrain boundaries are effective barriers to dislocation motion. The ability of a subgrain boundary to act effectively against glissile dislocations depends on the number and the arrangement of dislocations in the subboundary,8 which in turn is a function of the creep strain.' This paper compares the creep rate of specimens possessing a large vacancy concentration to that of annealed specimens and with specimens having a stable subboundary wall containing Lomer locks produced by prior fatigue damage. EXPERIMENTAL PROCEDURE Vacancy Loops. Aluminum creep specimens, 99.9995 pct pure, were machined from zone refined rod to a 2-in. gage length and a cross-sectional area of 0.025 in. These specimens were annealed at 538oC, 1000°F, for 10 min and quenched into ice water. They were then reheated to the creep test temperature of 260°C, 500°F, and held at this temperature for 18 min to al-low for vacancy condensation and loop formation. The specimens were subsequently creep tested at 1000 psi constant stress at 260°C. The temperature along the gage length was monitored by three Pt-Pt, 13 pct Rh thermocouples, two embedded in the upper and lower shoulders of the specimen and one attached to the mid-point of the gage length- The temperature gradient was held to ± l°C or less throughout the Creep test. Elongations were measured by two LVDT's 'connected in an averaging 'On-figuration, with quartz extension arms placed against the specimen at the extremes of the 2-in. gage length. i K 3J?--' * *H - >^> - ' ^L r 77 . *y J OJL, J Fig. 1—Transmission electron micrograph of pure aluminum strained 8.2 pct at 260°C and 1000 psi constant stress. Note the presence of numerous loops near the light dislocation tangle

Jan 1, 1970

By M. R. Tek, F. F. Craig, J. O. Wilkes, C. S. Goddin

The waterflood performance of a water-wet, stratified system with crossflow is computed by a finite difference procedure. The effects of five dimensionless parameters on tile oil displacement efficiency, water saturation con-tour.7 and crossflow rates are evaluated in the absence of gravity forces. Crossflow due to viscous and capillary forces is shown to exert a significant effect on oil recovery in a field-scale model of a two-layered. water-wet sandstone reservoir. The crossflow is at a maximum in the vicinity of the front advancing in the more permeable layer. Under favorable mobility ratio conditions, the comparted oil recovery with crossflow always is interrnerliate between that predicted for a uniform reservoir and that for a layered reservoir with no crossflow. INTRODUCTION The important erects of reservoir heterogeneity on waterflood performance are commanding increased attention in the technical literature. Much of this attention is centered on two categories of layered reservoirs: those in which layers are non-communicating and those in which crossflow of fluids occurs between the layers. In the first category, the reservoir is assumed to consist of discrete layers, each uniform within itself and differing from the others only in such properties as thickness, porosity and absolute permeability. The performance within each layer is calculated by one-dimensional flow theory, and the performance of the total reservoir is obtained by summing individual layer performances. Capillary and gravity effects usually are not considered. Representative publications dealing with thi5 type of reservoir are those of Stiles,' Dykstra and Parsons,' Hiatt,3 Warren and Cosgrove' and Higgins and Leighton." Prediction of performance for reservoirs in the second category is considerably more difficult since viscous, capillary and gravitational forces all play important roles in causing crossflow between layers. A number of authors have investigated the simpler problem of two-dimensional displacement flow in a stratified system with a mobility ratio of unity and negligible capillary and gravity effects.'; Others have considered two-dimensional, non-steady-state flow of a single, slightly compressible fluid in a stratified reservoir. A limited number of laboratory oil displacement tests in layered models with crossflow have been reported. Miscible floods (with resultant zero capillary forces) in layered five-spot models were conducted by Dyes and Braun," who studied the effect of mobility ratio with zero gravity forces, and by Craig et al. 12 who studied the effect of gravity forces at constant mobility ratio. Waterfloods in layered five-spot models (with cross-tlow due to capillary, viscous and gravity forces) were conducted by Gaucher and Lindley,"' who showed the effect of gravity forces in causing underrunning of the injected water and by Carpenter, Bail and Bobek, 14 who demonstrated the reliability of Rapoport's" dimension-less parameters for scaling layered systems. Waterfloods in rectangular layered models were conducted by Richardson and Perkins."' who investigated the effect of velocity at constant mobility ratio and with zero gravity forces, and by Hutchinson," who studied the effects of varying mobility, layer permeability and layer thickness ratios. The differential equations which rigorously describe waterflooding in a heterogenous porous medium are non-linear and do not facilitate analytical solution. By using finite difference approximations it is possible to obtain a solution to any desired degree of accuracy. Such a solution, using an alternating direction implicit procedure (ADJP), is described by Douglas, Peaceman and Rachford.18 In the present study, a computer program using ADJP explores systematically the effects of important parameters on waterflood performance of a two-dimensional, two-layered, field-scale model of a water-wet sandstone system. Particular attention is given to evaluation of the water saturation contours and crossflow rates at the interface between layers to gain improved understanding of the crossflow mechanism. PROCEDURE BASIC: FLOW EQUATIONS The basic flow equations for two-dimensional, two-phase, immiscible, incon~pressible flow in a porous medium are:

Jan 1, 1967

By E. L. Claridge, K. S. Lee

Areal sweep efficiency of oil displacement by enhanced-viscosity water exhibiting pseudoplastic behavior was measured in a Hele-Shaw model representing one-quarter of a five-spot pattern. The pseudoplasticity of polymer solutions and the velocity distribution in the five-spot pattern produced a condition under which the mobility ratio between the displacing and the displaced fluid could not be assigned a single value. Instead, the movement of the displacement front is governed by local mobility ratios which are also time dependent. The areal sweep at breakthrough with polymer solutions was poorer than the sweep obtained with Newtonian fluids of comparable viscosity. However, the areal sweep at 1 PV throughput was greatly improved as compared to flood water without polymer. It was also demonstrated that, even after the oil-cut had declined to a low value during a regular waterflood, switching to polymer flood efficiently swept out the oil remaining in the model. INTRODUCTION The behavior of fluid displacements in isotropic porous media for various patterns of injection and production wells has been extensively investigated.1-7 These investigations all concerned Newtonian fluids, i.e., the viscosity of each fluid was constant regardless of flow rate. The generally unfavorable influence on areal sweep efficiency of higher mobility of the displacing fluid as compared to the mobility of the displaced fluid has been established for both miscible and immiscible fluids. The principle was also established that a close correspondence exists between miscible and immiscible flood front behavior,2, 3 although oil recovery in a waterflood at unfavorable mobility ratio may be less than that observed in a miscible displacement at the same mobility ratio. This is true even when oil recovery is expressed on the basis of movable oil. The reason is that oil saturation only slowly achieves its final value behind the water-flood front in accordance with the Buckley-Leverett8 simultaneous flow relations. It is convenient to use miscible displacements for laboratory simulation of waterflood frontal advance since the interfacial tension forces which are negligible in proportion to viscous forces on a reservoir scale9 are thus made nonoperative in the laboratory model. For miscible displacements, the Hele-Shaw type of model adequately represents a porous medium so long as the appropriate scaling rules5-10 are observed in its design and operation. During simulation of waterflood front behavior in the laboratory by using miscible displacements, the behavior of connate water may ordinarily be disregarded since it is usually indistinguishable from flood water in this process. However, when the flood water is deliberately thickened to improve the mobility ratio between water and oil. the effect on the sweep efficiency due to generation of a connate water bank during the process must be considered. In a uniform porous medium, such a bank is generated11 and efficiently displaced by injection of thickened water.l2 The oil originally in-place at the start of the waterflood is then displaced by connate water followed by thickened water. If the flood water must be thickened to obtain a favorable mobility ratio, the mobility of the oil phase is appreciably less than that of the connate water. Hence, the oil phase is inefficiently displaced by the connate water bank, and a considerable proportion of the oil comes in contact with and is displaced by the thickened waterflood front' A recent study 12 of five-spot displacements using miscible Newtonian fluids in which connate water banks of various sizes were simulated shows that, for connate water banks as large as 40 percent of 1 PV, major improvements in sweep efficiency resulted from use of a thickened

Jan 1, 1969

By C. B. Manula, H. Gezik

In planning modem-day mining operations, management needs to pass from the area of subjective decision-making to an area of objective decision-makirlg. Planning procedures currently being practiced by many mine managers depend primarily on common sense and trial-and-error methods which oftentimes never reach optimum solutions. For this reason, an attempt is made to provide an understanding of linear programming to solve production planning problems. More specifically, a linear programming model is proposed to solve the production scheduling problem as it was found to exist in the crushed stone industry. Optimal production scheduling is of prime concern here, since the industry has to cope with seasonal demands for its products. The solution, as determined through the application of the model, provides management not only with an optimal production schedule but also with the following information regarding the improvement of existing facilities: (I) the location of bottlenecks which limit production and profit; (2) the directions in which the company i fixed equipment should be expanded or modified; and (3) the evaluation of unused resources. Future trends in stone production will quite likely parallel those in the construction industry. Continued record expenditures for structures such as dams, highways, and buildings are the reasons for this accelerated growth. Production forecasts for the next decade are predicted to be around one billion tons with an approximate 7% annual increase. The problem of increasing production from 570 million tons in 1966 to 1 billion tons by 1976 offers no small challenge. This situation is further compounded by other problems and economic factors. Population expansion, government regulation, labor agreement, and changes in technology wd cause managers and operators to be confronted by many more complexities than characterize present-day operations. Some of these. problems and economic factors are: (1) In some metropolitan and other urban areas, easily accessible sources of stone are almost depleted. Transportation costs will be a prime concern of the industry. (2) Population growth has forced residential expansion to a point where many quarry operators must conform to urban and suburban regulations. (3) Government regulation in regard to air and water pollution, noise, and reclamation as well as zoning developments necessitate that crushed stone producers keep abreast of changing conditions and requirements. (4) The crushed stone industry operates under pressure of a highly competitive situation. Every company's objective is to retain a favorable marketing position by exercising great attention to production and marketing costs. This is rather difficult, since the crushed stone industry is one segment of modern mining which has to cope with a seasonal demand for its products. Among the problems posed here, the one of exercising great attention to production and marketing costs is of current interest to management. Under seasonal demand conditions, this problem creates further difficulties for managers who need to plan and forecast. STATEMENT OF THE PROBLEM Various methods of Operations Research have been developed to handle a problem whose genesis is explained above. It is very seldom, however, that these quantitative means are used in practice by the mining industry. This paper describes the use of Operations Research methods to solve a production planning problem when demand is seasonal. The seasonal production planning problem is an important one because of the substantial costs associated with changing production, employment, and inventory levels to adjust to seasonal sales patterns. Most companies try to operate within reasonable levels of employment changes and hold inventories within storage capacities. In general, little or no attempts were made to develop minimum cosl schedules or to provide objective means to compare various production policies. When demand varies within a fixed time period, 11:-oduction requirements may be met by following one

Jan 1, 1970

By John C. Sawyer

A brief review is given to show that catastrophic oxidation of stainless steel in the presence of lead oxide is similar in many respects to catastrophic oxidation as caused by MoO3 and V2O5, Experimental data are presented to show that lead oxide accelerates the normal oxidation process of stainless steel by chemically altering the normally protective oxide film formed on the base metal. Data are given to show that similar reactions occur during oxidation of stainless steel in the presence of MOO3,and V2O5. THE mechanism by which pure metals oxidize at high temperatures has been the subject of intensive research for many years and much has been accomplished towards understanding the processes involved. The mechanisms by which alloys oxidize have also been investigated, but due partly to the fewer number of investigations and partly to the complexity of the oxide systems present, available knowledge is not as extensive as with pure metals. When oxidation processes are further complicated by atmospheres contaminated with sulfide, halide, or oxide vapors, the knowledge of fundamentals becomes meager or altogether lacking. Thus, it is not surprising to find that the catastrophic oxidation of stainless steels in metal oxide vapors of MOO, and V2O5 is not a well understood phenomenon despite the fact that a number of investigations have been undertaken. The literature in dealing with the subject of catastrophic oxidation is primarily concerned with the effects of MOO3 and V2O5, but does mention that other metal oxide vapors, one of which is lead oxide, appear to promote a similar type of accelerated oxidation. While this latter effect is virtually unknown, the deleterious effect of lead oxide is of considerable commercial importance, as it is believed to be one of the fundamental reasons for exhaust-valve failure in gasoline engines. For this reason studies have been in progress for many years attempting to better understand the mechanism involved. Modern exhaust-valve alloys"2containing 20 to 25 pct Cr plus nickel and/or manganese are in effect high-grade austenitic stainless steels exhibiting ex- cellent oxidation resistance in air at exhaust valve operating temperatures of 1100° to 1700°F; yet these same alloys are prone to corrode during engine service because of the conditions attending their operation. These conditions are such that the valves become coated with deposits of complex lead compounds mainly derived from gasoline additives of tetraethyllead, ethylene dibromide, and ethylene dichloride. Literature searches3-7 and discussions with various investigators show that considerable difference of opinion exists as to the nature of the lead-bearing deposits found on the surface of exhaust valves and the atmospheres to which they are subjected; yet, all seem to generally agree that corrosion, a gradual transformation of alloy to nonmetallic products of corrosion, is a consequence of the presence of lead compounds. This majority opinion seems to be substantiated by laboratory engine tests which show that corrosion does not occur to any appreciable extent when lead-free fuels are used. The first observations of catastrophic oxidation reported in the literature were probably made by pfeil,8 but it was Leslie and Fontana9 who undertook the first detailed study of the phenomenon. Their work was chiefly centered on the acceleration of oxidation which took place when stainless steels containing more than 3 pct Mo were heated in air with restricted circulation. These investigators reported that the presence of other oxides near the surface would further accelerate the rapid oxidation which occurred. While they recognized that the effect was attributable to the presence of molybdenum, it was Rathenau and Meijering" who emphasized the fact that the presence of molybdenum in the alloy was incidental to the presence of MOO3 vapor in the atmosphere surrounding the specimens. Evans,11 as well as Cunningham and Brasunas 12 studied similar

Jan 1, 1962

By C. H. Li, R. J. Stokes, T. L. Johnston

The phenomenon of solid-solution strengthening has previously been studied in a number of binary alloy systems.1-8 There has been, however, very little information published concerning the strengthening effects of specific solutes in magnesium alloys. Schmid and his coworkers9,10 investigated the behavior of aluminum and zinc in binary and ternary magnesium-alloy single crystals, but this work was done before high-purity materials were readily available. Workers at the Dow Chemical companyH reported increases in the critical resolved shear stress upon alloying magnesium single crystals with zinc and indium. Recently, Hardie and parkins* investigated the extent of hardening produced by the addition of various solutes to poly crystalline magnesium. The present work was undertaken to provide data on the strengthening effects of a number of solutes in high-purity magnesium, for the purpose of supplying further experimental information which might be useful in describing a general mechanism for solid-solution strengthening. To simplify interpretation of the results, the mode of deformation was limited to basal slip, through the use of single crystals. EXPERIMENTAL PROCEDURE Magnesium and magnesium binary-alloy single crystals containing indium, cadmium, thallium, aluminum, and zinc were prepared in the form of l/2-in.-diam rods, approximately 8 in. long, using a modified Bridgman technique in which the crystals were slowly cooled from the melt in a steep temperature gradient, in an atmosphere of helium. It was found that a solid-liquid interface travel speed of about 0.7 in. per hr consistently produced single crystals. The poly crystalline starting materials for the preparation of single crystals consisted of 1/2-in.-diam extruded rods, prealloyed to the desired composition. The unalloyed magnesium and magnesium-indium extrusions were furnished by the Dow Chemical Co. The other alloys were prepared at the Rensselaer Polytechnic Institute Metallurgical Laboratories. Some preliminary studies of the strengthening effect of indium were conducted at room temperature using a constant rate of tensile loading of 1040 g per min. Stress was applied to the crystals by means of a simple 5:l lever arm. The lever was loaded by a constant rate of flow of shot into a receptacle suspended from the long arm of the lever. Plastic flow was detected by an electrical resistance strain gage cemented to the surface of the crystal normal to the minor axis of the predicted glide ellipse. The crystals were acid machined to produce tensile specimens with a reduced section between the grip ends, with an additional step in the center of the reduced section to compensate approximately for the slip-retarding effect of the strain gage. The balance of the single crystals were extended in an Instron tensile tester. This machine employs a synchronously driven screw and automatically records applied load as a function of cross-head separation. It became unnecessary therefore, to employ a strain gage to detect deformation or to acid machine a reduced section. Preparation for testing reduced merely to chemical cleaning of the surface, thus minimizing the amount of handling and probability of accidental damage to the crystals. The rate of elongation employed in all the tests was 0.002 in. per min, corresponding to a strain rate of approximately 0.0004 per min. EXPERIMENTAL RESULTS The investigation confirmed the finding of many workers that the yield strength of hexagonal metal crystals is strongly orientation-dependent. This is illustrated in Fig. 1, in which o, the tensile stress required to produce gross plastic flow in eleven unalloyed magnesium crystals stressed at a constant rate of loading is plotted as a function of sin X, cos Ao where X, is the angle between the

Jan 1, 1960

By J. Szekely, V. Stanek

Steady state heat transfer and fluid flow rates are predicted for the counter-current contacting of gas and solid streams. For a .fixed pressure drop across the bed the convective heat transfer rate is affected by the gas flow rate, which in turn is affected by the temperature distribution in the system. Variations in bed voidage and in particle size were found to change the predicted temperature profiles quite markedly, in particular when the thermal capacities of the gas and solid streams were equal. IN Part I we examined the effect of step and pulse changes in the inlet gas temperature on transients in the flow of gases through fixed packed beds. Here we shall discuss heat flow and fluid flow interactions in the counter-current heat exchange between packed solids and a gas stream. Counter-current heat exchange in packed beds has been extensively studied. In his now classical paper, Furnas' assumed constant property values and developed a simple analytical. expression for the temperature profiles in the gas and in the solid stream. A highly sophisticated analysis of counterflow systems, allowing for transients and interparticle conduction was developed by Amundsen. In preparing a very readable review of this topic Elliott' also presented a number of elegant solutions, including the case involving heat generation. A great deal of work has been done on the application of these concepts to the modelling of the blast furnace. However, in all this previous work the mass flow of the gas was regarded as a quantity independent of the gas-solid heat exchange process. In the following we shall consider counter-current heat exchange between a gas stream and a packed bed of solids, under conditions where the pressure drop across the bed is specified, rather than the mass flow of the gas. FORMULATION Consider a packed bed, extending from y' = 0 to y' = L, in which moving solids and a moving gas are contacted in a counterflow arrangement. Let cold solids, with an initial temperature Tsi and mass flow GA, enter at y' = L, and hot gas with an initial temperature of Tgi and mass flow Gg enter at y' = 0. In contrast to the problems considered in Part I this system has a nontrivial steady state and we shall restrict our attention to the study of the pressure drop flow relationships and the resultant temperature distributions for steady state conditions. Under these conditions the formulation of the problem requires the statement of the following: 1) heat balance on the gas 2) heat balance on the solids 3) the mechanical energy balance equation and in addition the knowledge of the equation of state, and a functional relationship between viscosity and temperature and an expression for the heat transfer coefficient in terms of the other parameters. On neglecting temperature gradients within the solids and for a constant mass flow of the gas, the solution to this problem is well known.'" The expressions for the temperature profiles of the gas and solid streams may be put in the following form: where Tg and T, denote the dimensionless gas and solid temperatures respectively, y is a dimensionless distance, and R is the ratio of the thermal masses. These quantities are defined as follows: For the particular case of equal thermal mass capacities by the gas and the solid, that is R = 1, the Eqs. [1]and[2] simplify to Eqs. (71 and[8], respectively: i.e., the resultant profiles are linear. As for most practical cases , solutions given by Eqs. [7] and [8] will contain only the inlet temperatures as parameters. For a fixed, predetermined value of C;b the above

Jan 1, 1970

By A. N. Crownover, H. T. Kennedy, E. P. Miesch, C. H. Bowman

The paper presents correlations of (l) molar volume of gaseous hydrocarbon mixtures with pressure, temperature, composition and properties of the C,-plus fraction; (2) shrinkage of oils during flash and differential liberation of gas, including the calculation of formation volume factor under various conditions; and (3) bubble-point pressure with temperature, composition and characteristics of C, plus. The data on which the correlations are based comprise 1,615 measurements on 900 hydrocarbon systems, including numerous systems containing nitrogen, hydrogen sul-fide and carbon dioxide. In each correlation, the number of data points covered and the accuracy is substantially greater than in previously available work. Thus, the equation yielding molar volumes of gases has an average deviation of 2.04 per cent, applied to mixtures having temperatures up to 313F and pressures up to 9,800 psis, compared to 2.37 per cent for the Benedict-Webb-Rubin equation applied to the same data, and 4.53 per cent for the method based on the law of corresponding states. he equations resented are all explicit in the dependent variable, and require no iteration on the digital cow puter. INTRODUCTION The ease and accuracy of determining the composition of hydrocarbon mixtures, compared to the difficulty of measuring their properties under reservoir conditions, makes it desirable to utilize composition as the key to physical behavior to the greatest possible extent. As a result, there are available correlations between composition, or easily measured characteristics dependent on composition, and practically every important engineering property of reservoir fluids. The task confronting us is one of finding more exact relationships between important variables rather than extending correlations to new properties. This paper describes new correlations of molar volumes of gases, formation volume factors, and bubble-point pressures with composition, temperature and where possible, pressure. Each correlation is obtained by employing a sufficiently large amount of data so the calculated properties are probably as least as accurate as the measurement on which they are based. MOLAR VOLUME OF GASES Although many equations of state have been proposed for pure gases, only a few methods are applicable to hydrocarbons at conditions comparable to those in petroleum reservoirs. Still fewer arc useful in describing the behavior of mixtures, with which the petroleum engineer is largely concerned. The correlation presented here involves procedures similar to those of Alani and Kennedy The van der Waals equation is modified to make a and b functions of temperature instead of being constants for each material.This is a cubic equation, which may have either one or three real roots. The lowest root corresponds to liquid volume, while the highest applies to gas. When the procedure was applied to 703 pressure-temperature points of 164 gases of known composition and volume, the average deviation was 12.08 per cent, and the standard deviation 8.15. The above calculations were made using the constants derived for liquids by Alani and Kennedy and the mixture equation developed by them. A closer approximation to measured molar volumes is obtained by employing different sets of constants for different areas on the pressure-temperature chart, and by changing the relationship between a,,, and b,,, for mixtures and the a, and b, for individual hydrocarbons. The various sets are shown in Table 1, and the areas for which they are recommended are plotted in Figs. 1 through 7. In these figures, the area designated by zero is in the critical region for pure materials, and the values obtained for them may be unreliable. For normal mixtures, however, the area may be added to any of the bounding areas without substantial error. The new relationships employed are

Jan 1, 1966

By A. W. Schlechten, C. M. Hsiao

The apparent vapor pressures of a number of metal sulphides were determined by measuring their rate of weight loss when they were heated under vacuum. The calculated pressures are due in some instances to the volatility of the sulphide and in others to the decomposition pressure. FIE direct reduction of metallic sulphides, as in the precipitation method for lead or antimony, is not a widely used process today. However, it might be employed for the production of a number of metals if the operation were carried out under reduced pressures. Two general procedures are possible; the reducing agent could form a stable nonvolatile sulphide from which the valuable metal could be removed by volatilization or liquation. The other approach would be to choose a reducing agent which would form a sulphide that could be volatilized away from the valuable metal. A vacuum operation would serve to protect the sulphides and metals from oxidation and would increase the evaporation rates. A related process is the separation of metallic sulphides by selective volatilization either at atmospheric pressure or under vacuum. This procedure is already used commercially for the elimination of lead and cadmium as PbS and CdS from zinc calcines; it might well have other applications. To predict such processes it is necessary to know the relative volatilities of the metallic sulphides and to know whether they volatilize as sulphides or whether the predominant weight loss is due to decomposition. Such data would also be of value in studying problems such as the use of sulphides for refractories, as described by Brewer,] or in the evaluation of the theory of ore deposition from the vapor state, as proposed by Brown.' Qualitative statements concerning the volatility of metallic sulphides can be found frequently in the literature and are well summarized by Mellor, but accurate determinations of rates of evaporation or vapor pressures are very scarce. Kelly" lists vapor pressure data only for PbS; his reference being to the work of Schenck and Albers,4 who measured the vapor pressure of PbS by a dynamic method between 850" and 995°C. More recently Veselovskii5 has determined the vapor pressures of the sulphides of antimony, lead, cadmium, and zinc at temperatures ranging from 360" to 1000°C using Knudsen's method. A second Russian worker, Pogorelyi,6 determined the vapor pressures of CdS and ZnS at temperatures from 900" to 1250°C. The results obtained by these workers are shown in subsequent graphs in comparison with the data of this paper. Brewer and coworkers1 have studied some of the metallic sulphides to determine their suitability as refractories in vacuum. Veselovskii's determinations on PbS checked well with those of Schenck and Albers, and his determinations of CdS check almost exactly with Pogo-relyi. For the vapor pressure of ZnS the data of Veselovskii and that of Pogorelyi are in fair agreement, as shown in Fig. 1, but it should be pointed out that Lumsden7 has discovered that the equation given by Veselovskii for the vapor pressure of ZnS: - 14,200 log P9.495 is in error by a factor of ten when compared to the experimental data. In contrast with the scarcity of data on vapor pressures, the decomposition pressures of many metal sulphides have been determined and are summarized by Smithells.8 Theoretical Considerations The experiments described in this paper consisted of determinations of the rate of weight loss of certain metallic sulphides under very low pressures at different temperatures. These data plus observation and analyses of the products reveal the mechanism of the weight loss and also provide a basis for the calculation of the vapor pressure or apparent vapor pressures of the sulphides.

Jan 1, 1953

By S. S. Marsden, S. H. Raza

An experimental study of the flow of line-textured. aqueous foams through Pyrex tubes is described. The foams range in quality F (ratio of gas volume to total volume) from 0.70 to 0.96 and behave like pseudoplastic fluids. At lower flow rates they exhibit laminar flow ad have apparent viscosities which increase with quality, and which cover a range of 15 cp to 255 poise for tubes of 0.25- to 1.50-mm radius rj. At higher flow rates a plug-like type of flow is developed, the extent of which increases with both T and ri. When the same foams flow through either open or packed Pyrex tubes, remarkably high streaming potentials ?E are often generated. These can easily reach 50v if nonionic foaming agents are used, but are at least an order of magnitude less for ionic foaming agents. A linear relationship between ?E and the pressure differential ?p is observed; this usually extrapolates to positive values of ?p at ?E of zero. The slope of the line increases with both F and ri. An equation was derived to describe the streaming potential of non- Newtonian fluids in circular tubes and was used to correlate experimental results. The calculated < potentials are of the right order of magnitude. INTRODUCTION Foams are both unusual and inniguing in their physical properties, and have been the subject of many scientific studies. However, present knowledge of foams is still fragmentary, specific and often contradictory. Apparent viscosity of foam is the physical property of greatest interest to both rheo1ogists and engineers. sibreel reported that the apparent viscosity decreased with increasing shear rate in a manner similar to some non-Newtonian fluids. Penny and Blackman2 reported that fire-fighting foams had both a limiting shear stress and a tensile yield stress. There is little doubt that some foams at least behave like non-Newtonian fluids, and have apparent viscosities considerably higher than those of either constituent phase. The high apparent viscosity of foam with its concomittant effect on mobility ratio and sweep efficiency no doubt prompted several attempts by research groups to use foam as a displacing agent in porous media. Based on recent experience, most of these groups probably succeeded in completely blocking fluid flow in the porous media and then abandoned their efforts. Two groups apparently found the successful combination of experimental parameters at about the same time.3,7 Others have recently added to our knowledge of foam flow in porous media and its use as a displacing agent.8-13 An experimental problem encountered by Fried3 was a transient blockage of foam flow in porous media when distilled water was used to prepare the foam-producing solution. Fried surmised that this was due to an electrokinetic effect and he eliminated it by using electrolytes in preparing foaming solutions. He also measured the streaming potential of a number of foams in capillary tubes which he found to be appreciably higher than those obtained when the constituent liquid flowed under comparable conditions. This paper presents results of a more comprehensive study of the streaming potential generated by aqueous foam flowing in both open and packed Pyrex tubes. It also adds to knowledge of the rheology of these foams as deduced from their flow behavior in open tubes. APPARATUS AND PROCEDURE A diagram of the apparatus used is shown in Fig. 1. Details of its construction, testing and use are described elsewhere.14 Careful selection of materials, extreme cleanliness and rather elaborate electrical insulation and shielding were necessary to obtain reproducible results (±5 percent). Both streaming potential and streaming current were measured with an electrometer. The design of the foam generator developed for this work is novel (Fig. 2). It will produce uniform fine-textured foams at pressures sufficiently

By G. R. Speich, D. J. Mack

The transformation of was studied by isothermal methods. At all temperatures, the ß transforms quickly to fine grained ß" which develops silver-rich striations. At higher temperatures the striations disappear, the final structure being Widmanstatten a in ß&apos;. At lower temperatures, the striated ß" is consumed by pearlite nodules of a + ß&apos; which in turn form the Widmanstatten a in ß&apos;. The transformation is unlike any previously reported for a eutectoid. AS part of a general program of study on the eutectoid reaction, the Ag-Cd eutectoid seemed particularly attractive because the ß phase which undergoes the reaction is another of the typical "electron compounds," AgCd, having the disordered body-centered cubic structure. This work deals with the eutectoid which occurs at 50.5 wt pct Ag, and not with the eutectoid which has been reported at 42.9 pct Ag.&apos; There has been no work reported specifically on the Ag-Cd eutectoid which was studied. The only information available on this eutectoid is that work done in the determination of the phase diagram. For the most part, the terminal portions of the diagram, Fig. 1, are well established; however, the region from 40 to 60 pct Cd is still somewhat in doubt. The X-ray evidence&apos;..&apos; in this region is in conflict with the results of thermal analysis. X-ray analysis3 does not show a polymorphic transformation at about 240°C as indicated by thermal analysis.1,5 Even the X-ray evidence, within itself, is confusing in this region of the diagram. One investigator2 reports the ß phase as hexagonal close-packed and has found three phases in the one-phase ß field. Other X-ray investigators"&apos; &apos; report the ß phase as body-centered cubic. The authors believe that much of the confusion resulted from incomplete homogenization of the alloys. It was found that it takes almost a week to produce a completely homogeneous two-phase structure at 370°C. Another possible source of error is the volatilization of cadmium from the surface and cracks of the specimen during annealing with the possible formation of a cadmium-poor phase. The only metallographic data available on this eutectoid alloy are those of Durrant1 and Fraenkel and Wolff.5 his information is very limited. The phase diagram chosen as most nearly correct is that in the Metals Handbook," which is almost identical with the phase diagram suggested by Durrant.1 Materials and Procedure The alloy was prepared from high purity silver (999.5 plus fine) and high purity cadmium (0.02 Pb, 0.005 Cu, trace Fe). The metals were melted in carbon under a nitrogen atmosphere. Some cadmium was lost by volatilization but by remelting several times and adding cadmium a 200 g ingot of approximately eutectoid composition was obtained. Because the CdO fumes are quite toxic, the exit gas from the furnace was led through a water trap to condense the CdO. The ingot was homogenized at 650°C for 24 hr, no attempt was made to prevent decadmiumization. This resulted in a narrow decadmiumized rim on the ingot which was removed. Filings for analysis were obtained from cross sections of the ingot. Silver was determined by the Volhard thiocyanate method,? cadmium by difference. The ingot analyzed 50.8 pct Ag at the top of the ingot and 50.7 pct at the middle and at the bottom of the ingot. This alloy was only slightly hypo-eutectoid compared to the accepted value of 50.5 pct Ag and was used for the isothermal studies. It was felt that further attempts to adjust the composition might be futile because of the volatility of the cadmium. The ingot was sectioned to specimens approximately 3/4 x 3/4 x 1/8 in. These specimens were "austeni-tized" for 3 hr in air at 650°C to yield homogeneous ß and then quenched into a salt bath held at constant subeutectoid temperatures where the transformation

Jan 1, 1954

By M. R. J. Wyllie, Walter D. Rose

The use of electrical well logs for the quantitative determination of such reservoir parameters as connate water saturation, formation permeability and connate water salinity has recently been attracting the attention of a number of workers1,2,3,4,5. While the theory of the determination of connate water salinity from the self potential S.P. log has received rather detailed treatment, relatively little attention has been paid to the theoretical aspects involved in the quantitative interpretation of resistivity data. It is clear that if electrical well logs can be used for the quantitative evaluation of physical characteristics of reservoir rock, they will provide a valuable tool to supplement cheaply information obtained by more laborious core analysis. In certain cases it is conceivable that the coring program could be considerably curtailed if the electric log could be relied upon to give reasonably accurate quantitative information. It is our object in this paper to examine the theoretical basis of quantitative log interpretation as expressed in such well-established logging concepts as formation factor and cementation factor. It is also our object to investigate the physical aspects of the rela- tionship which is presumed to exist between resistivity index and brine saturation in reservoir rock. In particular, we will endeavor to draw attention to the fact that it is possible to express these logging concepts in terms of capillary pressure-saturation relationships, permeability and tortuosity, parameters which we will consider in this paper to be fundamentally indicative of rock texture. The probability of being able to obtain from log recordings alone the data theoretically essential to permit quantitative log interpretatiou will be examined, and consideration will also be given to the problem of formulating simple semi-empirical relationships for use in the field. THEORETICAL CONSIDERATIONS The concept of formation resistivity factor, or as it is now commonly called, formation factor, appears to have been introduced by G. E. Archie",". Formation resistivity factor as defined by Archie is the resistivity of a rock 100 per cent brine-saturated divided by the resistivity of the brine. This relationship had previously been used by physical chemists and the concept is, for example, implicit in an early treatment by Fricke7 of the conductivity of aqueous slurries. Archie was able to establish the fact that although the relationship between formation factor and porosity did not change too greatly for different rock textures, the relationship between formation factor and permeability was markedly dependent on rock texture. The relationship found by Archie between formation factor and porosity was expanded by H. Guyod8. Archie had suggested that the formation factor, F, could be expressed in terms of the porosity, F, as follows: F =F-m .... (11 where, F = effective porosity fraction of the sand, and, m = an exponent greater than about 1.3. Guyod introduced the term cementation factor for the exponent, m, since it may be qualitatively shown that formation factor for any given porosity tends to increase as the sand becomes more cemented. It would appear that the concept of m as an exponent representing the cementation in a sand is strictly only of qualitative value, since the quantitative estimation of m by means other than the determination of F and F has not been possible. Guyod suggests that m can be estimated from a visual examination of a sample of rock, but this method appears to be both arbitrary and unreliable. It can be said, however, that in intergranular reservoirs of practical interest, the range of m is from about 1.4 to 3.0. A possible explanation for the fact that m is found to lie between these limits will be offered below. The determination of formation factor is an essential preliminary to the

Jan 1, 1950