Search Documents

By L. J. Ingvalson, Roy H. Miller

IN 1948, as part of a program to expand the copper tube mill facilities of the American Brass Co. plant at Kenosha, Wisconsin, plans were formulated to convert the 100 ton capacity anode casting furnace at the Great Falls, Montana plant of the Anaconda Co. to the production of 3 in. phosphor -ized billets. The furnace was rebuilt completely and construction was started on an addition to the building. Members of the Great Falls staff visited various billet casting plants to draw upon the experience of other producers of this type of casting. Many features of other plants were incorporated into the Great Falls plant but, in addition, many new designs were made and perfected to make the plant the most advanced of its type. Construction was begun in 1949, with the intention of having the plant in production in the early part of 1951. In January 1951, as the plant was being rushed to completion, the government placed restrictions on the production of copper tubing. As a result, work was halted and the plant lay idle for two years. Interest was reactivated in May 1953 with the lifting of restrictions on copper. Plans were laid to go into production in July and the first charge was taken out July 1, 1953. Basically the plant is designed to produce from cathode copper a 3 in. diam billet, 50 in. long. weighing 110 lb, and containing 0.013 to 0.036 pct P, at a rate of 150,000 Ib per operating day. These billets are to be pierced and drawn into tubing. Furnace The furnace is a 100 ton capacity reverberatory type, gas fired furnace; 21 ft long and 9 ft wide. Gas is introduced to the furnace through two 38 hole Mettler gas burners. These burners are of chrome nickel alloy with four jets per hole. The furnace is charged with cathodes from the Electrolytic Copper Refinery plus return scrap from the casting operation by means of a Wellman Seaver Morgan Co. 7,000 lb capacity charging machine. The charge crane has a 20 ft boom and fork peel provided with an air ram to push material off the peel. The machine travels on overhead tracks and can be moved in any direction. The boom can be raised and lowered and swung through an arc of 180". With well stacked units of 6,000 lb, an experienced operator can put the first charge into the furnace in 40 min. After charging, the steel framed silica brick charge door is put in place and sealed with clay. After melting and skimming off the slag, the copper is blown until an oxygen content of 0.40 pct is reached. The metal is then heated to 2180°F, covered with coke, and poled. Poling is continued until a 3 in. diam test billet, 8 in. long, shows a flat set. The oxygen content is usually 0.030 pct at this time. At the end of poling, a temperature of 2160" to 2180°F is desired. In the side of the furnace opposite the charging door is the taphole slot, extending 36 in. upward from the bottom of the furnace. It is 4½ in. wide at the outside face and widens to 9 in. at the inside face of the furnace. It is filled with a wet mixture of sand, fireclay, and crushed coal. Flow of metal

Jan 1, 1957

By J. L. Thompson

There are two basic forms of the LPG-slug process: the gas driven and the water driven. The pressure required for miscibility between gas and LPG prohibits the use of the gas-driven LPG process in shallow reservoirs. The water-driven LPG slug process normally exhibits good sweep efficiency. However, displacement of the LPG by water is poor. An improvement in this process appears possible by injecting a slug of carbon dioxide between the LPG slug and the water drive. Laboratory experiments were conducted in homogeneous linear systems to determine the effect of pressure on the various displacement zones. In the carbon dioxide-LPG displacement zone, the effect of pressure on the displacement efficiency is negligible above the miscibility pressure. However, the displacement efficiency decreases linearly with decreasing pressure below the miscibility pressure. With water displacing carbon dioxide, it was found that the recovery of carbon dioxide after water breakthrough can be predicted from its solubility in water. A displacement test was conducted with LPG and carbon dioxide slugs large enough to avoid interference between the oil-LPG, LPG - carbon dioxide and carbon dioxide - water displacement zones. Under these conditions, essentially complete oil and LPG recovery was obtained. However, a substantial amount of carbon dioxide was left in the core at water breakthrough. INTRODUCTION It has been shown in homogeneous linear systems that essentially complete oil recovery can be obtained by displacement with a miscible fluid such as LPG. Economic considerations dictate that the valuable LPG also be recovered. Methods proposed for recovering the LPG include (1) displacement with natural gas at a pressure high enough to make the two fluids miscible1 (the disadvantage of this method is that the required pressure of about 1,500 psi is too high for shallow reservoirs) and (2) displacement of the LPG with water2 (principal disadvantage is that a large amount of LPG remains in the formation because of poor displacement efficiency). This paper treats an improvement in the water-driven LPG slug process3,4 where carbon dioxide is injected between the LPG and water. Carbon dioxide is completely or at least partially miscible with propane and soluble in water; thus, displacement efficiency is improved at both ends of the carbon dioxide slug. To minimize effects of other variables, all experiments were conducted in linear core systems. GENERAL PROCESS DESCRIPTION The displacements of propane by carbon dioxide and carbon dioxide by water are considered separately, assuming enough of each is injected to prevent mutual interference. Displacement descriptions are confined only to the part of the pore volume containing mobile fluids unless otherwise specified. The schematic representation is given in Fig. 1. Starting at the producing well, a zone containing only oil is followed by a zone containing a single-phase mixture of propane and oil. The concentration of propane increases in the direction of the injection well and ultimately grades to pure propane. Then follows a zone containing a mixture of propane and carbon dioxide in continuously varying proportions until 100 percent carbon dioxide is reached. In this zone, the concentration of carbon dioxide in the immobile water phase also varies from zero to full saturation. Whether the displacement of propane by carbon dioxide is of complete or partial miscibility depends on pressure and temperature (Fig. 2). If the horizontal line representing the prevailing pressure intercepts the boundary curves of the phase envelopes for the prevailing temperature (i.e., if the pressure is lower than the critical

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 Philip R. Strand, O. F. Stewart

Vermiculite is the name generally applied to the group of hydrated ferromagnesian aluminum silicates that are characterized by the ability to expand when heated. This process, called exfoliation, results in a lightweight product of commercial value. Most uses of vermiculite are for the expanded material. The chief markets are in construction, agriculture, and horticulture, with lesser uses in general industry. Composition and Properties Vermiculite, in its natural state, has the characteristic micaceous habit, a perfect basal cleavage which causes splitting into thin laminae that are soft, pliable, and inelastic. The structure of vermiculite is basically that of a talc. The prominent monoclinic crystal faces are often marked by lines at 60° and 120°. Hardness varies from 1.5 to 2 or more; specific gravities are between 2.1 and 2.8; color varies from almost clear to amber, bronze, brown, green, or black. Vermiculite feels like talc, especially when wet. Although much research has been performed on the chemical and structural composition, there is not yet complete agreement on the exact formula. This is to be expected when different workers have examined the many different varieties. Vermiculite is not considered to be a single mineral species but a family of related minerals. The structural formula for a trioctahedral vermiculite may be written: (H20)-(Mg,Ca,K)-(A12,Fe,Mg)-(Si,A1,Fe)4O10 (OH) 2 Hydrobiotite also occurs with vermiculite and is usually considered a vermiculite for commercial uses. When heated quickly to elevated temperatures, vermiculite expands by exfoliating at right angles to the cleavage into wormlike particles. The name vermiculite is derived from the Latin vermiculare, to breed worms. This characteristic of expansion is the result of the mechanical separation of the layers by the rapid conversion of contained water to steam. The decrease in bulk density of commercial grades is usually approximately 10 times, from 50 to 5 lb/ft3, but varies depending on the quality, size, and furnace efficiency. Individual flakes may expand up to 30 times. Vermiculite may also be expanded by soaking in chemicals such as hydrogen peroxide, weak acids, and other electrolytes. Color change during expansion is dependent upon the type of vermiculite and furnace conditions. Heating in an oxidizing atmosphere produces a dull gray or tan color, whereas a reducing atmosphere can produce a bronze or gold color. The expansion of the vermiculite crystal results in large pores being formed between the platelets. Thus, exfoliation makes available a large increase in void volume which is important in the application of vermiculite as a chemical carrier. Mathieson and Walker (1954) state vermiculite must be regarded as a true clay mineral. The characteristic properties of the mineral, such as high cation exchange capacity, organic complexing ability, and variable interlamellar distance are very similar to those of montmorillonite. The cation exchange capacity of vermiculite is one of the highest of all the clay minerals. The interplatelet space is accessible to penetration by some electrostatically neutral molecules. Water and glycerine are two common substances whose molecules may be so imbibed. In the natural state and under normal atmospheric conditions, water occupies the spaces between the silicate layers. The crystal d-spacing is near 14.2A. By differential thermal analysis, it has been determined that the water

Jan 1, 1983

By J. S. Aronofsky, J. P. Heller

This paper presents a mathematical analysis of the fluid mixing which occurs during flow through porous media. The analysis is based on the well-known diffusion equation with mass transfer term. It is pointed out that the use of this equation is justified by a single general assumption which does not specify any particular mechanism for the mixing process. Formulas are given for tracer fluid concentrations for two different boundary and initial conditions. Calculated numerical values compare closely with some published results of experiments in which no viscosity or density difference existed between displaced and displacing fluids. INTRODUCTION A phenomenological theory of the mixing and diffusion process is presented for the flow in a porous medium of two miscible liquids of equal viscosity and density. This is based on the classical diffusion equation which has been used to explain such processes as Brownian motion and heat conduction. More recently, this same differential equation has been derived (in a manner analogous to the work in this paper) for the problem of heat transfer during fluid flow in porous media.' The purpose of this paper is to show that the available experimental evidence justifies the use of this simple diffusion equation. The recent data of Koch and Slobod' and of von Rosenberg' are used to test this diffusion model by comparing concentration curve shapes and calculating "effective diffusion coefficients". The same differential equation for miscible displacement, as well as some experimental confirmation, has been published in Japan by K. Yuhara,4 and recently in this country by Day.5 Yuhara arrives at his con-clusions through analogy with turbulent diffusion, de- scribing the microscopic velocities in terms of "eddy motion" — although unlike the ordinary eddies of hydraulics. Day proceeds directly from Scheidegger's6 concept of the dispersion of velocities. After a very clear description of Scheidegger's theory, he applies it to calculate the motion first of a drop and then of a thin layer of brine into a bed of sand. For the latter case, Day derives and solves the same differential equation as Eq. 1. In our presentation no detailed assumptions are made as to the mechanism of displacement. In fact, several microscopic processes may well be involved—the only general assumption made is that their cumulative microscopic effect may be described by the addition of a second-order term to the differentirl equation of continuity. DERIVATION OF EQUATIONS Let us consider the displacement in one spatial dimension of Liquid A from a sample of porous media by Liquid B. A and B are assumed completely miscible and of the same viscosity and density. Both fluids are also assumed to be incompressible and it is specified that the total flow is the same for any cross section. It is evident that the rate at which Liquid B would be carried across such a section would be given by the product of the fluid velocity, u, and the concentration, C, of B in the fluid mixture assuming the constancy of u and C over the cross section. Since the volume of Fluid B is conserved, this leads to the one-way wave equation - u ?C/?x = ? ?C/?l where x increases in the direction of u—the "Darcy velocity", Q/A (or velocity flux per unit area), and 9 is the fractional porosity. Consider the case where fluid B is displacing Fluid A and the concentration distribution is represented by a step funcLon (with a vertical discontinuity). The solution of Eq. 1 predicts no mixing

Jan 1, 1958

By A. S. Emanuel

A correlation has been developed for calculating the phase behavior of, gas-condensate systems at reservoir conditions. The correlation is based on the principle of corresponding states and has been coded for an IBM 7094. Experimental K-values were determined for several gas-condensate systems at reservoir conditions to evaluate various semi-empirical parameters of the correlation. The approximate range of application of the correlation is 150 to 300F and 1,500 to 6,000 psi. INTRODUCTION The rapid development of digital computers during the past several years has made feasible the calculation of hydrocarbon phase behavior by methods based on rigorous thermodynamic principles. Good correlationsl have been developed for low to moderate pressures, but these techniques have not yet been extended successfully to reservoir fluids at high pressures. Consequently, the determination of phase behavior of oil and gas systems at reservoir conditions is still based almost entirely on generalized data correlations or on experimental studies of the fluid in question. While these methods have been used successfully many times, they do have inherent limitations that restrict their applicability. Generalized correlations, such as the NGShlA K-charts, are limited to the range of pressure, temperature and components for which the data were determined. The accuracy of these correlations is often questionable because the effect of total system composition is not well defined. Experimental studies offer a reliable method for determining phase behavior, but usually the studies are costly and time consuming. Recently, Leland and co-workers presented a new approach to calculating phase behavior from the principles of corresponding states.6-8 Corresponding states methods determine the thermodynamics properties of a given system by comparison with a reference substance whose properties are known. The accuracy of this approach depends on close chemical and structural similarity between the reference substance and the system in question and between components within the system itself. For high accuracy, it is usually necessary to correct for chemical and structural dissimilarities. In principle, however, the corresponding states method should be no less accurate at high pressures than at low pressures, provided reference substance properties are known. This paper describes an empirical modification of the basic correlation proposed by Leland, et al.6 for the specific purpose of calculating the phase behavior of gas-condensate fluids at reservoir conditions. The modified correlation, which has been programmed for an IBM 7094, may be used for either approximate or precise determination of fluid behavior depending on the amount of analytical and, experimental data available for the system. BASIC THEORY The basic theory of the corresponding states phase equilibria correlation was first published by Leland, Chappelear and Gamson.6 Subsequently, Leland, Chappelear and Leach published methods for improving the accuracy of the original theory.7-9 The aim of the correlation is to calculate the K-value of each component of a given system as a function of pressure, temperature, and over-all composition, where Once the K-values are known, the phase behavior may be determined directly by an appropriate flash calculation. The basic equation for calculating component K-values was taken from the work of Joffe. 4 For any component i of a mixture, the K-value is given by

Jan 1, 1969

By D. E. Peacock, B. Harris

Thermal-expansion and electrical-resistivity measurements have been carried out below 400°K on niobium and two niobium alloys containing tungsten. For anonaly in the expansion us temperature curve below aboct 300°K. This reslt is consistent with an earlier obsercation that Young's modilus joy Cb-132 lms an unusiully strong -tenperati-e dependence below roon temperature. No anonulous behavior was observed in the resistiuity vs temperature curve for either alloy or in the magnetic susceptibility us tenperature curve for Cb-132. These results are disclssed in terms of various electronic and bonding theories of the transition metals and their alloys, but it has not been possible to reach any definite conclusion as to the cause of the results obtained. ALONG with investigations of the low-temperature mechanical behavior of niobium and its alloys, some of the physical properties of these materials have also been studied. In the work described here, thermal expansion and electrical resistivity were measured below 400°K for niobium and two of its alloys, Cb-132 and D-43. In addition, the magnetic susceptibility of Cb-132 was measured over the same temperature range. The alloys are representative of two distinct tion-hardened alloy such as might find service as a turbine blade, while D-43 (Nb, 5.3 at. pct W, 1.1 at. pct Zr. and 0.08 at. pct C) is a solution- and dispersion-hardened alloy, of considerably lower strength than Cb-132. which is a potential vane material. Their significant common feature is the presence of tungsten as the principal source of solution hardening. EXPERIMENTAL TECHNIQUES Thermal-expansion measurements were made with a quartz-tube dilatometer using baths of liquid nitrogen, liquid bromoethane, and a liquid polyalkylene glycol derivative to cover the range 80 to 400 K. Results are recorded as expansion relative to that of quartz, and the probable maximum error in any value of relative expansion. All, is about il pct. The metallurgical condition of the specimens, which were in rod form. will be discussed in the section dealing with experimental results. Resistivity was measured on recrystallized wires 0.030 in. in diam using a potentiometric technique. The same baths were used for temperatures between 80' and 400°K as were used for the expansion measurements. and in addition the resistance of D-43 and Cb-132 was measured at 4.2 K by immersion in liquid helium. The maximum expected experimental error in any resistivity value is about 3 pct and that in any change in resistivity about i0.05 pct. Temperatures in both sets of experiments were measured with liquid-in-glass thermometers with an accuracy usually better than 1°K. EXPERIMENTAL RESULTS 1) Thermal Expansion. In Fig. 1 is shown the relative linear expansion of recrystallized niobium. D-43. and Cb-132. The curve for Cb-132 is composed of results from several runs on three different samples. Two of the samples were annealed for an hour at 1500°C to cause recrystallization and the third was annealed for an hour at 2000'C to remove some mi-crosegregation of tungsten which was known to be present in the wrought alloy received from the manu-

Jan 1, 1967

By M. K. Yen, W. R. Hibbard

Previous studies of plastic deformation of metals have emphasized the important role of bending and constraints during strain under relatively pure stresses.1"5 Some new phenomena such as early conjugate slip6 and polygonization7 are intimately concerned with the relief of bending stresses, the former by slip and the latter by a process analogous to re-crystallization. However, few analyses of the bending deformation of single crystals were found in the literature, and those8,9 are sufficiently previous to the modern interpretations of slip to invite further investigation in this area. The transverse bending mode of testing, using the two-point loading method, was selected because the center portion of the specimen may be considered as under a theoretically pure bending stress with a uniform bending moment. Considerable attention was given to both the lattice distortion and flow characteristics of the crystals during deformation. Experimental Proeedare A specimen under the action of equal and opposite couples at both ends is said to undergo pure bending. The shear and moment diagrams of a specimen under the two-point loading method of transverse bending are shown in Fig 1. The moment over the middle-third of the span is equal to PL/6, and the vertical shear is zero. Fig 2 (A) illustrates the stress distribution in the elastic range on a plane across the middle of the specimen. It can be seen that the normal stress has a maximum value at the extreme surfaces on the tension and compression sides, and decreases proportionately to zero as it approaches the neutral plane. However, when the plastic range is reached, the stress distribution is no longer a straight line relation, but changes in a manner which can be best illustrated as shown in Fig 2(B). The following calculation and stress analysis are mainly in accordance with the works reported by Timoshenkol0 and Kochendörfer.9 From the simple beam formula, the maximum normal stress will occur at the surface of the specimen. For a round specimen with a radius r and for a rectangular specimen with a height 2h and width 2b, the maximum normal stress, Sp, will be given by the following formulas: For round specimen: 2PL Sp = 3pr3 [1] For rectangular specimen: Sn = 8bh2 [2] Since it is reasonable to consider that slip takes place under the same conditions as in uniaxial loading, the resolved shear stress S?, along the operative slip direction will be as follows: For round specimen: 2PL sin x cos ? Sp = 3pr3 sin x cos ? [3] For rectangular specimen: PL S8 = 8bh2 sin x cos ? [4] where X is the angle between the specimen axis and the slip direction and X, the angle between the specimen axis and the major axis of the glide ellipse. It was also reported by Kochenclorfer9 that the critical bending moment, that is, the bending moment exerted on the specimen to initiate slip, should be higher than the value given by the equations stated above. Based

Jan 1, 1950

By John McGhee, Jack Newcombe, M. J. Rzasa

A series of water floods was made on laboratory prepared unconsoli-dated sand columns to study the effects on oil recovery of the solid-water-oil contact angle, the oil-water interfacial tension, flood rate, and oil viscosity. A procedure was developed for treating silica sand with silicone polymer to produce surfaces of varying degrees of wettability as evaluated by measuring the contact angle on a flat silica plate that had been simultaneously treated with the sand. Fluid interfacial tensions were changed by the use of surface active agents in the flood water. Water-wet, oil-wet, and intermediate wettability systems were studied at high and at low values of interfacial tension, at several flood rates, and with oils of different viscosities. For both oil-wet and water-wet systems and a low viscosity oil, recoveries were functions of the oil-water interfacial tension; also, increase in flood rate resulted in increased oil recoveries. High interfacial tension floods were more efi- cient than low interfacial tension floods on water-wet systems, while low interfacial tension floods were more eficient on oil-wet systems. Intermediate or neutral wettability systems were less sensitive to rate of flood advance and interfacial tension than either oil-wet or water-wet systems. The eflects of surface forces on oil recovery for high viscosity oils were not so well defined as for low viscosity oils. INTRODUCTION Reservoir rock surfaces vary in their wettability1, some being water-wet while others are apparently oil-wet. The factor of wettability enters into reservoir performance in both primary and secondary recovery, and into laboratory measurements of capillary pressure, connate water saturation, and relative permeabilities.' In theory, the degree of wettability is measured by the contact angle for the system solid-water-oil. A water-wet system is defined as one in which the advancing water contact angle is less than 90". An oil-wet system is one in which the advancing water contact angle is greater than 90°.3 The surface forces in a solid-water-oil system are related by the Young-Dupree equation as follows4: These surface forces are commonly expressed as the difference in adhesion tension, or the product of the water-oil interfacial tension and the contact angle. If the contact angle is less than 90°, the adhesion tension difference is positive and oil displacement by water from a solid surface will be spontaneous and favored by a high oil-water interfacial tension. If the angle is greater than 90" the displacement will not be spontaneous but will require less energy the lower the oil-water interfacial tension. The contact angle and interfacial tension also appear in the capillary pressure equation5: PC = 2ocosejr .... (2) The quantity cos was used by Rapoport and Leas6 s a factor in the generalized form of their "sealing coefficient." Warren and Cal-houn' presented the results of cos on oil recovery in oil-wet systems and correlated their data using a "scaling coefficient" suggested by Jones-Parra and Calhoun.8 The "sealing coefficients" were advanced to help correlate the effects of the variables influencing oil recovery by

Jan 1, 1956

By B. Fazan, O. D. Sherby, J. E. Dorn

McLean's technique was employed to determine the effect of temperature on the contribution of grain boundary shearing to the total creep strain in pure aluminum over the range of 610° to 747°K. The contribution of grain boundary shearing at the same stress was found to depend upon a temperature-compensated time, where the activation energy AH was practically identical with that for total high temperature creep. SEVERAL outstanding investigations by McLeanl-' have demonstrated that the total creep strain in polycrystalline aluminum at 473°K arises from microscopically detectable slip, subgrain tilting, and grain boundary shearing. Whereas microscopically detectable slip arises from dislocations that move out of the grains, subgrain tilting is due to entrapment of dislocations of like sign in subgrain bound-arles. Consequently, creep in polycrystalline aggregates appears to be the resultant of two mechanisms, namely, migration of dislocations and grain boundary shearing. This suggests that an appropriate law for creep should contain two additive terms ?1 = ?d + ?yb where the total creep strain, ?1, depends on the sums of the creep strains arising from motion of dislocations, ?d, and from grain boundary shearing, ?yb. But extensive investigations5 over wide ranges of alloys, grain sizes, temperatures, and stresses have shown that the total creep strain is quite accurately dependent on the single term functional relationship ?1 = f (?), s = constant where s is the stress (constant),* ? = te-?H/RT, t is the time under stress, AH is the activation energy, R is gas constant, and T is absolute temperature. Since Eq. 2 was found to be valid over those ranges where appreciable grain boundary shearing is observed, it appears probable that the activation energy for grain boundary shearing might be identical with the activation energy for migration of dislocations. Thus, under a given stress, it could be thought that where fd and fgb represent the appropriate functions for each mechanism of creep. The primary objective of this investigation was to obtain a direct check on the accuracy of this hypothesis. Experimental Procedure and Results Rolled aluminum sheet (99.987 pct Al), 0.100 in. thick, was used in this investigation. Creep specimens having a 0.25 in. wide, 3 in. long gage section were so machined that their axes were in the rolling direction. All specimens were annealed for 9 hr at 894°K before testing to produce a mean grain diameter of about 0.1 in. One face of each specimen was polished electrolytically in 10 pct HBF4. In order to permit accurate measurements of grain boundary shearing, the polished surface was lightly scribed with a ruling machine accurate to ±0.0001 in. to give a uniform grid of lines which were 0.01 in. apart. Although it was thought that surface working might modify the creep behavior in the vicinity of the scribed grid lines, this technique nonetheless was employed because such local differences in creep properties should not modify seriously the fractions of the total creep strain due to grain boundary shearing or to migration of dislocations. All creep tests were conducted under constant load with an initial stress of 250 psi. The specimen tem-

Jan 1, 1955

By C. E. Morris, S. R. Pollack

Current -vo1tage -temperature characteristics were studied for Al-A12O3-Al, Au structures. The oxide film was grown by gaseous oxidation in an 0, glow discharge. The electron-transfer mechanism was identified as tunneling and the barrier heights and thickness were determined. It was found that the barrier height at the parent aluminum-oxide interface is 1.5 ev and the built-in vloltage with an aluminum counter electrode is 0.35 v. The growth mechanism of the oxide is also discussed and it is suggested that Mott&apos;s theory is not applicable in describing the growth. Discrepancies were observed in the detailed comparison of tunneling theory, based on a trapezoidal model, with the data, and several possible explanations are discussed. THEORETICAL analyses of electron tunneling from one nonsuperconducting metal to another through an extremely thin insulating film have been presented.1-6 These calculations relate the current density, J, to the voltage, V, between the two metals, and to the temperature, T. The results of these calculations indicate that the magnitude of J and the detailed shape of the J-V characteristic will in general depend most strongly upon such properties of the insulating film as the thickness, the dielectric constant, and the structure of the electron energy barrier profile. The experimental results of many7-11 workers have corroborated the main features of these calcu- lations and have in fact been used to determine properties of the thin insulating film. It has been demonstrated, for example, that the tunneling current formulation of simmons4 is adequate to describe the observed J-V characteristic over approximately 10 orders of magnitude of current. Furthermore, from these observations and analyses, insight into the nature of the A1-A12O3 (thermally grown) interface region has been obtained.&apos;&apos; The temperature dependence of the tunnel current has also been studied and, although a self-consistent set of tunneling parameters can be chosen to describe these data at low voltages (i.e., <0.5 v), the tunneling parameters are different from those obtained from higher voltage (i.e., of the order of 1 to 5 v) measurements." We have measured the effect of temperature and voltage upon the tunnel current at higher applied voltages and find that the theory is inadequate to describe all of the effects reported here. A discussion of the origin of this inadequacy is presented following a description of the samples used for these experiments, and the experimental results. SAMPLE PREPARATION Metal-insulator-metal structures were fabricated in the following way. A 99.999 pct A1 charge was used to deposit, in vacuo, thirty films, 2000Å thick, onto a 1 by 2 in. glass substrate (7059 Corning glass). Prior to deposition the substrate was glow-discharged for final cleanup and then baked at 200°C for 2 hr at a pressure of 5 x 10-7 torr. The pressure during deposition of the metal films was 2 X 10-6 torr and the deposition rate was 50 to 100Å per

Jan 1, 1965

By J. R. Kyte, L. A. Rapoport, R. J. Stanclift, S. C. Stephan

Experimental studies covering a wide range of core materials and fluid properties have been conducted to determine the mechanism of oil displacement by water in a partially gas-saturated porous medium. In all instances, the presence of a gas phase was found to have a beneficial effect in reducing residual oil saturations. The practical significance of this benefit is discussed, and a simplified procedure is outlined for evaluating the effects of free gas on water flooding by means of short core tests. INTRODUCTION In addition to oil and water, reservoirs subjected to water flooding frequently contain, also, a gas phase. Common engineering procedures account for the presence of this "free gas" only from the viewpoint of volumetric balance, implying that the only role of the gas consists of providing "fill up" space. It is usually visualized that during the initial stages of the water invasion, the oil, moving ahead of the water, displaces part of the gas and that subsequently the remaining portion of the gas phase is totally compressed and dissolved in the advancing oil bank. Thus, consideration of a two-phase, water-oil flooding mechanism supplies an adequate basis for predictions of oil recovery if the pressure build-up caused by the flood is sufficiently great, so as to reduce the free gas saturation to a negligible value in any portion of the reservoir by the time that portion is reached by the advancing flood water. In many in- stances, however, waterflooding operations are carried out when the reservoir pressure is still relatively high SO that the pressure build-up associated with the water flood may not result in complete dissipation of the gas phase. Under such circumstances it is necessary to ascertain whether or not the presence of free gas has an effect on waterflood behavior and, if so, to account for any such effect in the evaluation of field operations. The difficulties encountered in attempting to arrive at a comprehensive picture of the displacement mechanism of oil by water in the presence of free gas stem from the fact that most of the published data on this subject were obtained under conditions that were either not clearly representative of those existing in a reservoir or not sufficiently controlled to permit evaluation of the effects specifically caused by the gas phase. Some of the most precise data found in the literature were obtained in tests using a "static" gas phase created prior to each water flood by circulating oil through partly gas-saturated cores until a "trapped," residual gas saturation was established throughout the system.&apos; Other investigations were conducted under more representative conditions where the gas phase was free to move during the water flood: but in most cases the results of these tests were obscured by auxiliary phenomena which arise from the compressibility and solubility of the gas and preclude quantitative conclusions.2,3,4 Prior to water flooding, a reservoir usually contains oil and gas, distributed in a manner in which both are free to move. Thus, the principal objective of the present studies was to determine the type of displacement mechanism occurring when water invades a porous medium containing oil and a "mobile" gas phase. Further-

Jan 1, 1957

By M. C. Flemings, T. S. Piwonka

Pore formation is examined analytically in cellu-larly solidified single crystals, in cylindrical castings of pure metals and alloys, and in unidirectionally solidified castings. Effects are considered of fluid flow. dissolved gas, and surface tension. Experimental results show that flow through a partially solid "mushy" aluminum alloy obeys relations similar to those applicable for flow through other types of porous beds. Flow rate is proportional to the square of the fraction liquid (fL) for fL less than 0.3. Results indicate models chosen for mathematical description of pore formation in alloy solidification are valid. ONE of the solidification variables which is influential in determining mechanical properties of castings and of wrought material produced from cast ingots is porosity. This porosity may result from shrinkage, dissolved gases, or a combination of both, and may take the form of either macroscopic pores in localized areas of the casting or microporosity evenly distributed throughout large areas of the casting. In alloys which freeze over a range of temperatures some microporosity is generally found after solidification. There has been much discussion as to whether the micropores are caused primarily by shrinkage, by gas, or by a combination of both.&apos;-3 There has, however, been little experimental or analytical work performed to delineate the relative contributions of shrinkage and gas in causing microporosity or to quantitatively describe the influence of processing variables on microporosity. Porosity other than microporosity (e.g., center line shrinkage) has also received only limited quantitative treatment.495 The aim of this work is to examine analytically the general problem of pore formation in solidification and to perform limited confirmatory experiments. Consideration is first given to shrinkage porosity in pure and nearly pure metals (shrinkage in cellular growth, centerline shrinkage in cylinders). Next, a discussion is given of formation of shrinkage-caused microporosity in unidirectional solidification and in "mushy" solidification. Finally, influences of dissolved gas and surface tension on microporosity are examined. Numerical results are given for aluminum and an aluminum alloy. Typical examples of microporosity are shown in Figs. 1 to 3, for 1) cellularly solidified aluminum, 2) sand-cast aluminum alloy, and 3) low alloy steel. SHRINKAGE POROSITY IN PURE AND NEARLY PURE METALS 1) Cellular Growth. In cellular freezing of an impure liquid, liquid grooves may exist over a considerable length of the growing crystal.6 The cross-sectional area of the grooves decreases with increasing distance from the cell tips; the root of the groove is blunt if the liquid at this point reaches eutectic composition. As a result of volume contraction during solidification, flow occurs down the groove, parallel to the growth direction; hence, there must be a pressure drop from the groove entrance to points within the groove. The geometry of the groove is simplified to that of a circular cylinder of radius r and length L, Fig. 4, L » r, and a simple and useful limiting case is obtained. This is that, in the circular cylinder freezing at constant velocity, the pressure at the root of the groove is lower than that pressure at radius r in grooves of other geometries where r decreases continuously to a lower limit at L. As discussed below, a pore resulting from a given pressure difference

Jan 1, 1967

By B. H. Bergstram, C. L. Sollenberger

When glass spheres are crushed by slow compression loading, the outer lune-shaped fragments resulting from the fracture consistently fly outward at high velocity. About 45 pct of the strain energy fed into the sphere-platen system reappears in the form of the kinetic energy possessed by these fragments. These measurements were obtained from high-speed photographs. This kinetic energy can be utilized for addtional comminution if the fragments are allowed to impact against a suitable structure. Comminution theory has advanced considerably in recent years by the study of this unit operation under laboratory conditions. Some of the parameters have been isolated, but the experimental measurements of the pertinent quantities, such as energy, particle size, etc., are still only made prior to and after the conclusion of a fracture sequence. The catastrophic fracture process itself has not been adequately studied from a comminution point of view. It has recently been suggested by Gilvarry that the fragments which result from the free crushing of a single glass sphere by slow compression contain a considerable amount of kinetic energy after fracture. This kinetic energy, if properly directed, possibly can be partially utilized to cause additional fracture. In this paper a hypothesis will be presented which suggests that this kinetic energy is necessarily present and that its dissipation is the process by which a large part of the thermal energy imparted to the crushed product is generated. Evidence supporting the existence of kinetic energy and its utilization for secondary fracture will also be presented. THEORETICAL A previous paper4 has shown that the load-deformation curve of a system composed of an elastic sphere gripped by a pair of parallel platens, can be predicted by a consideration of the theories of Hertz as summarized by Timoshenko.&apos; The deformation A, of one half of the system is given by which states that the deformation is a function of the mutually applied load, P, the radius of the sphere, R; and the respective physical properties of Poisson&apos;s ratio, v, and Young&apos;s modulus of elasticity, Y. The shape of experimental load-deformation curves conform to this relationship. The energy input to the sphere-platen system, can be predicted from the integral of the load acting through the appropriate distance, dA. Since dA can be evaluated in terms of P from the derivative of Eq. 1, this integral can be integrated between limits of 0 and P. The test sphere, however, is loaded at not one but at two points of contact, i.e., at both ends of a diameter, and hence the total energy input, E, is twice the energy input at each end, or The energy input to the system calculated from this equation corresponds to the energy computed from the measured area under the load-deformation curve. Timoshenko, considering the impact of two perfectly elastic spheres, has shown that the inter-penetration distance A is given by - m f (D (^ * ^t w as modified for the special case where one sphere is of infinite radius, i.e., its surface is a plane; where M is the mass of the finite sphere of radius R, and V is the velocity of approach at the start of an impact and also equals the final velocity of separation for perfectly elastic bodies. Because A is also a function of the load and the elastic properties, as shown in Eq. 1, the velocity, V, can also be calculated in terms of the elastic properties and maximum load, P, or

Jan 1, 1961

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge&apos;s method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey&apos;s theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey&apos;s method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge&apos;s6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

J. B. Cohen (AVCO Mfg. Co.)-The model presented by the authors is only one of several possible mechanisms for low-temperature aging. It is true that the motion of dislocations to the copper atoms, followed by diffusion along the dislocation lines, is one such mechanism. However, the details of such a process which would :result in copper-rich platelets on (100) planes in the f.c.c. matrix are not clear. Furthermore, there is at least one other mechanism that may be used to e.xplain the authors&apos; results, so that the data do not clearly favor the Turnbull model. The zones may form by "normal" diffusion. The diffusion constants would be higher at room temperature than extrapolated data would suggest due to vacancies trapped in the specimens in the quench. The diffusivity would then be multiplied by the ratio of the number of vacancies present at the solution temperature to that present at room temperature. These vacancies could remain in the alloy specimen for long times at room temerature, due to interaction with the copper atoms. &apos; The activation energy for diffusion is smaller than that for bulk diffusion at elevated temperatures for it is only that for the motion of vacancies existing after the quench; after these are used up, the aging process effectively ceases. It is somewhat questionable if the laws governing transformation kinetics can be directly applied to changes in elastic moduli, especially in the case of aging where the structure and the composition of the zones may be changing during aging, If such laws are assumed, it is possible to explain the exponential factor of % with this vacancy model. It is probable that the structure of the zones and the matrix match only on (100) planes. Certainly the precipitates that form in this system after various aging treatments show orientation relationships similar to those of the zones. Strains on the edges of the zone may set up fields which do not allow the zone to grow in an edgewise manner. The exponent in the growth of a thin plate in a direction perpendicular to its face would be %. It may also be that copper atoms are attracted to dislocations. If this is the rate-controllin factor in the growth process Ifthisthen an exponent of % or % may arise. The writer would like to know if the authors have tried plotting l/a; proportional to as might be expected for "dislocation drainage." R. H. Doremus (General Electric Research Laboratory)—An important new method for following the growth of Guinier-Preston zones is presented in this paper, and I wish to congratulate the authors on their experiments. Their interpretation of the results, however, is doubtful in certain respects. To determine the exponent m in Eq. [4] the authors plotted -log at against t to various powers (their Fig. 6). The initial value of frequency used for these plots was found by extrapolating the data to zero time. In Fig. A the data from their Fig. 2 are plotted against a linear, rather than logarithmic, time scale. To obtain graphs similar to their Fig. 6 for these data one must use an initial frequency of about 52720; such a plot for m = % is shown in Fig. B. The straight line in Fig. B fits the data in an intermediate range only. The curve drawn in Fig. A was mapped from this line; comparison of

Jan 1, 1960

By Mary Barrett, Walter Rose

INTRODUCTION The conventional way to program reservoir problems* for solution on the electric analyzer1 is one which concentrates the block units towards the inner radius of the prototype system being represented by the analog model. This is achieved by the device of considering block units of equal value radial resistance, leading to the consequences: (1) the ratio of the external to the internal radius of each block unit is a constant (i.e. r1/ri-l = constant a); (2) the external radius of the nth block unit is equal to the innermost radius of the system times the constant a raised to the nth power: (3) the capacitance of the block units increase by a factor of a2 going out radially from the innermost to the outermost boundaries of the system; (4) given n block units in the radial direction, the capacitance ratio of the outermost to the innermost block unit is a (2n-2); and (5) the vertical resistance of the block units decrease by a factor of a going out radially from the innermost to the outermost boundaries of the system, such that the ratio of the maximum to minimum vertical resistance is also given by the factor a(2n-2) PRESENT PROGRAMMING METHOD While this conventional way of programming problems is quite adequate for routine uses of the electric analyzer, it leads to difficulties in those cases where the inner radius of the prototype system (e.g., the well-bore) is small compared to the external radius (e.g.. the external reservoir boundary). To handle such a case by the conventional method would require either that a large value be selected for the constant a or that a large number n, of radial block units be chosen. In consequence, the capacitance and vertical resistance ratios become impossibly large, as will be evident from the relations given above. In any event, the size of the smallest and largest capacitors will always be limited by the size of commercially available components, and the analyzer itself will determine how many mesh points are available for radial block units. For example, it can be shown that if the maximum capacitance ratio which can be electrically handled is a factor of l04, the maximum external radius of the system which can be considered (assuming an internal wellbore radius of 0.25 ft is: It is thus seen that in direct consequence of having the inner wellbore radius small compared to the external radius, models extending beyond 30 to 40 ft are not possible when the capacitance ratio is limited to l04. (NB. A value of 10-9 farads can perhaps be taken as the minimum electrical capacitance possible if corrections for the capacitance of the lead wires are to be avoided. Values greater than 10-5 farads would require the use of unstable and inaccurate electrolytic capacitors). On the other hand, models extending out to 500 to 5,000 ft may be desired in single-well problems, and hence the need for a revised method of problem programming. * REVISED PROGRAMMING METHOD The revised method of programming is one which attempts to retain the desirable features of the conventional method (viz. concentration of block units towards the inner boundary, simplicity in computation, etc.) without having as a consequence that impossible capacitance and vertical resistance ratios will be required. Other advantages will also be achieved by the revised

Jan 1, 1957

By M. D. Hassialis

[Surface Areas of Concentrates and Collector Coatings (T.P. 2002, by A. M. GAUDIN and G. S. PRELLER, Min. Tech., May. Discussion by M. D. HASSIALIS and the authors) . I Activation of Minerals and Adsorption of Collectors (T.P. 2082 by J. ROGERS and K. 1&apos;. SUTHERLAND, Min. Tech., Nov. Discussion by A. F. TAGGART and the authors). . . . 3 Correlation between Mineral Behavior in Cataphoresis and in Flotation (T.P. 2005 by A. If. GAUDIN and S. C. SUN, Min. Tech., May. Discussion by M. D. HASSIALIS). . . . . . 8] Surface Areas of Concentrates and Collector Coatings M. D. HASSIALIS.*-In the current flotation theories there exists a common gap spanned by a common assumption pertaining to the arrangement and orientation of the collector ions at the solid surface. Also common is the belief that contact angle and flotation performance are dependent upon the arrangement and orientation of collector ions. This stimulating and informative paper sheds some light upon this general problem. One of the indices of arrangement is the area available per collector ion; in this connection, Gaudin and Preller have applied the B.E.T. method for surface determinations. In a surprising number of cases, their results show a greater specific surface for tailing than for concentrate. This cannot be accounted for by the differences in specific gravity and is attributed to the presence of clay minerals of "primary" sluices in the tailing. This is not a satisfactory explanation, for experience shows that when free slimed clay minerals are present in a flotation cell they will be carried mechanically into the concentrate. A part of the answer to the decreased specific surface of the concentrate would appear to be its usual flocculated condition, yet this should not entirely account for the large spreads reported. The possibility occurs that the nitrogen adsorption isotherms for clean mineral and collector-coated mineral are different. If that were true, and if the coated mineral surface exhibited decreased adsorptive powers, the discrepancy might vanish. In discussing the disposition of collector ions at mineral surfaces, the authors adopt the more consistent viewpoint that there exists some agreement between the ions of the substrate and those of the superstrate. In the case of xanthate ions on galena, a possible arrangement is shown in Fig. 2. This arrangement shows every xanthate ion being shared by two lead ions, whereas the chemical composition calls for two xanthate ions for each lead ion (p. g). This whole discussion is misleading. In the first place, it is known that xanthate ions are not abstracted by unoxidized galena, hence Fig. 2 should show as the substrate not Pbs but PbS„ Om. It would be equally erroneous to use the lattice arrangement of PGSO4, since the surface properties of xanthate-conditioned anglesite differ from those of xanthated galena. It is highly probable that the state of the galena surface exploited in its flotation is ephemeral, since weathered galena is only poorly, if at all, collected by xanthate. Just how a substrate of oxidized galena is to be pictured depends upon one&apos;s conception of the transformation from galena to anglesite. This writer is inclined to the following sequence of events. Surface sulphide ions undergo progressive oxidation. As more oxygen is taken up by the sulphide ions they become displaced from their equilibrium positions in the galena lattice, thus accommodating for the larger volume of the S„ Om ions. Concomittant with this change the forces acting on the lead ions

Jan 1, 1946

By T. D. Murphy

The element silicon, with its usual partner, oxygen, plays the same role on this planet relative to inorganic materials as carbon and hydrogen play with respect to living organisms. The crystallographic structure of silicon dioxide (SiO,) consists of one atom of silicon well-nigh inseparably bonded to four contiguous atoms of oxygen, a gas of but slightly more than half its atomic weight, forming a three-dimensional network of SiO, tetrahedra (Dietz, 1968). This wispy arrangement is the fabric of the mineral quartz which is one of the harder, more abrasive, and chemically stable raw materials to be found in Mother Nature&apos;s cupboard. Most of us are vaguely familiar with the observation that nearly 60% of the lithosphere to an approximate depth of 10 miles is composed of these two elements, which proportion actually rises to about 75% if the atmosphere, hydrosphere, and biosphere are included (Parker, 1967). Elemental silicon is a brittle steel-gray metalloid with a density of 2.42 and does not occur in nature. In its limited commercial form, which is still brittle and metallic in appearance, its density drops to 2.33. It combines with oxygen to form the gaseous oxide SiO, or the solid dioxide SiO, (analogous to its congener carbon in CO and CO,); a tetrachloride SiCl, (analogous to CC1,); or, combining with oxygen and one or more metal ions, forms the largest rock-making family of all, the silicates. Further, next to the feldspars, quartz is the most abundant known mineral and, in some form or another, accounts for roughly 12% of surficial terrestrial rock. It is this so-called "free silica" to which your attention is chiefly directed in this chapter. In any such consideration of silica, its several crystalline phases with their attendant differing physical properties appear worthy of delineation. Alpha quartz is the only one of this interesting polymorphic family that is thermodynamically stable at normal pressures and temperatures up to about 573°C. The low phases of tridymite and cristobalite can exist in the so-called metastable state below 117°- 200°C but seldom occur in nature in these two different crystal forms (Sosman, 1965). Opaline silica, now considered by most researchers to be a variety of cristobalite, occurs in abundance and is the only exception. The three remaining polymorphs, keatite, coesite, and stishovite, also are metastable under ambient conditions. Keatite, however, is not known to occur naturally and has only been synthesized in the laboratory. The other two have been identified in submicron sizes and in microscopic amounts, occurring naturally as a result of high order shock-wave pressures generating quite elevated temperatures-such as might accompany the impact of a large meteorite on a hard sandstone terrain (Chao, et al., 1962). An intriguing characteristic of these silica minerals, and one which has very practical applications, is their extremely variable densities. The classic 2.65 specific gravity of alpha quartz drops to 2.20 in tridymite, cristobalite, and lechatelierite (natural silica glass), as well as in many man-made glasses of commerce which are merely noncrystalline or amorphous solids in the form of super-cooled liquids which can no longer change their shapes (Pauling, 1950). It rises again to 2.50 and 2.93 for keatite and coesite, respectively, and finally becomes a reported 4.35 for stishovite-nearly a 50% range in density for the known physical combinations that silica naturally assumes (Frondel, 1962; Konnert, et al., 1973; Pecora, 1960).

Jan 1, 1975

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior