Search Documents

By W. C. Winegard, J. R. Carruthers, A. Plumtree, L. R. Morris

The unidirectional solidification of dendrites containing central twin planes was studied in A1-Ti alloys. Once nucleated, the twinned dendrites are a Twore ejficient form for solute redistribution and therefore grow in preference to the normal columnar dendvites. Comparison of these twinned dendrites to adjacent colunmar dendrites by means of decanting experinzents and electron-probe rnicvoarmlysis indicates that these special dendrites grow with less undercooling than normal dendrites. These findings are further supported by the effect of forced convection on the dendrite morphologies. COMMERCIAL semicontinuous cast ingots of most aluminum alloys frequently exhibit large grains which appear to be composed of hundreds of parallel, continuous, thin lamellae. This structure has been termed "basaltic",' "fiederkristall", or, commonly, "feathery grain". The lamellae are 'about 100 p thick, several inches long, and each lamella contains a central (111) coherent twin boundary. The feathery grain has been reported to have a (112) direction2 and a (110) direction4 in the twin plane parallel to the casting direction, in contrast with the usual columnar structure where a (100) direction predominates. Aust et uZ.~ proposed that the twin boundaries were growth twins nucleated by stacking faults on the octahedral planes. chalmers6 has suggested that feathery grain may grow by a re-entrant edge mechanism, as proposed by wagner7 for twinned dendritic growth in germanium. Cahn et ~1.~ have concluded that the occurrence of feathery grain is evidence of some form of lateral layer growth rather than the atomically continuous growth normally observed in metals. The postulates by Chalmers and Cahn would seem to be contradicted by the work of Nakao~"' who showed that feathery grain only occurs when growth rates are high and the aluminum contains some solute. Specifically, Nakao found that in order to obtain feathery grain, in small castings solidified unidirec-tionally upward, the rate of growth must be above 2.4 cm per min and a critical solute concentration must be present. Below this solute concentration the grain structure was totally columnar. The critical solute concentration was found to be approximately: 0.04 wt pct Ti, 2 wt pct Cu or Mg, and 8 wt pct Zn. As pointed out by Chalmers it is not obvious why a twin-plane re-entrant edge mechanism would occur in aluminum which is thought to have a diffuse solid-liquid interface. The present experiments were undertaken to determine the growth mechanisms in- volved and to study the solute segregation in more detail. EXPERIMENTAL PROCEDURES Alloys ranging in composition from 0.05 to 0.23 wt pct Ti* were prepared from 99.9 pct pure Ti and 99.993 pct pure Al. Two-hundred-gram samples of A1-Ti binary alloys were solidified unidirectionally and vertically upward from a water-cooled copper base, in a heated insulated mold. The 1-in.-diam, 6-in.-length mold was made of Marinite (manufactured by Johns-Manville Co.). The mold was attached to a 24-in. pivoted arm such that by dropping a weight the mold was rotated 180 deg, throwing the liquid metal from the solid. In this way, the solid-liquid interface was revealed by decanting. A sketch of the decanting mold is shown in Fig. 1. The alloy was poured into the mold at temperatures ranging from 680" to 750°C, partially solidified by water cooling from the base, and decanted after a measured time interval. Growth rates for each metal-pouring temperature were calculated from solidification-time vs length solidified curves. Temperature gradients in the melt were measured using four No. 34 gage thermocouples which protruded into the mold cavity. Grain orientations were determined by X-ray diffraction using Laue back-reflection techniques. Grain substructures were examined metallographically, using polarized light, by applying a thin epitaxial anodic film to polished sections after the method of Hone and pearson." Titanium micro segregation was measured by electron-probe microanalysis using a NORELCO AMR/~ with a mica crystal and proportional-flow counter. Several of the cast samples exhibited feathery and columnar dendrites growing in the same direction and

Jan 1, 1967

By K. H. Kilgren

At high formation pressures the distillate produced from a gas-condensate reservoir may be black in color. In this event the dense gas phase existing above the dew point is correspondingly dark. Volumetric phase data and an analysis of a reservoir fluid exhibiting these characteristics, together with a description of the visual equilibrium cell in which these observations were made, are presented in this paper. INTRODUCTION Previously the author, like many others in the oil and gas industry perhaps, tacitly assumed that the expressions black or dark oil system, crude oil system and bubble-point system were synonymous. Crude oil reservoir fluids are bubble-point systems and yield a black or dark stock-tank oil of relatively low API gravity. Conversely, a clear or amber colored trap distillate of high API gravity is assumed indicative of a dew point system or a gas-condensate reservoir fluid. This broad classification appears satisfactory for shallow reservoirs, but as the following study demonstrates, may be misleading when applied to deep reservoirs. Theoretically, there is no reason to exclude the possibility of producing a dark, low-gravity distillate from a gas-condensate reservoir. At sufficiently high values of pressure and temperature, heavy, dark-colored hydrocarbons may exist in the vapor state of a multi-component system. If enough dark-colored components are present in the reservoir vapor phase, the resulting condensate will be dark. The reservoir fluid investigated in the present study supports this contention. Stock-tank production was black in color and measured 29" APT gravity. From outward appearances, the liquid closely resembled a medium gravity crude oil. Experimental measurements proved the reservoir fluid was in reality a gas-condensate system. Volumetric phase data for this high-pressure system and a description of the visual cell in which the study was conducted successfully are presented in this paper. THEORY Phase behavior of a reservoir fluid can be predicted accurately with reference to a pressure-temperature phase diagram. If the reservoir temperature is lower than the critical temperature of the hydrocarbon fluid in place, bubble-point behavior will be observed. If the reservoir temperature lies between the critical and cricondentherm temperature, dew point behavior and retrograde condensation will occur. For reservoir temperatures above the cri- condentherm, only a single gas phase can exist in the reservoir regardless of pressure. Providing the composition of the reservoir fluid were known, it would be possible to predict the critical temperature and estimate the phase behavior from equilibrium relationships. However, the usual practice is to obtain a sample of reservoir fluid, subject it to varying pressures at the reservoir temperature and observe the phase behavior experimentally. The latter method was used to obtain the data reported here. WELL AND TRAPPING INFORMATION A summary of pertinent data for the well from which the reservoir fluid was sampled is presented in Table 1. This well is located offshore Louisiana. Except for the pressure which substantially exceeds hydrostatic pressure, the information does not appear unusual. Prior to the sampling program, the well was produced for 22 hours at a stock tank oil rate of 139 B/D. Average trapping conditions and gauging data for the six-hour test period that followed are summarized in Table 2. Samples of the first-stage trap gas and liquid were obtained during the latter portion of the test period. Ambient temperature remained 5 to 10F below trap temperature and presented no problem for sampling. Surface wind and moderate foaming of stock-tank oil presented some difficulty in obtaining accurate stock-tank gauges. SAMPLE ANALYSIS Compositions of the gas and liquid samples are shown in Table 3. The trap gas was analyzed by isothermal chromatography which revealed only a trace of heptane in the stream. The trap liquid was initially analyzed by low-temperature fractional distillation, yielding a bottom product of heptane and heavier components. Specific gravity of this fraction was measured and the mol weight was determined by freezing point depression. The hexane and lighter overhead gas collected during distillation was

Jan 1, 1967

By Carl F. Lutz, Robert F. Mehl

A layer of brass was formed on 0 brass using a vapor-solid reaction technique. The variation in composition with distance within the phase layer and across the a -ß interface was determined by an electron probe. The diffusivities were calculated as a function of concentration; the diffusion coefficient in brass was found to change by a factor of twenty-five over a compositional range of 8 at. pct. An explanation is proposed to account for the large change in the diffusivity, based on possible stresses in the diffusion zone; it is suggested that the concentration dependence of the thermodynamic factor might decrease the activation energy with increasing zinc content. ALTHOUGH the body of data on diffusion coefficients (D- values) in terminal solid solutions and in isomorphous systems is quite large, there is but little quantitative information for intermediate alloy phases. This paper recounts measurements of D-values for the ? phase in the Cu-Zn system. The experimental method employed consists in the forming of alloy layers and determining and analyzing the concentration-distance (c-x) curve. The experimental method employed could have been applied to several systems; the Cu-Zn system was chosen because the D-values in the a and ß phases are well known, and can thus be used for purposes of comparison. The rates of growth, of intermediate alloy layers are known for a number of systems. In nearly all cases the thickness varies parabolically with time, as to be expected if the rate of growth is controlled by the D-values. In a few cases, nonparabolic behavior has been noted.1 Nonparabolic growth may be the result of a) a variation of the composition at the phase interface with time, b) interface-controlled kinetics, c) short diffusion times coupled with long vacancy lifetimes (in vacancy diffusion),2 and d) crystal reorientation during diffusion in anisotropic systems.3 In the Al-Ni system4 and in the Al-u5 system parabolic growth kinetics are slowly approached after an initial transient period. In general those phases stable at a given temperature in a system A-B appear as separate phase layers when A and B are put in contact and given a diffusion he at-treatment. In most cases the compositions at the interface of adjacent phase layers are those read from the phase diagrams for the termini of the respective phases. This discontinuity in concentration is taken as independent of time, and the growth of one phase at the expense of another is assumed to be independent of interface kinetics, i.e., the rate of interface movement is controlled by volume diffusion in the phases concerned. wagner6 has given many solutions for multi-phase diffusion processes, assuming that the chemical diffusion coefficient, D, is independent of concentration, as have others. Jost7 has pointed out that the familiar Boltzmann-Matano solution is as valid for a (c-x) curve exhibiting concentration discontinuities at phase boundaries as for the usual single-phase case. Heumann and associates8,9 have applied this solution to the multi-phase case, but lacking concentration data within the several layers were forced to assume linear concentration gradients, thus calculating only average D-values. The purposes of the present study were: 1) to determine the dependence of D on concentration, 2) to calculate the intrinsic diffusion coefficients D?Cu and D?Zn, 3) to establish the time-law for the movement of the interface, 4) to determine the concentration limits at the ß-? interface, and, 5) to attempt to clarify the mechanism of diffusion. EXPERIMENTAL PROCEDURE The ? phase is exceedingly brittle; conventional solid couples and conventional sectioning methods were thus inapplicable. Vapor-solid couples were accordingly employed. Unfortunately the equilibrium vapor pressures of Zn for the ? phase are unknown; to escape this awkwardness, samples were enclosed in a capsule containing powder of an alloy of known composition, sufficiently large in quantity as to be effectively an infinite source, maintaining the concentration of Zn at the sample surface constant with time. The sources of Zn-vapor which can be employed in reaction with Cu, or a brass, or ß brass, to form a layer of ? brass satisfactorily wide in concentration range, are alloys of the phases ?, or ? + E or ? + d (depending on the temperature).

Jan 1, 1962

By J. B. Cohen, B. W. Howlett, M. B. Bever

The partial molar heats of solution at infinite dilution in tin of aluminum at 300° and 350°C and of gallium, indium, and thallium at 240°, 300°, and 350°C have been measured by tin solution calori-metry. Aluminum, gallium, and thallium are endo-thermic on solution; indium is exothermic. Any temperature dependence of the heats of solution lies within the experimental scatter. Over the dilute ranges investigated, only aluminum has a measurable change in its heat of solution with composition. HEATS of solution of one element in another reflect the interaction between them. The investigation of partial molar heats of solution in dilute alloys is of particular interest as the properties of the solvent are altered to only a limited extent by the presence of the solute and also as the interaction between solute atoms is small. When the heats of solution of a related series of elements in a solvent are known, a systematic comparison may be made. In the investigation reported here, the partial molar heats of solution of the Group III elements aluminum, gallium, indium, and thallium in dilute solution in tin were measured. This work follows an investigation of the heats of solution in tin of the Group IB elements.' EXPERIMENTAL PROCEDURES Materials. Samples of gallium, indium, and thallium were obtained from Johnson, Mathey and Co., Ltd. Indium and thallium were supplied as wire, 1.6 mm in diam; gallium was in the form of irregular pieces. The supplier reported the following minimum purities: gallium—99.95 pct; indium—99.99 pct; thallium—99.99 pct. The aluminum, obtained from Alcoa Research Laboratories, was reported to be 99.995 pct pure. The tin was supplied by Baker and Co., Inc.; the reported analysis indicated a tin content of at least 99.96 pct with lead as principal impurity. Calorimeter. this description will cover only the essential features of the calorimeter with special attention to modifications made since an earlier description was published.' A Dewar flask containing the tin bath was held in a constant-temperature bath of a near-eutectic mixture of lithium, sodium, and potassium nitrates. This bath, which was stirred vigorously, was heated by a primary resistance winding in the container wall and by a secondary winding immersed in the salt. The voltage supplied to both windings was stabilized. The temperature of the salt bath was controlled by means of a platinum resistance thermometer in one arm of a Wheatstone bridge. The light from a mirror galvanometer in the bridge circuit fell on a photocell which controlled the current in a saturable reactor in series with the secondary winding. In this manner, the temperature of the salt bath was controlled to ±0.003°C and that of the tin bath to at least ±0.002°C. Each of these temperatures was measured by two iron-constantan thermocouples in series, coiled in a helix to minimize heat loss and immersed in the salt and tin baths in protective sheaths. The temperature of the laboratory was kept constant to ± 1°C during a run. Specimens were dropped into the tin bath from an addition arm held at O°C which was part of the cal-orimetric system. The system was evacuated to about 0.02 1 to minimize oxidation and to reduce transfer of heat. The bath was stirred by a glass stirrer introduced through a double Wilson seal. The samples were scraped clean before weighing, which was carried out as rapidly as possible. Each sample was immediately placed in the evacuated addition arm to minimize contamination. These precautions were especially necessary with aluminum. After the runs with gallium and thallium at all temperatures and with indium at 240° and 350°C (Series I) were completed and before the runs with indium at 300°C and aluminum at 300° and 350°C (Series 11) were begun the following changes were made. The shape of the Dewar flask was changed so as to result in a lower surface to volume ratio of the bath and at the same time the amount of tin was reduced from 500 to 400 g. The paddle type stirrer was replaced by a helical screw and the rate of stirring was increased to about 150 to 200 rpm.

Jan 1, 1962

By S. L. Lopata, E. B. Kula

A variation of yield and tensile strength with direction has been noted in heat-treated 4340 steel which had been warn-worked by rolling in the austenitic condition prior to quenching. Measurements by X-rays showed a preferred orientation of martensite crystals. Various ideal preferred orientations of martensite were calculated by assuming a preferred orientation for austenite and transforming this to martensite by the Kurdjumov -Sachs and Nishiyama relationships. Comparison with the measured preferred orientation showed that the ideal rolling texture for austenite in 4340 steel could be expressed as a (112) [111] texture with a component of (110) [112]. No evidence was found for any fiber axis in the rolling direction, or fop- a (123) [412] texture. The observed martensite pole figures can be adequately accounted for by these ideal austenite textures and by the Kurdjumou-Sachs or Nishiyama transformation relationships between austenite and martensite. ALTHOUGH preferred orientations are found in cold-rolled steel sheet and in many nonferrous metals, they are not generally encountered in wrought, heat-treated high-strength steels, since such steels are generally not cold worked sufficiently to develop a strong texture, and any texture that is present is removed during heat treatment by the phase transformations taking place. This report gives details on a preferred orientation found in warm-rolled and hardened 4340 steel. In a recent study, a special combination heat-treating and working technique was used in which 4340 steel was deformed in the austenitic condition, and before re crystallization of the austenite could occur, quenched to form martensite.' Since the austenite was unrecrystallized, it was cold worked and hence could have contained a preferred orientation which could be transmitted to the martensite on quenching, and this in turn could be manifested as a directionality of mechanical properties. Such a directionality of properties was found as a result of working, and this prompted an X-ray study of the preferred orientation. Literature Survey— Rolling Texture in Face-Cen-tered-Cubic Austenite—M a preferred orientation is found in a rolled and hardened steel, it is first appropriate to consider what texture existed in the austenite as a result of the rolling but prior to the transformation to martensite. The changes taking place as a result of the austenite-martensite transformation are then considered. In view of the instability of austenite in 4340 steel at low temperatures, no direct determination of the austenite texture was made. However, an estimate of the texture for austenite in 4340 steel may be made by considering the textures found for other face-centered-cubic metals. Barrett' has reviewed the rolling textures for many materials. The most common orientation in face-centered-cubic metals has the (110) planes parallel to the rolling plane and the [112] directions parallel to the rolling direction, or (110) [112]. Other orientations found are (112) [111], as well as (100) [001], (110) [001], and others. Beck and his co-workers, using quantitative pole figures, report that (123) [121],3 later corrected to (123) [412],4 better describes the ideal orientation. smallman5 reported that the ideal preferred orientation for high-purity face-centered-cubic metals is (123) [121], but that there is a transition to (110) [112] as solid solution elements are added. The latter texture tends to occur more readily at lower temperatures. This transition between the two textures was explained by differences in the slip characteristics with concentration or temperature. More recently, Jones and Fell6 concluded that the (123) [412] texture for face-centered-cubic metals could equally well be represented by a mixture of (110) [112] and (112) [111] textures, and hence it had no physical reality. From the above results it can be seen that there is no agreement on the "ideal" texture for fcc metals and, hence, no firm basis for assuming a texture for austenite in 4340 steel. Therefore, the orientations (110) [112], (123) [412], and (112) [111] have been considered in the following.* *The (123) [4121 is identical to the (123) [121] rotated 90 deg. Orientation Relationship between Austenite and Martensite—When martensite forms from austenite, certain crystallographic planes of the martensite are parallel to planes in the austenite. This orientation relationship was first determined by Kurdjumov and

Jan 1, 1960

By F. E. Jaumot, R. L. Smith

Self-diffusion in zinc at temperatures below 200°C has been studied using both single crystal and polycrystal samples. Anomalous results were obtained for single crystal samples, the data indicating that in some cases grain boundary diffusion predominated. These anomalous results are presumed to be due to low angle lineage boundaries in the single crystals. When volume diffusion occurred in either single or polycrystal samples, the values for the diffusion coefficient were as much as several orders of magnitude larger than would be expected from high temperature data. ONSIDERABLE work has been done on self-^-*diffusion at temperatures in the neighborhood of two thirds of the melting temperature (in OK) and higher, but very little has been done at lower temperatures. This is particularly true of self-diffusion in single crystals. This paper reports work on self-diffusion in zinc at temperatures below 200 °C where both polycrystalline and single crystal samples have been used. The measurements were made using both the sectioning and absorption techniques. The experiments on diffusion at low temperatures were initiated primarily for three reasons: to obtain data from single crystal diffusion covering a wide range of temperatures, to obtain values of the activation energy for grain boundary diffusion, and to obtain, if possible, a lower limit to the temperatures at which the absorption technique may be used. It is possible to determine whether volume or grain boundary diffusion predominates at low temperatures by analyzing the manner in which the concentration of the diffusing atoms varies with the depth of penetration. The analysis is possible only for the sectioning technique. For volume diffusion and for the boundary conditions of the present experiment, the curve of the logarithm of the concentration, C, vs the square of the penetration depth, x, is a straight line with the slope proportional to — 1/D. If grain boundary diffusion predominates, Fisher' has shown that a curve of In C vs the first power of the penetration depth should be a straight line. This linear dependence of concentration on depth has been observed experimentally by several investigators. (See, for example, refs. 3 and 4.) Following Fisher's analysis, the grain boundary diffusion coefficient can be written as where 8 is the grain boundary width, D, is the volume diffusion coefficient at the temperature involved, t is the time of diffusion, C is the concentration of diffusing material, and x is the depth of penetration. A more general analysis employing the same model as Fisher's but with fewer assumptions has been given by Whipple.' In contrast to Fisher's analysis the more general analysis of Whipple does not always yield a straight line relation between In C and x. Recently Turnbull and Hoffman have discussed both analyses and have applied the Fisher-Whipple model to a more refined dislocation model of the grain boundary. They have also discussed the possibility of a large effect of grain boundary direction on diffusion. In the case of the absorption technique, for which results are also given, one simply gets a value for the diffusion coefficient with no indication of the predominating mechanism. Experimental Techniques Single crystals of 99.99+ pct Zn were grown using a modified Bridgman method. Polycrystal samples were prepared in three grain sizes (small, about 4 mm2; medium, about 9 mm'; and large, about 25 mmz) by quenching, air cooling, and furnace cooling. The grains of the polycrystalline samples tended to be columnar, with some preferred texture such that the c-axis was inclined at about 20" to 25" from the normal to the diffusion face, particularly in the medium and large grain size samples. The samples were cut from the rod, polished, and heavily etched, and faces cut accurately parallel using a jeweler's lathe. They were again etched to remove the lathe smear and annealed for 3 hr at 250" to 300°C. After annealing, they were etched, examined, and the orientation determined by Laue back reflection photographs. They were then electroplated with Zn. The plated surfaces were examined for uniformity; nonuniform coatings were removed by etching and the samples replated without additional preparation. For the diffusion anneals, the samples were placed in pairs, with active faces together, in evacuated Pyrex tubes. Constant temperature baths were used, and the variation in temperature did not exceed ±l°C, Diffusion runs were made at l00°,

Jan 1, 1957

By John M. Tinsley, Calvin D. Saunders, H. K. van Poollen

In the past few years much con-sideration has been given to the evaluation of the effect of hydraulic fracturing on the productivity of wells. Generally, these studies included the evaluation of fracturing materials, fracture extension and formation damage due to the use of various fracturing fluids. Only little consideration has been given to the characteristics, and in particular the flow capacity, of the fracture itself and its effect on well productivity. This paper presents the results of laboratory investigations pointed toward fie evaluation of the efficiency of various fracrures with special emphasis on the flow capacity of these fractures. Data presented in this paper are the results of both an electrical model study and physical testing. Under consideration are (I) effect of overflush, (2) premature production of well after treatment, (3) "tailing-in" with coarse sand near the end of the treatment, (4) effect of propping agent size and concentration, (5) reduction in effective frac-ture permeability caused by formation caused by formation fines, silt and clays, and (6) effect of various fluids on formation strength and competency. The results of this investigation indicate that the flow capacity of a fracture is affected by any or all of the various parameters mentioned above. The authors believe that a better understanding arid utilization of these factors should result in more efficient formation fracturing. INTRODUCTION Hydraulic fracturing has become almost a standard practice of many companies for stimulating production from old and new wells. Although most companies utilize this service, techniques of application vary widely between companies and areas. Probably too often when a well in an area responds favorably to a particular technique all future wells in the same area are treated in a similar manner. Possibly a modification of the technique would result in a further production increase. Variables, of which many are extremely difficult to evaluate from field results, hamper the selection of procedure changes. Attempts are being made by a number of organizations to analyze statistically treating techniques from production data. This is a very worthy and necessary approach but very possibly laboratory investigations may aid in evaluating some of the variables which tend to affect the results of a fracturing treatment. Some of the factors cannot be studied from practical field experience and only laboratory tests can show the possibilities which might exist. One of the factors which appears to be of major concern today is the flow capacity of the created fracture and how it can be changed. Papers on this general subject have been primarily concerned with the size of propping agent and the extent of fracture. Papers have also been written on the possible permeability damage to formations by fracturing fluids. In addition, it might be possible that another type of flow restriction is prevalent. This would be a restriction of flow through a sand-packed frac- ture caused by foreign materials integrating within the propping agent. This paper presents preliminary data obtained in an attempt to evaluate the effect of some factors affecting flow through sand-packed fractures. No attempt is being made to offer a fracturing technique adaptable to all areas and conditions, but to furnish data tending to show the possible effects which might be caused by variations in procedures and materials. PROPPING AGENT PERMEABLLITIES A hydraulically induced fracture containing sand as a propping agent may theoretically be classified as a packed-sand system. The flow of fluids through such packed systems has been the subject of much research. Although there have been numerous methods proposed for the evaluation of such systems, most writers agree on the general properties affecting their flow capacity or permeability. These properties include porosity, particle size, sphericity and the roughness of the particle. In some methods of evaluation . the particle size and sphericity terms are combined to produce an equation which is a function of the surface area of the particles. In this study the permeability of various fracturing sands was both measured and calculated. The apparatus for the permeability measurements consisted of a 52-in. Lucite tube with a 2.5-in. ID. A screen and drain plug were fitted in the bottom of this tube to retain and hold the sand in place while allowing fluid flow. Two pressure taps consisting of thin, highly perforated

By G. W. Tracy, R. D. Carter

Numerical calculations were made to determine the behavior of reservoirs with high-pressure drawdown and wide well spacing where the initial productivity is low and the wells are completed by hydraulic fracturing. The two-phase flow equations were solved for the flow into a single well. This well was assumed to be producing from a reservoir with hydraulically created horizontal fractures (four different systems with fractures were studied). For comparison purposes, additional two-phase flow calculations were made assuming a reservoir with uniform rock properties. The two-phase flow results were also compared with the conventional calculation methods, which do not include the effect of saturation gradients resulting from a simultaneous flow of oil and gas which are normal to this type reservoir. It was found that the conventional methods predict (1) a high and too optimistic value of ultimate recovery, (2) a high producing rate and a high reservoir pressure at a given oil recovery and (3) a low trend of gas-oil ratio with oil recovery. Included in the two-phase flow calculations were provisions to control the oil production rate by an allowable rate and, also, by a gas-oil ratio penalty rule. For the systems with hydraulic fractures, the producing rate was controlled by the gas-oil ratio penalty rule for most of the life. This is in contrast to the system with uniform rock properties which went "on decline" almost immediately. An unexpected characteristic of the systems which included fractures was the early rise in producing gas-oil ratio from 730 cu ft/bbl to approximately 1,200 cu ft/bbl, followed by a "leveling off" before the normally expected gas-oil ratio rise began. Additional features which are a result of hydraulic fracturing are (I) greater ultimate recovery, (2) higher average producing rates and (3) a lower average reservoir pressure at a given oil recovery. INTRODUCTION Some oil fields discovered during the past few years are producing from certain volumetric ally controlled reservoirs (often referred to as solution or internal gas- drive reservoirs) which are characterized by high-pressure drawdown at the wells. Since the available pressure drawdown at a well is limited by the static reservoir pressure and the producing rate is controlled by the available drawdown, wells completed in this type of reservoir usually produce at a rate less than the allowable from the time of completion. Because of this, this type of reservoir is referred to as a low productivity reservoir. Economic considerations require the use of wide well spacing and well stimulation by hydraulic fracturing to make commercial wells in this type of reservoir. Performance predictions for volumetrically controlled reservoirs have been made using a combination of two standard equations. 1. The "Schilthuis" or "Muskat" type material balance equation is used to relate the average reservoir pressure and the cumulative oil recovery. 2. The results from the material balance equation and the productivity factor as described by Pirson' are used to relate the cumulative recovery with producing rate and time. The material balance equation assumed uniform pressure and liquid saturation conditions throughout a reservoir. The steady-state radial flow formula allows for a pressure gradient toward a well but assumes uniform liquid saturation. These calculation methods are adequate for application to reservoirs wherein the drawdown at the well to realize satisfactory producing rates is small compared to the total pressure. In low productivity, volumetrically controlled reservoirs, the pressure drawdown at the well is large cornpared to the total pressure. Although a precise number cannot be given for the magnitude of a large pressure drawdown, values in excess of 1,000 psi would definitely be included. For practical considerations, this usually occurs when the formation flow capacity is less than about 100 md-ft. However, this limit of formation flow capacity will vary with the well producing rate. The low pressure in the neighborhood of the well which results from a high drawdown causes evolution of large volumes of gas. This causes the gas saturation to be higher near the well than at a greater distance—-hence, a non-uniform gas saturation. Also, the relationship between the relative permeability to oil (K/K) and gas saturation is nonlinear but decreases approximately in an exponential way with increases in gas saturation. Because of this, the following chain reaction is established.

By P. B. Baxendell

The performance of a water-drive reservoir having a gas cap depends primarily on the movement of the gar-oil and oil-water contacts. The movement of the contacts during production depends in turn on fluid withdrawals and how the reservoir pressure changes as fluids are produced from the reservoir; that is, on how effectively the aquifer maintains pressure by replacing withdrawals. inasmuch as pressure changes, fluid withdrawals, contact positions, and produced gas-oil and oil-water ratios are interdependent, the analysis and prediction of the performance of a reservoir produced in a given way must take into account this interdependence throughout depletion. This paper presents an analysis which, within the limitations of the assumptions made, yields an engineering approach to predicting future performance based on reservoir pressure and production history. The most significant assumption is the method of extrapolation of future gas-oil and water-oil ratios. The extrapolation procedure smoothly increases both ratios to preselected values as the remaining oil column undergoes a specified decrease in thickness. This preselection is made on the busis of previous field experience in depletion of similar reservoirs under similar conditions. For computing future performance a vdumebic balance L combined with the differential equation defining pressure distribution in the aquifer to obtain positions of water-oil and gas-oil contacts. From these positions are extrapolated produced water-oil and gas-oil ratios. Reservoir performance can be investigated when oil production rates are dependent upon various factors including the performance of the reservoir itself. Examples of practical application of the procedure are included. INTRODUCTION To predict performance of water-drive reservoirs with gas caps and thin oil columns, it is necessary to describe the motions of the fluids within the reservoir during the entire production period to depletion. These motions depend on the pressure changes and on the withdrawal of oil, gas and water. Production of oil from the oil zone primarily causes water to move in to take its place; production of gas from the cap tends to cause oil to migrate into its place. In addition, the volumes of oil and gas remaining in the reservoir depend on changes in the pressure, since a decrease in pressure causes fluid expansion, gas liberation and oil shrinkage. A method of relating the future pressure to total withdrawals is used to describe the motions of fluids within the reservoir under conditions arising in possible modes of production. The analysis is based on two relationships and will reduce the problem to one amenable to digital computation. The first of these concerns the dependence of the pressure distribution in the entire aquifer furnishing the water drive upon the total reservoir withdrawals. This dependence is dictated by the permeability distribution and the extent of the aquifer; at present, such information is most readily obtained from the performance history by means of the resistancecapacitance reservoir analyzer. The second relation involves withdrawals, pressure in the reservoir, and movement of oil, gas, and water within the reservoir. The analysis is subject to certain simplifying assumptions that are necessary to permit solution of the problem. A comparison of the methods of this paper with those presently practiced is pertinent. One method of analysis is to use the reservoir analyzer to predict reservoir behavior based on aquifer characteristics determined from production history. Inasmuch as the total withdrawal rate depends upon the gas-oil and water-oil ratios, which depend in turn upon, among other things, the positions of the gas-oil and water-oil contacts, these ratios must be assumed in advance. In this paper these ratios are, instead, related to the computed positions of gas-oil and water-oil contacts. Thus, our method can be applied to problems in which the oil production rate is limited by produced gas-oil ratio, or to problems in which it is desired to determine the variable gas injection rate that will maintain the gas-oil or water-oil contact stationary. These problems cannot be worked satisfactorily on the analyzer. Another method presently used is based on a paper by Hurst.' Since his procedure is dependent upon using the solution of the heat flow equation, which requires constant permeability within the aquifer, and our procedure recognizes variations of permeability in

By W. T. Weller

In Part I, a study of pressure build-up curves calculated for conditions under which both oil and gas flow led to the conclusion that the presence of a dispersed free gas phase in an oil reservoir must be taken into consideration to estimate accurately average reservoir pressure and permeability from build-up curves. Familiar methods based on the assumption of no free gas can be extended to the two-phase case by using total compressibility and mobility in place of single-fluid compressibility and mobility. These methods give correct values for average pressure and permeability when gas saturation is small. Errors become larger as the gas saturation increases. However, for the use to which the results will be put, the methods are satisfactory for reservoir engineering purposes. An improved method of calculating the performance of depletion-type reservoirs is presented in Part 2. Because the mathematical relationships describing simultaneous flow of oil and gas are quite involved. simplifying assumptions are made to provide means of obtaining approximate solutions of reasonable accuracy. One such approximate method now in use is the constant-GOR solution. It involves the assumption that at any instant, the ratio of total gas flow rate (both free and disolved) to oil flow rate is the same at all points in the reservoir. The approach is not applicable unless the free gar saturation in the reservoir is everywhere greater than the critical gas saturation. This paper presents a modification which, by avoiding the constant-GOR assumption, makes the method applicable to all reservoir conditions, and so far appears to be more accurate than the constant-GOR solution and to he comparable in required compuctation time. PART I— BUILD-UP CURVE ANALYSIS INTRODUCTION Pressure build-up characteristics of shut-in wells have been used for many years by engineers to evaluate average reservoir pressure, effective permeability thickness of the pay section, effectiveness of well completion (skin effect) and reservoir size. A number of methods of analysis have appeared from time to time."" Without exception, these methods are based on the assumption that the reservoir contains but one fluid of constant small compressibility and constant mobility. It has been suggesteda" hat these single-fluid methods may be applied to data from reservoirs containing both oil and gas by substitution of some effective total properties of the multiphase system in place of the corresponding single-phase properties. The present investigation was undertaken primarily to evaluate this approach. METHOD A number of theoretical build-up curves were calculated for conditions of two-phase flow, under the assumption of certain reservoir and fluid properties, and were analyzed by single-fluid methods with appropriate total compressibility and total mobility values for the corresponding single-fluid properties. Results of the analyses were compared with the assumed conditions. The theoretical build-up curves were completed by procedures similar to those of West, Garvin and Sheldon." Since these calculations require a considerable amount of computer time, an attempt was made to derive an approximate calculation method. The attempt was unsuccessful for calculating build-up curves, but the effort did result in a new approximate method of calculating the performance of solution gas drive reservoirs, which appears to be an improvement over the constant-GOR method" used previously (see Part 2). The West, Garvin and Sheldon calculations involved the following assumptions: (1) the reservoir is circular and completely bounded, with a completely penetrating well at its center; (2) the porous medium is uniform and iso-tropic, with a constant water saturation at all points; (3) gravity effects can be neglected; (4) compressibility of rock and water can be neglected; (5) the composition and equilibrium are constant for oil and gas; (6) the same pressure exists in both the oil phase and the gas phase; and (7) no afterproduction occurs, i.e., the well is shut in at the sand face for build-up. These assumptions make it possible to describe two-phase flow of oil and gas by the partial differential equations:

Jan 1, 1967

By J. F. Martin, J. E. Friedline, L. M. Melnick, G. E. Pellissier

A method has been developed for determining gases in steel in which the gases are extracted by vacuum fusion and analyzed by mass spectrometry. This method is especially applicable for determining small amounts of oxygen and nitrogen. The lower limit of detection for both oxygen and nitrogen is 0.0001 pct with an average deviation of about 0.0001 pct for oxygen and 0.00008 pct for nitrogen. Hydrogen analyses by this method are in good agreement with analyses by vacuum tin-fusion. A complete analysis of all gases extracted from a steel sample can be accomplished in 40 min. It has long been known that certain of the gaseous elements exert particular effects on the physical and mechanical properties of steel. With the recent advent of new experimental steelmaking practices it is now possible to obtain steels with extremely low contents of such elements as oxygen and nitrogen. Present methods for determining these two gases are of insufficient sensitivity as applied to steels made by these newer practices, and therefore a new technique was considered necessary. A review of the literature showed that the conventional vacuum-fusion method1 was the procedure most suitable for modification to obtain a sensitive method for determining the gaseous elements present in steel simultaneously. However, for trace amounts of gases it is impractical to increase the sample size because of the physical limitations of the apparatus and the increased degassing time which would be necessary. Also, when only one of the desired gases is present in trace amounts, quantitative separation by conventional vacuum fusion is almost impossible. Since small analytical mass spectrometers are now commercially available, it was decided to use this instrument in conjunction with conventional vacuum fusion for both qualitative examination and quantitative determinations of the gaseous elements in steel. This technique has already been used with success for the analysis of copper and lanthanum.2"4 APPARATUS A vacuum-fusion unit was first designed and built, Fig. 1. This unit was designed especially for the determination of small amounts of nitrogen and/or oxygen in steel since an acceptable method (vacuum tin-fusion)5 was already in use for the determination of hydrogen in steel. The vacuum system was con- structed of Pyrex glass throughout. Mercury check valves were used in place of stopcocks in the high-vacuum section of the system. Ferromagnetic samples are normally introduced into the vacuum system through a mercury lift. In case of a break in the lift, the mercury is prevented from entering the furnace by a self-locking, mercury safety check valve, Fig. 2. This check valve, developed at our laboratory, is a modified T/S ball-and-socket joint placed in the sample inlet system between the mercury and the furnace. A sudden entry of mercury or air pushes this valve into the closed position, completely blocking the top of the mercury lift. The check valve is held in the closed position by the vacuum of the apparatus. A separate entrance tube also was placed in the system, so that samples which would amalgamate with mercury, or which are nonmagnetic, could be introduced into the furnace. A slightly modified Guldner-Beach furnace is used in this system.8 The optical flat was connected to the system by a ground-glass joint to facilitate cleaning, a quartz delivery tube was substituted for Pyrex, and the crucible was used uncovered. Fusion takes place in a high-purity graphite crucible insulated by — 200-mesh high-purity graphite powder. The crucible set (crucible and funnel) is suspended in a quartz crucible. The furnace envelope and the sample inlet system are exhausted by 2 single-stage, mercury-diffusion pumps in series. The first of these is a high-speed pump, which will work only with a very low backing pressure. This pump is used to exhaust the gases from the furnace as rapidly as possible. The second mercury-diffusion pump is used as a backing pump for the first; the capacity of this pump is much lower. However, it will work efficiently against a back pressure of 3 mm of mercury. To assure positive pumping, a modified Naughton and Uhlig condenser7 is used. With the use of proper valving the second pump is also used as a circulation pump in the catalyst system.

Jan 1, 1959

By G. L. F. Powell, L. M. Hogan

Samples of silver and Ag-O alloy, 0.12 wl pet, have been undercooled to a maximum of 250°C by melting in a slag of commercial soda-lime glass. Grain refinement occurred in undercooled silver samples whet the degree of undercooling exceeded a critical value which varied from 153° to 175°C, depending on the melting conditions. In under cooled Ag-0 samples, the grain structure was fine and equiaxed at degrees of undercooling larger than 50oC. The grain refinement in both silver and Ag-0 alloy samples was the result of recrystallization during or immediately following freezing. When the oxygen content was increased, the degree of undercooling at which recrys-tallization occurred decreased. INVESTIGATIONS of the grain structure of strongly undercooled melts of nickel and cobalt1-4 have revealed a change in room-temperature grain size as a function of undercooling. At large degrees of undercooling the grain size was fine, but samples undercooled by smaller amounts exhibited a much coarser grain structure. The reported undercooling at which this grain size transition occurred in nickel varied between 140° and 175'C undercooling. Walker1 attributed this grain size change to an enhanced nuclea-tion rate. Cavitation at a melt-solid interface was con sidered to generate a high-pressure wave. This, in turn, produced nucleation in the liquid as a result of an elevated freezing point increasing the undercooling in the melt. A similar grain size transition was observed in undercooled silver between 133° and 153°C undercooling.5 The grain size transition in undercooled silver was found to be dependent on oxygen con tent. More recently, the authors5 showed that a simila transition occurred in an undercooled Cu-O alloy, =0.08 pet, but not in oxygen-free copper undercooled more than 200°C. The transition was attributed to recrystallization of the solidified grains after freezing; the greater the degree of undercooling, the more complete the recrystallization. In view of this effect of oxygen content on grain size transition, additional undercooling experiments were carried out with silver. This paper reports the results of these studies on silver and Ag-O alloys undercooled more than 50°C EXPERIMENTAL 1) Materials. The silver was supplied by Matthey Garrett Pty., Ltd., as fine granulate, silver mini' mum 99.9 pet, and contained the trace impurities shown in Table I. The silver originated from two different refineries resulting in the difference in impurity content shown in Table I. However, the variation in composition had no discernible effect on the undercooling behavior. Samples of the silver were melted in an open-ended, vertical, cylindrical resistance furnace, wound with Kantha] wire, the temperature of which was adjusted by means of a variable transformer, 2) Melting and Undercooling. Bulk samples of silver were undercooled by preoxidation and removal of oxidized impurities in a glass slag, as previously described,5 but special precautions were necessary to allow comparison between oxygen-free silver and a Ag-O alloy. Samples free of oxygen were prepared by first melting weights of granulate varying from 200 to 400 g in clean vitreous silica crucibles in air for approximately 15 min. A sample was frozen and a complete glass cover formed over the top surface by adding granulated glass rod. The sample was then melted and frozen several times until no oxygen bubbles appeared in the glass slag. Subsequent thermal analysis indicated that the oxygen content was negligible. The glass acted as an efficient barrier to absorption of oxygen from the air. Ag-O alloys were prepared by melting silver in silica crucibles with only a partial cover of glass, so that the top surface of the melt was in contact with air. The silver was held molten in air for several hours at 1000° to 1050°C, in the expectation that the oxygen in solution would reach equilibrium with atmospheric oxygen. At T - 1000°C the equilibrium oxygen content was calculated as 0.12 wt pet O, using Sieverts law and Eq, [l], derived by Mizikar et a!.7 from the data of Sieverts and Hagenacker8 for 1 atm of oxygen:

Jan 1, 1970

By M. E. Wadsworth, T. L. MacKay

Sintered UO, samples were leached in sulfuric acid solutions of various concentrations. A pressurized system was used so that it was possible to investigate the kinetics of the reaction to 270°C with oxygen overpressures as high as 900 psi. The rate was observed to be a function of the concentration of hydrogen ions and directly proportional to the partial pressure of oxygen. Evidences are presented which indicate that a UO2 surface site reacts with a molecule of water to form a hydroxyl complex which in turn can dissociate with the characteristics of a weak acid. A rate determining step has been proposed which involves the reaction between an oxygen molecule and the hydroxyl complex on the surface of UO,. ThE 2 principal methods for uranium dissolution are carbonate and acid leaching. The sulfuric acid leach is the more popular and is used for the treatment of the majority of the ores of Africa, Canada, and the United States. Low recoveries in basic leach circuits led investigators at the University of British columbia1 to study leaching of uranium ores in pressurized vessels. Early in the study of dissolution of uranium it was found that only uranium in the hexava-lent state could be leached in acid or basic circuits. Therefore, the use of oxygen over-pressure in an autoclave offered an interesting approach to solving the low recovery problems of carbonate leaching. UO, was used in this study because it could be obtained in high purity and also because it is representative of the most refractory of the primary uranium minerals. A pressurized system was used to provide a means whereby important temperature and pressure parameters could be varied for the evaluation of the kinetic processes. The mechanism for the dissolution of UO2 should be similar for any of the uraninite type minerals. At the present all of the kinetic studies that have been conducted have been in carbonate media. Peters and Halpern2 carried out a kinetic study of the leaching of pitchblende. Their specimens were pitchblende ores selected on the basis of high grade and homogeneity. Pearson and wadsworth3 conducted a kinetic study of the dissolution of UO2 in carbonate solutions with results very similar to those obtained by Peters and Halpern. EXPERIMENTAL The UO2 as received* was found by spectroscopic *Mallinckrodt reagent, supplied by the Atomic Energy Commission. analysis to have a purity of 99.94 pct 0.03. It was ground in a mechanical agate mortar and screened through a 400-mesh sieve. Thin flat disks of UO2 were prepared by pressing the sized powder in a specially constructed die at a total pressure of approximately 25 tons per sq in. The pressed disks were approximately 0.25 cm thick, 1.6 cm diam, and 3.5 g in weight. As pressed, the samples were approximately 65 pct of theoretical density. It was essential that a binder such as poly-vinyl alcohol be added before pressing to prevent formation of cracks when fired. It was found by trial and error that pressed disks sintered close to 1870°C in a hydrogen atmosphere, reach densities between 91 and 91.5 pct of the theoretical density. At this density the samples had no measurable porosity based upon a 2-hr emersion in boiling water. An X-ray examination of these sintered UO2 disks showed them to be identical with the original UO2 received. Reaction of UO2 with the alundum furnace core was prevented by placing molybdenum between the samples and the core according to the method of Corwin and Eyerly.4 Leaching studies were carried out in a specially designed autoclave, the details of which have been presented e1sewhere.6 The main features of the equipment as applied to a kinetic study are: 1) solution samples may be removed during the course of a single run; and 2) temperature and agitation are carefully controlled. Two-ml samples removed from the autoclave during the course of a run were analyzed for U3O8 content with a Beckman model DK-2 spectrophotometer by the method of

Jan 1, 1959

By J. J. Gilvarry, B. H. Bergstorm

The first part of this paper describes a new approach to the problem of energy relationships in fracture and comminution. The basic theoretical method used (as contrasted to previous empirical or semi-empirical approaches) is an attempt to verify or disprove the hypothesis of von Rittinger, In the second part of the paper, theoretical conclusions for the three headings of single fracture, plural fracture, and comminution are tested. Agreement between theory and experimental results is demonstrated. The classical hypotheses on the energy relationship in brittle fracture are those of von Rittinger-that the energy required is proportional to the new surface formed— and the alternative one of Kick.' Over the years, repeated attempts have been made to discriminate experimentally between the two assumptions, or to replace both of them, as has been done by Bond? However, no semblance of agreement seems to exist among different investigators. This paper is the first in a series essaying a new approach to the entire problem of energy relationships in fracture and comminution. The purpose of this program of research is to verify (or disprove) the hypothesis of von Rittinger. In contrast to previous empirical or semi-empirical approaches to the problem, the basic method adopted is theoretical. The prime aim is to determine the distribution function for fragment size theoretically and to verify the predicted distribution experimentally. This approach should make possible determination of the total surface of the particles by integration over the corresponding distribution function; one may hope that this direct procedure will yield more definitive results than use of gas adsorption or other methods. The surface area so obtained can then be compared with experimental energies, as obtained by Bergstrom, Sollenberger, and Mitchell,4 for example. It is convenient to distinguish between single fratture, plural fracture, and comminution. Single fracture of a brittle solid is defined as fracture by an external stress system which is removed instantly and permanently when fracture is initiated. Plural fracture consists of single fracture of the original specimen, followed by a sequence of only a few secondary fractures; comminution differs in that the sequence of secondary fractures consists of a large number of repetitive steps. For single fracture, Gilvarrys has given a rigorous derivation of the proper distribution function for fragment size, based on a closely defined physical model and deduced strictly by the laws of probability. Gilvarry and Bergstrome have compared the predictions of the theory with experiment, and have found excellent agreement. The prior work of Gilvarry and of Gilvarry and Bergstrom considers only single and plural fracture. The purpose of the present paper is to extend the discussion to the case of comminution. For prefatory purposes, the theory of Gilvarry for single fracture will be outlined. The considerations will then be applied to the cases of plural fracture and comminution. SINGLE FRACTURE The derivation is based on the Poisson law.7 For Points distributed with mean density y over a domain, this law states that the probability p(t)dt of one point lying in the range t to t + dt with none in the interval 0 to Ms p(t)dt = ye-yt dt. For this distribution to be valid, it is necessary and sufficient that the points be distributed at random, individually and collectively in the sense of ~~~.7 The former qualification requires that the position of One Point be independent of another, and the latter implies that the probability of a region containing a

Jan 1, 1961

By R. J. Blackwell, J. R. Rayne, W. M. Terry

This paper presents results of an experimental investigation of factors that control the efficiency with which oil is displaced from porous media by a miscible fluid. The study was made to elucidate the relevant processes both on microscopic level (within individual or between neighboring pore spaces) and on macroscopic level (within a large sand body). Mixing of miscible fluids on the microscopic level was studied in sand-packed tubes. It was found that molecular diffusion is the dominant dispersion mechanism for reservoir conditions of rate, length and pore sizes. Macroscopic channeling was studied for various mobility ratios in reservoir mode1s scaled to relate viscous, gravitational, and diffusional forces. The formation of channels was due to viscous firzgering, gravity segregation and variations in permeability. With adverse mobility ratios, it was found for reservoirs of realistic widths that diffusion will not be effective in preventing the formation and growth of fingers, even in homogeneous sands. At sufficiently low rates channeling was eliminated by gravity segregation in tilted reservoirs. The dependence of recovery on mobility ratio, length-to-width ratio, flow rate and angle of dip is presented. INTRODUCTION Oil recovery by solvent flooding is finding increasing application in the field. While the process promises high recoveries from the region swept by solvent, under adverse conditions only a small fraction of the reservoir volume may be swept at the time solvent breaks through to the producing well. Further, the high cost of the solvent encourages its use only as a bank whose size must be kept at a minimum. Thus, two important questions arise: (1) what fraction of the reservoir can be swept, practically, by solvent? and (2) what is the minimum size solvent bank that can be used to carry out the displacement? The answers to these questions require knowledge of both macroscopic channeling processes and microscopic mixing processes. The studies described here were carried out to gain this knowledge. Microscopic mechanisms which cause mixing will be discussed first, because an understanding of these mechanisms is necessary for proper interpretation of the experimental work on channeling described later. INVESTIGATION OF MICROSCOPIC DISPLACEMENT PHENOMENA FOR MISCIBLE FLUIDS The questions to be resolved concerning microscopic displacement behavior are: (1) how much mixing occurs between oil and solvent in the direction of flow? and (2) does the solvent completely flush the oil from the pore spaces in the region invaded by solvent? The classic work of Sir Geoffrey Taylor' and its extension by Aris2 have shown that it is possible from theory to describe the amount of mixing in single straight capillaries when solvent displaces a fluid of equal viscosity and equal density. Aris showed that the lengths of the mixing zones, 81, corresponding to the 10 per cent and 90 per cent concentration levels of the solvent can be calculated from the formula, 81 = 3.62 \/Kt, .........(1) where 8l is the length of the mixing zone, cm; t is the time, sec; and K is the effective dispersion coefficient, cm2/sec, given by K/D = 1 + (1/48) (au/D)2......(2) where a is the radius of the capillary, cm; D is the molecular diffusivity, cm2/sec; u is the average velocity

By David A. Thomas, John F. Peck

Fibrous microstructures and their development have been studied by metallography and by hardness and quantitative metallographic measurements. Thin, curved grains were observed in transverse sections of commercial tungsten wire and were developed in laboratory-drawn iron wire at reductions greater than about 70 pct. The unusual grain shapes are attributed to the nonuniform mode of deformation of the grains associated with the characteristic [I101 preferred orientation that develops in bcc metals during drawing. Transverse sections of heavily swaged iron wire showed a ospiral nebulao structure that appears to result from the elliptical shape of the cross section of the swaging dies and their rotary motion. Hardness and quantitative metallographic results show unusual trends that parallel the changes in microstructure. Other phenomena are considwed in light of the results, including the strength of fine wires of bcc metals, ocylindricalo textures, and the ductility of fibrous tungsten. 1 HE microstructure of tungsten wire in a longitudinal section shows grains highly elongated in the working direction, giving the appearance of a bundle of parallel fibers, Fig. l(a). Wires of tungsten and other bcc metals show a pronounced preferred orientation in which a [I101 direction is parallel to the wire axis.l The microstructure and texture of tungsten develop during warm working (working above room temperature but below the recrystallization temperature) into wire by swaging and drawing. Fibrous tungsten so produced has limited but useful room-temperature ductility, whereas annealed tungsten is brittle at room temperature. During an investigation of the effect of recovery annealing on the substructure of commerical 0.030 in. diam tungsten wire, it was found that transverse sections of the as-received material showed thin, curved grains, Fig. 1(b), instead of the expected equiaxed grains. The fact that Figs. l(a) and (b) show different sections of the same material at the same magnification is at first disconcerting, because the long dimension of the curved grains in the transverse section is several times the average fiber width in the longitudinal section. However, a longitudinal section cuts most of the curved grains across their narrow dimension, giving the appearance of a multitude of narrow fibers. On the other hand, the plane of a longitudinal section occasionally includes a long dimension of a curved grain, resulting in a broad ofibero, such as A in Fig. l(a). EXPERIMENTS AND RESULTS To investigate the origin of fibrous microstructures and particularly the unusual transverse-section structure of tungsten, rods of iron (also bcc) were drawn and swaged in the laboratory. One rod of annealed Ferrovac E was drawn in six steps to a total reduction in area of 87 pct (0.10 to 0.036 in. diam), and a second rod was swaged to the same reduction, using approximately the same reduction steps. Longitudinal and transverse sections of the swaged and drawn iron wires are shown in Fig. 2, for reductions in area of 39, 73, and 87 pct. A fibrous microstructure in longitudinal section is first evident at 73 pct reduction, and it is further developed at 87 pct reduction. The grains in transverse section remain equiaxed until about 73 pct reduction, when they become somewhat elongated and curved. At 87 pct reduction, the drawn wire shows thin, curved grains in the transverse section, similar to those observed in tungsten. The transverse sections of the iron wires swaged 83 and 87 pct reduction in area show the striking ospiral nebulao structure in Fig. 2(d). This struc-

Jan 1, 1962

By Wieland D. R., Hurst R. E., A. R. Hendrickson

Analysis of acidizing techniques, in correlation with reservoir data and a backlog of past treatments, has resulted in the development of a valuable engineering guide for planning acidizing treatments. Such treatments fall into three categories: (I) acid injection into the pores of the matrix; (2) acid injection into natural formation fractures at less than parting pressure; and (3) combination acidizing-fracturing treatments in which acid solutions (without propping agents) are injected at treating pressures sufficient to open and extend fractures through which the acid flows. Because the spending tirnze of acid during a specific well treatment does not change appreciably, maximunl penetration is attained when the first increment of injected acid is completely spent. Additional acid injection cannot be expected to further extend the benefits of the treatment. Depth of penetration will depend upon the reaction rate of the acid under treatment conditions, the injection rate of the acid into the matrix or fractures and the area-volume relationship existing in the flow channels. Based on Darcy's flow formula, extremely low injection rates must be used in order to keep bottom-hole injection pressures below formation fracturing pressure. As a result, only limited penetration of unspent acid will occur. Treatment records indicate that, in most acidizing treatments, formation parting pressures are exceeded, greatly extending acid penetration. Under these conditions, stimulation benefits are limited to the fracture area produced during the spending time of the first increment of acid injected into the formation. This area may be calculated from laboratory and well data to estimate depth of penetration. This, in turn, may be correlated with productivity data to assist The art of gas and oil well acidizing has been characterized by many changes in treating materials and techniques since its inception. These developments have been designed to provide greater production increases. prolong production declines and shorten payout time. Such improvements have been based primarily on data derived from laboratory research and field experience. As more of the variables influencing these treatments have been recognized and evaluated, acidizing has become less of an art and more of a science. Recent studies of fracturing treatments,' in light of individual well conditions and the results of thousands of fracturing treatments, made possible the formulation of an engineering guide that is now being used to select optimum treating techniques and to forecast probable results of such treatments. A similar analysis of the factors controlling acidizing treatments has been made and is the basis for this paper, The findings herein can be used as a guide in the selection of acidizing solutions and techniques, tailored to fit specific well conditions and to provide optimum stimulation per dollar cost. Acidizing treatments may be classified into three basic categories—(1) treatments in which the acid is injected uniformly into the pores and flow channels of the matrix, (2) treatments in which the acid enters natural fissures and fractures in the formation at less than fracturing pressures and (3) injection of acid into the formation at a pressure sufficient to open and extend fractures into the rock through which the acid penetrates (without the inclusion of a propping agent). TYPE 1 — MATRIX ACIDIZING This category consists of treatments in which acid solutions are injected into a homogeneous carbonate

By V. O. Struebing, K. A. Johnson, F. H. Ellinger, C. C. Land

The system Pu-Ce-(Pu,Ce)Co2 has been investigated by thermal, micrographic, and X-ray diffraction methods. The principal characteristics of this system were found to be: 1) a "eutectic" valley extending across the phase diagram from the 6 Pu/Pu3 Co eutectic at 408°C to the 7 Ce/Ce24Co11 eutectic at 42Z°C, the maximum temperature in the valley being 443°C; 2) three ternary inter mediate phases, each having an appreciable mutual solid solubility for plutonium and cerium; 3) eleven four-phase equilibria comprising one eutectic, four peritectics, and six Class II reactions; 4) a pseudobinary section between the ternary phase A and the y Ce solid solution, and another such section between the 6 Pu solid solution and the 77 phase (Pu3 CO) containing cerium in solid solution. ALLOYS of plutonium that are completely molten at temperatures in the vicinity of 500°C are of interest for use in the liquid state as fast reactor fuels. Since exploratory work1,2 on Pu-Ce-Co alloys had indicated that a "eutectic" valley extended across the ternary system between the Pu-co[6 Pu/n (Pu3co)] and the Ce-CO[? Ce/? (ce24co11)] binary eutectics, thus allowing alloys of any plutonium concentration to be selected within the temperature limitation, the present investigation was begun for the purpose of accurately defining the "eutectic" valley. However, when the existence of three ternary intermediate phases and the consequent complexity of the liquid-solid equilibria were discovered, the investigation was expanded to cover the complete system, except for the corner between the isomorphous series PuCo2-CeCo2 and pure cobalt. 1) EXPERIMENTAL The purities of the alloy components are given in Table I. Alloy buttons, each weighing about 5 g, were prepared by arc melting the accurately weighed components (±0.1 mg) under a purified He-Ar atmosphere. Weight loss as a result of melting was usually less than 1 part in 1000; therefore, the buttons were not chemically analyzed. One half of each button was re-melted by induction heating in vacuo in a magnesia crucible and furnace-cooled, so as to produce a somewhat coarser microstructure than that of the rapidly cooled arc-melted button. About 20 g of each alloy to be examined by thermal analysis were remelted in vacuo in a tantalum crucible of the inverted thermocouple well type. The thermal, X-ray diffraction, heat treating, and micrographic methods used in this investigation were essentially the same as have been described in earlier reports.3,4 Two etchants for these alloys were found to be satisfactory: A) 8 (parts by volume) orthophosphoric acid, 5 ethanol, and 5 glycol; B) 7 tetraphosphoric acid, 36 water, and 57 2-ethoxyethanol. Etchant A was used at 5 to 6 v for 15 to 20 sec, and etchant B was used at 10 v for 10 sec. 2) THE LIMITING BINARY SYSTEMS 2.1) Plutonium-Cerium. The phase relationships of this system have been investigated by Ellinger, Land, and Cramer,5 and by Selle and Etter.6 The over-all agreement of their phase diagrams is good, but there are some discrepancies in the positions of some phase

Jan 1, 1967

By J. Szekely, J. W. Evans

A formulation is given for radiant heat loss from the surface of molten steel held in a ladle and numerical solutions are presented for the resultant integrodif-ferential equations. The results are shown on plots of the net rate of heat loss against time, with the geometlry, emissivities, thermal properties of the walls, and Preheat as parameters. It was found that the rate of the vadiant heat loss from the exposed metal surface was comparable to the conductive losses from the Portion of the metal in direct contact with the walls. The radiant heat loss was markedly affected by the geometry and by the thermal properties of the walls the net radiant heat loss from the steel could be drastically reduced by providing a top cover and also by the provision of preheat. The results of this work are thought to be relevant to heat loss calculatious concerned with vacuum degassing and with operatzons in which molten metals are held in ladles. WhEN molten steel is held or transferred in ladles heat losses inevitably occur. The accurate assessment of these heat losses is of importance in process design considerations, as the temperature of the metal may have to meet rigid specifications in certain stages of the processing sequence. As shown in Fig. 1, the principal mechanisms responsible for heat loss are: i) Conduction into the walls and ii) Thermal radation from the top metal surface (possibly augmented by natural convection) The conductive heat loss may be readily calculated as much information is available in the literature on both the general theory1'' and on its particular applications to ladle heat loss calculations. Evaluation of the radiant component of the heat loss is less straight-forward. In two recent papers it has been shown that the presence of a slag layer 2 to 3 in. thick would prevent any significant heat loss for periods up to an hour,5'" however, should bare metal be exposed, constituting the radiant surface, substantial heat losses could occur. This latter case is of considerable practical interest in various degassing operations and in the vacuum processing of certain specialty steels, where the top surface of the metal is not covered by a protective slag layer. Under these conditions a substantial radiant flux is emitted by the top metal surface; part of this flux leaves the system, part of it falls in the ladle walls. The radiant flux incident on the ladle walls is in part absorbed and conducted away, and in part reflected and re radiated. The purpose of the paper is to present a general formulation of the problem, the resultant solution of which should enable computation of the net radiant heat loss from the metal surface. The principal variables, the effect of which will be investigated, include initial metal temperature, geometry, thermal and radiative properties of the wall, and various modes of preheat. FORMULATION Consider a cylindrical enclosure, shown in Fig. 2, with a diameter D, height H, and wall thickness W, the cover at the top of the cylinder having a wall thickness L. If the temperature of the base (representing the top surface of the metal) is kept constant, say, Tb, the problem is to calculate the net rate of heat loss from the base as a function of time, geometry, the initial temperature distribution, and the thermal properties of the system. As customary, the formulation will consist of the mathematical statement of the conservation of thermal energy within the system.

Jan 1, 1970

By T. F. Bower, M. C. Flemings, H. D. Brody

A series of experiments are reported which show thai reasonable assumptions for analysis of solute redistribution in solidification of castings and ingols illclr(de: negligible undercooling before nuclcation of solid phases, negligible increase of solute in advance of the tips of growing dendrites, complete diffusion within the Liquid oer distances the order, of dentrite arm spacings, platelike dentrite t/o?$holog?.. Qzmli/ali,e and periment is obtained JOY Al-Cu alloys front measurements of amount of nonequilibrium second phase, minimum solute content at the centers dentrites arms (after solidification), and 3) change in minimum solute content in dentrite arms during solidification. In a previous paper' analyses were given describing solute redistribution in dendritic solidification. The analyses were undertaken primarily to permit prediction of "severity of microsegregation" in a casting or ingot after solidification and cooling to room temperature, but can also be used to predict "local fraction solid' as a function of temperature in the solidifying alloy. Numerical results were given for Al-Cu alloys. The aims of this paper are to examine the validity of some of the assumptions of that analysis, and to compare predictions of the theory with experiments. All experiments reported are on Al-Cu alloys, the bulk on Al-4.5 pct Cu alloy. EXPERIMENTS Experiments conducted were on small ingots, 41/2 in. in diam by about 10 in. high. solidified in a controlled thermal environment (i.e., in a resistance-heated, re-circulating air furnace, whose temperature could be controlled to 1°C). In one series of experiments, the ingots were solidified very slowly (solidification time was varied from 8 to 1000 hr) and at constant rate of heat extraction. In a second series of experiments the ingots were solidified ( 'unidirectionally". Fig. 1 is a sketch of the method of making the two different types of ingots. The unidirectionally solidified ingots were fully columnar, and the slowly solidified ingots had equiaxed structures, Fig. 2. In making the slowly solidified ingots, the mold employed was carbon steel, coated with inert refractory. Typically the mold and its environment were heated to 700°C. The molten alloy was then poured into the mold and the temperature of the system equilibrated above the liquidus. The environment of the mold was next cooled with a controller responding to a differential thermocouple that measured the difference between the temperature of the melt and the temperature of the mold environment. The differential thermocouple consisted of one chromel-alumel junction in the melt connected in opposition with a second junction in the furnace. The potential across the differential thermocouple was input to the controller which acted to keep this potential at a set. constant value; hence the temperature difference between the two junctions was maintained nearly constant. Major resistance to heat flow in the above apparatus was between the metal surface and the furnace atmosphere, not in the metal itself. By maintaining a constant temperature difference between the metal and the furnace atmosphere, the controller maintained an essentially constant rate of heat extraction during solidification, regardless of amount of heat released by the metal per unit temperature change. When the metal temperature dropped just below the eutectic temperature. the ingot was removed from the furnace and water-quenched. The unidirectionally solidified ingots were cast in a similar mold to the foregoing except that provision was made to extract heat through the bottom face of the mold by passing air or water through the mold base, Fig. 1. For these experiments, metal and mold temperature were again equilibrated above the liquidus. Temperature of the mold environment was then held

Jan 1, 1967