Search Documents

By L. A. Wilson, P. J. Root

The relative costs of heating a reservoir by steam injection and by combustion have been examined. The comparison was based on a model similar to that proposed by Chu.' The cost of boiler feed water, the price of fuel, pressure and plant capacity were parameters in determin-ing the costs of air compression and steam generation. The analyses compare the cost of heating to the same radius by the two methods. Results suggest that the two primary factors for comparison are the price of fuel and the amount of crude burned during underground combustion. The cost of fuel has a greater effect on the cost of heat from steam than it does on its cost by combustion. As a result, analyses indicate that when the price of fuel is low, steam may be unequivocally cheaper than air. The influence of heat loss is such, however, that as the heated radius increases combustion becomes relatively more competitive depending upon the amount of crude burned. This implies that steam may be cheaper for small stimulation jobs (huff and puff) but combustion may be more economically attractive for heating large areas (flooding). INTRODUCTION Use of thermal methods of recovery is an accepted fact today. After an induction period of several years, processes are being widely used that involve reservoir heating to augment recovery. Of the several techniques, steam injection and forward combustion appear to be destined to dominate the field. Although the objectives of both are the same, the basic differences between generating heat in situ and injecting heat after surface generation influence the cost in different ways. This study compares the cost of heating a reservoir either by steam injection or by forward combustion. There has been no consideration of recovery. Presumably, recovery from the swept region would be high in either case. The sole consideration was the cost of heating to the same radial distance by either process. PROCEDURE THE MODEL The basis for comparison was a mathematical model similar to that used by Chu' for combustion. The model simulates a radial heat wave in two-dimensional cylindricaI coordinates. It includes heat generation, conduction and convection within the reservoir and conduction in the bounding formations. Thus, heat losses from the formation are considered. Three significant modifications were made. 1. Equal logarithmic increments rather than equal increments were used for the mesh spacing in the r direction. By this technique large distances were simulated with relatively few mesh spaces. 2. A backward difference approximation to the convection term was used to avoid troublesome oscillations which result from a central difference approximation when the convection term is large. 3. The radial increments of the combustion zone motion were not necessarily uniquely related to the mesh configuration. The cumbersome step function introduced by the heat of vaporization of steam was circumvented by assuming the enthalpy of the steam to be a linear function of temperature between reservoir temperature and steam temperature. This is equivalent to assuming an average heat capacity numerically equal to the difference between the enthalpy of saturated steam and the enthalpy of water at reservoir temperature divided by the difference between the two temperatures. Heat losses obtained by this model are in essential agreement with those obtained by the analytical solution of Rubenshtein.' A detailed description of the model is presented in the Appendix. Using the model, the times required to heat to particular radial distances were obtained as a function of injection rate and other physical parameters. For the steam case, injected fluid was assumed to be saturated steam at pressures of either 500, 1,000 or 1,500 psia. The corresponding temperatures are 467, 544 and 596F, respectively. Thickness ranged from 10 to 50 ft and injection rate ranged from 100,000 to 1 million Ib/D. Reservoir and overburden temperatures at the injection well were assumed to be that of saturated steam at the injection pressure. The effect of maintaining the overburden temperature at the well at a different temperature (initial reservoir temperature) was examined with no significant change in behavior. The influence of wellbore heat losses for the steam case was determined in the following manner. The rates of heat loss as a function of time were estimated using an approach similar to that suggested by Ramey." he data were based on injection through 2%-in. tubing in 7-in. casing. Integration of these data over the entire iniection period yielded the total heat loss. Total heat losses were then corrected to their equivalents in steam (this number resulted from dividing the total heat loss by the latent heat). This was considered additional steam required to accomplish the reservoir heating and the total cost was increased accordingly.

Jan 1, 1967

By P. E. Rider, K. A. Gschneidner, O. D. McMasters

The solid solubilities for thirteen rare-earth metals in gold were determined by using the X-ray parametric method. Solubilities ranged from 0.1 at. pct for lanthanum in gold up to 8.8 at. pct for scandium in gold. The solubilities from lanthanum to gadolinium were very small and essentially constant, but a sharp increase occurred from gadolinium to scandium. The large solubilities for the heavy rare-earth metals were not expected because of the large size and electrochemical differences between rare-earth atoms and the gold atom. Contributions from first- and second-order elasticity theory plus an electronic contribution were found to reasonably account for a more favorable size factor. Electron transfer from the rare-earth metal to the gold Is thought to occur such that the resultant rare-earth and gold electronegativities are favorable for solid-solution formation. It was also found that this mutual adjustment of size and electronegutivity does not occur if the pure-metal size factors are greater than a critical value of 25 pet. The eutectic temperatures for ten systems were determined and these remained fairly constant at approximately 809 "C for the lighter lanthanide metal-gold systems until the Er-Au system was reached, at which Point the eutectic temperature successively increased reaching a maximum of 1040°C in the Sc-Au system. This rise was correlated to the size factor becoming more favorable for solid-solution formation at erbium. The valence state of ytterbium was found to change from two in the pure metal to three when ytterbium is dissolved in the gold matrix. RECENT results1 reported concerning the solubility of holmium in copper, silver, and gold, showed that the solubility of holmium in gold was quite large, 4.0 at. pct, compared with 1.6 in silver and 0.02 in copper. The small solubilities of holmium in silver and copper are quite reasonable in view of the large size difference (22.2 and 38.2 pct, respectively), large electronegativity difference (0.59 for both systems), and possible unfavorable valency factor (assuming one for silver and copper and three for holmium). The large solubility in gold, however, is unexpected because these same factors are also unfavorable for holmium and gold (22.5 pct size difference and 0.69 electronegativity difference), and because the light rare-earth metals, lanthanum, cerium, and praseodymium, have negligible solid solubilities in gold.2 In view of this unexpected behavior, it was felt that a study of the solid solubilities of most of the rare-earth metals in gold would be desirable to better understand the factors involved in the formation of solid solutions. Of the rare-earth metals added to gold in this study, only ytterbium is divalent in the pure metallic state (the other rare-earth metals are all trivalent) and many of its physical properties (such as the metallic radius, electronegativity, and so forth) are much different from those of the normal trivalent rare-earth metals.' The properties of ytterbium are such that one would expect solid-solution formation to be less favorable for ytterbium in gold than for any of the normal trivalent rare-earth metals. But chemically ytterbium is known to possess a stable trivalent state, and it is quite possible that ytterbium may alloy as a trivalent metal under certain conditions rather than as a divalent metal. Because of the dual valency nature and because so little is known about the alloying behavior of ytterbium, the gold-rich ytterbium-gold alloys are of special interest. EXPERIMENTAL PROCEDURE Materials. The gold used in this investigation had a purity of 99.99 pct with respect to nongaseous impurities. In general the rare-earth metals were prepared by reduction of the corresponding fluoride by calcium metal.3 The impurity contents of the metals used in this study are given in Table I. Preparation of Alloys. Two- or 3-g alloy samples were prepared by arc melting. The samples, with the exception of some of the Er-Au alloys, had weight losses of 0.5 pct or less. All alloy concentrations noted in this paper are nominal compositions. After arc melting, the alloys were wrapped in tantalum foil, sealed off in quartz tubing under a partial atmosphere of argon, homogenized for approximately 200 hr at 780°C, and then quenched in cold water. X-Ray Methods. The X-ray parametric method was used in determining the solubility of the rare-earth metals in gold. filings were sealed in small tantalum tubes by welding under a helium atmosphere. The tantalum tubes were then sealed in quartz tubing under a partial argon atmosphere, and annealed for times ranging from 1/2 to 3 hr (the length of time was inversely proportional to the an-

Jan 1, 1965

By R. C. Ku, A. J. McEvily, T. L. Johnston

The microplastic response of a series ofas-quenched Fe-Ni-C martensites has been measured at 77°K. At strains less than JO'3 the flow stress is governed primarily by the transformation-induced dislocation structure of the martensite. Only at strains in excess of 10-3 is the influence of carbon manifested in the flow stress. At these macroscopic strains, typically 10-2, the solid-solution hardening is proportional to (wt pct C)1/3, and, in an alloy containing 0.39 wt pct C, amounts to 50 pct of the flow stress. THE technological significance of high-strength ferrous martensite has stimulated many investigations of its structure and properties. Although our knowledge of the characteristics of martensite has increased immensely, especially with the advent of high-resolution techniques, an understanding of the basic strengthening mechanism still remains elusive. The purpose of the present paper is to consider certain aspects of micro-plastic behavior of Fe-Ni-C martensite which we feel can help to resolve this important problem. Such alloys are particularly suitable for experimental investigation because their compositions can be adjusted to reduce the M, to a temperature low enough essentially to eliminate the diffusion of carbon in the freshly formed martensite.1 The mechanical properties in this condition are of interest inasmuch as they reflect a state that is free of the important but complicating influence of precipitation processes. In this virgin martensite the carbon is distributed as it was inherited from the parent austenite; i.e., it is present interstitially, and gives rise to tetragonality through strain-induced ordering.' In order to determine the source of strength of such alloys, Winchell and Cohen1 investigated the low-temperature macroscopic stress-strain behavior of a series of virgin martensites of increasing carbon content but of common M, temperature (-35°C). They found that the flow stress increased rapidly with carbon content up to 0.4 wt pct; beyond this point the flow stress increased at a much slower rate. It was concluded that martensite is inherently strong. To account quantitatively for the strength of virgin or as- quenched martensite in terms of the role of carbon, Winchell and cohen3 suggested that the carbon atoms, trapped in their original positions by the diffusionless martensite transformation, interfere with dislocation motion according to a model akin to that of Mott and Nabarro. 4 In this treatment, individual carbon atoms are considered to constitute centers of elastic strain and thereby generate an average stress resisting the motion of dislocations throughout the lattice. The additional stress necessary to move dislocations, over and above that necessary for motion in a carbon-free martensite, is given by where L is an effective length of dislocation capable of motion. L was assumed to be limited to the spacing between the twins that are an essential structural element of Fe-Ni-C martensites. They assumtd the spacing to be invariant and of the order of 100A. However, recent work5 has shown that L is variable and can be in excess of 1000Å, so that the assignment of an appropriate value of L is not straightforward. In contrast to the above conclusion that there is an intrinsically high resistance to plastic flow, it has been suggested by Polakowski6 that freshly quenched martensite is in fact "soft" in the sense that dislocations are initially free to move upon application of stress. The high indentation hardness and macroscopic yield stress of ferrous martensites are then a consequence of rapid strain hardening that depends upon carbon in solution. Consistent with this point of view are the results of Beau lieu and Dubé who measured the rate of recovery of internal friction as a function of aging (tempering) temperature in a freshly quenched steel containing 0.90 wt pct C, 0.37 wt pct Mn, 0.1 wt pct Cr, and 0.07 wt pct Ni. The kinetics were clearly consistent with the idea that many dislocations are unpinned in the as-quenched state and that during aging they become progressively pinned by carbon at a rate controlled by carbon diffusion in the body-centered martensite lattice. In order to provide a basis upon which to distinguish between the "hard" and "soft" interpretations indicated above, we have made studies of the initial stages of plastic deformation in Fe-Ni-C martensites similar to those'used by Winchell and Cohen. It will be shown that the results support the contention that dislocation segments in as-quenched material are indeed

Jan 1, 1967

By T. F. Kassner

The dissolution kinetics of type 304 stainless steel in the Bi-Sn eutectic alloy have been investigated under the well-defined hydrodynamic conditions produced by the rotating-disc sample geometry. In addition, the mutual solubilities of iron, chromium, nickel, and manganese from 304 stainless steel in the eutectic alloy were determined over the temperature range 450" to 985°C. The convective -diffusion model for mass transport from a rotating disc was used to interpret the experinlental dissolution data. The dissolution process was found to be liquid-diffusion-controlled under specific conditions of temperature and Reynolds number. Liquid penetration into the 304 stainless steel resulted in a reduction of the di,ffusion-controlled mass flux and thus precluded the calculation of the diffusion coeficients of the four components from 304 stainless steel in the Bi-Sn eutectic alloy. The convective-diffusion model for diffusional limitations of electrode reactions and mass transport at the tationssurface of a rotating disc set forth by Levich 1,2 has found wide applicability in the investigation of electrochemical and dissolution phenomena in aqueous systems. Riddiford 3 and Rosner have reviewed the model and also include numerous references on work of this nature. More recently the rotating-disc system has been applied to the investigation of hetereogeneous reactions in liquid-metal systems. Shurygin and Kryuk 5 have measured the dissolution rates of carbon discs in molten Fe-C, Fe-Si, Fe-P, and Fe-Ni alloys. Shurygin and shantarin6 also studied the dissolution kinetics of iron, molybdenum, chromium, and tungsten, and the carbides of chromium and tungsten in Fe-C solutions with a rotating-disc sample geometry. In these systems it was possible to distinguish between diffusion and reaction control mainly through experimental confirmation of the velocity dependence of the dissolution rate predicted by the model. However in the absence of dependable solubility data and the virtual lack of diffusion data in these systems, a quantitative check of the magnitude and the temperature dependence of the rate was not possible. In many instances, estimates of the activation energy for solute diffusion and the diffusion coefficient based upon the experimental dissolution data are not credible. A recent study by this author7 has resulted in a critical test of the model in a liquid-metal system. The solution rates of tantalum discs in liquid tin were measured over a wide range of temperature and velocity conditions. In addition, the solubility and diffusion coefficient of tantalum in liquid tin were determined as a function of temperature. The latter data were used with the model to predict both the magnitude and the temperature dependence of the dissolution flux. In that work it was also deemed necessary to reevaluate the solution to the convective diffusion equation to incorporate the effect of the lower range of Schmidt numbers encountered in liquid-metal systems. Good agreement between the model and the experimental dissolution data in the region of diffusion control was obtained in the Ta-Sn system. The Bi-Sn eutectic alloy is used as a seal between the reactor head and the reactor vessel in the Experimental Breeder Reactor-11. The alloy is fused periodically prior to fuel-handling operations. In that connection, it was necessary to investigate the compatibility of the liquid alloy with the type 304 stainless-steel containment material. The results of a rotating-disc study in this multicomponent system are presented. EXPERIMENTAL METHOD The 5.08-cm-diam discs were machined from 0.317-cm-thick plate. Chemical analysis information for the type 304 SS material is given in Table I. The discs were ground flat on metallographic paper and given a final polish on Linde B abrasive. A thin support rod was threaded into the disc and the region around the threads was fused under an inert gas. The support rod was fitted with a quartz protection tube and then was attached to a supporting shaft which passed through a rotary push-pull vacuum seal. The disc and supporting shafts were dynamically balanced prior to insertion into the furnace tube. The apparatus is shown schematically in Fig. 1. The 58 pct Bi-42 pct Sn eutectic alloy melts were prepared from 99.995 pct pure Bi and Sn by fusing the components in a 7-cm-ID Pyrex crucible. The system in which the melts were made was evacuated to a pressure of 1 x 10-6 Torr and back-filled with purified argon several times before melting the charge. The ingot was reweighed and placed in a slightly larger-diameter Vycor crucible used in the dissolution runs. A run was started by lowering the disc into the liquid

Jan 1, 1968

By W. C. Maurer

A study of bit-tooth penetration, or crater forniation. under simulated borehole condirions has been made. Pressure conditions existing when drilling with air, water and mud have been sirnulated for depths of 0 to 20.000 ft. These crater tests showed that a threshold bit-tooth force must he exceeded before a crater is .formed. This thresh old force increased with both tooth dullness and diflerenrial pressure between the borehole and formalion fluids. At low differential pressures, the craters formed in a brittle manner and the cuttings were easily removed. At high differenlial pressures, the cunings were firmly held in the craters and the craters were formed by a pseudoplas-tic mechanism. With constant farce of 6,500 16 applied to the bit reeth, an increase in differential pressure (sitnulated mud drilling) from 0 to 5,000 psi reduced the crater volumes by 90 per cent. A comparable increase in hydrostatic fluid pressure (simulated water drilling) produced only a 50 per cent decrease in volutne while changes in overburden pressure (simulated air drilling) had no detectable effect on crater volume. Crater tests in unconsolidated sand subjected to differential pressure showed that high friction was present in the sand at high pressures. Similar friction belween the cuttings in craters produces the transition from brittle to pseudo plastic craters. INTRODUCTION The number of wells drilled below 15,000 ft increased from 5 in 1950 to 308 in 1964. Associated with these deep wells are low drilling rates and high costs. High bottom-hole pressures produce low drilling rates by increasing rock strength and by creating bottom-hole cleaning problems. This paper describes an experimental investigation of crater formation under bottom-hole conditions simulating air, water and mud drilling. Although numerous investigators have studied bit-tooth penetration (cratering) at atmospheric pressure conditions, only limited work has been done on cratering in rocks subjected to pressures existing in oil wells. Payne and Chippendale2 have studied cratering in rocks subjected to hydrostatic pressure using spherical penetrators. Garner et aLJ conducted crater tests in dry limestone by varying overburden pressure and borehole fluid pressure independently and using atmospheric formation-fluid pressure Gnirk and Cheathem4,5 have studied crater formation in several dry rocks subjected to equal overburden and borehole pressure and atmospheric formotion pressure. Podio and Gray studied the effect of pore fluid viscosity on crater formation using atmospheric borehole and formation-fluid pressurc and varying overburden pressure. Although these studies have provided useful information on crater formation under pressure, they were limited in that the three bottom-hole pressures could not be varied independently and, therefore, that many drilling conditions could not be simulated. The prersure chamber used in this study allowed visual observation of the cratering mechanism and independent control of the borehole, formation and confining pressures. By using different fluids in the chamber, pressure conditions existing in air, water and mud drilling to depths of 20,000 ft were simulated. The mechanisms involved in cratering at these different pressure conditions were studied for teeth of varying dullness and at different loadins rates. High-speed movies (8,000 frames/sec) and closed-circuit television were used to visually study the crater mechanism under pressure. EXPERIMENTAT PROCEDURE PRESSURE CHAMBER The Pressure chamber in Fig. I was used to simulate bottom-hole pressure conditions. This chamber has been pressure-tested to 22,500 psi and is normally operated at pressures up to 15,000 psi. The chamber contains four lucite windows' used for illuminating and observing the crater mechanism under pressurc. A closed-circuit television and a Fastax camera (8,000 frames/sec) have been used in these studies. Cylindrical rock specimens (8-in. diameter X 6-in. long) were subjected to three independently controlled pressures simulating overburden, borehole fluid and formation-fluid pressures. Overburden pressure, which corresponds to the stress induced by the overlying earth mass, was applied by exerting fluid pressure against a rubber sleeve surrounding the rock. Borehole pressure, which is the pressure exerted by the column of mud in the wellbore, was simulated by applying pressure to the fluid overlying the rock in the chamber. Formation pressure was simulated by applying pressure to the water saturating the rock. The borehole and formation pressures were equal except when mud was used in the chamber, in which case the differential pressure between these fluids acted across the mud filter cake.

Jan 1, 1966

By T. C. Boberg, R. B. Lantz

This paper presents a method for calculating the producing rate of a well as a function of time following steam stimulation. The calculations have proved valuable in both selecting wells for stimulation and in determining optimum treatment sizes. The heat transfer model accounts for cooling of the oil sand by both vertical and radial conduction. Heat losses for any number of productive sands separated by unproductive rock are calculated for the injection. shut-in and production phases of the cycle. The oil rate increase caused by viscosity reduction due to heating is calculated by steady-state radial flow equations. The response of successive cycles of steam injection can also be calculated with this method. Excellent agreement is shown between calculated and actual field results. Also included are the results of several reservoir and process variable studies. The method is best suited for wells producing from a multiplicity of thin sands where the bulk of the stimulated production comes from the unheated reservoir. The flow equations used neglect gravity drainage and saturation changes within the heated region. INTRODUCTION This paper presents a calculation method which can be used to predict the field performance of the cyclic steam stimulation process. The calculation method enables the engineer to select reservoirs that have favorable characteristics for steam stimulation and permits him to determine how much steam must be injected to achieve favorable stimulation. While the calculation represents a considerable simplification of physical reality and the results are subject to numerous assumptions which must be made about the reservoir, it has been found that realistic calculations can be made of individual well performance following steam injection. The duration of the stimulation effect will depend primarily on the rate at which the heated oil sand cools which, in turn, is determined by the rate at which energy is removed from the formation with the produced fluids and conducted from the heated oil sand to unproductive rock. A complete mathematical solution to this problem is a formidable task, and finite difference techniques would undoubtedly have to be used. The calculation method pre- sented here utilizes analytic solutions of simple related heat transfer and fluid flow problems. The method is sufficiently simplified that it can be used as a hand calculation, although the calculations are somewhat lengthy and laborious. For that reason, the analysis was programmed for an IBM 7044 digital computer. Well responses observed at the Quiriquire field in eastern Venezuela&apos; have been matched using this program after making suitable approximations for reservoir and wellbore conditions. One of the most valuable uses of this calculation method is to assess the effect of reservoir and proc-cess variables on the stimulation response. This paper contains results of several studies made of key reservoir and process parameters. Among the most important of these is the influence of prior wellbore permeability damage. If a well is severely damaged prior to stimulation, a higher stimulation response will be observed than if it is undamaged. If a portion of this damage is removed, a permanent rate improvement will occur. THEORY I)ES(:KJJ&apos;TION OF CALCULATION METHOD The process of cyclic steam stimulation is essentially one of reducing oil viscosity around the wellbore by heating for a limited distance out into the formation through the injection of steam. Suitable modifications of the calculation technique presented here can be made so that stimulation of wells by hot gas injection or in situ combustion can also be calculated. A schematic drawing of the heat transfer and fluid flow considerations included in the calculation method is shown in Fig. 1. In brief, the calculation assumes that the oil sand is uniformly and radially invaded by injected steam. For wells producing from several sands, each sand is assumed to be invaded to the same distance radially. In calculating the radius heated rn energy losses from the wellbore and conduction to impermeable rock adjacent to the producing sands are taken into account. After steam injection is stopped, heat conduction continues and oil sands with r < ra cool as previously unheated shale and oil sand at r > r, begin to warm. The effect of warming of oil sand out beyond r, has little effect on the oil production rate compared to the effect of cooling of the oil sand nearer the wellbore than ra. Thus, in computing the oil production rate, an idealized step function temperature distribution in the reservoir is assumed where the original temperature exists for r > rn and where an average elevated temperature exists for r < rn. The average temperature in the oil sand for the region r < rn is computed as a function of

Jan 1, 1967

By W. M. Wojcik, R. F. Kowal

Oxygen and sulfur distributions in commercial, 5-ton ingots of killed, medium carbon steel are described. Oxygen distribution is found to vary with deoxidation practice. Irregular distribution of oxygen within ingots makes necessary special precautions in sampling of rolled products for analysis of oxygen. Oxygen distribution is discussed in terms of recently published solidification concepts which had been successfully applied to simpler cases of segregation. These concepts have been found inadequate to explain observed oxygen distributions. Convective movements of the liquid metal, as determined by tracer elements, are shown to be capable of accounting for the observed distributions of oxygen. IN an effort to explore the origin of surface and subsurface imperfections in pierced steel products, a study of oxygen and sulfur segregation was made on ingots cast in open-top and hot-top molds. The results of our previous investigations1"3 have indicated the importance of the location and amount of oxide inclusions in an ingot. Inclusions close to the surface of the ingot have been found to contribute greatly to the formation of imperfections in the surface of finished products. This study of the effects of deoxidation and casting practice on segregation and the resulting oxygen distribution in ingots was initiated to determine the parameters controlling the location of inclusions in an ingot. Segregation of solute elements during solidification of low-melting binary alloys has been studied in the past.1, 5 Formation and growth of inclusions in iron melts have been studied under specific conditions."- In spite of these and other recent studies,10-12 segregation during solidification of commercial, killed steel ingots is not well understood. Consideration of solidification rates, of segregation during solidification of the chill, dendritic, and central zones, and of material balances for the segregated elements has indicated that a simplified theoretical solidification model is not adequate. However, the observed high oxygen contents in localized volumes of the dendritic zone can be rationalized if additional effects of convection currents in the ingots, precipitation, and rapid growth of new phases are considered. EXPERIMENTAL PROCEDURE Steelmaking and Processing. A group of nine killed. medium carbon steel heats having compositions listed in Table I have been studied. The deoxidation and mold practices used were varied to give a wide range of steel oxygen contents. The amounts of aluminum added to the ladle and the ingot casting practices (hot top and open top) were the main variables. The steel was made by a duplex practice in 160-ton tilting basic open-hearth furnaces. All nine heats were top-cast into 24 by 24 in. big end down, fluted molds, to a height between 60 and 76 in., using both open tops and exothermic hot tops. The deoxidation practice and the tapping and teeming details for each heat and ingot studied are given in Tables II and III, respectively. Hot-top practice is indicated by the letter H following the heat designation. Furnace and ladle temperatures were measured by standard disposable-tip, Pt/10 pet PtRh thermocouples. Teeming-stream temperatures were obtained as described by Samways et al.,13 by immersing a Pt/10 pet PtRh thermocouple, covered by a silica sheath, into the teeming stream under the nozzle. The output of this thermocouple was recorded with Leeds & Northrup Speedomax potentiometer. Calibration of the latter thermocouples was based on the freezing point of a pure iron/oxygen alloy (2795°F). The accumulated errors of measurements were within ±10°F. The thermocouple measurements were supplemented in this investigation by continuous recording of a ratioing, two-color pyrometer (Shawmeter), protected from smoke by a blast of clean air within the sighting tube, and calibrated to read with better than ±10°F accuracy. Following teeming of three heats, P, R, and T, tracer elements were added to the steel in the molds to obtain a record of the progress of solidification. As soon as the teeming stream was shut off, a 0.010-in.-thick steel can containing a mixture of crushed standard ferro-titanium and ferro-vanadium (0.05 pet of each alloy element) was plunged into the middle of the steel pool to a depth of 6 in. In about 30 sec no indication of the can or its contents remained. The surface of the open-top ingots solidified in 20 to 30 sec. A study of liquid metal movement and the precipitation of oxides was facilitated materially by use of the tracer technique as titanium has a low distribution coefficient between solid and liquid steel while vanadium has a high distribution coefficient.

Jan 1, 1965

By A. L. Mincher, W. F. Sheely

Mechanical properties of columbium have been studied over the temperature range of -196 to 1093oC. The decreased strengthening influence of cold-work at temperatures below ambient has been interpreted in terms of the Peierls-Nabarro effect. Maxima in the rate of strain hardening observed during tensile testing in the range 250-600°C. have been correlated with interstitial impurities to indicate the temperature ranges at which carbon, oxygen, and nitrogen, respectively, are responsible for strain aging. THE growing need for structural materials for use above the useful service temperatures of the iron-, nickel-, or cobalt-base alloys has caused the refractory metals to be considered as potential engineering materials. These metals, which include columbium, tantalum, molybdenum, and tungsten, are called refractory because the lowest melting point among them,that of columbium, is about 1000°C higher than the average melting temperatures of conventional high-temperature alloys. They are all body-centered cubic transition metals and, as such, their mechanical properties have basic characteristics which distinguish them from the face-centered cubic metals. For example, all show a much steeper rise in strength with decreasing temperature below room temperature than do the face-centered cubic metals, and their mechanical properties are strongly influenced by interstitially dissolved impurities. In order that these new metals may be used efficiently, it is necessary that their characteristics of behavior be fully known. In this paper, the mechanical properties of columbium will be examined over a wide range of temperatures. In particular, the influences of cold-work and individual species of interstitial impurity atoms on mechanical properties will be described, and basic mechanisms which may control the observed characteristics will be explored. EXPERIMENTAL The material used in this investigation was Union Carbide Metals Co. columbium roundels consolidated to four 4-in. diam ingots, three by consumable-electrode arc melting and one ingot by electron beam melting. Impurity contents of the ingots and methods of ingot conversion and treatment are summarized in Table I. The only metallic impurity occurring in any significant quantity was tantalum at about 0.1 pct. Iron, silicon, titanium, and zirconium were each less than 0.015 pct; boron was 1 ppm or less. This should have no appreciable influence on properties. The electron beam melted material, being the purest, will be used as the basis for comparison in the discussions to follow. Tensile tests were conducted from-196 to 1093oC, on both cold-worked and fully recrystallized arc-melted and electron-beam melted columbium using standard 1/4-in. diam, 1-in. long gage length test specimens. A strain-rate of 0.005 in. per in. per min was employed until the 0.2 pct yield strength was achieved and then the strain-rate was increased to 0.05 in. per in. per min for the balance of the test. Samples were protected in an inert atmosphere at tests above 300°C. The tensile properties obtained on the electron-beam melted columbium, E, in both the cold-swaged and recrystallized conditions are given in Fig. 1. The yield strength data of Dyson, et al.,&apos; obtained on recrystallized electron beam melted columbium and the tensile strength data reported by Tottle2 on powder metallurgy columbium are included in Fig. 1. The material used by Tottle had been purified by vacuum sintering. There is excellent agreement between Dyson&apos;s data and those obtained in the present investigation. The tensile strengths obtained by Tottle were slightly greater than those obtained in this investigation on electron-beam melted columbium but varied with temperature in a similar manner. Tottle&apos;s data showed a maximum in tensile strength near 500°C, as did our data on electron-beam melted material, and also showed a small maximum at 300°C. The significance of these maxima will become evident later in the discussion. The tensile properties of cold-swaged and recrys-tallized arc melted columbium are plotted in Fig. 2. It was found that the properties of the recrystallized arc-melted columbium from all three heats showed very close agreement except at temperatures between about 500" and 800°C. A reason for this range of disagreement will be suggested in the discussion. The generally good agreement, however, attests to the ability of cold-working and subsequent recrystal-lization to erase the effects of the three different primary breakdown procedures and to produce nearly equivalent structures in the samples derived from the three different heats. wesse13 reported tensile data on columbium having interstitial impurity contents between those of the

Jan 1, 1962

By R. H. Jacoby, V. J. Berry

When dry gas is injected into a reservoir containing a volatile crude oil, a significant amount of the reservoir liquid phase may become vaporized. The resultant rich gas phase, when subsequently produced, con-tributes to tank oil production. This contribution assumes greater importance, the more volatile the oil in-olved. Oil recovery may be sub-stantially greater than that predicted by conventional frontal-drive methods, which do not consider the vaporization equilibrium between the reservoir oil phase and the injected gas. A calculation method has been (developed to account for vaporiza-. tion of the reservoir liquid phase during gas-injectiorz operations, and for the tank oil production which results from this factor. Recovery performance calculations are presented for a reservoir containing a highly volatile oil. Tank oil recovery is calculated to be about twice that predicted by the use of the conven-rional frontal-drive equations. In contrastto usual pressure maintenance performance results, in which the gas-oil ratio rises at an increasing rate after gas breakthrough, the pre-dicated gas-oil ratio rises rapidly to about 12,000 scf/bhl and then rises much less rapidly. During gas inje-tion, most of the reservoir liquid phase contacted is evaporated by the dry injection gar. The gas-oil ratio during this period is dependent upon reservoir pressure. The higher the operating pressure, the lower the gas-oil ratio. The predicted behavior i., in accordance with laboratory PVT tests made to .simulate the vaporization behavior. In addition to recovery performance predictions, results of the calculation procedure include complete wellstrearn composition data of value in the design of gasoline plant facilities often used in con-rrction with gas-injection operations. INTRODUCTION In the cycling of gas-condensate reservoirs, dry gas is injected to maintain reservoir pressure during wet gas production and to thereby eliminate or reduce ultimate loss of liquids due to retrograde condensation within the reservoir. Gas injected into crude oil reservoirs has a dual function. It displaces oil to the producing wells and at the same time serves to partially or fully maintain reservoir pressure. Oil shrinkage which would occur upon pressure reduction is thereby minimized or eliminated. Accepted calculation methods are available&apos;= for predicting recovery performance of either gas-condensate reservoirs or crude oil reservoirs which are being subjected to gas injection. In gas-condensate reservoirs, any retrograde liquid formed does not flow, and it is necessary to account only for the vaporization equilibrium between this liquid and the injected gas. Conversely, in normal crude oil reservoirs, both the oil and gas phases flow, but it has not been considered necessary to account for any vaporization of the reservoir liquid which might occur upon contact with the dry injection gas. Recently, high shrinkage reservoir fluids known as "volatile oils" have been found in increasing amounts. These oils are characterized by tank oil gravities above 4.5" API, solution gas-oil ratios above 1,000 scf/bbl, and reservoir volume factors above two. Special techniques have been devised for predicting depletion performance of reservoirs containing such oils.74.5.1; One of the characteristics of reservoirs producing volatile oils is that the reservoir gas phase carries a significant amount of oil which is recoverable as stock-tank liquid. This unusual vaporization be-havior implies that an appreciable amount of reservoir liquid would be vaporized upon contact of the oil with dry injection gas. In a gas-injection operation, tank oil recovery would he obtained not only through frontal displacement of the reservoir liquid by the injected gas, but also through production of the rich gas phase. This means that not only are improved methods needed to predict recovery of such oils from reservoirs undergoing gas injection, but it would also be expected that high oil recoveries might be obtained by such operations. When relatively dry gas is injected in to a volatile oil reservoir, phase equilibrium between the injected gas and the reservoir oil will tend to bc established. Initially the most volatile components, such as propane, the butanes and the pentanes, will account for most of the material transferred from the oil phase to the gas phase. As the partially stripped oil phase is contacted with additional dry injection gas, the heavier intermediate components, such as the hexanes, the heptanes and the octanes, will gradually transfer to the gas phase in increasing amounts. This is because the supply of lighter components in the oil phase dwindles due to the stripping action of the injected gas. This stripping action will usually continue to be effective down

By Peter R. Morris

The effecl of preferred orienlation on X-vay dzffvaction measurements of retained austenzte was investigated for four precipitation-hardening staznless steels in sheet form. A method is preserzted for estimating the ervor in measurement associated with a given samplirig direction. The method was used to select an "optimum" sampling direclion in order to minimize errors in measurement due to preferred orientation. hleasuremenls of retained austenite content employing lhe proposed sampling direction are conzpaved to measuretnents enzploying the more commonly used normal direclion for a series of sawzples. THE first application of X-ray diffraction to the measurement of retained austenite in steels is due to Sekito, 1 who employed a photographic technique in which the (111) reflection from a thin strip of gold affixed to a cylindrical sample was employed as a standard. Averbach 2 introduced the "direct comparison" method in which the ratios of observed to calculated random intensity are assumed proportional to the austenite and/or martensite contents. Averbach&apos;s work forms the basis of most subsequent X-ray diffraction methods for the determination of retained austenite. Subsequent improvements are due to: Averbach and Cohen,3 who employed a sodium chloride crystal to monochromate cobalt radiation; Averbach et a1.,4 who introduced a bent sodium chloride monochromator; Mager,&apos; who used a bent quartz crystal to monochromate chromium radiation ; Littmam, who first used a geiger counter diffractometer for this purpose; Beu and Beu and Koistinen, 11,12 who studied effects of absorption factor, surface preparation, sample geometry, integrated intensity vs peak height, choice of radiation, monochromator, and filter. The possibility of errors in measured values due to orientation effects was noted by Miller,13 who suggested examination of a surface other than the plane of rolling. Lopata and Kula 14 have developed an experimental technique in which the preferred orientation is measured in each sample. They illustrated the method for a sample containing 42 pct retained austenite. Application of their technique to the 1 to 15 pct range typical for the precipitation-hardening stainless steels does not appear feasible. EXPERIMENTAL PROCEDURE The nominal compositions of the precipitation-hardening stainless steels investigated are listed in Table I. Ingots were solution-treated, hot-rolled to approximately 0.2 in., and reduced to 0.050 in. by a suc- cession of cold rolling and annealing operations. After this treatment the 17-4PH sample was in the marten-sitic condition, while the 17-7PH, PH 14-8Mo, and PH 15-7Mo samples were in the austenitic condition. Samples of 17-7PH and PH 15-7Mo steels in the mar-tensitic condition were obtained by heating to 1750&apos;F for 10 min and holding at -100°F for 8 hr. A sample of PH 14-8Mo steel in the martensitic condition was obtained by heating to 1700°F for 1 hr and holding at -100°F for 8 hr, followed by aging at 950" for 1 hr. POLE FIGURE DETERMINATIONS Samples were thinned to 0.003 to 0.005 in. by etching in a solution containing 250 ml reagent-grade phosphoric acid (85 to 87 pct H3PO4), 250 ml technical-grade hydrogen peroxide (30 to 35 pct H 2 O 2), and 50 to 100 ml reagent-grade hydrochloric acid (37 to 38 pct HCl). The specimens were placed in an "integrating" sample holder which provided a 1-in. oscillation in the plane of the sample. The diffractometer was aligned to measure the intensity diffracted by planes of the particular {hkl} type being studied. The sample was Set for a given latitude angle, a, measured from the plane of the sheet, and diffracted intensity recorded as the longitude angle, 0, measured in the plane of the sheet from the rolling direction, was increased from 0 to 360 deg. After a 360-deg scan of B, a was incremented by 5 deg, and the process repeated. Random standards obtained by spraying suspensions of powdered iron (bcc structure) and nickel (fcc structure) in lacquer were used to correct observed intensities for absorption and geometrical effects. Zirconium-filtered molybdenum radiation was used to determine the transmission regions of the (111) (0to 45 deg), (200) (0 to 60 deg), and (220) (0 to 45 deg) austenite and (110) (0 to 45 deg), (200) (0 to 50 deg), and (211) (0 to 35 deg) martensite pole figures. Vanadium-filtered chromium radiation was used to

Jan 1, 1968

By R. E. Schilson, F. H. Poettmann

The direct determination of the stabilized performance behavior of low capacity, slowly stabilizing gas wells is extremely time-consuming and wasteful of gas. From both field experience and theoretical considerations, a test procedure has been evolved by which the stabilized hack pressure behavior of such gas wells can he predicted without having to revert to long time flow tests. The method consists of using the isochrona1 test procedure to establish the slope of the back pressure curve, "n", and the short time variation of the performance coeficient, "C", with time. From this short time transient flow data and theoretical considerations, the value of C at large times can he established. By assuming the radius of drainage of a well to be half the distance between wells, one can calculate the stabilization time for various well spacing patterns. Once the stabilization time for a given .spacing has been determined, the value of C can be calculated and the stabilized back-pressure curve can he. establbrhed. The calculated perfortnance coefficient as a function of tinze was cotnpared to the experintentally measured values for a number of gas wel1.e. The deviation of the calc~*lated from the cxperimental res1tlt.s vary depending on the set of short time experimental points used to evabrate the parcrttzeters of the equation. The longer the tirne far the flow test data user1 in the calculations, the better was the agreement with the experimental results. The time necessary to obtain this data from well tests varies consitlerably, depending on the physical natrcre oj the rt3scrvoir under consideratiot~. INTRODUCTION For many years, the U. S. Bureau of Mines Monograph 7&apos; has served as a guide for testing and evaluating the performance of gas wells by means of the back-pressure method. The back-pressure performance of a gas well is expressed by the following equation: Q~ CW-W.........(I) where the characteristics of the back-pressure equation are determined by C, the performance coefficient, and n, the exponent which corresponds to the slope of the straight line when Q and {P* - PN2) are plotted on logarithmic paper. Q is gas flow rate at standard conditions, and P, and P, are equalized and flowing bottom-hole pressures, respectively. Prior to the development of the back-pressure tcst, the "open flow" capacity method of testing a well was common. By this method, the wells were flowed wide open to the air and the flow rate measured. Such procedure was wasteful of gas and did not provide information on the deliverability of the gas to the pipe line. Monograph 7 Procedure The back-pressure method of testing wells was developed to overcome these shortcomings. Although much has been learned regarding the laws of the flow of gas through porous formations, the original development of the back-pressure relationship was based entirely on empirical methods. The back-pressure behavior provides the engineer with information essential in predicting the future development of a field. It permits him to calculate the deliverability of gas into a pipe line at predetermined line pressures, to design and analyze gas gathering lines, to determine the spacing and number of wells to he drilled during the development of a field to meet gas purchasers&apos; requirements, and to solve many other technical and economic problems. As described in Monograph 7, the flow-after-flow method of back-pressure testing, when applied to fast stabilizing and usually high capacity wells, correctly characterized the behavior of the well. However, as the value of the gas at the wellhead increased, small capacity gas wells having slow rates of stabilization became economically operable. The flow-after-flow method of testing could not be used to describe the behavior of these slowly stabilizing wells. The procedure of Rawlins and Shellhardtl for establishing the back-pressure behavior of a gas well was based on the rcquirement that the data be obtained under stabilized flow conditions; that is, that C is constant and does not vary with time. C depends on the physical properties of the reservoir, the location, extent and geometry of the drainage radius, and the properties of the flowing fluid. In a highly permeable formation, only a very short period of time is required for the well to reach a stabilized condition, and, consequently, the requirements for the test procedure described in Monograph 7 are met. For a given well, n is also constant

By W. A. Oates, J. A. Lambert, P. T. GaIIagher

Monte Carlo methods have been used to compute the arrangements of interstitial atoms dissolved in tetrahedral sites in bcc lattices. It is assumed that the presence of an interstitial atom "blocks " a certain number of neighboring sites and prevents their occupancy. Sites "blocked" by more than one filled site are allowed for. The computed values of. the mean occupation number (defined as the ratio of the total number of sites blocked to the number of solute atoms are used to calculate the configurational entropies of the solutions. These entropies are compared with those resulting from previous theoretical studies of this problem and also with available experin~ental data for the p Zr-H, Nb-H, V-H, and Ta-H systems. Evidence is also given that the "blocking" explanation of low limiting compositions in these systems, rather than this being due to initial limitations on the number of sites available, is probably correct. THE ideal partial configurational entropy of mixing of an interstitial solute in a metal is given by: where p is the number of interstitial sites per metal atom and Xi is the atomic fraction of the interstitial. For the bcc lattice. which we shall be concerned with in this paper, the interstitial positions are shown in Fig. 1. It can be seen that for the tetrahedral sites, p=6. whereas for the octahedral sites, p = 3. Different emphasis has been placed on the relative importance of energy and entropy effects in determining deviations from ideality in interstitial solid solutions. In some cases the same system, e.g., Fe-C, has been described by the contradictory regular and athermal solution models indicating that the enthalpy and entropy functions, derived from equilibrium data, are frequently not accu.rate enough to differentiate between these treatments. However, for certain metal-hydrogen solutions the equilibrium data is available over sufficiently wide ranges of temperature and composition to permit a reasonably accurate determination of the compositional variation of the heats and entropies. Hoch&apos; has attempted to interpret the results of interstitial solid solutions in terms of a regular solution model. In the case of the Ta-H system where 13 = 6, this model entails fitting the experimental relative partial entropies of solution, asH, to the equation: where ASgs is the relative partial excess entropy of solution of hydrogen. Hoch found that the results of Mallett and Koeh1 could be fitted to this equation with an approximately constant value of AF up to XH = 0.25. However, it is apparent from the solubility isotherms in this system which become asymptotic to the composition TaH that, since (Xh /6 - ~Xh ) becomes infinite only at TaH6, it is necessary that AS<&apos; tends to infinity at TaH. In other words, the low saturation composition of TaH, instead of the anticipated TaH,, eliminates the possibility of applying regular solution theory to such systems. Rather large negative excess configurational entropies must exist at higher hydrogen concentrations in order to explain the lower saturation values. To account for these low limiting compositions and excess entropies two distinctly different approaches have been followed. Rees and many others1-l2 have assumed that not all interstitial sites are crystallographically equivalent with respect to the interstitial addition; that is, in Eq. [I] p is less than the value anticipated from geometrical considerations. To describe, say, a bcc metal-hydrogen system with a limiting composition of MH by this approach one would consider that p = 1 in the first instance instead of p = 6.&apos;j3 In some cases, nonintegral values of B have been taken in order to improve the fit with the experimental data over limited ranges of composition. The other approach which has been used to explain the low saturation compositions is to assume that, although all sites are available for occupancy, strong repulsive interactions exist between the neighboring interstitial atoms, and hence occupancy of any site excludes or blocks a certain number of neighboring sites from being occupied. Earliest treatments of this concept considered the exclusion of an integral number, of nearest-neighbor sites from being occupied at all concentrations. In this case, the partial configurational entropy is given by: These early treatments failed to allow for the overlap of the blocked sites which will arise at all but the very lowest concentrations. More recently attempts have been made to calculate the effect of this decrease in the number of blocked sites on the configurational entropy. Using the quasichemical treatment of interstitial solid solutions as given by Lacher and assuming that an infinite repulsive interaction energy existed between the solute atoms. atom obtained an approximate configurational entropy applicable to the blocking with overlap case:

Jan 1, 1970

By J. Szekely, G. P. Martins

A formulation is given for describing the rate of expansion of spherical bubbles rising in liquids the freeboard of which is evacuated. The computer solution of the resultant differential equations has shown that, for low freeboard pressures (less than about 5 mm Hg), on approaching the free surface the bubbles expand much less rapidly than predictable from equilibrium considerations. In other words, in this region the pressure inside the bubbles will be significantly larger than the static pressure in the liquid corresponding to the position of the bubble. These theoretical predictions were confirmed by experimental work, using two-dimensional air bubbles rising in mercury. The important consequence of these findings is that the expansion of gas bubbles in vacuum degassing operations will be a great deal less than expected from hydrostatic considerations. This would lead to a significant reduction in the available interfacial area and may explain the apparent poor efficiency of many vacuum degassing units. VACUUM degassing as a treatment for liquid steel has gained widespread popularity in recent years; the number of known installations exceeds several hundred at the present time.&apos; Although much information is available on both the thermodynamics of the system and the overall performance of various industrial units, much less is known about the fundamental aspects of the process kinetics.2-4 The basic physical situation common to virtually all vacuum degassing operations is the interaction between gas bubbles (swarms of bubbles) and the surrounding molten metal, held in a container, the freeborad of which is at a low absolute pressure. Once formed (or introduced from an external source) the bubbles will ascend, due to the buoyancy forces, and, during this ascent, a significant increase in their volume will occur. This progressive increase in the bubble volume is due to two factors: a) the continuous reduction in the static pressure acting on the bubble during its rise; and b) mass transfer into the bubble from the surrounding molten steel. In a recent paper Richardson and Bradshaw developed equations5 for describing mass transfer into gas bubbles from molten metals at reduced pressures. However, in deriving these expressions it was assumed that the pressure inside the bubble was identical to the static pressure in the adjacent liquid. In other words, the volume of the bubble was considered to be in equilibrium with the pressure of the fluid adjacent to it. This assumption, thus their analysis presented, was thought to be reasonably accurate for most of the bubble&apos;s ascent; however, it was unlikely to be valid in the immediate vicinity of the free surface, held at a low pressure. It was pointed out in the discussion6 that the region close to the surface may be of considerable importance as both the driving force and the interfacial area available for mass transfer are at their highest value here. The &apos; &apos;anomalous" behavior of gas bubbles when approaching a free surface at low pressures was recently confirmed in a preliminary investigation by Szekely and Martins. ?1 Here high-speed motion photography was used to study air bubbles rising in a column of n-tetradecane with a freeboard pressure of 0.1 mm Hg. It was found that significant distortion of the bubbles occurred on approaching the free surface; furthermore, the expansion observed was much less than what one could expect from hydrostatic considerations, i.e., factor a previously discussed. It follows from the foregoing that a detailed study of these phenomena would be justified both from fundamental considerations and because of their potential relevance to technology. The purpose of the paper is to present a more realistic formulation for the expansion of a gas bubble approaching a low-pressure region, together with a comparison of the theoretical prediction with experimental measurement. An inert bubble will be considered in the first instance; it is thought that the understanding of the fluid mechanics is an essential first step toward the realistic formulation of the mass transfer process. This latter problem will be the subject of a subsequent publication. FORMULATION The Physical Model. Consider a spherical bubble, of initial radius Ro, rising in a fluid, having a density pL. Initially let the bubble be at a distance H from the free surface, and at a pressure Pgo, as illustrated in Fig. 1. Pgo, the initial pressure in the bubble, is given by the following equation: pgo = Po = Ptp +pLgH [ la] where Po is the pressure in the liquid corresponding to the initial position of the bubble and Ptp is the pressure at the free surface. The fluid pressure at

Jan 1, 1970

By V. B. Kurfman

Grain refinement of Mg-Al melts by carbonaceous additions has been attributed to nucleation by aluminum carbide. The effects of process and alloy variables are interpreted and predicted in terms of the dispersion and chemistry of this phase. The grain coarsening action of Be, Zr, Ti, R.E., chlorination, temperature extremes, and prolonged holding times is described. Measures necessary to insure an adequate dispersion of the catalyst are discussed. CARBON inoculation treatments have become fairly well known and used for grain refinement of magnesium alloys containing Al. Although there is general agreement that a nucleation process occurs, the process is not understood and the inoculants are used in a rather empirical fashion. The treatment is applied to the class of alloys containing 3 to 10 pct Al, i.e., AZ31A to AM100A. Typical methods involve melting, alloying, and adjusting the temperature to 1400° to 1450°F. Then 0.01 to 0.5 pct C as CaC2, C6C16, or lampblack is added by any convenient means, and the melt poured within 10 to 30 min. Investigators generally have been impressed by an assumed similarity of this refinement process to superheat grain refinement, which depends on heating approximately the same alloys to a temperature in the range of 1550" to 1650°F, then pouring promptly after the melt is cooled to the pouring temperature. Various predictions have been made that carbon refinement would replace superheating in commercial practice due to reduced process costs, but this replacement has not fully taken place because of production difficulties and conflicting observations. Davis, Eastwood, and DeHaven1 agree with Nelson2 and wood3 in suggesting that an excess of inoculant may be harmful. Wood however says that overtreat-ment is not a problem in production use of hexa-chlorobenzene inoculation, and Hultgren and Mitchell4 claim no evidence of harm from excess additions. Various grain coarsening reactions are known to occur, including the possibility of overtreatment mentioned above. Trace amounts of Be,2 Zr, and Ti may prevent refinement by either a carbon treatment or a superheat. Occasionally treatment with cl25 may cause coarsening, although the Battelle refinement process&apos; uses a CC14-C12 blend. Grain coarsening also tends to occur on holding at temperatures below 1350°to 1400°F, especially after a superheat treatment, and for this reason Nelson2 stresses the desirability of a refinement method useful at lower temperatures for open pot melting practice. Since a carbon treatment can be made to work at temperatures below 1400°F, it seems desirable to investigate the mechanism of the refinement and the mechanisms of the coarsening reactions in order to establish control conditions for use in commercial production. The identity of the nucleating phase must first be established and then the factors affecting its chemistry and physical dispersion must be determined. THE IDENTITY OF THE NUCLEATING PHASE Davis, Eastwood, and DeHaven suggested that the nucleating phase in this system is Al4c3,1 but Mahoney, Tarr, and LeGrand8 disagree, largely because they found no evidence of the compound in alloys after carbon treatment and because there is no indication that aluminum carbide should be unstable over the temperature range used in the superheat treatment. This latter objection is based on the assumption that both the carbon treatment and the superheat treatment introduce the same nuclei. Electron diffraction studies have been made to identify the nucleating phase. Samples of grain refined A292 have been selectively etched SO that clean surfaces are obtained and so that secondary phases are in relief. Electron diffraction patterns from these surfaces have established that the carbon treatment of A292 introduces into the metal a large number of small, plate-like particles with a structure very similar to Al4C3. In most cases, the plate-like nature of the particles prevented positive identification but in the cases where the identification could be made the particles proved to be AIN A14C3. However, enough variation in lattice constants was observed so that all compositions from pure A14C3 to the 50:50 solid solution A1N.Al4C3 were probably present.14 In A14C3 and especially AlN.Al4C3 the A1 atoms occur in layers within which they have the same hexagonal symmetry and spacing as the Mg atoms in a single basal plane of a magnesium crystal. The solid solution spacing lies between the 3.16 of AIN and the 3.3? for Al4C3, in satisfactory agree-

Jan 1, 1962

By Marius S. Darrow, William B. White, Rustum Roy

The PbS-PbTe systen has been studied by quench-ing and D.T.A. techniques f?om 400&apos; to 1150°C. Runs were made in evacuated silica tubes so that all equilibria are at the vapor pressure of the system. Lattice parameters of the quenched salnples , measured by X-ray diffraction, show a complete crystalline-solution series existing over a narrow temperature range between approximately 805" and 871°C. An exsolution dome extends from a maximum of about 805"C (approximately 30 mole pct PbTe) to 1 and 96.5 pet PbTe at 400°C. A narrow melting region, deternined by D.T.A., extends form 918c (mp PbTe), The shapes of the liquides and solidus curves imply the existence of a minimum at 871°C at approximately 65 pct PbTe. THe exact composition of the minimum could not be established due to the very narrow two-phase region. At compositions containing less than 50 pet PbTe, liquidus temperatures begin to increase, while the solidus remains almost flat to about 15 mole pet PbTe before beginning to vise toward the mp of PbS (1075 C). LEAD sulfide and lead telluride are isostructural (NaC1 type) semiconductors whose electrical and optical properties have been extensively studied and used in recent years. If appreciable crystalline solution exists between these compounds, the variation of physical properties with composition could be of interest. The purpose of this investigation was to determine the extent, if any. of crystalline solution, and to obtain the phase diagram for the system. To the knowledge of the authors, only three studies of the system PbS-PbTe have been reported, and, in chronological order, each investigation found an increasing amount of crystalline solution. In 1956, Yamamoto reported finding no evidence of crystalline solution between the compounds. Sindeyeva and Godov-ikov,&apos; in 1959, found very limited crystalline solution. but only under conditions of excess tellurium concentration. Finally Melevski s3 investigation in 1963 indicated that one solid phase exists in the region from PbS to 7 pct PbTe and from 82 pct PbTe to PbTe at 886&apos;C, with an eutectic at 55 pct PbTe at that temperature. Detailed data on the solvus boundary were not given. EXPERIMENTAL EQUIPMENT AND MATERIALS Commercially produced PbTe and PbS powders were used as starting materials. Batches of specific mole percent composition were accurately weighed and mixed in a plastic bottle, in a shaker mill. An analy- sis of impurity content is given in Table I for pure PbS and PbTe and for two randomly selected batches after the powders were mixed. Individual samples, ranging in weight from 0.2 to 0.5 g, were sealed in evacuated silica tubes which had been thoroughly washed and rinsed with acetone and distilled water. Thus all data taken were at the pressure of the system. Subsolidus relations were studied down to 400°C by heating the samples in a vertical tube furnace for 24 hr. The sealed tubes were quenched in water with quench time from the hot zone not exceeding 1 sec. Temperatures were measured by a chromel-alumel thermocouple and controlled to 53°C for most runs. The number and composition of phases present were determined from powder X-ray diffraction patterns taken at room temperature on a Norelco diffractome-ter, using silicon as an external standard. Above 850°C quenching techniques were, in general, found to be unsatisfactory, and differential thermal analysis (D.T.A.) was used to determine melting relations. The evacuated tubes were recessed about 1 cm at one end to accommodate the differential thermocouple. Al203 was used as the reference material in a similar tube containing the other side of the differential couple. For temperature measurements, a separate thermocouple was placed in the recess of the tube containing the sample to be measured, thus providing an opportunity to obtain thermal, as well as differential, analysis. All thermocouples for these measurements were Pt-Pt 10 pct Rh. Temperature and differential curves were recorded separately on synchronized strip-chart recorders. Thermocouples and recording equipment were calibrated using NaCl and gold standards, using the melting points 801" and 1063 C, respectively, which span most of the temperature range of interest. Heating and cooling rates generally were from 4 to 7°C per min. It was found, in fact. that rates ranging from 1.5 to 25°C per min did not significantly change the data obtained.

Jan 1, 1967

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner&apos; makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,&apos; is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors&apos;&apos; derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel&apos; terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones&apos;" initiates steady-state

Jan 1, 1967

By J. L. Lehman, J. D. Sudbury, J. E. Landers, W. D. Greathouse

A full scale field experiment on cathodic protection of casing answers questions concerning (1) the proper criteria for determining current requirments, (2) the amount of protection provided by different currents, and (3) the transfer of current at the base of the surface pipe. Three dry holes in the Trico pool in Rooks County, Kans., were selected for cathodic protection tests. The three holes were in an area where casing failures opposite the Dakota water sand often accur in less than a year. Examination of the electric togs showed the wells to be similar to other wells in the field where casing in four of seven producing wells has failed. The three holes were cleaned out and cased with 75 joints of new 51/2-in. 14-tb J-55. Each joint was visually inspected and marked before it as run. The casing was bull plugged and floated in the hole 50 that the inside might remain dry and free of excessive attack. Also, if a leak occurred, a pressure increase could be observed on gawge at the surface. Extensive testing was done, including potential profiles, log current-potentid curves and electrode measurements from both surface and downhole connections. Based on these data, a current of 12 amps was applied to one well and 4 amps to mother. The third well was left to corrode. During the two-year period when the casing was in the ground, [he applied current was checked weekly, and reference electrode measurements were made about every two months. Three sets of casing potential profi1e.c were run. When the three strings were pulled, each joint was examined for type of scale formed, presence of sulfate-reducing bacteria, extent of corrosion nttnck and pit depth. Since the pipe was new when run, quantitative determination of the protection provided by current was possible. This is the first concrete field evidence to help resolve the many arguments about the proper method for selecting adequate current for cathodic protection of oilwell (-using. INTRODUCTION A casing string is run when a well is drilled. This pipe is supposed to protect this valuable "hole in the ground" for the life of the well. Often the casing does not last the life of the well; it is with these casing failures that this work is concerned. The cost of repairing a casing failure varies from field to field—from as much as a $30,000 per leak average in California to $5,000 per leak in Kansas. Additional costs other than actual repairs are also important. These include formation damage, lost production, etc. Casing damage caused by internal corrosion is important in some areas. Treatment normally consists of flushing inhibitor down the annulus, but further research is being done on control measures. The test described in this paper is concerned only with external corrosion. The problem of casing failure from external attack has appeared in several areas including western Kansas, California, Montana, Wyoming, Texas, Arkansas and Mississippi. Cathodic protection is currently being used in an attempt to control external corrosion. From reports in the NACE there are thousands of wells currently under cathodic protection. The quantity of current being applied ranges from 27 amps on some deep California wells to a few tenths of an amp being supplied from magnesium anodes on wells in Texas and Kansas. Considerable field and laboratory effort1,9,5,6 was exented on the problem of cathodic prctection of casing, and it became fairly obvious that this method could be used to protect wells. Early workers showed that current applied to a well distributed itself over the length of the casing and was not concentrated on the upper few hundred feet. Basic cathodic protection theory had shown that corrosion attack could be stopped by applying sufficient current. The problem resolved itself, then, into one of trying to decide just how much current was necessary. Various criteria were utilized in installing the many existing cathodic protection installations. These methods included the following. 1. Applying sufficient current to remove the anodic slope as shown by the potential profile." 7. Applying enough current to maintain all areas of the casing at a pipe-to-soil potential of .85 v.&apos; 3. Applying the current indicated by a log current-potential (or E log I) curve." 4. Supplying the current necessary to shift the pipe to-soil potential .3 v." 5. Applying 2 or 3 milliamps of current per sq ft of casing."

By E. J. Rapperport, E. D. Levine

The boundaries of the (a +ß) field in the Al-Li system were determined between 150°and 550°C utilizing quantitative metallography and lattice-parameter measurements. The solubility of lithium in aluminum decreases from 12.0at. pct Li at 550°C to 5.5 at. pct Li at 150°C. P Al-Li is saturated with aluminum at 45.8 at. pct Li and has this boundary value constant over the temperature range 150°to 550°C. THE solid solubility of lithium in aluminum has been determined by several investigators, 1-6 but, as shown in Fig. 1, there is little agreement among the various determinations. The earliest investiga-tions&apos;-&apos; are suspect because of the use of impure materials. Although high-purity materials were employed in more recent work,4&apos;5 the experimental techniques may have led to contamination of the specimens. Probably the best work has been that of Costas and Marshall,6 who obtained close agreement between results obtained by two independent phase-boundary techniques: electrical resistivity and mi-crohardness. No detailed studies of the solubility of aluminum in the bcc ß phase, Al-Li, have been reported. Cursory investigations1,2,6 have indicated only that the (a+ß) -p boundary lies between 40 and 50 at. pct Li and is relatively independent of temperature. The present work was undertaken in order to provide an independent check on Costas and Marshall&apos;s determination of the solubility of lithium in aluminum, to extend knowledge of this solubility limit to temperatures below 225°C, and to make an accurate determination of the solubility of aluminum in Al-Li. EXPEFUMENTAL Alloy Preparation. In view of the difficulties encountered in previous investigations of the A1-Li system, close attention was paid to the use of methods of alloy preparation and treatment that would minimize contamination. Aluminum sheet (99.99 + pct Al) was vacuum-induction melted in a beryllia crucible to remove hydrogen. Lithium (99.9 pct Li) was charged with pre-melted aluminum into a beryllia crucible, in a helium-filled drybox. The crucible was sealed in a Vycor tube and transferred from the drybox to an induction furnace. Melting of alloys was performed by induction heating in a helium atmosphere. Solidification was accomplished by means of a suction apparatus, shown in Fig. 2, in which the alloy was forced by changes of pressure into a 3/16-in. inside diam closed-end beryllia tube. This technique produced rapid solidification of a small portion of the melt, resulting in alloys with a high degree of homogeneity. Typical lithium distributions are presented in Table I. Transverse sections 1/8 in. long were cut from the alloy rods, and each section was split in half longitudinally. One half of each section was analyzed for lithium, and the opposing halves were employed for phase-boundary determinations. Lithium contents were determined by flame photometry with an accuracy of 1 pct of the amount of lithium present. Thermal Treatments. Homogenization and equilibration heat treatments were performed in electrical-resistance furnaces with temperatures controlled to ± 2OC. Calibrated chromel-alumel thermocouples were employed to measure temperature. Homogenization was performed in helium-filled l?yrex tubes for 1 hr at 565°C. The encapsulated specimens were then transferred directly to furnaces maintained at lower temperatures for equilibration. Equilibration times were 2 hr at 550°C, 8 hr at 450°C, 27 hr at 350°c, 90 hr at 250°c, and 285 hr at 150"~. These times were chosen on the basis of conditions employed by previous investigators. Alloys were quenched from the equilibration temperatures by breaking the capsules into a silicone oil bath. By performing all possible operations either in sealed capsules or in a helium-filled drybox, the specimens were given minimum exposure to the atmosphere. Quantitative Metallography. Metallography of Al-Li alloys is difficult because of the atmospheric reactivity of the ß phase. It was found possible, however, to prepare surfaces of good metallographic quality by preventing contact with moisture during preparation. Grinding through 4/0 paper was performed in the drybox. The specimens were then transferred under kerosene to the polishing wheel. Three polishing stages were employed: 25-p alundum with kerosene lubricant on billiard cloth, 1-µ diamond paste on Microcloth, and 1/4-p diamond paste on Microcloth. Between stages the samples were cleaned by rinsing in trichloroethylene and buffing

Jan 1, 1963

By R. G. Aspden

Crystals with an (001) [loo] initial orientation of an iron-base alloy containing 3 pct Si were cold rolled with and without the use of constraints. A major difference in the rolling and annealing textures was observed between crystals rolled with and without constraints. These data show that the contribution of constraints at grain boundaries in a poly crystalline sheet should be considered in applying textural data on single crystals to grains in an aggregate. SILICON-iron alloys with a cube texture have been recently developed and their magnetic characteristics reported.1-4 Of interest in the development of this texture were the textural changes of single crystals accompanying rolling and annealing and the influence of constraints at grain boundaries in an aggregate on the behavior of individual grains. The present study was primarily concerned with the effect of constraints during rolling on the textures of 3 pct Si-Fe crystals initially (001)[100]. Barrett and Levenson5 were among the first to observe an influence of constraints at grain boundaries on the textural changes of individual grains during deformation. They tested Taylor&apos;s6 theory of plastic deformation of face-centered-cubic metals in which deformation textures were predicted. About one-third of the grains in poly crystalline aluminum did not rotate as predicted. Grains of the same initial orientation were observed to rotate in different directions under the influence of applied stress and anisotropic flow of neighboring grains. Recently, the various inhomogeneities of flow of crystals in an aggregate have been studied7&apos;8 and reviewed.9-11 Barrett and Levenson" rolled (001) [loo] iron single crystals inserted in close-fitting holes in copper to limit lateral flow and to simulate rolling of grains in an aggregate. Deformation bands were formed after a 90 pct reduction in thickness, and the cold-rolling texture contained two components described by rotating the (001)[100] about 35 deg in both directions around the normal of the rolling plane. No annealing textures were reported. Chen and Maddin13 rolled molybdenum single crystals initially (001) [loo]. The crystals were mounted between two hardened silicon-iron plates and 96 pct reduced in thickness by rolling at a low rate of reduction, about 0.0001 in. per pass. The deformation texture had the mean orientation of (001) [loo], and the azimuthal spread included orientations described by rotating (001) [loo] about 35 deg in both directions about the pole of the rolling plane. The presence of deformation bands were not reported by Chen and Maddin or detected in subsequent work of Ujiiye and Maddin.14 The ideal orientation of the annealing texture was (001) [loo]. Recently, Walter and Hibbard 15 reported on the textures of 3 pct Si-Fe alloy crystals initially near (001) [loo]. Each crystal was in an aggregate cut from a columnar ingot. After 66 pct reduction by rolling, the texture consisted of two symmetrical components which had the orientations described by rotating (001) [loo] about 30 deg in both directions about the pole of the rolling plane. Annealing texture was near (001) [loo]. In the above work, the textures of body-centered-cubic crystals were studied after rolling under the influence of constraints. The deformation textures varied from (001) [loo] to near the (001) [110] type and appeared sensitive to the manner in which the crystals were rolled. No textural data were available on the effect of rolling (001) [loo] crystals with and without constraints. The purpose of the present work was to evaluate the influence of constraints during rolling on the textures of 3 pct Si-Fe crystals initially (001) [loo]. Rolling and annealing textures were studied for a) crystals rolled with no constraints at different rates of reduction, and b) crystals rolled with constraints imposed by neighboring grains and by plates between which a crystal was "sandwiched". PROCEDURES AND EXPERIMENTAL TECHNIQUES Data are presented on four crystals which are representative of several crystals studied. The orientation of each crystal prior to rolling was (001) [loo] as determined by the Laue X-ray back-reflection method," i.e., each crystal had an (001) within 3 deg of the rolling plane and [100] within 3 deg of the rolling direction. These crystals were obtained from two iron-base alloys containing 3 pct Si by weight which were prepared by vacuum melting electrolytic iron and a commercial grade of silicon. Crystals 1, 2, and S-1 were cut from a large single crystal grown from the melt of one alloy by the Bridgman technique17 in an apparatus described by

Jan 1, 1960

By J. R. Elenbaas, J. A. Vary, D. L. Katz

The development of gas fields, oil fields and aquifers for storing natural gas is treated from two main vieu.-points: (I) the volumetric storage capacity for gas in a given situation and (2) the prediction of the number of wells required for the delivery of gas. Other experiences in the design and operation of storagc fields are incluclerl INTRODUC TION Storage of natural gas in underground reservoirs near the terminus of long distance pipelines has been the prime factor in opening the space heating market to the natural gas industry. Storagc has permitted a major. increase in both the load and the load factor of pipe-lines; some are now operating at steady load throughout the year. Thus, underground storage has been responsible for the rapid increase in demand for natural gas in recent years. Three types of reservoirs have been used for gas storage: natural gas reservoirs, oil reservoirs, and waterbearing sands or aquifers. This paper presents the factors to be considered when developing gas storage reservoirs of these three categories. There are two prime considerations tor any storage reservoir: (1) the volume of gas which a given reservoir will store advantageously and (2) the number 01 wells needed to provide the required peak deliverability. These two problems will be considered for the three types of reservoirs just noted STORAGE IN PARTIALLY DEPLETED GAS FIELDS Early storage operations consisted of replenishing the natural gas in a depleted gas field situated adjacent to the market. Today, newly discovered fields near the market may be considered for storage, and this discussion applies equally to both types of reservoirs. For reservoirs originally containing gas or oil, the question of the impermeability of the cap rock nor-mally does not arise. However, such fields are likely to have many wells drilled either to or through the reservoir under consideration. Positive assurance must be obtained that such wells are or can be made mechanically tight. Corroded casings may need to be lined or permanently plugged. Abandoned wells should bc reopened and properly cemented. The volumetric capacity for gas storage depends upon space available in the porous rock as well as pressure and temperature of the gas in the reservoir. The production-pressure decline data on partially depleted gas reservoirs without water drive permit calculation of the reservoir space for gas. Isopachous maps of sand volume and porosity data for the reservoir rock provide an alternate method of calculating the pore volume for water-drive reservoirs. The pressure range selected for the storage cycle depends upon ()) the safe upper limit of pressure. 2) the flow capacity of wells and (3) compression requirements when injecting gas into the reservoir or delivering to market. Normally, gas and oil fields have pressures at discovery in the range of 0.43 to 0.52 psi/ft of depth. Pressures of around 1.0 to 1.2 psi/ft of depth appear to lift the overburden1-3 and invite uncontrolled movement of fluids in the porous rock. Some top pressure is normally selected for a storage reservoir ranging from below discovery pressure for deeper reservoirs to 0.65 psi/ft of depth for shallower reservoirs. Pressures to 0.66 psi/ft have been experienced without difficulty. The lower pressure limit is set by water intrusion accompanying low pressures, reduced flow capacity for wells at lower pressures and compression requirements. Depletion-type gas reservoirs often encounter water problems in the later stages of gas production. Such water intrusion may be due to movement from the surrounding aquifer. Accordingly, displacement of this water back into the aquifer by gas pressure and subsequent surges of water corresponding to the gas storage pressure cycle must be considered. Storage fields often produce in four months a volume of gas equal to its initial content. Rapid decreases in reservoir pressure occur, such as 20 psi/day. Accordingly, closed-in pressure observation wells which reflect the pressure in the bulk of the reservoir are required for following the operation of the reservoir. It has been found that a plot of observation wellhead pressures against gas content, Fig. 1, is very useful in observing operation of the field, checking the inventory and predicting future behavior. The plot is based on a given quantity of base or cushion gas in place. The injection and withdrawal curves may spread depending upon the homogeneity of the reservoir rock. permeability of the rock, well spacing and flow rates.