Search Documents

By J. W. Cunningham, W. A. Tiller, D. H. Lane

The effect of ultrasonic vibrations on ingot solidification has been considered both theoretically and experimentally. The theoretical section elucidates the mechanisms by which the ultrasonic vibrations might refine the ingot structure. The experimental section describes a very efficient means of introducing the ultrasonic vibrations to the freezing interface. The application of this technique to the consumable-electrode arc melting process shows that the columnar mode of freezing may be effectively suppressed yielding an ingot structure consisting predominantly of equiaxed grains. The results support the contention that the ultrasonic vibrations increase the nucleation frequency in the layer of liquid adjacent to the freezing interface. CERTAIN structural characteristics of as-cast ingots are undesirable since they create considerable difficulty for subsequent processing and often produce defects in the finished product. These structural defects are illustrated in Fig. 1(a), and are: 1) planes of weakness where columnar grains growing from the different mold faces meet, 2) strength anisotropy due to the preferred orientation of the columnar grains and 3) inhomogeneity due to segregation at dendrite boundaries within the columnar grains and at the columnar grain boundaries. The type of ingot structure that eliminates 1) and 2) and diminishes 3) is the fine-grained equiaxed structure illustrated in Fig. 1(c). Many attempts have been made to produce the desirable ingot structure illustrated in Fig. 1(c) by inoculating, casting with varying superheats, vibrating or stirring, or by irradiation with ultrasonic vibrations. These attempts have met with various degrees of success, the ingots generally having the structure illustrated in Fig. 1(b). All four methods have been the subject of numerous investigations,1"4 the latter producing the most erratic behavior and inconsistent success. The method of ultrasonic irradiation during ingot solidification has shown considerable promise on small-scale ingots (diam 2 in.) but has not yet been successfully adapted to large-scale operation. The present paper considers the application of ultrasonic radiation to ingot solidification both theoretically and experimentally. The purpose of the paper is to elucidate the mechanism by which ultrasonic energy refines the ingot structure and to describe a method for producing the desirable ingot structure in both laboratory scale experiments and in ingots of commercial size. THEORY Conventional Casting—Before discussing the possible effects of ultrasonic radiation upon ingot structure let us first briefly review the reasons for the breakdown of columnar crystallization to equiaxed crystallization during ingot solidification. The following is a slightly oversimplified description. In the conventional method of casting, the alloy having a certain degree of superheat, T, is poured to fill a cold mold which is at room temperature, Tc,. Due to conduction of heat through the mold wall the outer rim of liquid metal cools rapidly to the liquidus temperature and begins to supercool. At some degree of supercooling, nuclei of solid form at random in this supercooled layer at the mold surface. The grains in this chill layer begin to grow inwards producing a solute build-up, thus lowering the freezing temperature at the interface as illustrated in Fig. 2. As growth proceeds, the temperature gradient in the solid conducts away not only the latent heat of fusion evolved during growth but it also conducts away some of the superheat of the melt so that the temperature

Jan 1, 1961

By E. H. Edwards, J. Washburn, E. R. Parker

The mechanism of strain hardening was discussed in connection with some recent observations on the stress-induced motion of dislocation boundaries and on the simple shear deformation of zinc, cadmium, and copper crystals. THE nature of work hardening accompanying plastic deformation of metals has long been the subject of speculation. As early as 1925, Becker' suggested a thermodynamic process to explain slip based on the statistical probability for the occurrence of local glide steps. The theory that plastic glide was due to the presence of lattice defects known as dislocations was first introduced in 1934 by Taylor: and Orowan, and Polanyi.4 Taylor presented a rather complete analysis of the behavior and interaction of dislocations. He assumed a certain type of stress field to be associated with a single dislocation, and postulated that work hardening would result from plastic flow in two ways. First, dislocations following one another across a slip plane might encounter a barrier which would impede their movement, thereby gradually diminishing and eventually stopping plastic flow on that plane. The second mechanism involved the concept of dislocations traveling in opposite directions on nearby planes attracting one another to form a meta-stable lattice. The interaction of dislocations in such an array was shown to be large enough to require an increase in stress before plastic flow would continue. Kochendorfer5,6 extended the concept of work hardening involving a back stress due to an accumulation of like dislocations on a slip plane. He suggested that the formation of a new dislocation is hindered by those dislocations already present. The hardening resulting from the interaction of newly forming dislocations with bound dislocations was designated as "formative hardening," and indicated to be the only type of hardening necessarily connected with the slip process. In another detailed picture of work hardening, Mott and Nabarro7 assumed that sufficient dislocations are primarily present in the crystal as a consequence of growth irregularities (mosaic structure) to initiate slip. Work hardening resulted from the interaction of migrating dislocations with localized internal stresses produced during deformation. More recently, however, Mott8 as modified his treatment to include the generation of dislocations at Frank-Read sources and he accounts for the hardening by the formation of sessile dislocations. It has been widely accepted that metals of hexagonal structure, usually possessing only one set of glide planes, exhibit mechanical behavior markedly different from those having cubic structures. However, the differences in behavior of the two classes may be less fundamental than suggested. Rohm and Kochendorfer9 have shown that a single crystal of aluminum subjected to approximately simple shear gives the type of stress-strain curve associated with the hexagonal metals. Similar behavior has been demonstrated for gold and silver crystals by Andrade and Henderson.'" Thus when glide in cubic metals is restricted to one plane, the strain hardening characteristics closely resemble those of the hexagonal metals. It would appear, therefore, that an investigation of some of the strain hardening properties of single crystals sheared in simple glide might provide a more complete understanding of the phenomenon of work hardening. This report relates a number of recent observations made on single crystals of copper, zinc, and cadmium tested in simple shear. In addition, pertinent observations of effects accompanying stress-induced movement of dislocation boundaries are reported. Experimental Procedure and Results Details of the production and advantages of the type of single-crystal shear specimen employed for the following tests have been presented in a previous publication.n The method of testing makes

Jan 1, 1954

By J. F. Evers, S. Sadiq, F. W. Preston, G. W. Swift

Proper scaling of unsteady-state flow in laboratory gas reservoir models can be achieved if the permeability/fractional porosity ratio is less than 0.5 md. Porous media meeting this specification can be obtained from mixtures of sand and Portland cement provided (1) the water/cement ratio is in the range of 0.35 to 0.40, (2) the sand/cement ratio is in the range of 1/1 to 2/1, (3) the sand grain size is in the range of 48-100 to 18-28 mesh, (4) the mortar curing time is 7 days or more and (5) the cured mortar is dried under vacuum to reduce the temperature level employed during the drying process. The areal homogeneity of the mortars tested is reasonably good, although significant trends in vertical permeability were noted. The Klinkenberg and non-Darcy coefficients of mortar are such that proper scaling of these effects in the laboratory model could be expected. INTRODUCTION A long term deliverability prediction for a gas well is generally based on data determined from flow tests that are of short duration relative to the producing life of the well. These data, together with an estimate of gas reserves in place, form the basis for establishing the parameters of the mathematical model used to predict performance at future times. Results obtained from a mathematical model represent only an interpretation of the physical process being modeled. Experimental verification of the production process from discovery to abandonment is needed where production occurs under carefully controlled conditions. Ideally, one might consider such an experiment to be conducted in the field. Even if the cost of such an experiment could be justified, the time required would make it impractical to consider. However, if the time scaling problem could be handled properly, an experimental model of the gas reservoir production process could provide the necessary data to study the mechanisms of the production process in a reasonable time period. An experimental model has been designed and constructed using the data of this study as a basis. The model reservoir system consists of one producing well located at the center of a right-circular cylinder of porous material. Ten to 15 observation wells are located at strategic offset positions from the producing well. These wells, and the producing well, are fitted with pressure transducers so that the pressure-time history can be recorded at various radial and azimuthal positions in the drainage area during the entire production life of the system. Time is scaled in the model so that the time required to produce gas from discovery pressure to abandonment pressure is about 2 hours. This compares with a 20-year producing life in the model's field counterpart and gives a time scale factor of approximately 1:100,000. DESIGN CRITERIA To establish design criteria for the experimental model, the reservoir is considered to be a finite right-circular cylinder of radius rr, and height h with a producing well of radius rw located at the center and completed through the producing formation. This reservoir consists of a porous medium which has a uniform permeability k and gas porosity and the gas contained in the reservoir has density p, viscosity µ and molecular weight M. Initially, the pressure at all points within the reservoir is uniform at po. and the reservoir temperature is T. Since the reservoir production process has been demonstrated to occur in two distinct transient modes for either Darcy or non-Darcy flow,l it is

By C. H. Samans, G. F. Tisinai

HARDENING and embrittlement of the ferritic chromium stainless steels at temperatures near 885 °F have been known for a long time.' However, no satisfactory explanation has been given. Both ordering and precipitation hardening have been suggested. Fisher, Dulis, and Carroll2 discovered a chromium-rich body-centered-cubic precipitate in embrittled alloys. Subsequently, Williams' and Wright4 also have found evidence of this body-centered-cubic constituent. William's hardness, electrical resistance, and magnetic studies led to his proposal of a new constitutional diagram for the Fe-Cr alloys. The new feature of this diagram is the eutec-toidal decomposition, below about 970°F, of the s phase into two stable body-centered-cubic solid solutions, one rich in iron and the other rich in chromium. An observation of Heger5 that the amount of s phase has increased and the amount of chromium-rich body-centered-cubic phase decreased in the interval 5000 to 32,600 hr at 900°F, even though the alloy increased in hardness, casts some doubt on this diagram. If precipitation hardening causes the embrittlement, the precipitate probably is always coherent with the matrix, since no evidence of overaging has been reported. It cannot be stated positively, however, that Fisher's body-centered-cubic constituent is the embrittling phase. According to Preston's data,6 such a chromium-rich phase should have a parameter less than about 0.5 pct larger than that of the iron-rich matrix lattice. Moreover, the constitutional relationships whereby Fisher's constituent forms, and whether it is stable, metastable, or a transition phase, are not yet completely clear. The recent studies of Wright4 and of Lommel,7 as well as the earlier studies of Becket8 and of Hochmann,9 have shown that 885°F hardening occurs in relatively pure Fe-Cr alloys. No evidence of its occurrence or absence in chromium itself has been reported. The rate of hardening in high purity alloys seems appreciably slower than that in alloys containing larger amounts of impurities (chiefly carbon and nitrogen). This may be due merely to the absence of hardening and straining effects from precipitated carbides and/or nitrides. These observations on high purity alloys make it seem unlikely that 885 °F hardening is caused by a constituent containing other than iron and chromium atoms, e.g., carbides and/or nitrides. Wright,4 Lommel,7 and Lena and Hawkes10 emonstrated the accelerating effects of a finely dispersed precipitate of (probably) chromium nitride on the hardening reaction. This appears to be evidence that lattice straining increases the rate of nucleation of the hardening constituent. Some new observations on 885°F hardening are reported here. They show an effect of preheat treatment that should contribute to a better understanding of the 885°F hardening phenomenon. They also clarify a discrepancy which exists in the literature. Some investigators have reported most rapid hardening at 885", whereas others have reported it at 930°F. The present work also shows that 885°F hardening does not occur in electrolytic chromium, and that any long range ordering type of reaction affecting any sizeable part of the volume is unlikely. Materials Studied Table I shows the currlposition of the alloys used in this investigation. In these alloys carbon and nitrogen are present in roughly similar amounts. Presumably both elements precipitate, as carbides and/or nitrides, during the heat treatments given.

Jan 1, 1958

By C. R. Johnson, R. A. Greenkorn, L. K. Shallenbarger

This paper describes a study, based on core data, of the directional permeability of a sandstone reservoir. Directional air permeabilities are explained and correlated with lithology by the tensor theory of permeability, which is extended to the more general case of heterogeneous anisotzopic porous media. In this case, the permeability tensor is made up of two components: (I) anisotropy (variation around a point) which correlates with bedding, and (21 heterogeneity (variation from point-to-point) which correlates with grain size. INTRODUCTION Previous studies of directional permeabilityl.5-8 have concerned themselves only with anisotropy at a point, rather than with both point-to-point differences and local anisotropies. In this paper, we present the existing tensor theory of permeability for anisotropic porous media and then extend it to the more general case of heterogeneous anisotropic porous media. The laboratory measurements of directional permeability are explained in view of this extended theory and the data are discussed in terms of the lithologic factors that correlate with it. The results show what the heterogeneity and anisotropy of the reservoir element are, and that the directional permeability correlates with lithology. Data used in this paper are measurements of air permeability in eight directions, spaced at 45° intervals, on 142 2-in. vertical plugs from 30 cores. The core data are meaningful in terms of directional permeability because the cores were oriented to within 45O during drilling and coring.3 Although we initially thought that only 60 per cent of the core material was reliable, subsequent study showed that almost all of it was reliably oriented. After determining the air permeabilities, the data were reduced to three independent variables for each plug: the major and minor permeability axes, and the direction of the major axis. These were obtained by converting the permeability data to the reciprocal square root of permeability and fitting the transformed data with ellipses according to the tensor theory of permeability.8 The point-to-point areal variation of the minimum permeability axis is related to grain size. The direction of the permeability axes, where the major axis exceeds the minor axis by at least 5 per cent (measured variation is about 4 per cent), correlates with the bedding. The permeability tensor used in this study must be considered as the sum of a scalar and a tensor, with the scalar being the minor axis permeability as a function of position, and the tensor the directional effect, which is additive permeability over the minor axis. In this case, the point-to-point variation or heterogeneity (minor axis) is substantially larger than the variation at a point or anisotropy. It may be important that this separation of "directional permeabilities" be recognized when considering migration of fluids due to permeability variation. Local migration may be due to anisotropy and point-to-point migration may be due to heterogeneity, but the direction and magnitude of these may not be the same. Furthermore, in truly heterogeneous systems, one would expect that anisotropy would be the smaller of the two effects. THEORY TENSOR THEORY OF PERMEABILITY FOR ANISOTROPIC POROUS MEDIA8 Darcy's law for flow in porous media in its usual form is where qis the flow rate vector, k is the permeability, is the viscosity, is the vector differential operator, and P is the Pressure. Consider flow, q, through a unit area in the direction of the unit vector n from elementary flow tubes parallel

Jan 1, 1965

By J. R. Kyte, C. C. Mattax

Previous workers have developed differential equations to describe oil displacement by water imbibition, but have not explicitly defined the relationship between recovery behavior for a single reservoir matrix block and its size. In the present work, imbibition theory is extended to show that the time required to recover a given fraction of the oil from a matrix block is proportional to the square of the distance between fractures. Using this relationship, recovery behavior for a large reservoir matrix block is predicted from an imbibition test on a small reservoir core sample. The prediction is then extended to analyze recovery behavior for fractured- matrix, water- drive reservoirs in which imbibition is the dominant oil-producing mechanism. Experimental data are presented to support the basic imbibition theory relating matrix block size, fluid viscosity level and permeability to recovery behavior. INTRODUCTION Imbibition has long been recognized as an important factor in recovering oil from water-wet, fractured-matrix reservoirs subjected to water flood or water drive. l-3 Recently, two approaches have been published which might be used to predict imbibition oil-recovery behavior for reservoir-sized matrix blocks. Graham and Richardson4 used a synthetic model to scale a single element of a fractured-matrix reservoir. lair,' on the other hand, used numerical techniques to solve the differential equations describing imbibition in linear and radial systems. This latter method requires auxiliary experimental data in the form of capillary pressure and relative permeability functions. These two approaches, i.e., synthetic models and numerical techniques, have been used to study a variety of reservoir fluid-flow problems. One pur- pose of this work is to present, with experimental verification, a third method for predicting imbibition oil recovery for large reservoir matrix blocks. This method uses scaled imbibition tests on small reservoir core samples to predict field performance. The imbibition tests are easier to perform than the capillary pressure and relative permeability tests required to apply the numerical method. Furthermore, when suitably preserved reservoir-rock samples are used, the properties of the laboratory system are the same as those of the field. This offers an important advantage over the use of synthetic models because there is usually some question as to how accurately reservoir-rock properties can be duplicated in such models. Based on the recovery behavior for a unit matrix block, an analysis is presented to predict oil recovery for a fractured, water-drive reservoir made up of many such unit blocks. In the analysis, it is assumed that the flow resistance and the volume of the reservoir fracture system are negligible compared with that of the porous matrix. These assumptions are generally consistent with observed characteristics of many fractured-matrix reservoirs, and have been employed in previous studies.2,3,6-8 further assumed that the effect of gravity on flow in the matrix blocks is negligible. On first thought, the latter assumption might appear to seriously limit application of the method. However, in a fractured reservoir, the effect of gravity on flow in a matrix block will be restricted by the height of the block. Furthermore, matrix permeabilities are often very low (10 md or less) in such reservoirs. This means that capillary or imbibition forces will be large, thus tending to minimize the relative importance of gravity. For these reasons, imbibition should be the dominant oil-recovery mechanism in many fractured-matrix, water-drive reservoirs. The predictive method presented in this report is applicable to such reservoirs. SCALING IMBIBITION THEORY The "scaling laws" applicable to the operation.

By N. J. Grant, C. F. Flo, F. B. Cuff

An investigation of the Ti-Cr system has shown the presence of a complete series of solid solutions in the ß phase, with a minimum in the solid us near 50 pct Cr. An intermetallic compound, TiCr2, forms during cooling from the ß solid solution. There is a eutectoid reaction at the low chromium side of the system. IN view of the recognition of the potentialities of titanium and its alloys as important structural materials there has arisen a need for a systematic investigation of various titanium binary diagrams. Of these, the Ti-Cr system was of particular interest because of the improved properties imparted to titanium by small chromium additions. A literature survey disclosed a partial constitution diagram that had been suggested by Vogel and Wenderott,¹ This diagram, shown in Fig. 1, is the result of work done on the Fe-Cr-Ti ternary system. Alloys containing less than 43 pct Cr were not investigated. In addition to this, a study of alloys in the Ti-Cr system was reported by McPherson and Fontana.² Adenstedt3 furnished a number of dilato-meter curves obtained from specimens containing 5, 10, and 15 pct Cr. Very recently, McQuillan4 submitted a proposed diagram covering the complete range of compositions between titanium and chromium. This diagram was in general agreement with the one described here. Experimental Procedure Sponge titanium (99.7 pct) and electrolytic chromium (Bureau of Mines analysis: 99.03 pct Cr, 0.40 pct Fe, 0.53 pct O2) were used for the initial determination of the system. Subsequent refinement of points, particularly for the low chromium side of the diagram, was accomplished by using iodide titanium and high purity chromium (National Research Corp. analysis: 0.050 pct C, 0.045 pct O2). The Ti-Cr alloys were prepared by melting 25 g charges in an enclosed, water-cooled copper crucible using a movable tungsten electrode. The atmosphere was helium, prepurified by first passing it through a drying agent and then through sponge titanium heated to 850°C. Particular care was taken to avoid any large scale segregation resulting from insufficient mixing of the components. Each charge was melted and kept molten for approximately 1 min, allowed to cool, inverted, remelted, and kept molten for an additional minute. In order to ascertain the amount of impurity pickup (oxygen and nitrogen) during the melting procedure, three iodide titanium specimens were melted in sequence and then checked for increased hardness. The specimens in the as-cast condition showed a relatively large hardness increase over the hardness of the as-received iodide titanium. The scattering on these readings was large. However, after the specimens had been homogenized at 840 °C for 2 hr and slowly cooled, the measured increase in hardness over the as-received titanium was small with a corresponding decrease in scatter. The reason for this hardness change is that in heat treating these specimens below the transformation temperature (885°C) there was a conversion from the as-cast structure to the more nearly equilibrium structure of the annealed specimen. The alloys were analyzed for chromium content and although the actual percentage was invariably low the deviation from the nominal composition never exceeded 0.6 pct Cr.

Jan 1, 1953

By D. J. DeLazaro, W. Rostoker, R. E. Riley, M. Hansen

KNOWLEDGE of the isothermal transformation behavior and the TTT chart method of graphically summarizing such information has been of invaluable aid to the ferrous metallurgist in understanding and developing heat treatments intended to provide specific mechanical properties. It is to be expected that similar lines of study directed to titanium-base alloys will be as profitable. This paper summarizes, on conventional TTT charts, the isothermal transformation products and pertinent reaction rate data for binary alloys of titanium containing 1, 3, 5, 7, 9, and 11 pct Mo, respectively. The results so presented were derived from metal-lographic examination. Some X-ray diffraction work was done to clarify certain points of question. Ti-Mo System To the authors' k.nowledge, TTT charts have only been used to describe the reaction rates of eutectoidal decompositions. It would seem that the kinetics of decomposition of any unstable solid solution might be conveniently described in this manner. The equilibrium diagram of the Ti-Mo system' is shown in Fig. 1. The similarity with the Fe-Ni system will be recognized. Additions of molybdenum serve to stabilize the /3 phase to increasingly lower temperatures. The solubility of molybdenum in the a phase is restricted to less than 1 pct at 600°C. Under normal conditions. a homogeneous ,8 structure quenched to a, temperature in the a f & field will reject a of invariant composition. The rate at which this occurs depends on the degree of undercooling and the degree of supersaturation of /3. The structural character of the rejected phase may also be expected to be related to the undercooling. Materials and Experimental Methods A sponge titanium was used which has a nominal impurity content as reported in Table I. The nominal composition of the molybdenum used is also included. The alloys were prepared in 150 g melts by arc melting in an evacuated water-cooled copper crucible of the type reported in an earlier paper.' The pancake-type ingots were forged to 1/4x1/4 in. bars and homogenized for 24 hr at 1000°C in evacuated Vycor bulbs. Specimens Vi in. thick were cut from these bars for heat treatment. The heat treatment cycles followed an invariant procedure. The specimens were soaked at 1000°C for 20 min under a helium protective atmosphere prior to quenching into a lead bath held at a temperature corresponding to a chosen degree of undercooling. After a prescribed period of time at this temperature, the specimen was removed and water quenched. A series of specimens so treated for progressively longer periods served to establish the times for initiation and completion of the decomposition of the /3 phase at a particular temperature. Although heat treatment of unprotected specimens in a lead bath does permit surface contamination by lead diffusion, this method is necessary to insure sufficiently rapid quenching rates. As long as the bath treatments were of less than 1 hr duration. the contaminated surface was not unduly deep and could be removed without excessive grinding. Metal-lographically, the contaminated zone could be readily distinguished from the pure binary alloy. Transformation Characteristics of Ti-Mo Alloys Metallographic examination showed that a supersaturated p phase may decompose by one of two processes. Alloys containing 1 to 9 pct Mo, when quenched in water from the homogeneous /3 field,

Jan 1, 1953

By Noel Jarrett, Allen S. Russell, Bernard M. Starner, Stanley C. Jacobs

A process has been developed to refine high-grade commercial aluminum to 99.99 pct purity. This enzploys precipitating titanium, vanadium, and zirconiu~ as borides. The upgraded liquid is partially crystallized under conditions that prevent massive freezing. The resulting crystals are compacted and washed free of adhernig impure liquid aluminum by melting the top crystals while draining off the liquid. Finally, the bottorn portion of the crystal bed is melted and removed as the pure product. Leakage from the crystallizat~on vessel is prevented by a layer of alumina granules that are not wetted by the aluminum. ALUMINUM of 99.99 pct purity is desired in a variety of applications. These include electrolytic capacitors, catalyst carriers, decorative trim for automobiles, aluminum powder, and high-strength alloys. This market was for many years filled by refining commercial-purity aluminum in a three-liquid-layer electrolytic cell of the type originally proposed by Hoopes. To meet the expanding need for lower-cost high-purity aluminum, a refining process based on chemical precipitation and fractional crystallization was developed. The availability of this super-purity metal was announced in November 1961, and three patents1-3 relating to the process were issued in 1965. An alternate process, comprising continuous fractional crystallization of a portion of a liquid-aluminum feed stream under conditions of strong agitation at the surface of the growing crystals, was published by Dewey4 in 1965. Earlier proposals by Larsen 5 and by Regner6 also emphasized the advantage of rapid agitation of the melt in the region of the crystallization zone. johnson 7 more recently proposed extracting purified aluminum from an alloy by mixing crystals of pure aluminum into the molten alloy, allowing them to grow, and removing the pure aluminum as by casting. Aamot's 8 recent patent describes a freeze refining process that involves controlled partial freezing against cooled moving retort walls. The Alcoa process to produce high-purity aluminum is schematically described in Fig. 1. Relatively pure aluminum is treated with boron in a furnace (a) to precipitate high melting point impurities, the supernatant liquid is transferred to a crystallization vessel (b), the liquid is partially crystallized with mechanized "tamping" (c), the remaining impure liquid is tapped out and the crystal bed partially and directionally re-melted to wash the underlying crystals (d), and the balance of the crystals is melted and withdrawn as high-purity aluminum (e). PRINCIPLES OF FRACTIONAL CRYSTALLIZATION Commonly encountered impurity elements form two types of systems with aluminum. Titanium, vanadium, and some other transition elements in periodic groups IV, V, and VI are the high melting point impurities. Initial solidification of aluminum containing small percentages of these elements will produce solid markedly enriched in them at temperatures slightly above the freezing point of high-purity aluminum. Most other elements, such as silicon and iron, produce eutectic type reactions with aluminum. In this case the first solid formed on initial solidification at temperatures slightly below the freezing point for high-purity aluminum is markedly depleted while the remaining liquid is enriched in these elements. The equilibrium change in composition during solidification of aluminum contaminated with silicon, iron, or titanium is illustrated in Fig. 2. Here, for clarity, the temperature scale is expanded beyond the

Jan 1, 1968

By Tell Ertl

The term oil shale is defined. Foreign oil-shale developments are outlined briefly. The richest and most extensive oil-shale deposits in the United States lie in western Colorado. The Bureau of Mines is opening two adjacent areas for underground mining. in western Colorado. From the first area the oil shale is mined selectively as required by the retorting plant. The second area is being developed as a unit of a full scale oil-shale mine to ascertain the most practical methods of mining oil shale and to determine costs. In both of these mining areas, research will be carried on to develop new and improved practices for mining oil shale. Introduction Oil shale is a sedimentary rock containing a solid mineraloid of indefinite composition known as kerogen. On heating oil shale, the kerogen undergoes thermal decomposition and some of the vapor products of the decomposition can be condensed as shale oil. Oil from shale can be a major source of liquid fuels since extensive deposits of oil shale are known to exist throughout the world. The United States considers the establishment of an oil-shale industry only when doubt exists of the adequacy of the domestic supply of crude petroleum. On the other hand, oil shale is not strange to foreign oil men. The industry is an old one in other parts of the world. The French industry has been in existence since 1838, the Scottish since 1859, the Australian since 1860, the Estonian and Swedish The since World War 1, and the Manchurian industry since 1929. The most highly exploited oil-shale deposits are the Carboniferous shales of Scotland. About 20 seams ranging from 4 to I2 ft in thickness are worked. The yield during recent years in scotland has been 18 to 25 gal of shale oil per ton, from 3&"30'000 tons per year' The French and Swedish oil-shale deposits are of lower grade than the Scottish. The Estonian oil-shale deposits are 7 ft thick, and the oil yield averages 50 gal per ton. Over 1,000,-ooo tons are mined each year, half of which is used directly as fuel to substitute for coal and the other half is used for the production of oil. Immense deposits of oil shale that will yield 10 to 15 gal per ton exist in Germany. Because of the impending loss of the Rumanian and Polish oil fields and the bombing of synthetic liquid fuel plants, the Nazis in 1943 started a large-scale oil-shale program in Germany. By V-E Day, it is said, oil-shale production had reached the rate of 2500 tons per day and rapid expansion of the industry was imminent. Australian oil shale, known as torbanite, occurs in beds 2 to 6 ft thick. It is of a very high grade and on retorting yields 80 to 160 gal per ton. The oil-shale industry of Manchuria was the largest in the world during the past decade. The well-known Fushun coal deposit in Manchuria has an overburden of 500 million tons of oil shale about 450 ft thick that will yield 10 to 15 gal of shale oil per ton. The oil shale is

Jan 1, 1949

By P. Hazebroek, M. Prats, E. E. Allen

Well locations in water floods frequently do not coincide with the ideal well locations associated with the flood pattern which is used. A question of importance is. "How much oil is lost because of this?". This is the problem treated in this paper. It is assumed that each pattern in a five-spot flood has the proper number of production and injection wells but that the wells may be off-pattern. The distance by which a well is off-patern is given in terms of the pattern spacing, defined as the length of one side of a five-spot. The influence of off-pattern wells on the performance of a five-spot flood has been investigated mathematically for unit mobility ratio. All production and injection rates are assumed to remain equal. The loss in recovery resulting from pattern irregularitie: has been calculated at a cumulative water injection corresponding to a 98 per cent water-cut in a regular five-spot. Displacing a single production or injection well by about one-third the pattern distance reduces the total field recovery by 4.4/N per cent of the recovery obtainable with regular spacing, where N is the number of five-spots in the flood (assumed large). When off-pattern wells are separated by at least one regular five-spot and are displaced by less than one-third the pattern distance, the total loss in recovery in a large-scale flood is, at most, 1 per cent of the recovery obtainable from that field with a regular five-spot pattern. However, much higher losses may occur if the off-pattern wells are not separated by a regular five-spot, if the rates are imbalanced, or if the mobility of the injected water is higher than that of the crude. Although the loss in oil recovery due to pattern irregularity often may be low, the off-pattern case will always require injection of more water to reach the same final water-cut. The volume of water injected per unit volume of oil displaced has been calculated for each five-spot adjacent to a five-spot containing a single 08-pattern well. The cumulative water injection at the time of 98 per cent water-cut may be increased by as much as 37 per cent when a production well is displaced one-third the pattern distance toward an injection well. The existence of differences between actual subsurface locations and the idealized well patterns used in secondary recovery operations is well recognized. However, the methods commonly used for predicting waterflood behavior ignore such differences, even when subsurface locations are known. The reason for this is simple. The idealized well patterns are highly symmetrical and, as such, the production history of the pattern (and even of the field) can be determined from that of a few production wells—usually from only one. Consideration of the many irregular well spacings found in a specific field would entail a large amount of work, and then the results would be applicable only to that specific field. In an attempt to provide an answer to the more general problem, this paper considers a water flood having an idealized well pattern, except for one well which is not at its regular location. When several off-pattern wells are present, their contributions are assumed to be additive unless the off-pattern wells are close to each other. Only four positions of the off-center well (at x, y. of Fig. 1) are considered in this paper. These are shown

By R. E. Carter, C. Wagner, F. D. Richardson

By means of inert markers of radio-platinum, it has been shown that cobalt metal oxidizes by out-ward diffusion of cobalt atoms through the oxide. Oxidation rates have been measured at various temperatures and oxygen pressures and have been found to agree with the rates calculated from the Wagner equation and the authors' values for the diffusion coefficient of cobalt in the oxide. The distribution of radio-cobalt in growing oxide loyers has been accurately measured and found to be different from that predicted from the Wagner oxidation theory. Attempts have been made to measure the change of lattice parameter of the oxide with composition. SEVERAL investigators have shown that cobalt metal oxidizes according to the parabolic law,1-3 but no detailed study has been made of the oxidation mechanism. Cobaltous oxide has a metal-deficient lattice, but the defect concentration which is dependent on the oxygen pressure is only approximately known.' From the similarities in the structures of "CoO" and "FeO" and the rate laws governing the oxidation of cobalt and iron, Gulbransen- as suggested that cobalt, like iron," oxidizes by diffusion of cobaltous ions through the oxide layer. To substantiate this hypothesis, he has attempted to measure the oxidation rate of cobalt as a function of the oxygen pressure but has been able to report only that the rate depends on a "small" power of the pressure. On the other hand, Preece," Valensi,' and Arkharov' all have suggested from indirect evidence that oxygen diffuses through the COO to the metal-oxide interface. The present authors, in a previous paper,' have shown that the rate of diffusion of cobalt in cobaltous oxide is dependent on the oxygen pressure. In that investigation, they also found that cobaltous oxide disks made by complete oxidation of cobalt metal were porous at their centers (Fig. 5L, ref. 4). Both these facts strongly suggest that oxidation of the metal proceeds via outward diffusion of cations. In order to establish the mechanism unequivocally, oxidation experiments were made with inert radioactive platinum markers, and the dependence of oxidation rate on oxygen pressure was measured. In addition, an attempt was made to determine the form of the cation vacancy gradient, which the marker experiment showed must exist across a growing oxide layer. This was done by depositing radioactive Co on the surface of a cobalt disk and then measuring the tracer distribution across the oxide layer produced by partial oxidation of the metal. Finally, an X-ray examination was made in an attempt to find the expected change of lattice parameter of COO with oxygen pressure. Experimental The apparatus and materials used were substantially the same as those previously described by the authors.4 Marker Studies: The use of inert markers to follow oxidation reactions was first attempted by Pfeil" in 1929, and radioactive markers have recently been used by Davies et al.' In the investigation reported here, the markers were applied as an (NH4)2PtCI6 slurry in several thin streaks across both faces of four disks of cobalt metal 1 mm thick and 14 mm in diameter. The disks were then oxidized almost to completion at 1200°C in pure oxygen. After oxidation, they were mounted and sectioned and auto-radiographs were taken to locate the markers. Oxidation Rate: The oxidation rate of cobalt was measured at 1148°C at oxygen pressures ranging from 0.0055 to 1.0 atm and at 1000°, 1148", and 1350°C under 1 atm of oxygen. The rate was deter-

Jan 1, 1956

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

By T. M. Geffen, H. N. Hall

One-third to one-half of the original oil is left in most reservoirs even after thorough gas driving or waterflooding. In contrast, essentially all of the oil contacted con be recovered by flooding with a solvent (a fluid miscible with the reservoir oil). It is not economical to fill a reservoir with solvent since the value of the solvent would very likely exceed that of the oil. A practical method would be to inject a small bank of solvent to displace the oil. This solvent bank in turn could be driven through the rock by a less valuable scavenging fluid, such as natural gas which is also miscible with the solvent. Such a method of oil recovery was studied in the laboratory using long (up to 95 ft) core systems. It was found that, as the solvent moves through the reservoir, the front end of the solvent bank becomes mixed with the oil and forms a relatively short mixing zone between undiluted crude oil and pure solvent. An engineering correlation was obtained which can be used to predict the length of this zone. Displacement of the solvent by natural gas can be achieved with complete miscibility (above the critical pressure of the gas-solvent system at reservoir temperature) or with partial miscibility at lower pressures. Under conditions of complete miscibility, it was found that the length of the gas-solvent mixing zone could be predicted from the correlation developed for the solvent-oil mixing zone. With partial miscibility, two phases exist and the solvent is displaced by a modified gas drive in which vaporization and condensation occur. Equations are developed to predict the performance of this type of displacement. Two plans are suggested for application of miscible fluid displacement. one for miscible displacement of the solvent and the other applying to lower pressures where partial miscibility occurs during solvent displacement. Calculations based on the laboratory tests indicated that a bank of solvent with a volume of less than 10 per cent of the hydrocarbon volume under consideration will generally be adequate to obtain high oil recovery from the portion of the reservoir contacted in miscible fluid displacement. INTRODUCTION The conventional methods by which oil is displaced from reservoir rock—solution gas drive, conventional gas drive, and waterflooding—are all processes in which the oil is expelled from the rock by direct displacement with a fluid immiscible with the oil. Because of the immiscible type displacement, complete recovery of the reservoir oil is not achieved. In economically successful flooding operations, as much as 50 per cent of the original oil may be abandoned in the reservoir as unrecoverable. The use of a solvent (miscible with oil, such as LPG) as an injection fluid does not have the undesirable displacement features common to the use of gas or water (immiscible with oil). It can be readily appreciated, however, that it would not be economically feasible to displace the oil and entirely replace it with a solvent. A process was therefore devised for use in field operations in which only a small volume of solvent would be required to obtain essentially complete oil recovery in the portion of the reservoir contacted. In this process, the solvent acts as a buffer between the reservoir oil and injected gas, allowing the natural gas to displace indirectly the oil with much greater efficiency than if the gas displaced the oil directly. DESCRIPTION OF OIL RECOVERY PROCESS The solvent oil recovery process will be described in

Jan 1, 1958

By Thomas L. Kesler

Present Raw Materials The lithium industry has had exceptional growth since publication of the last edition of this volume, and the present scale of mining and consumption of raw materials is a great many times that of 1949. The materials of present economic importance consist of four pegmatite minerals and lithium-bearing brine. Spodumene, a lithium pyroxene, occurs in greatest abundance, and is the chief source. It has the formula LiAlSi2O6 and theoretically may contain 8.03 pct Li2O, but in the more important deposits the content is 6.5 to 7.5 pct. The hardness is 6.5, the specific gravity 3.1 to 3.2, and the color mostly pale green to white although deep green, lilac, and yellow varieties exist. The crystals are prisms with strongly developed longitudinal cleavages, and crushing characteristically produces lath-shaped particles. Giant crystals up to 47 ft long are known, but in most occurrences the predominant size range is 0.5 to 12 in. long and 0.1 to 2 in. wide. Hypogene* alteration has formed lithia-bearing sericite and coarser muscovite and rarely eucryptite, while supergene* alteration yields kaolinite and hydrous mica. Lepidolite, of second importance, is a mica of varying composition approximating K2Li3-A14Si7O21(OH,F)3. It has lamellate cleavage, book structure, and a specific gravity of 2.8 to 3. The books range in size from microscopic to a diameter of several inches. The color is mostly violet, but other colors include pink, gray, yellow, and white. Most commercial lepidolite contains 3 to 4 pct Li20 and chemically is in the higher grade range of a series of lithium micas that includes polylithionite, zinnwaldite, and protolithionite. In this series, the content of Li2O is exceptionally as high as 7 pct. Petalite and amblygonite are also of economic importance, but are used in smaller amounts. Petalite has the formula LiA1Si4O10 and is white, gray, or pink. The hardness is 6, the specific gravity 2.4, the cleavage in two planes 114° apart, and the Li2O content commonly 3 to 4.7 pct with a maximum of 4.88. Amblygonite has the formula LiA1FPO4 in which OH may substitute for F to form montebrasite. It is white and has polysynthetic twinning in two directions at 90°. Amblygonite superficially resembles potash feldspar, but has a bluish or grayish tint unlike the cream or salmon tint of the feldspar.[t] The hardness is 6 and the specific gravity 3 to 3.1, and the content of Li2O is the highest among lithium minerals that occur in significant quantity, ranging mostly from 7.5 to 9 pct with a maximum of 10.2. The brine of Searles Lake in California30 has been a commercial source of lithium since 1938. See Chap. 7, p. 108. The "lake" is a saline evaporite deposit, and the exposed part has an area of 12 sq mi. To an average depth of 71 ft, this part of the deposit consists of crystalline salts with 50 pct interstitial voids. The voids are filled with brine that contains mainly sodium and potassium compounds, but the brine also contains 0.015 pct Li2O, which is recovered as dilithium sodium phosphate containing 21 pct Li20. This deposit overlies a

Jan 1, 1960

By B. Wilhelmi, W. H. Somerton

Simultaneous measurement of pore and elastic properties of socks under a wide range of triaxial stress conditions may be made with equipment and methods developed in this paper. Tests were run on three outcrop sandstones at several confining pressures and axial stresses up to approximately 80 percent of failure stress. The greater reduction in porosities and permeabilities occurred in the hydrostatic loading portion of the stress cycle, but substantial reductions in both properties did occur upon application of deviator (triaxial) stress. Triaxial elastic moduli were calculated from strain gauge measurements. Both Young's moduli and shear moduli were found to increase with increasing confining pressure. Poisson's ratios increased with axial stress but showed no consistent change with confining pressure. Only very general relations were observed between pore and elastic properties. The lesser decrease in porosity and permeability at higher stress values is related to increase in rigidity of rocks at higher stresses. INTRODUCTION Many investigators have studied the effects of hydrostatic stress on the properties and behavior of rocks. A few have studied the effects of triaxial stress fields on individual physical properties of rocks. To our knowledge no one has measured simultaneously pore and elastic properties of rocks under a wide range of triaxial stress conditions. This is the purpose of this work. Rocks exist in the subsurface under what Jaegerl refers to as "polyaxial" stress conditions (all principal stresses unequal). When two of the three principal stresses are equal, this is referred to as triaxial loading. The triaxial method of loading is most commonly used in modern testing work because experimental difficulties in measuring physical properties of rocks are greatly reduced over polyaxial testing. Although triaxial loading may not reproduce subsurface conditions precisely, subsurface stress conditions are never known precisely anyway. Triaxial testing would seem to be more realistic than testing under hydrostatic loading (all principal stresses equal) for application of results to subsurface problems. The objectives of this work were threefold: (1) to develop methods of measuring simultaneously pore and elastic properties of rocks under triaxial stress conditions, (2) to determine whether changes in these properties are of sufficient magnitude to be considered in subsurface calculations and (3) to attempt to relate changes in pore properties and elastic properties since they will have been measured at the same time under identical loading conditions. REVIEW OF EARLIER WORK PORE PROPERTIES Several investigators have studied the effects of stress on the physical properties of rocks. Fatt and Davis2 showed the reduction of permeability of rock with hydrostatic loading. McLatchie et al.3 and Knutson and Bohor4 related permeability reduction to compessibility of reservoir rock, again under hydrostatic loading. Wble,5 and later Redmond,6 evaluated the effects of applied radial pressure on conductivity, porosity and permeability of sandstones. In this latter work, no axial stress was applied. Dobrynin7 reviewed much of the previous work done . on the effects of overburden pressure on physical properties of rocks. He concluded that changes in physical properties are controlled to a large extent by pore compressibility of the rock. Several sets of curves developed from empirical relations are presented. These curves permit estimation of change of properties with net overburden pressure, knowing the pore compressibility and porosity of the rock. Net overburden pressure is defined as (Pe - 0.85 Pi ) where Pe is the

By M. H. Gaskell, D. M. Kehn, G. T. Pyndus

A correlation of the bubble point pressure for black oil systems is developed using the standard physical-chemical equations of soiutions. The correlation is based on 158 experimentally measured bubble point pressures of 137 independent systems and is expressed in terins of the usually measured field pnrameters— flash separation gm-oil ratio, rank oil gravity, total gas gravity, and reservoir temperature. The data were obtained on systems produced in Canada, Western and Mid-Continent01 United States, and South America. The average error (algebraic) in the representation is 3.8 per cent, and the maximum error encountered is 14.7 per cent. INTRODUCTION In the absence of experimentally measured properties of reservoir fluids, it is often necessary for the field engineer to make estimates regarding the fluid properties based on the usually measured producing parameters. To aid in these estimations, various correlations have appeared in the literature in recent years. Among the pertinent properties of interest is the bubble point pressure. A correlation for this parameter has been reported by Standing'. However, this correlation was based essentially on California produced crudes and this limitation was pointed out with its presentation. The correlation presented in this paper utilized data on crude oil systems from Canada, Wes- tern and Mid-Continental United States, and South America. CORRELATION DEVELOPMENT The basic assumption used in this development is the same as employed by Standing', There is a wide variety of ways to combine these parameters; however, in this instance the combination was made on the basis of Henry's law2. Accordingly, the relationship proposed is Although Eq. 2 defines an individual system, it is of limited value since H' is a function of gas-phase composition and the system temperature. It was observed that for the systems where the bubble point was measured at several temperatures that the ratio of the bubble point pressures and the ratio of the corresponding absolute temperatures (R) were practically identical. Thus, for correlation purposes the bubble point pressure may be taken as a direct function of the absolute temperature. This relation is valid only for those systems that are not near the critical point. Accordingly, this correlation will be inadequate for systems in the region of the critical point. The solubility of the various hydrocarbons found in the gas phase increases with the molecular weight. Thus, the saturation pressure should be inversely related to the gas gravity. Applying these principles to Eq. 2 and rearranging terms gives: The variables 02 the left side of Eq. 3 were designated as the "hubble point pressure factor". The number of mols of tank oil per barrel is a function of the "mclecular weight" of the tank oil. Although the tank oil is a complex mixture, it was assumed for the purposes of this correlation that a unique molecular weight could be assigned to a given oil. This was designated as the "effective molecular weight", and was related to the oil gravity, This empirical relationship was developed simultaneously with the correlation by assuming values of M,. and working to obtain a smooth curve for both the correlation and the effective molecular weight. The relationship between the oil gravity and the effective molecular weight used in this correlation is shown in Fig. 1. The effective molecular weight is somewhat higher than the molecular weight of the C fraction. The difference between these values is largest for the low-gravity systems. It is noted that this effective molecular

By Nicholas J. Grant, Hanna Ibach, Bill C. Giessen

The system Ta-Rh was investigated over the entire comnposition range using metallogvaphic and X-ray techniques as well as thermal analysis. Terminal solubility limits, solidus temperatures, and the crystal structures of two new intermediate phases were determined. There are five intermediate phases: o, tetragonal, isostructural with o(Fe-Cr); a1, orthorhomnbic : a2, orthorhombic, isostructural with CozSi; 0 3, a high -temperature phase of unknozun structure; and a TaRh3, cubic, AuCu3 structure. a, aL, a,, and a melt peritectically; and a TaRh, has a maximum melting point. Five peritectic and one eutectic reactions occur. THE Ta-Rh phase diagram has not been fully treated in the literature; however, two inter metallic phases have been described.&apos;-&apos; This study attempts to work out the diagram covering the details of phase boundaries as well as the crystallography of the intermediate phases. EXPERIMENTAL METHODS Starting materials were tantalum powder of 99.8 pct purity, supplied by National Research Corp., Cambridge, Mass., and rhodium powder of 99.9 pct purity, supplied by Bishop Co., Malvern, Pa. Analyses as given by the manufacturers are presented in Table I. In addition to the spectroscopic analysis, an oxygen analysis has been made of the tantalum powder. After preliminary trials, the standard procedure employed in making the alloys was as follows (see Ref. 4 for further details). The powders were weighed, mixed manually, compacted under 16,000 psi, and arc-melted at not more than 500 amp under thoroughly titanium gettered welding-grade argon of 400 mm Hg pressure. The gettering process was repeated after each melt. It had been shown previously5 that under these conditions no oxygen pick-up exceeding 0.001 pct in a hafnium button could be detected. The buttons, ranging from 5 to 10 g in weight, were inverted or, if possible, crushed and remelted three times. Melting presented no difficulty, except for tantalum-rich al- loys, which showed considerable spashing during melting. The average weight loss of 1 wt pct rose to 2.5 wt pct in these cases. Subsequent melting, however, resulted in losses of less than 0.1 wt pct. The melts obtained frequently showed small in-homogeneities, mostly in the tantalum-rich part of the system below 25 at. pct Rh. This was not due to insufficient blending of the elements, but to crystallization taking place from the cooled bottom. This condition was not improved by further re-melting; therefore the following precautions were taken. Cross sections of as-cast buttons were made. For the investigation only buttons were used where these inhomogeneities did not exceed 1 at. pct (estimated from the amounts of haSeS resent). In addition, only the center portions of such buttdns were utilized to minimize errors. Smaller buttons showed smaller concentration gradients and were used as a check. Finally, all alloys were homogenized to eliminate coring (see below). After this, the error in alloy composition due to inhomogeneity was cl at. pct. In order to confirm the concentrations of the alloys and to provide an impurity check, wet-chemical analyses and oxygen-fusion analyses of four key alloys were carried out. The results, together with the concentrations as calculated with the weight balance method, are presented in Table 11. It was found that for three alloys the agreement is very good (< 0.2

Jan 1, 1964

By D. A. Flanagan, B. L. Landrum, P. B. Crawford, B. C. Norwood

Herern is described a heat conduction model which has been developed to predict the pressure distribution within a reservoir being subjected to rrnsteady-state depletion. Although potentiometric models have been widely used to study reservoir pressure distribution, their application is limited to the study of steady-state cases. Electrical analyzers may be used to study unsteady-state phenomena, but their application and simulation have certain limitations. Much of the now accepted theory of unsteady-state fluid flow in porous media has evolved from previous heat conduction knowledge. An analogy between the flow of fluid in an unsymmetrical reservoir and con-duct ion of heat in a similarly shaped unsymmetrical metal plate is made to show adaptability of the hear conduction principle to transient reservoir cases. A description of the model and necessary instrumenrarion is included. Data show influence of well location on pressure distribution during unsteady-state deple[ion. Application of results is made to a hypothetical reservoir and to a reservoir located in an irregular aquifer. INTRODUCTION Development of methods for predicting behavior of petroleum reservoirs has resulted partially from adaptation of electrical and thermal flow concepts to the flow of fluids through porous media. Extensive investigation of electrical and thermal flow problems has resulted in mathematical solutions applicable to many reservoir problems. The relative ease of measuring electrical flow has brought wide acceptance of the use of the potentiometric model to study certain field problems. Potentiometric models are restricted to the study of steady-state conditions because of the instantaneous response of electrical flow to changes in the applied potential. The electrical reservoir analyzer has been one of our most useful instruments for studying certain phases of unsteady-state behavior1. Other than approximation techniques, no method exists for determining pressure distribution in an irregular reservoir during primary depletion. Heat flow, like fluid flow, is inherently transient. The model described here was developed to use this transient nature of heat flow, measurable in the laboratory, to simulate transient flow of fluids in the porous medium of the reservoir above the bubble point. A metal plate shaped like the reservoir to be studied was heated by a "well" source. Temperature increase in the plate is analogous to pressure decline in the field if the field is being depleted, or pressure increase if fluid is being injected. This study is intended to determine the feasibility of using the transient nature of the heat flow to study the transient behavior of fluid flow. It is hoped that this new technique will aid in the understanding of one of the more complex phases of reservoir engineering. THE FLUID FLOW-HEAT CONDUCTION ANALOGY The similarity of the equations which describe the flow of fluids through a porous medium and the conduction of heat through a solid is well known2. Rate of flow of fluids is given by Darcy&apos;s law3. while the rate of flow of heat is given by Fourier&apos;, (I,, - KA 4T........(2) Combining Eq. 1 with the equation of continuity for a compressible liquid", i.e.,results in Expressing y in terms of the compressibility, i.e., pr = peey ,..........(5) and by expanding the exponential and neglecting higher orders of cp one can obtain

By George D. Graffin.

The first use of graphite is lost in the mists of time. It was used by primitive man to make drawings on the walls of caves and by the Egyptians to decorate pottery. As early as 1400 A.D. graphite crucibles were being made in the Haffnerzell district of Bavaria. Through the Middle Ages graphite was confused with other minerals, especially galena and molybdenite. Two common names that are still used for the mineral are "Plumbago," meaning leadlike, and "Black Lead." The latter name implied that graphite is either composed of lead or at least contains a large percentage of it. Gessner is credited with having recognized it as a separate mineral in 1565, but its composition was not determined until 1779, when Scheele demonstrated that graphite oxidized to carbon dioxide, thus proving its carbon constitution. In 1789, Werner named it "graphite" from the Greek word "graphein," meaning to write. The United States has long relied on foreign countries for its graphite supplies. This reliance has been due more to the lack of suitable domestic graphites than to any particular preference for buying minerals overseas. Ceylon (now Sri Lanka) graphites were imported as early as 1820, and Madagascar graphites have been imported since the properties on that island were opened up about a hundred years later. Mexico has long supplied the United States with most of its requirements for amorphous graphite. Physical Properties and Classification Graphite is found in laminated, flaky aggregates disseminated in schistose rocks. It also occurs in veins and exhibits a foliated or fibrous structure. Graphite is a black lustrous mineral that crystallizes in the hexagonal system, with rhombohedral symmetry. The crystals have tabular form and are six-sided; the faces are commonly striated. The flakes have perfect basal cleavage and are opaque. When well crystallized, the flakes have a black metallic luster while the amorphous material is black and earthy, with a microcrystalline compactness. The flakes feel greasy. Graphite is an excellent conductor of heat and electricity. It melts at approximately 3550°C at a triple point under 1250 psi. It sublimes between 3300-3500°C at 14.7 psi. At 3726°C a second triple point occurs at approximately 100,000 atmospheres. Thermal oxidation in the presence of oxygen begins at 300°C and the rate increases with temperature. The three forms of carbon (charcoal, graphite, and diamond) are distinguished by chemical and physical tests. The specific gravity of charcoal is 1.3 to 1.9, of graphite 2.266 g/ cm3 (crystal density), and of diamond 3.5. Graphite has a hardness of 1 to 2 (Mohs&apos; scale). Natural graphite can be divided into three classes: disseminated flake, crystalline vein (fibrous or columnar), and amorphous. Flake graphite is a lamellar form found in metamorphic rocks, such as marble, gneiss, and schist. Each flake is separate, having crystallized as such in the rock. Crystalline vein graphite is found in the form of well defined veins or as pocket accumulations along the intrusive contacts of pegmatites with limestones. The graphite of such deposits is of two types, foliated and columnar. The important Sri Lanka (Ceylon) deposits are vein type. Amorphous graphite is commonly found as minute microcrystalline particles more or less uniformly distributed in feebly metamorphic rocks, such as slates and shales, or in beds consisting practically entirely of graphite. The latter usually are metamorphosed coal seams and carry as much as 80-85% graphitic carbon, while the former, being altered carbonaceous sediments, commonly range from 25 to 60%

Jan 1, 1983