Search Documents

By H. Dykstra

Calculated wellbore pressures were obtained for parameters of radilcs ratio and permeability. In all cases bur two, after-production was allowed to occur for one day. The calculated pressure build-up data were compared with actual pressure build-up data from a condensate well. The paper discusses the comparison, gives reasons for, inability to match calculated (2nd actual data, and presents conc1usions derived from the comparison. INTRODUCTION Considerable information has been published on the decline in wellbore pressure for oil and gas fields resulting from production of fluids. Considerable information also has been published on pressure build-up for oil fields, but relatively little has been published for gas and condensate fields. Perrine' summarized the equations derived for pressure build-up in oil wells and showed how the different methods could be applied. On the other hand, very little theoretical information has been published on pressure build-up in gas or gas-condensate wells.2-4 Tracy,' by pointing out the similarity between the equation describing oil wellbore pressure decline and the equation describing gas wellbore pressure decline, showed how methods of pressure build-up analysis could be applied to a gas well. In analyzing pressure build-up data for gas or gas-condensate wells, it would be very desirable to compare actual pressure build-up data with calculated pressure build-up data. Such a comparison would lead to a better feeling for the quantitative picture of gas-well build-up curves and would help to establish a degree of confidence in a given build-up curve. This paper discusses the comparison and presents conclusions made from it. Calculated wellbore pressures are given for parameters of radius ratio and permeability. In all cases but two, after-production was allowed to occur for one day. The work was done in an attempt to evaluate the reserves of a gas-condensate field. As will be discussed later, it was not possible to make an evaluation. METHOD OF CALCULATION The method used for solving the depletion and pressure build-up behavior for a gas-condensate well was that developed by West, Garvin and Sheldon.5 or the IBM 704 computer program used, it was assumed that liquid condensing out of solution would have only a minor effect on the flow behavior. Therefore, the field was treated as having single-phase flow, with the gas saturation S, remaining constant at unity minus the connate-water saturation Sw. As will be discussed later, this may not have been a good assumption. It also was assumed that oil production could be considered as equivalent gas production by using a conversion factor based on an analysis of the liquid produced. An outer boundary condition of no-flow was assumed. BASIC DATA The data required for the study were reservoir properties, fluid properties and production data. Reservoir properties included the following: thickness h, 179 ft; porosity, 0.194; connate-water saturation, 0.43; gas saturation, 0.57; permeabilities k, 1.0, 0.5, 0.25 and 0.15 md; radius ratios re/rw, 800, 1,000 and 2,000; wellbore radius rw, 5 in.; initial reservoir pressure P, 6,529 psia; and reservoir temperature, 272°F. The selection of permeabilities was based on core-analysis data showing an average air permeability of about 2 md. With a relatively high connate-water saturation, the average effective gas permeability would be about one-half, or less, of the air permeability. The selection of radius ratios was based partly on the observed rapid decline in wellbore pressure during production and partly on the fact that the well was believed to be located in a relatively small fault block. Reservoir fluid properties were obtained from a reservoir fluid study on a recombined sample of gas and condensate. Gas formation-volume factors were calculated from pressure-volume measurements made at reservoir temperature. Gas viscosities were calculated by the method of Carr, et al,G from an analysis of the re-combined sample. These data are shown as a function of pressure in Table 1. The production schedule from the gas-condensate

By G. T. Hahn, C. N. Reid, A. Gilbert

An evaluation is made of the possible cautsal relationship between twinning and fracture in molybdenum. For both single and poly crystalline material no instance of twin-induced fracture was observed. Instead brittle fracture was found to be slip-induced at heterogeneities. For single-crystal material, the yield stvess in compression and the fracture stvess in tensiow obey a similar angular- relationship which follows approximately a 1/cos2 ? law. where 0 is the angle between the specimen axis and the nearest (100) plane. This sittlilarity between yield and fracture behavior casts doubt on the interpretations made previoltsly that a 1/cos2? relationship supports a critical normal fracture stress criterion. SINCE fracture in brittle materials takes place at stresses an order of magnitude lower than the theoretical strength of the lattice, it must be postulated that some stress-concentrating effect is operative during the fracture-initiation process. Several mechanisms have been proposed whereby the necessary stress concentrations could be produced. The well-known Cottrell mechanism1 describes a dislocation interaction which can lead to the formation of an incipient cleavage crack on a cleavage plane, and both zener2 and stroh3 have discussed models which predict the stress concentration at the intersection of a slip band with a grain boundary. When twinning occurs twin/twin intersections and twin/grain boundary intersections represent another possible means of fracture initiation. Although still controversial, the concept of twin-induced fracture is supported by a weighty mass of evidence, which has been reviewed in a recent paper.4 This investigation was conducted in order to assess the extent to which the brittleness of molybdenum can be ascribed to this cause. Emphasis was placed on seeking direct metallographic evidence for twin-induced fracture such as arrangements of twins located at the fracture origin. EXPERIMENTAL PROGRAM Materials. The experimental materials are described inTable I. From the chemical analyses and estimates of the amounts of interstitials likely to be retained in solution,5 it is concluded that all were mul-tiphased systems. X-ray diffraction experiments6 showed no evidence of preferred orientation in the polycrystalline materials. Single crystals were grown from Molybdenum X and Y using the floating-zone technique,' in a vacuum of better than 3 x 10-5 Torr. However difficulties were experienced with Molybdenum X due to violent gas evolution from the molten zone, and additional crystals (Crystals 3 to 6) were produced by annealing for 5 1/2 hr at 2300°C under a vacuum of 10-4 Torr. Techniques. A series of tension, compression, and bend tests was conducted. Tension tests were used to demonstrate brittle behavior, to measure brittle-fracture stresses, and to provide fractured specimens for metallographic inspection. Compression tests were employed in order to obtain some ductility and a measure of the yield stress. Furthermore, it was considered that compression loading would permit the study of crack nucleation in the absence of propagation. The bend tests were conducted to facilitate identification of the fracture origin, which would be expected to be at, or near, the tensile surface of the sample. Tests were carried out at temperatures between 78" and 298°K attained by means of a liquid-nitrogen evaporator of modified Wessel design.' Specimens were fashioned by mechanical grinding, and, prior to heat treatment, were electropolished at 10 v in a 3:l mixture of ethyl alcohol and sulfuric acid, using a stainless-steel cathode. Except where noted, a strain rate of 4 pct per min was used. All specimens had a gage section 0.75 in. long; diameters of 0.10 and 0.20 in. were adopted for the tension and compression specimens, respectively, the former of single shoulder design, the ends being gripped in split collets with a conical bearing surface. In the bend test, electropolished specimens measuring 0.080 by 0.25 by 1.25 in. were deformed at 78°K by four-point loading. The bending device had spans of 1 in. between the outer and 4 in between the inner fulcra; it was stressed between the compression anvils of the Instron machine at a constant deflection rate. EXPERIMENTAL RESULTS Polycrystals. Tension and Compression Tcsls. It is apparent from Fig. 1(a) that above 170°K fracture of Molybdenum X takes place at stresses equal to or greater than the compressive yield stress whereas

Jan 1, 1967

By A. V. Bradshaw, R. I. L. Guthrie

The behavior of large bubbles in the size range 4 to 25 cm3, rising through molten silver, has been studied. It was found that rising velocities were equivalent to those in aqueous systems of low viscosity. Mass transfer coefficients for oxygen bubbles dissolving in silver were found to be 0.036 ± 0.007 cm sec-1, being close to those predicted for transfer through the front surface of the spherical cap bubble only. It is suggested that the surface active nature of oxygen in silver could account for the relatively low coefficients obtained. MANY metallurgical processes involve interactions between gas bubbles and liquids. Examples include the removal of carbon monoxide in Open Hearth Steelmak-ing, the removal of sulfur by blowing air through copper matte during converting, and the removal of hydrogen from steel during vacuum degassing or inert gas flushing. The steps involved in such refining processes include; transport of the dissolved species to the bubble interface, adsorption and chemical reaction of the species at the interface, desorption of product molecules from the interface, and transport of product gas into the bulk gas phase of the bubble. It has been concluded1 that all the interfacial steps involved proceed so rapidly at steelmaking temperatures that transport of the solutes, present in the metal, become the important rate controlling factors provided nucleation phenomena are not restrictive. The O-Ag system was chosen for the investigation into gas bubble-molten metal interactions due to the relatively high solubility of oxygen that enables rates of oxygen transfer to be measured from changes in bubble volume. Other advantages of this system include the absence of a stable oxide phase at an oxygen pressure of 1 atm and the relatively low melting point of the metal which permits the use of a metallic container, providing that it is resistant to oxidation. In those metallurgical processes where bubbles have an important influence, bubble volumes are usually greater than 5 cm3. For this reason the present study relates specifically to single large bubbles of oxygen rising in silver. These bubbles adopt the characteristic spherical cap shape similar to that shown in Fig. 1 for a 30 cc bubble rising in water. After an initial investigation to determine the velocities of inert (nitrogen) bubbles rising in molten silver, experiments were carried out with oxygen and the rates of mass transfer between the oxygen bubbles and the silver were measured. EXPERIMENTAL Apparatus. The apparatus, Fig. 2, for containing molten silver, was constructed from "Nimonic 75" Alloy (75 pet Ni, 20 pet Cr, 5 pet Fe, Mn) and provided for the release of single bubbles from an hemispherical cup, situated at the bottom of the column. The cup was turned by translating the rotation of the drive shaft through 90 deg. This was accomplished by the use of a bevelled gear system, and a smooth drive was provided by the lubricating action of the silver on the gears. Since reliable high temperature seals at 1000°C were found to be impracticable, the filling and drive shaft tubes were extended outside the 3.5 kw resistance wire tube furnace, where connections were made using easily accessible O-ring seals. The apparatus remained gas tight to the atmosphere at pressure differentials far in excess of those used. The filling tube was connected via a small bore tube to the differential pressure transducer. Gas could be bubbled into the inverted cup from two i-in. tubes which passed down the inside of the column to preheat the gas. The temperature of the silver was maintained at 1020°C during all experiments. Measurement of Bubble Volume. In order to calculate mass transfer rates, it was necessary to obtain a continuous record of the bubble's volume during its passage through the column of molten silver. The method adopted for measuring the bubble volume involved closing off the top gas space to the atmosphere prior to each experiment, and recording the variation in gage pressure of this space during the formation and rise of the bubble. Since any change in bubble volume results in an equal change in top space volume, Boyles Gas Law may be applied (for isothermal con-

Jan 1, 1970

By W. A. Tiller

The conventional techniques1 for determining the liquidus and solidus surfaces of an alloy system containing more than two components are extremely tedious to use and do not provide a complete picture of the equilibrium relations between solid and liquid alloys. These techniques are unable to yield "tie-line" information concerning solid and liquid phase equilibrium, a very important parameter in the solidification description of the liquid alloy and a very necessary parameter in the preparation of "zone-levelled" alloy crystals. The tie-line in a polycomponent system is analogous to the partition coefficient, k, in a binary system, it gives the composition of the solid, Cs, in equilibrium with a liquid of composition CL. That is, CS = kCL, where CL = [CO, c1,...Cn] denotes the concentrations of the n + 1 constituents, and k = [ko, kb ,...kon] denotes the gross partition coefficient for the elements between the two phases; thus, k is the tie-line in this system. To provide complete information concerning the two-phase equilibrium in a polycomponent system it is necessary to know both the liquidus surface and the gross partition coefficient. From these two the solidus surface is obtainable. The conventional techniques are unable to provide this information in other than a binary system so we must look elsewhere. In recent years considerable insight has been gained into the correct description of the liquid-solid transformation, 2,3 and controlled solidification experiments may now be designed to both facilitate and enhance equilibrium diagram studies. In the present paper consideration is given to two methods for obtaining the liquidus surface and the tie-lines in polycomponent systems. The methods to be described below deal with the solidification of an n + 1 constituent liquid alloy of initial composition [co, Co ,...Con], where the superscripts refer to the elements present and the Oth element is considered as the solvent in which the n solutes are dissolved. The general assumptions made in the treatment are the following: (i) the solid and the liquid at their interface are in equilibrium during the growth of the solid phase, (ii) there is no diffusion in the solid phase, and (iii) the liquid phase is completely mixed and is therefore homogeneous in concentration. METHOD I In this section a general experimental method will be described which, in principle, is capable of giving an exact description of the liquidus surface and tie-lines. Consider the unidirectional solidification of the liquid alloy specimens L cm in length (freezing either horizontally or vertically). Allow the sample to be frozen very slowly from one end, as indicated in Fig. 1, with complete mixing in the liquid, and then analyze the solid bar to determine its chemical constitution as a function of position, x, along the bar. Let Fig. 2 represent a possible distribution of the ith constituent along the bar. At the point x' the concentration of the ith constituent is C1/5(x'), and the average concentration of the rest of the bar between x' and L is given by C:(X' - L) where [ A(x) c1s(x)dx c1/s(x'-L) = f A(x) dx x * and A(X) is the cross-sectional area at x. In a similar manner. all the Cf(x' - L) may be determined. Thus, the gross artition coefficient k for liquid of composition [Cs(x - L),...Cs(X' - L)] is given by -ko= CUs')/Cos(x' - L)" k0 = csn{x')/cs(x'- L) During the freezing of a charge of length L, the liquid composition may vary over a wide segment of the phase diagram and the gross partition coefficient over this segment of the phase diagram may be determined from one, bar. Fig. 3 illustrates the magnitude of Cs(g)/C,' as a function of the fraction g of the bar which has Solidified.2 We can see that the concentration in the bar will vary over a range of about 15 wt pct for C,' = 10 wt pct and k0 = 0.5. It appears that 4 or 5 specimens would be adequate to study a simple binary eutectic phase diagram.

Jan 1, 1960

By H. Groeneveld, C. A. Connally, P. J. Hoenmans, J. J. Justen, W. L. Mason

A miscible displacement pilot using a slug of LPG driven by separator gas was conducted in the Cardiurn reservoir of the Pembina field. The injection pattern was a 10-acre, inverted, isolated five-spot. Upon completion of the LPG-gar phase, an experiment was conducted using a slug of water followed by gas. Calculated performance of the pilot is compared with actual performance. Equations are developed to calculate the distribution of LPG into zones of varying permeability, to estimate the progress of the flood at different times in the various zones and to estimate gas rates after breakthrough. The analysis indicates that permeability stratification was a dominant factor in controlling oil recovery and that oil was completely displaced from the swept pore volume. The results of the pilot indicated that miscible flooding is a practical means of pressure maintenance in this reservoir. The total recovery from the pilot area was good in spite of the early breakthrough of LPG. The effects of stratification were reduced by injecting a slug of water into the partially swept reservoir. INTRODUCTION The Pembina field,' located in Alberta, is the largest oil field in Canada and one of the largest in the North American continent. The reservoir is a stratigraphic trap producing from the Cardium sand. Neither bottom water nor free gas has been found. The recovery of oil by the natural depletion mechanism has been estimated at 12.5 per cent. Pressure maintenance studies of various areas have indicated that the recovery can be increased 21/2 times by water flooding, and a large area of the field is presently under water flood. However, reservoir studies of the North Pembina area indicated that miscible flooding might be competitive with water flooding. A pilot test was conducted to evaluate the performance of a miscible flood. A 10-acre, inverted, isolated, five-spot pattern was selected for the pilot. The pattern area was large enough to minimize wellbore fracturing effects and contained sufficient oil to provide significant working numbers. The performance of each of the four producers could be evaluated individually and compared. In the event of breakthrough in one direction, the effect would be isolated from the other producers. The use of a single injector minimized the volume of LPG required, and, because of the high mobility of gas, one well was sufficient to inject the necessary daily volume to replace the high rate of production. With four producers, the test could be completed in time for results to be evaluated, additional engineering studies to be made and a unit to be formed before the reservoir pressure in the North Pembina area declined below the bubble point. The pilot was located in an area developed on staggered, 80-acre spacing. The injection well was drilled at a regular location, while the four producers were drilled 467-ft north, east, south and west of the injector. Each quadrant and its associated producer were identified according to their direction from the injector— that is, north, east, south or west. The eight surrounding producers on 80-acre spacing were shut in to isolate the pilot area and provide for reservoir pressure observation. The pilot wells were completed using permanent-type completion techniques. After coring, casing was run through the pay section and cemented. Inside 51/2-in. casing, 2 1/2-in. tubing was hung. The wells were perforated opposite the Upper Cardium sand and lightly fractured. Fracturing volumes, rates and pressures were low to minimize the extent of the fractures. The fracturing treatments average 1,000 lb of 20-40 mesh sand in 700 gal of a low fluid-loss sand-carrying agent. Feed rates and wellhead fracturing pressures averaged 5.5 bbl/min at 2,535 psig, respectively. After fracturing, the productivity index was measured in each of the five pilot wells. The average PI of the four producers was 0.41 BOPD/psig drawdown. The measured PI'S were approximately the same as PI'S calculated from core analysis data, indicating that the fracturing treatments were just sufficient to overcome

By M. S. Rashid, C. J. Altstetter

Allotropic transformations in cerium have been studied by dilatometric, resistometric, X-ray diffraction, and metallographic techniques. The dilatometric study indicated that, on cooling below O°C, the high-temperature fcc phase, y, transforms partly to the hexagonal phase, ß, and, on further cooling, to the collapsed fcc phase, a. The amount of $ phase present at room temperature is increased by repeated cycling through the a-y transformation. It has been shown metallograPhically that the y-ß transformation has many characteristics of a martensitic transformation. In contrast to the y-ß transformation the ?-a transformation does not give the manifestation of a shear transformation. Small cellular ? domains of random shape and size collapse to a in a short time with no apparent coordination with neighboring domains. The considerable confusion in the literature over the existence of more than one high-temperature fcc phase is discussed. Two such phases have been reported in the literature and an attempt is made in this study to clarify the situation. Twelve fcc and two hcp structures have been shown to be easily reproduced or eliminated. It is proposed that the two "additional" allotropes reported in the literature and fourteen of the phases detected here are not allotropes of cerium but are due to contamination. CERIUM exists in several allotropic forms, but there is some disagreement over what the forms are. Furthermore, the conditions favoring the presence of a particular allotrope and the nature of the transformations from one form to another are uncertain. The objectives of this research were 1) to ascertain the allotropic forms of cerium, 2) to establish the conditions under which the allotropes exist, 3) to study the effects of annealing and thermal cycling on the allotropic transformations, and 4) to study the transformation mechanisms. Dilatometric, resistometric, metallographic, and X-ray diffraction techniques were employed. The form of cerium commonly found at room temperature is fcc and is designated ?. A complex hexagonal phase, 8, forms when y is cooled to slightly below room temperature. At still lower temperatures the y fcc structure transforms to an fcc form with a much smaller lattice parameter, termed a cerium. A bcc form, 6, which exists just below the melting point (800°C), will not be considered further in this work. There is a substantial body of experimental evidence (reviewed by Gschneidnerl) which favors the acceptance of these four allotropes, though some investigators have tried unsuccessfully to observe the ß hexagonal form.'-' There is disagreement, however, over the phase-transformation temperatures, due, in part, to broad hysteresis and overlapping of the transformations between the a, ß, and ? forms. The transformations are also sensitive to prior thermal and mechanical treatment. The differing purity of cerium used by different investigators is undoubtedly a factor. Cerium is difficult to separate from other elements and is quite reactive, igniting spontaneously when it is filed in air. The highest purity of cerium to date is reported to contain several hundred parts per million by weight of impurities, and early investigations were carried out on cerium containing several percent of impurities. There have been reports of more than one fcc allotrope at room temperature. Gschneidner, Elliott, and McDonald5 obtained diffraction patterns of an fcc phase with a lattice parameter about 1 pct less than that of the ? phase, instead of the y phase, on slowly cooling cerium filings from 23° to -198°C and warming them back to room temperature. However, when the sample was heated to 447°C and cooled to room temperature it consisted of only the ? phase. They have designated this new fcc phase "a-? intermediate", and say it is quite sensitive to impurities. After prolonged high-temperature treatment of a powder specimen, Weiner and Raynor2 obtained a diffraction pattern of an fcc phase of lattice parameter about 1 pct less than the ? phase. This they called the y' phase. It could not be reconverted to the y phase and is claimed to be different from the a-? intermediate phase.5,6 Dialer and Rothe3 reported two fcc phases* after cycling their powder specimens between room temperature and -192°C. Gschneidner, Elliott, and McDonald5 suggested that one of the fcc structures obtained by Dialer and Rothe was equivalent to their "a-? intermediate" phase. Table I presents some pertinent data on the proposed allotropes. For the ?(fcc)-ß(hexagonal) transformation

Jan 1, 1967

By David Cornell, D. L. Katz

Procedures for computing turbulent flow of gas in steady state near the well bore and a graphical method for predicting unsteady state laminar flow at distances from the well have been combined to compute pressure gradients in gas reservoirs. Methods are discussed for predicting single and multiple transients at constant flow rate in an infinite reservoir, predicting constant How rate. in a finite reservoir, reproducing and interpreting back pressure test data, and prediction of the behavior of a closed-in gas well. An example of the graphical method is given for a single transient and the results are compared to a published analytical solution. A typical calculation of the pressure gradient in a reservoir and the production of gas is made .starting with data from a hack pressure test. INTRODUCTION The calculation of pressure gradients throughouse a ga-reservoir at any point in its production history is a complex problem involving turbulent flow near the well bore and unsteady state flow of the gas from the reservoir. The problem is complicated further by the variety of boundary conditions that may be imposed upon the flow. Such boundary conditions include finite or infinite reservoirs. constant or variable rates of production, constant or varying bottom hole pressures, complex initial pressure distributions throughout the reservoir, and pressure maintenance through cycling. In addition to these problems there are the practical difficulties associated with natural gas reservoir analysis These might include: (1) variation. in the permeability and porosity throughout the formation. (2) variations in the thickness of the formation, (3) communication of the producing formation with other producing zones. (4) radial variations in permeability due to influx of drilling fluid. acidizing. or the buildup of liquid in the formation around the well, (5) partial penetration of the producing zone. (6) water flood, and many others. These problems will be considered eventually. In the mean- time, it is necessary to present methods of handling the case of radial flow through a homogeneous, regular, producing stratum to a single well uncomplicated by other factors. The idealized case of steady state, radial flow from a natural gas reservoir was studied by Elenbaas and Katz5 and by Mac-Koberts.' An analysis of the unsteady flow of gases through porous media has been made by Aronofsky and Jenkins' for the one dimensional laminar case with the boundary conditions of constant downstream pressure and uniform initial pressure. Solutions for the partial differential equation for unsteady state radial flow of liquids through porous media have been given by Van Everdingen and Hurst" for the constant bottom hole pressure and constant production rate cases for laminar flow in a reservoir initially at a constant pressure. The analogous heat transfer equation has been treated by Perry and Berggren9 for the case of quenching a cylindrical hole in an infinite medium to a constant temperature. The back pressure curve, which involves the variation of the bottom hole pressure as it is determined by the behavior of the reservoir as a whole. has been studied by Rawlins and Schellhardt,10 Binck-lev,' Baumel and Breitung,' and others. No complete. adequate treatment has been given for the unsteady state, radial flow of gases into a cylindrical hole from a porous medium with or without the presence of turbulent flow and for any set of boundary conditions. If one relies only on analytical solutions of the basic equations, a new solution must be obtained for each new boundary condition. Furthermore, the complexity of the mathematical expressions for boundary conditions other than the very simple cases limits the use of analytical solutions. Graphical methods exist, however. that are general in nature, accurate, and readily employed. The analysis of steady and unsteady, laminar and turbulent, radial flow for various boundary conditions with application to specific natural gas well problems by means of graphical procedures forms the content of this paper. STEADY STATE RADIAL FLOW EQUATION FOR LAMINAR AND TURBULENT FLOW A steady state radial flow equation for laminar and turbulent flow through unconsolidated sands has been given by Muskat as discussion of the work by Elenbaas and Katz.5 This equation is based on the properties of unconsolidated

Jan 1, 1953

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey' have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks' have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.' Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method

By M. B. Kasen

eralized strain equation appear quite different they are really identical. This identity can be shown in a simple mathematical rearrangement. Referring to Eq. [I], the substitution of ln(1 +?E) for ? (where ?E is the engineering strain) yields: a =s8- (s8?so)e?1/?c) ln(1 +s8??E) [4] or s =s8?(s8?so)(l +?E)-1 [5] Replacing (1 + cE) with its equivalent, (l/lo): s=s8?(s8?so)l/€c [6] Adding and subtracting the quantity (s8 - so) yields: s=s8?(s8?so)+ (s8-so) l-(^)l/£c [7] = s0+(s8?so)ri-(^0)1/^l [8] Letting n = 1/?c it follows that: This equation will be seen to be identical to the generalized strain expression in Eq. [3], the coefficients a and b being equal to so and (s8?so)/?c or n(s8?so), respectively. Confirmation of the relationship noted above was provided by fitting experimental stress-strain data5 for 304 stainless steel (total strain from 0.072 to 1.15) at room temperature at a true total strain rate of 4 x 10-3 sec-1 to Eqs. [I] and [3]. In this study nonlinear regression analyses were used to evaluate the equation constants. The results of this evaluation are presented in Table I. When these constants are used in Eqs. [I] and [3] the curves obtained are coincident. It will also be noted that the value of n is equal to YE,, the value of a is equal to so, and the value of b is equal to (s8 - so)/Ec. These equalities are illustrated in Table I where the Voce constants in parentheses were calculated from a, b, and n values. Excellent agreement is seen to exist. It is concluded, therefore, that the Voce equation is, indeed, identical to the linear expression involving generalized strain. In applying the generalized strain concept to stress-strain data, it has been noted3 that data within what was termed the transition region (first few points in the plastic flow region) must be excluded from the analysis in order to obtain the linearity in Eq. [3]. A similar exclusion applies in the use of the Voce equation since an expression which would describe the strain-hardening region would usually be found to be not too effective in the elastic-plastic fillet2 which connects the elastic line to the strain-hardening region. These considerations further substantiate that Eqs. [I] and [3] do indeed describe identical behavior. Note added in proof. Since this manuscript was submitted for publication, correspondence with Professor S. R. Davies, University of Edinburgh, (of Ref. 3) indicated that his associates have also confirmed the identity between the Voce equation and the generalized strain concept. It is sometimes desirable to obtain a small grain size in super-purity metals. Unfortunately, the absence of impurities, the low recrystallization temperatures, and the low activation energies for boundary motion result in rapid growth of relatively few grains when conventional recrystallization practices are used. For example, Arajs et al.1 were unable to obtain a mean grain diameter less than about 0.4 mm in iron containing 38 at. ppm of impurity. The <1 at. ppm purity range in many modern superpurity metals results in yet larger limiting grain sizes. It is the purpose of this communication to describe a technique by which mean grain diameters as small as 30 µ have been produced in aluminum having an impurity content < 0.5 at. ppm. The technique should be applicable to superpurity metals other than aluminum and has potential applications to grain refinement of commercial alloys. The necessary and sufficient requirements for production of a small grain size are as follows: a) a large number of potential nucleating sites must be created randomly within the cold worked material, b) sufficient thermal energy must be provided at a rate that will simultaneously activate a large fraction of the potential nuclei, and c) the grain growth must be stopped shortly after grain impingement.&apos; A high density of potential grain nuclei may be cre-

Jan 1, 1970

By D. G. Westlake

Alloys of poly crystalline Nb-H have been tensile tested at 77" and 120°K after slow-cooling and after quenching from room temperature. A rationale has been developed to explain the effect of cooling rate on ductility in terms of microstructure. The inferior ductility of the slow-cooled alloys in this temperature range can be attributed to the precipitation of coarse hydride particles that fracture in the early stages of plastic deformation. Quenching Produced much smaller hydride particles and a decreased tendency for initia-tion and popagation of cracks. Prismatic dislocation loops, punched into the matrix of quenched alloys, eliminate the yield point and may Play some role in the low temperature ductility of the quenched alloys. The difference in ductility between quenched and slow-cooled V-H alloys was even greater than that observed in the N b-H alloys. It is well known that the ductility of Nb-H and V-H alloys decreases with increasing hydrogen concentra-tion and with decreasing temperature.1-l1 The follow-ing are some of the numerous mechanisms that have been adopted in attempts to explain these observations: 1) hydrogen-dislocation interarction,4&apos;6 2) lowering of the crack surface energy by adsorbed hydrogen,9&apos;10 3) localized lowering of the niobium bond strength,9 and 4) exertion of pressure inside voids by molecular hydrogen.6 A recent study1&apos; indicates that, in some way, the embrittlement is related to the precipitation of a hydrogen-rich phase, because embrittlement and hydrogen solubility have the same temperature dependence. Most of the investigators observed that at very low temperatures there was a return of the ductility. Thus, hydrogenated niobium and vanadium may be said to exhibit a low temperature ductility minimum."1&apos;8&apos;9 It has been reasoned by most observers that diffusion is required for the embrittlement mechanism to be operative, and that as hydrogen becomes immobile at very low temperatures, the ductility increases. This investigation of the mechanical behavior of Nb-H al-loys under various conditions is an attempt to gain a better understanding of the low temperature ductility minimum. EXPERIMENTAL Tensile specimens, with gage length dimensions 20 by 3 mm, were machined from sheet stock" and both faces were ground on wet 600 grit Sic paper to obtain a final thickness of 0.5 mm. Annealing as described earlier12 resulted in average grain diameters of 0.3 mm. The recrystallized niobium tensile specimens were electropolished in 90 parts H2S04, 10 parts HF (48.9 pct) by volume, at 287°K and 30 v. Hydrogen, which may have been present in the starting material and/or added during specimen preparation, was removed by annealing at 1073°K in a dynamic vacuum of 2 X 10-6 torr. Hydrogenated specimens (1.15 at. pct H) were prepared by reaction with a known quantity of pure hydrogen (obtained by thermal decomposition of UH3) at 1073oK, followed by cooling at the rate of 100" per hr. These were vacuum encapsulated, homogenized at 573°K for 3 to 5 hr, and furnace-cooled. Tensile tests were conducted on an Instron machine at a strain rate of 4 x l0-5 per sec. The gripping device was such that the specimen was positioned below the crosshead. Thus, subambient test temperatures could be obtained by immersion of the device in a cooling bath of either liquid nitrogen or Freon-12. The temperature of the Freon-12 was controlled within ± 0.1°K of the desired temperature by automatically monitoring the flow of liquid nitrogen through a copper coil in the bath. Slow-cooling to 77°K was accomplished as follows. The testing device was surrounded by a rubber jacket, open at the top and with small holes in the bottom. Liquid nitrogen, brought up to the bottom of the jacket, entered the holes and vaporized quickly. The temperature was monitored at the specimen with a copper-constantan thermocouple, and a constant cooling rate of 2 deg per min was easily achieved by controlling the level of the liquid nitrogen outside the jacket. RESULTS AND DISCUSSION Low Temperature Ductility Minimum. Resistometric studies of both Nb-H12 and V-H13 indicated that at tem-peratures less than 150°K, virtually all of the hydrogen

Jan 1, 1970

By M. M. Mebta, G. W. Dean, W. H. Somerton

With the advent of underground heating operations, interest has developed in the alteration of rock properties by high-temperature treatment. In the present work a number of sandstones were heated to temperatures in the range of 400°C to 800°C under both atmospheric and simulated reservoir pressures. Pertneabilities increased by at least 50 per cent and sonic velocities and breaking strerlgths decreased by an equivalent amount. Differential thermal expansion and other reactions of constituent min-era1 grains are the causes of these alterations. INTRODUCTION In the underground combustion of petroleum reservoirs, temperatures of the order of 600C are reported to have been reached in the combustion zone.&apos; At this temperature rocks are subject to extensive thermal alteration. Temperatures of this magnitude and higher may also occur in subsurface formations when subjected to bottom-hole heating, thermal drilling operations, and underground nuclear explosions. Temperatures of this magnitude might also be generated by conventional rock drilling methods at points of bit-tooth contact. In earlier work, the permanent deformation of rocks resulting from heating was reported. Major structural damage of rocks occurs due to differential thermal expansion of mineral constituents. A number of mineral alterations including crystal inversions, loss of water of crystallization and dissociation, may also contribute to changes in physical structure and properties of rock. In the present work, samples of three typical sandstones were heated to several temperatures up to a maximum of 800C and then allowed to cool to room temperature. Heating was done under both atmospheric pressure and simulated reservoir pressure conditions. Physical properties of the samples were measured before and after heating and comparisons made. Measured properties included permeability, sonic velocity, breaking strength and fracture index. Changes in physical properties were compared to changes in mineralogical characteristics as determined by thin-section. X-ray diffraction and chemical tests. EQUIPMENT AND PROCEDURE Two outcrop sandstones (Bandera and Berea) and one sub-surface sandstone (St. Peter) were selected for the tests. These samples have a wide ranee in composition and physical properties as shown in Table 1. The first series of tests was made on 2-in. diameter by 5-in. long test specimens. Test specimens used in all later work were 3/4-in. diameter by 1 1/8-in. long, this being the specimen size required for heating at simulated reservoir pressures. After careful washing, the cores were oven dried at 100 ± 5C for a minimum of 24 hours before the tests were run. Test specimens were heated in an electric furnace at a constant rate of temperature increase of 3C per minute. When maximum temperature of the run was reached, the sample was allowed to soak for one hour. The furnace was then cooled to room temperature at the same rate of 3C per minute. The entire heating operation was designed for reproducibility without subjecting the test specimens to excessive thermal shock. For samples heated under simulated reservoir pressures, a pressure cell designed by Dean was used (Fig. 1).3 The core sample was inserted into a thin-walled (0.006 to 0.01-in.) copper cup which was then mounted in a high-pressure cell. Provisions were made for the application of internal pore pressure as well as confining pressure. Tests showed that the thin-walled copper cup closed tightly around the core and satisfactorily transmitted confining pressure to the core. The core was heated by placing the entire cell into the electric furnace. The heating program was the same as that used in the atmospheric pressure runs: tempera-ture rise of 3C per minute to maximum temperature of the test, soaking at maximum temperature for one hour, and cooling at a rate of approximately 3C per minute. The cell was designed to withstand 5,000 psi at 1,000C. However, since it was considered likely that repeated heating and cooling would in time weaken the steel, 2,000 psi at 850C was set as a working limit. In the present series of tests, the pore pressure was held constant at 750 psi and the confining pressure at 1.500 psi. The pressure source was a high-pressure nitrogen tank. The two pressures were controlled manually but are accurate well within ± 50 psi.

Jan 1, 1966

By L. L. Handy

The development of a theory for miscible liquid displacement requires evaluation of the variables which affect growth of the mixing zone between solvent and displaced oil. Factors which appear to be important are individual fluid viscosities, viscosity ratios, flood rate, fluid densities, flow characteristics of the porous medium and molecular diffusion coefficients of the fluid components. The primary purpose of this paper is to evaluate diffusion effects. Theoretical treatments to date have been limited to floods for which the viscosity ratio is one. Two principal theories have been proposed. Von Rosenberg adapted for porous media a theory derived for capillary tubes by G. Taylor.&apos; ,&apos; In this theory molecular diffusion perpendicular to the direction of flow is a primary factor in sharpening the flood front. Slow floods give sharper fronts for a given distance traveled than fast floods. An alternative theory considers miscible liquid displacement as a statistical problem.3,4,5 Diffusion is not an important factor in this theory, but it leads to the same general type of equation as von Rosenberg&apos;s. Both theories predict S-shaped concentration profiles with the same dependence on distance traveled. The statistical or "dispersion" theory predicts rate independence, however. To supplement rate studies a direct measurement of a diffusion effect would be helpful in evaluating which of the two proposed mechanisms best describes miscible liquid displacement for one-to-one viscosity ratio systems. No quantitative theory has been proposed for floods in which a low-viscosity fluid displaces a high-viscosity fluid. It might be anticipated, however, that the extensive fingering observed in floods with adverse viscosity ratios would increase opportunities for an exchange of components between displaced and displacing liquids by a diffusional process. Even if molecular diffusion were not an important mixing mechanism for one-to-one viscosity ratio systems, it could be significant in systems with adverse viscosity ratios. METHOD FOR EVALUATING DIFFUSIONAL MIXING In solvent flooding the displacing liquid can become mixed with the displaced liquid by a number of processes; in general, either mechanical or diffusional in character. To evaluate the contribution of diffusion, a method was sought that would distinguish diffusional mixing from mechanical mixing. Suppose two soluble tracers are added to the displacing liquid and that one tracer has a diffusion coefficient much greater than the other. Then, if diffusion is important during miscible liquid displacement and if diffusional transport is primarily transverse to the direction of flow, substances with high diffusion coefficients will have sharper concentration profiles than those with comparatively low molecular diffusion coefficients. In order to have material balance the profiles for the two tracers must intersect. For one-to-one viscosity ratios the intersection would occur at the 50 per cent concentration point. It is generally assumed that diffusion in the direction of flow is negligible. This assumption appears reasonable because flow velocities are ordinarily orders of magnitude greater than diffusional velocities. If this is not a valid assumption, however, the shapes of the fronts will be further degraded by longitudinal diffusion. The effects of rate and molecular diffusion coefficients are the opposite for longitudinal diffusion to those for lateral diffusion. Higher rates or lower diffusion coefficients would tend to give sharper fronts. The double tracer method offers several unique advantages in evaluating effect of molecular diffusion. First, the method gives a direct measure of the effect of diffusion on mixing zone size. Second, the effect of a difference in diffusion coefficients is determined from a single experiment. Comparison of several experiments at different rates is not required. Third, the method is applicable for any viscosity ratio. Fourth, the method is unaffected by density differences between the two fluids which might result in rate dependent gravity effects. EXPERIMENTAL PROCEDURE Suitable tracers are substances whose diffusion coefficients differ as widely as possible, but whose adsorptions on rock surfaces is negligible. The two tracers selected were methanol and sucrose in water solutions. The approximate diffusion coefficients for methanol and sucrose, as given in the International Critical Tables, are 1.3 X 10-5 and 0.3 X 10-5 m2/sec, respectively. Concentrations in the water solutions were 200 gm/liter of solution for both methanol and sucrose. These high

By R. G. Ward, A. E. El-Mehairy

A technique was developed to calculate the density of liquid metals from the profile of a weighed levitated drop obtained by emitted light photography and calibration. The density of liquid copper was determined over the temperature range 1370" to 2100°K and expressed with the equation: The molar volumes and thermal coefficients of expansion are calculated. Copper, upon melting, was found to expand 3.34 pct. A knowledge of the density of liquid metals is essential as it is a parameter appearing in most theories of the liquid state. The liquid density is also necessary for the calculation of the commercially important contraction occurring on solidification. Unfortunately good density values are rarely available and the present work was undertaken to establish a new method of density determination free from the limitations of the existing methods. For liquid metals with high melting points, the variation and the limited number of density values reported in the literature have been caused primarily by difficulties in technique, a matter which has been adequately discussed by Urbain.&apos; In all the methods used for measuring the density of the high-melting liquid metals there are sources of error, and corrections have to be made. For example, in the method of balanced columns, the formation of gas bubbles in one or both arms was a great source of error. In the pyknometer3 and immersed-sinker methods, corrections of certain volumes for thermal expansion to the operating temperature exceed 4 pct. In the maximum-bubble method5 a correction for the thermal expansion of the bubble-tube is used (about 1.6 pct); a second correction takes into account the change in the volume of liquid displaced by the bubble-tube (about 3.5 pct). Another method used is based on the shape of a drop of liquid, at rest on an inert refractory plaque, obtained photographically or radiographi-cally. The volume is calculated from the shape, which is sensitive to the surface tension and supporting material, by some intricate empirical relation. The drop contour changes with 90 deg rotation of the support.7 The application of all methods is limited by reaction between liquid metal and the apparatus. The successful stabilization of a mass of liquid metal during levitation melting8 afforded a technique for measuring the density of liquid metals out of contact with solid materials. The technique was developed by using copper as the density of liquid copper has already been studied extensively, thus allowing comparison with the new method. It consists of calculating the volume from the dimensioned profile of the levitated drop obtained by emitted-light photography and careful calibration. The molten egg-shaped drop spins around the magnetic axis of the coil (lined vertically) at a speed which seems to increase with increasing temperature. The method does not involve corrections and the experimental setup allows the use of vacuum or special atmospheres. The maximum operating temperature is limited by excessive evaporation of metal from the droplet. EXPERIMENTAL Levitating System. A stabilizing ring was used in conjunction with a 3-turn levitating coil energized by a 450-kc per sec, 10-kw Toccotron with a high current-low voltage output stabilized to f 1 pet.&apos; APyrex glass tube with a Pyrex Optical flat window inserted through the coil and ring allowed the use of special atmospheres during levitation melting, as shown in Fig. 1. Optical System. Fig. 1 shows the optical arrangement used to photograph the profiles of the drop, parallel and perpendicular to the vertical spinning axis, simultaneously on the same frame of a high-speed 35-mm film. The optical system was adjusted to avoid any distortion of the profiles and to ensure that both profiles are in focus. The length of path from object to objective was about 30 cm, and a magnification factor close to 1 was obtained with the single-lens reflex-camera system. Temperature Measurement and Control. In operation, the specimen temperature was constant to 5 C for a fixed power input to the levitating coil and fixed rate of gas flow through the tube. The temperature was measured by an automatic two-color pyrometer (Coloratio) checked at the melting

Jan 1, 1963

By H. H. Klepfer, M. E. Snyder

UNTIL very recently, the commercial availability of the rare earths as metals has been very limited. Fabrication of mill products from these metals has not been studied in most cases. This note reports the results of the development of fabrication techniques for thulium. Thulium has a melting point of about 1550°C and a hexagonal close-packed crystal structure. It oxidizes in air to give a black oxide (TmzO,). The procedures for producing thin thulium foil were developed on one ingot weighing about 220 g and were subsequently applied in processing about a pound of metal to foil. Several alternates for the various fabrication steps were investigated and will be discussed. The metal fabricated was in the form of commercial chill-cast ingots 1 in. in diam and 2 in. long weighing approximately 220 g. Impurities in the ingots were reported by the vendor to be 4000 ppni tantalum, 2000 ppm calcium, 200 ppm nickel, 100 to 200 ppm iron, 100 ppm europium, and less than 100 ppm copper, lutetium, and ytterbium. In addition to these impurities, several salt-like inclusions as large as l/8 in. in diam were revealed along the center line of the one ingot sectioned. Preliminary tests indicated that small wafers cut from the as-cast ingot would not fabricate readily by rolling. Forging of copper-jacketed wafers was therefore attempted. At 1550"F forging was satisfactory but an apparent reaction of copper with thulium demanded investigation of lower temperatures. Therefore, the remainder of the test ingot was forged at 1450°F—with only minor cracking. All ingots forged were inserted into copper tubes of 1 in. 1D and 0.125-in. wall thickness. The jackets were sealed by flattening the ends of the tubes and welding under helium. Heating time at 1450°F was 30 min. Press forgings of 1/8 in. per pass were used, followed by 10 min reheats. When the ingots had been squared and reduced to 0.250 in. in thickness, the original copper jacket was stripped off and replaced by a new jacket in preparation for hot rolling. Hot rolling at 1450" F without edge cracking was readily accomplished after forging. Excellent results were obtained with 10 pct reductions of thickness followed by 5 to 10 min reheats. After reduction of thickness from 0.250 to 0.100 in. the copper jacket was removed. It was found, in fact, that hot rolling in air was possible. A tenacious black oxide similar to that seen on zirconium was formed during 3 min reheats at 1450°F. Reduction in air to 0.010 in. foil was possible taking 10 pct reductions per pass. The oxide coat formed during hot rolling in air could best be removed by sand blasting and pickling. Common pickling solutions containing polar solvents were found to attack the metal too rapidly and a concentrated nitric-hydrofluoric acid mixture attacked neither the oxide nor the metal. The most satisfactory pickling solution was 52 vol pct concentrated nitric acid-48 vol pct glacial acetic acid. After forging to 0.250 in., vacuum annealing and cold rolling was found to be another satisfactory alternate to hot rolling in copper jackets. After forging. a hardness of Rockwell B76 was found. Annealing in vacuum (2 X l0-5 mm of Hg) for 1 hr at 1200"F did not alter this value. Annealing at 1470"F for 1 hr brought the hardness down to R;]63. With cold rolling (5 pct per pass) the hardness returned to about R,1:76 after 20 pct reduction and edge cracking became noticeable. However, cold rolling to a total of 40 pct reduction in thickness (RI,,83) was possible before edge crack propagation became serious. Good surface finish was obtained, and the metal loss due to oxidation was minimized by cold rolling and vacuum annealing. Using this procedure the yield of 1.25 in. wide by 0.010 in. thick foil from a 1-in. diam ingot was about 40 pct.

Jan 1, 1961

By H. W. Hitzrot

*Original Pages Missing From Book

Jan 5, 1954

By P. Chiotti, J. T. Mason

Thermal, X-ray, metallographic, and resistivity data were obtained to establish the phase diagram for the Yb-Zn system. Six intermetallic compounds and five eutectics have been identified. The compounds position at 695°C. No evidence for the compound YbZn,, was observed. Whether this compound is stabilized or inhibited by small amounts of impurities is uncertain. SEVERAL aspects of the Yb-Zn system have been reported in the literature. Much of this work deals with the crystallographic identification of some of the compounds. The YbZn compound has been reported by Iandelli and palenzonal and the YbZn, structure and thermal data for the composition range between YbZn and Yb3Znll have been reported by Michel et al. A compound of uncertain stoichiometry in the range 60 to 68 wt pct Zn with a cubic structure has been reported by Kuz&apos;ma et aL4 The YbzZn17, YbZn,,, and YbZn13 compounds have been reported by both Kuz&apos;ma et a1 and Iandelli and Palenzona.&apos; YbnZn17 has also been identified by Laube. This investigation was undertaken to establish the phase diagram of this system. Ytterbium in many respects behaves as an alkaline earth rather than a rare earth metal6 and it is of interest to note the resemblance between the Yb-Zn system, the Ca-Zn,&apos; and typical rare earth-zinc systems.g&apos;10 MATERIALS AND EXPERIMENTAL PROCEDURES The materials used in this investigation were Ames Laboratory ytterbium and Bunker-Hill slab zinc, 99.99 pct pure. Analyses of the ytterbium showed 38 ppm Ta, 82 ppm 0, and 15 ppm N with Ca, Cu, C, Mg, Al, Si, Cr, Fe, Ni, Sc, Y, Er, Tm, and Lu all less than 50 ppm. The ytterbium was dehydrogenated by heating in tantalum under vacuum for 8 hr at 900°C. The metal was kept under argon except when being cut with shears and when being weighed. Tantalum was found to be a suitable container for alloys in the composition range between ytterbium and Yb2Zn17. With higher zinc concentrations the tantalum was penetrated by the alloy and failed after several days at temperatures of 600" to 700°C. Approximately 20 g of zinc and ytterbium metal were sealed under argon in *-in.-diam tantalum cruci- bles equipped with thermocouple wells. The crucibles were in turn sealed in stainless-steel jackets. The alloys were equilibrated in an oscillating furnace and given various heat treatments in a soaking furnace. The steel jackets were removed and differential thermal analysis data were obtained under an atmosphere of argon. The thermal analysis apparatus was calibrated against the melting points of pure zinc, aluminum, silver, and copper and found to give thermal arrests within 1°C of the accepted melting temperatures for these elements. The thermal analysis apparatus has been described in an earlier paper.g Alloys were further heat-treated as necessary and subjected to metallographic examination and X-ray diffraction analysis. Most alloys were etched with Nital (0.5 to 10 pct nitric acid in alcohol). EXPERIMENTAL RESULTS The results of the thermal analysis are summarized by the points on the phase diagram, Fig. 1. Six compounds and five eutectics were observed. The eutectic temperature and compositions in weight percent zinc are 491°C and 12 pct, 645°C and 32 pct, 687°C and 54 pct, 742°C and 67 pct, and 750°C and 77 pct. The compounds YbZn, YbZn2, Yb2Zng, Yb2Zn17, and YbZnll melt congruently at 650°, 751°, 75Z°, 754", and 755"C, respectively. The compound Yb3Znll undergoes a peritectic decomposition at 695°C. The transformation and melting temperatures for the ytterbium were observed to be 797" and 817°C. re-

Jan 1, 1969

By F. F. Poland

The purpose of this paper is to describe the recent applications and improvements made in the process and equipment for the recovery of metallic zinc from secondary metals by means of high temperature electric resistor furnace distillation. The process and high temperature furnace, as developed by Revere Copper and Brass Incorporated, Research and Development Department, was described in a previous paper. † Since the publication of that paper the process and equipment have been licensed to others. Consequently its field of application has been expanded, particularly for the treatment of scrap brass, galvanizers dross and scrap zinc. Other applications are being considered, particularly the melting of metallic titanium and special alloys requiring a controlled atmosphere. Consideration is also being given to modification in design that would allow the melting of copper cathodes for the production of oxygen-free and electrolytic tough pitch copper shapes without the necessity of blowing and poling. In the design of a furnace for melting copper cathodes it is intended that the fundamental principle of charging the solid material into a large pool of molten metal through a double door charging vestibule as shown in Fig 3 of the previous paper, † would be adhered to and, in addition, either dehydrated and desulphurized producer gas or commercial nitrogen, that is, approximately 96 pct nitrogen, 2 pct hydrogen 2 pct CO, also dehydrated, would be used to provide a protective atmosphere for the resistors and to flush out more soluble gases from the molten copper. The operation of a small unit producing approximately 1000 lb of copper per hr has shown that the commercial nitrogen will allow the production of tough pitch electrolytic copper when the metal is transferred from the meltdown unit contemplated to a special type induction unit where the oxygen absorption and temperature can be controlled. Standard practice in treating low grade scrap brass has been to oxidize the zinc by blowing air through the molten metal. This operation results in the recovery of approximately two-thirds or less of the zinc content of the brass in the form of an impure zinc oxide having only a fraction of the value of the zinc in metallic form. The balance of the zinc is lost to slag, fume and drosses. Fig 1 shows the revised flow sheet of the process for the treatment of scrap brass. For this application and as shown on the flow sheet a low frequency induction furnace has been substituted for the Wilkins-Poland type as a melting unit. This type of furnace has been found the most suitable for the purpose as will be shown by the following discussion. Serap Brass The melting of refinery scrap brass preparatory to distillation presents a melting problem quite different from the usual melting of brass mill scrap for recasting because refinery scrap is a heterogeneous mixture of alloys each having different melting and boiling points, and coated with a variety of corrosion products and foreign materials. In addition, brass and copper plated pieces of sheet iron, stainless steel, ceramic materials, and others are present that do not melt except at temperatures above the boiling point of the liquid bath. In the previous paper it is noted that difficulty was encountered in the melting of refinery grade scrap brass without excessive loss of zinc by volatilization and that 5 to 10 pct razorite flux was used to alleviate this condition. While this flux was successful in preventing the excessive loss of zinc by volatilization it interposed a low heat conductivity layer between the source of heat and liquid charge and resulted in lowering the melting rate of the furnace. This same problem presents itself in all types of furnaces in which the metal is heated by radiation and conduction from above the charge. It has been found in melting this type material that the melting rate per hour or per square foot of hearth area

Jan 1, 1950

By R. E. Goodman, M. A. Mahtab

The, state of stress around a vertical wellbore in rock following nonlinear stress-strain laws is examined by means of finite element analysis. The wellbore, is considered an axisymmetrier body with axisymmetric loading. The initial vertcal and borizontal stresses are "locked" in the rock elements around the wellbore and a new state of stress is generated by the displacements which occur around the borehole. A point-uise variation of the elastic moduli is made on the basis of the new stress state and the triaxial data. The initial stresses tire now reintroduced along with tile changed moduli and original boundary constraints. This procedure is repeatrd until convergent stresses are reached. The effect of nonlinearity on streases is examined for a 6,000-ft wellbore in a scbistose gnezss and Berea sandstone using results of laborratory triaxial compression tests. The results sbou that the effect is restrited to one well radius from the bottom periphery of the bole. Beyond a distance of one-quarter radius, the elfect of nonlinearity on stresscs IS almost always less than percent for the cases conszdered, The consideration of a static pressure inside the well does not magnify the effect of monlinearity on borehole stresscs. INTRODUCTION The terms " wellbore" and "borehole" here designate cylindrical openings in the ground with vertical axis and a circular cross-section. A knowledge of the stress redistribution that occurs on excavating a wellbore is important in understanding (1) the behavior of the lined or unlined hole, (2) hydraulic fracture response, and (3) the effect of stress redistribution on drillability; also it is important in predicting initial stresses in the virgin ground, and in analyzing the response of measuring instruments placed in the borehole. Our knowledge of the state of stress around a wellbore has been restricted to homogeneous, isotropic, elastic material and derives chiefly from the analysis by Miles and Toppingl and the photoelastic work of Galle and Wilhoit2 and Word and Wilhoit.3 In this investigation the state of stress is examined for a nonlinear elastic material by means of finite element analysis. Many rocks possess stress-strain curves that depart notably from straight lines in their initial or final portions. While the literature contains abundant stress-strain data from triaxial tests (axisymmetric loading) on cylindrical rock specimens, there is little information on rock deformability under nonaxisymmetric loading conditions such as occur at each point around the bottom of a wellbore. Although there is some knowledge of the effect of intermediate principal stress on rock strength, there is virtually nothing known about its effect on rock deformability; therefore, we have assumed here that the effect of intermediate principal stress can be ignored. A schistose gneiss4 and Berea sandstone5 were selected as representative rocks for this analysis. The traditional graphs of deviator stress (ol -o3) vs axial strain were reworked to give the tangent modulus as a function of the deviator stress for varying values of the minor principal stress. The result is a nesting family of skewed, bell-shaped curves for tile gneiss (Fig. 1A) and the sandstone (Fig. 2A). A similar replotting of the lateral strain data defines the variation of Poisson&apos;s ratio (v) with the deviator stress and confining pressure. These curves, shown in Fig. 1B for the gneiss and in Fig. 2B for the sandstone, are not so well ordered as the tangent modulus curves. However, all of these display an increase of v with deviator stress application, but the rate of increase diminishes with confinement. The ET and v curves for the two rock types are tabulated in Tables 1A and 1A for use in a digital computer so that material properties corresponding to a given state of stress can be assigned by interpolation. STATEMENT OF THE PROBLEM Given a state of stress in the ground prior to drilling a hole (initial stresses), we assume that

Jan 1, 1969

By Leonard L. Melton, Calvin D. Saunders

The design and performance of many operations common to the petroleum industry depend upon the unique properties of a class of materials known as non-Newtonian fluids. The art of measuring and describing these properties, called rheology, has fallen far behind the field applica-tion of these fluids in this industry. Empirical methods are still employed although theoretically sound methods of analysis are now available. This paper illustrates a completely general method of analysis of the rheological properties of non-Newtonian fluids which id independent of the type of viscometer employed and also the rheological type of fluid encountered except those of a thixo-tropic or rheopectic nature. Values of the rheological properties of a fluid obtained using this method of analysis are based upon an absolute scale of measurement, eliminating the necessity of employ-jng arbitrary concepts of the flow properties of non-Newtonian fluids, The fundamental nature of this method of analysis is further demonstrated by evaluating the precision with which measurements have been made of these rheological properties by employing a dimensionless criteric: —the Fanning friction factor-Reynolds number correlation. INTRODUCTION Non-Newtonian fluids are not new in the petroleum industry. Particular properties of non-Newtonian materials have been recognized and used to advantage in drilling fluids and have found excellent application in the hydraulic fracturing process.&apos;.&apos; While usage of these fluids is not new in this industry, the art of measuring and describing their rheological properties is still in its infancy. This is reflected in the accepted method of evaluating drilling fluids, API RP 29, (1950), which recommends a one point "viscosity" measurement. These shortcomings are recognized by various groups in the industry and more adequate methods have been prposed.3,4 The efforts Of these people are to be commended; but, at the same time, the limitations of their method of analysis and viscometer must be appreciated. This method Of analysis is applicable Only to Bingham plastic fluids as is pointed out early in their articles. The consistency of a Bingham plastic fluid must also lie within the range of the viscometer which is also well defined by the author. It is believed that the apparent simplicity of this method of analysis and the utility Of Operation Of their viscometer will soon make the terms non-Newtonian fluid and Bingham plastic fluid synonymous as used in this industry; a tendency which has already been encountered. Application of this method of analysis to a non-Newtonian fluid other than a Bingham plastic can lead to erroneous results. The field of rheology has long been hampered by the adherence to arbitrary concepts of Of of the flow properties Of non-Newtonian fluids. This has resulted in the design and operation of viscometers so as to fulfill these concepts rather than following the more logical method of analyzing experimental data so as to yield a fundamental truth. completely Bingham proposed the concept Of plastic flow and yield asu rement. data presented by Buckingham could hardly be interpreted so as to substantiate this theory; but rather succeeded in establishing the fact that compliance with this flow concept could never be realized in the capillary tube vis-cometer.5 Bingham, in reply, sug-gested "Rather than complicate our formula, it would seem to me a much better plan to so change the conditions of flow that the formula will be linear... Reiner0 stated .This plan has been followed by me in collaboration with Miss Riwlin. . . . From the present paper (Reiner and Riwlin, 1929) it will be seen that it is possible to choose the conditions of flow so that the formula based on Bingham&apos;s law of plastic flow remains linear. This is realized in the rotation viscometer," Thus, a viscometer, with design and operation based on the Reiner and Riwlin equation for plastic flow, is restricted to a specific range of flow

Jan 1, 1958

By H. D. Outmans

Present first-order theory for frontal stability and viscous fingering of immiscible liquids is improved by including the nonlinear terms in the equations describing conditions at the interface of the liquids. This leads to the revision of several conclusions which were based on first-order theory. They concern the growth rate and changing shape of a sinusoidal disturbance, particularly in relation to the wave number of this disturbance. Nonlinear aspects of effective interfacial tension are discussed and it is shown that this tension is not simply a positive proportionality constant in a linear relation between pressure difference and curvature at the interface. Scaling requirements are determined from the dimensionless groups which govern fingering. Gravity and interfacial tension invalidate a previously formulated conclusion that the shape of a finger after a given displacement is independent of the displacement velocity. Also, similarity of fingering (and, hence, of sweep efficiency in case of an unstable front) requires geometrical similarity of the initial disturbance in the model and the reservoir with a scale factor which is the same as the one for scaling down the dimensions of the reservoir. INTRODUCTION The critical velocity at which the transition zone in the vertical displacement of immiscible liquids becomes unstable was calculated by ill. More recently, the same was done for nonvertical displacements and it was also shown that the stability is affected by an effective interfacial tension between the liquids. 2-4 Although in the oil reservoir a stable front is always desirable from a recovery point of view, the necessary velocity may have to be so low that the corresponding production rate is not economical. In that case, a knowledge of the fingering unstable front is necessary for predicting recovery at break- through. Studies in this direction have not gone beyond the very moment at which fingering first occurs.5,6 Conclusions about stability and fingering in these references are all based on linear theory. In this theory the nonlinear terms in the equations describing conditions at the front are neglected. The usual justification for this lack of mathematical rigor is that the nonlinear terms can be made small relative to the linear terms and, supposedly, small causes produce only small effects. It has been recognized, however, that this is not necessarily true in the nonlinear problems of hydrodynamics.7 The hydro-dynamic stability and fingering of the front between two liquids, accelerated in a direction normal to the front, for instance, is strongly influenced by the nonlinearity of the boundary conditions.8-10 Since these boundary conditions are not unlike those arising in the problem of stability and fingering of the front during a slow immiscible displacement in porous media, it was thought that there, too, the nonlinear terms should be taken into consideration. A summary of the results obtained in linear theory and a comparison with experimental data precede the nonlinear theory developed in this paper. This summary serves to introduce some quantities which will be used in the further develop ment. It should also be pointed out that the method of higher-order approximation by which the nonlinear stability problem is solved has application in other reservoir studies.11 SUMMARY OF LINEAR THEORY 1. A plane interface remains stable if its velocity is smaller than a critical velocity. vc = (p2-p1)gn/(l/A2- l/A1 .........(1.1) This equation is derived for uniform flow in a thin layer inclined to the horizontal plane. The displacement takes place in a direction normal to the intersection of this layer and the horizontal plane. The undisturbed (line) interface is normal to the velocity. Fluid 1, the upper fluid (gas), is displacing Fluid 2 (oil). (VcQS would be the injection rate in cc/