Search Documents

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By A. Kenneth Schellinger

DURING the comminution of a brittle material in the presence of dry air, no known phase change or chemical reaction takes place. The energy changes associated with the comminution are those of the transformation of kinetic energy into heat energy and the energy absorption represented by the actual rupture of molecular and/or ionic bonds. The random nature of the actual grinding process in a tumbling-mill often precludes the effective application of a rupture force on each individual particle of brittle mineral present. The result must be a great deal of useless vibration on an atomic scale in both the brittle mineral and the grinding medium below the elastic limit, which vibration is eventually manifested as heat energy in the mill. The actual rupture of molecular, or ionic, bonds creates new surface on the brittle minerals; each square centimeter of this new surface then must possess a number of potential bonds, each of which could reform theoretically with the release of energy. The sum total of this potential energy of bonding per unit of new surface, minus bond satiation effects of the enclosing gas or liquid, must be the surface energy of the mineral. If the comminution of the brittle mineral were carried out within the chamber of a calorimeter, the energy released as heat could be determined rather accurately. Moreover, the total kinetic energy input to the comminuting mass (mineral and medium) conceivably could be derived by a suitable torquemeter and revolution counter if the whole calorimeter were revolved. The difference between these two energy figures should then represent the net energy (gross energy input—heat enerb output) used in creating new surface, gliding or slip on atomic planes, twinning, and possibly other lattice rearrangements within the mineral. If it be postulated further that with a hard, brittle crystalline material, such as quartz, the energy used to produce gliding, twinning, and other permanent lattice movements is a very small fraction of the total net energy absorbed, then this energy absorption within the grinding calorimeter can be ascribed to the creation of new surface energy alone. Such energy is, from the thermody-namic standpoint, the true effective work of the grinding device, in this case a tumbling-medium mill. The ratio of this effective work (in calories) to the total energy input (in calories) should express the thermodynamic efficiency of the comminution taking place under the conditions used in the grinding calorimeter. This approach to the measurement of the energetics of tumbling-mill grinding was suggested by some previous work by Fahrenwald and coworkers' with a laboratory ball mill and a thermometer. They measured the efficiency of grinding based on the formula: Pct Grinding Efficiency = Energy (calories) — Thermal Energy x 100 1npZ Output (calories) Energy input (calories) Efficiencies thus measured, while grinding quartz sand, ranged from 19.7 pct to 7.4 pct. Cook in 1915,' used this approach to the problem of grinding efficiencies in stamp batteries and obtained the energy input by stamp weights times drop distances and the energy output by means of thermometers and pulp flow rates. He calculated that 80 pct of the input energy to the stamp shoes was wasted as sensible heat in the ore pulp. No other studies of comminution by this method of energy differences came to the author's attention. Reasonable agreement as to the order-of-mag-nitude of comminution efficiencies in tumbling-media mills is lacking in the literature. There is a widespread impression among technical workers that the grinding efficiency of the industrial ball mill may be on the order of 1 pct, or even a fraction

Jan 1, 1952

By A. Gittens

In two recent papers stevens13,14 has compared the methods used by various investigators to calculate the contribution of grain boundary sliding to the total creep strain. In previous work Gittins15 has shown that the strain due to grain boundary sliding Egb may be estimated from the lateral, d, and vertical, v, displacements of a longitudinal marker line as it crosses a grain boundary. He obtained the formula, where nL is the number of grains per unit length measured on the specimens before creep, and \p and ?2 are, respectively, the angle between the grain boundary and the stress axis in the plane of the surface and normal to the surface (see Ref. 13, Fig. 1). The subscript L indicates that observations are made in a longitudinal direction parallel to the stress axis. It was generally found more convenient to measure the grain size of a specimen before creep; this was necessary to characterize the specimen. Stevens13,14 however prefers to define Egb in terms of NL, the number of grains per unit length after creep. Since nL - NL(l + €1) where cr is the total strain, Eq. 1 11 becomes The quantities d and ?, can be determined directly from the marker lines, while v and ?2 are found by sectioning normal to the surface. In the interior when symmetry considerations suggest that v = d and ?1 - ?2 the strain due to grain boundary sliding may be estimated from a suitable internal longitudinal marker line by the equation To determine cgb at the surface using either Eqs. |1| or [2] requires a large number of observations, particularly if many estimates of cgb are required, and also it is impossible to determine the variation of cgh with creep strain on the same specimen because the specimen must be sectioned to determine ?2. Gittins15 found that for a number of materials the sliding vectors d and v were related to each other so that provided one was measured the other could be estimated. Also during creep, grain boundaries tend to migrate to meet the surface at angles close to 90 deg,16 so that the term nL(v/tan?2)L was comparatively small when 4,- 90 deg. By determining the frequency distributions of ?2 on control specimens this term can be estimated without sectioning the creep specimen. In this way an empirical equation based on Eqs. |1| and |2| enables cgb to be determined with reasonable accuracy with fewer observations of the angles and sliding components. For B brass, the relationship between the sliding vectors indicated that ceb = 0.85nLVL . In Table I of Stevens paper this equation is stated to be incorrectly derived. However it should be pointed out that this is an empirical equation derived from Eq. |1|, based on a knowledge of the relationship between d and v and on the frequency distributions of ?1 and ?2 for a particular material under given experimental conditions. It was used only when these relationships were known, to reduce the number of observations required to estimate cb. Therefore the comment "derivation incorrect does not apply. Gittins did in fact make use of Eqs. |1| and 131. Author's Reply I must apologize to Dr. Gittins for the inadvertent misrepresentation of his work. Other slight errors in Table I have been pointed out to me. Langdon (No. 12 Table I), Gifkins (No. 7, Table I), and possibly other workers did not measure the number of boundaries per unit length after creep (A') but the corresponding quantity bejove creep (n). As Dr. Gittins points out, a completely analogous set of equations to those given in the paper can be derived by the same method: using the number of boundaries per unit length before instead of after creep. The basic equations corresponding to Eqs. [2] and 171 of the paper, from which all the others can be derived, are

Jan 1, 1968

By P. B. Mee

The secondary re crystallization behavior of two vacuum-melted, nominally 3 pct Si-Fe alloys has been studied, one alloy being produced from pure electrolytic iron and semiconductor-grade silicon and the other from commercial-purity hot band. The vacuum-melted and vacuum-remelted alloys were found to develop {110}(hhl) and {lll)(hkl), respectively, on annealing in high vacuum and {100}(0kl) and {11l}(hkl), respectively on annealing in pure hydrogen of dew point, < — 60°C. H2S doping of the hydrogen or annealing in commercial argon was found to effect a reversal for the vacuum-remelted alloy from (111) to (100) secondary recrys-tallization, while prolonged vacuum annealing at 2 x 10-6 Tow resulted in a tertiary re crystallization from (111} to (110). The growth selection of the {lll}(hkl) component was found to be governed predominantly by the {111}/planar misorientation, the maximum deviation of the resulting secondaries being 5 deg. It is considered that the driving force for {111}(hkl) secondary recrystallization arises from the preferential lowering of ?111, possibly by phosphorus or nitrogen diffusion, just as oxygen and sulfur under certain conditions have been found to preferentially lower ?100. When heavily cold-rolled metals are annealed at high temperatures, primary recrystallization usually occurs rapidly and may be followed by considerable grain growth. However, in the presence of a fine dispersion of inclusions1"3 or with a strong matrix preferred orientation4&apos;&apos; the matrix may be stabilized and a secondary recrystallization may then be effected. Also, in vacuum-melted Si-Fe it has been shown possible to develop {100}(001) and {110}(001)6 surface energy driven secondary recrystallization textures. Further, it has been found that, on annealing very pure metals, grain growth normally ceases when a grain size approximately twice the strip thickness is acquired, due to the formation of V-shaped thermal grooves which result in an equilibrium being set up between surface and grain boundary stresses. If, however, a surface energy difference ??s, exists between adjacent grains (of the order of a few percent of the grain boundary energy ?b), then according to Mullins8 a secondary recrystallization may be effected, since the grains having lowest ?s would grow at the expense of their higher-energy neighbors. For clean bee metal surfaces ?110 should be lowest; however, for 3 pct Fe-Si (100) secondary structures have frequently been developed and have been attributed by Walter and Dunn9 to an alteration of the relative surface energies due to the preferential chemisorption of oxygen on (100) planes. Kohler,10 however, indicated that this reversal may be achieved by employing minute quantities of highly polar compounds in the annealing atmosphere, i.e., sulfur compounds, in par- ticular H2S being beneficial in (100) growth. This work has been verified by Kramer" using {100)/{110} bicrystals, while Mager12 demonstrated the importance of sulfur diffusion within the strip. Detert13 has intimated that carbon may be responsible for the preferential lowering of ?110 although no direct evidence has appeared yet to support this supposition. It must , be appreciated that the above-mentioned examples of surface energy driven secondary recrystallization were obtained with very pure vacuum-melted materials with oxygen contents (which is now known to be a critical factor") less than 35 ppm. During the course of a textural study on Si-Fe it has been found that under certain conditions (111) secondary recrystallization can be obtained by what is believed to be a surface energy growth selection mechanism. The purpose of the present paper is to briefly describe this result and to try and postulate the possible causes. EXPERIMENTAL Two vacuum-induction-melted, nominally 3 pct Si-Fe alloys were produced, one from pure electrolytic iron and semiconductor-grade silicon designated A and the other from commercial-purity hot band, designated B.

Jan 1, 1969

By Avnulf Muan, P. V. Riboud

The effect of Cr2O3 on melting relations of iron oxide at oxygen pressures slightly above those prevailing in contact with metallic iron has been determined. Liquidus and solidus temperatures of wüstite are raised by the addition of Cr2O3, from 1385°C for pure iron oxide to 1420°C at the peri-ledie point where wüstite, iron chromite, liquid, and gas coexist in equilibrium under the atmospheric conditions used in the present investigation (CO2/H2 = 1). It has been shown in a previous investigation1 that the liquidus and solidus temperatures of iron oxide are drastically increased by the addition of Cr2O3 under the relatively oxidizing conditions prevailing in air. This is to be expected in view of the easy substitution of cr3+ for Fe3+ in crystalline oxide structures which are stable at high oxygen pressures, as demonstrated for instance by the completeness of the solid solution series Fe2O3-Cr2O3 at high temperatures (>1000°C). The stable phase of iron oxide at oxygen pressures approximating those prevailing when metallic iron is present is wüstite. This phase, although commonly represented by the simplified formula FeO for sake of convenience, in reality has a defect structure in which a part of the iron ions are present as Fe3+. It is to be expected that at least a limited amount of cr3 ions will substitute for Fe in this structure when Cr2O3 is added to wüstite. Depending on the ease of this substitution in comparison with that of cr3 substitution for Fe in iron oxide-Cr2O3 liquids under the same atmos- pheric conditions, the liquidus and solidus temperatures of wüstite may increase or decrease upon addition of Cr2O3. The methods of establishing these melting relations and the results obtained are presented in the following. I) EXPERIMENTAL METHOD 1) General Procedure. The quenching technique was used in this study. Oxide mixtures were heated at chosen constant temperatures in an atmosphere consisting of equal parts of CO2 and H2 until equilibrium was established among gas and condensed phases. The samples were then quenched in mercury to room temperature, and the phases present were determined by microscopic and X-ray investigation. 2) Preparation of the Mixtures. Starting mixtures were prepared by weighing and thoroughly mixing Fe2O3 and Cr2O3, or wüstite and Cr2O3. The Fe2O3 and Cr2O3 used were "Baker Analyzed" reagent grades, whereas wüstite of composition was synthesized by reducing Fe203 at 1325°C in an atmosphere consisting of equal parts of CG andHz. In order to avoid contact with crucibles, the equilibration runs were made with small pellets suspended from thin (-0.004-in.) platinum wires. Mixtures made from FOs and Cr203 were readily pelletized at a pressure of 500 atm in a 9.5-mm die, whereas mixtures made from wüstite and Cr203 when treated identically were too weak for handling. In order to obtain a stronger product, these pellets were sintered by heating for 2 min at -1200°C in the same atmosphere as was subsequently used for the equilibration runs (CG/H2 = 1:1). During the sintering, these pellets were wrapped in thin (0.0004-in.) platinum foil, which was stripped off before the equilibration runs. In order to reduce the time needed for equilibrium to be attained, preliminary equilibration of relatively large pellets of each mixture was carried out at -1200°C for 4 to 12 hr. Smaller pellets

Jan 1, 1964

By C. H. Li, J. A. Sarteli, W. Bonfield

The stress to cause a permanent micros train of 2 x 10-6 in. per in. (defined as the microscopic yield stress) in beryllium is found to be very sensitive to surface condition. The initiation of plastic flow in as-machined specimens, which contain a high density of twins and large residual compressive stresses to a depth of about 0.010 in. from the surface, occurs by the nucleation of slip from high stress fields around twin tips in the surface layers. Removal of the damaged surface layers by chemical polishing results in an appreciable increase in the microscopic yield stress, which is attributed to the remoual of stress raising twins rather than to the release of residual stresses. IT is well known that the macroscopic mechanical properties of beryllium are very sensitive to the nature of its surface finish. Matthews and coworkers1 have shown that the tensile strength of rolled beryllium sheet can be considerably enhanced by chemical polishing, presumably due to the removal of the high concentration of notches, twins, and large residual compressive stresses normally found in the surface layers of ground specimens. There was no similar effect of surface condition on the macroscopic yield stress (the stress to produce 0.2 pct strain), a point to which we shall return. To date there have been no extensive studies on the effect of surface condition on the microstrain region of beryllium. Recent measurements2 of the stress to cause a permanent microstrain of the order of 10-6 in. per in. on machined hot pressed QMV beryllium, gave values ranging from about 600 to 3500 psi. Reasons for the considerable scatter in the results were not discussed. In this paper, we wish to show that the microstrain region of polycrystalline beryllium is, in fact, very sensitive to surface condition. The stress to cause a permanent strain of 2 x 10-6 in. per in. which for convenience we will refer to as the microscopic yield stress, can be increased by almost a factor of 2 by minimizing the amount of surface damage. The variation of microscopic yield stress is attributed to the presence of twins in the surface layers, rather than to residual compressive stresses. EXPERIMENTAL PROCEDURE The starting material for this investigation was commercial purity hot pressed QMV beryllium, with an average grain size of about 20 . The effects of surface treatment were assessed by comparing the microstructure with residual stress measurements and microstrain properties. Specimens were prepared for metallographic examination by wet grinding with 240, 400, and 600 papers, followed by polishing with a suspension of fine alumina in warm 20 pct oxalic acid. After this treatment, the surface could be examined, either in the unetched condition with polarized light, or after etching in 2 pct hydrofluoric or boiling 20 pct oxalic acid with normal light. Thin films of beryllium suitable for electron transmission study were prepared from bulk specimens by chemical polishing with hydrofluoric and phosphoric acids, and were examined in an RCA-EMU-3 electron microscope. Residual stress measurements were made by a procedure similar to that of Treuting and Read.3 Specimens 4.0 in. by 0.625 in. by 0.080 in. were prepared by three different machining procedures, namely: turning, grinding, and sawing, in order to introduce different levels of residual stresses. An electrical resistance strain gage was attached to one side of the specimen, and that surface was then protected by an acid-resistant enamel, while the opposite surface was removed by etching approximately one thousandth of an inch at a time. The change in specimen curvature in the longitudinal direction was noted and the residual stress was calculated by the following formula:

Jan 1, 1963

By P. M. Blair

This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. ln the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. INTRODUCTION The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradients, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942.1 Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers.2"4 The influence of capillary forces in laboratory water floods has also been discussed by several authors.6 While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes7 pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured-matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been &apos;published on the experimental techniques involved?? Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. However, recent advances in the applica-

Jan 1, 1965

By F. Gordon Smith

ALTHOUGH several geological indicators of the critical type are known, including quartz inversions and decomposition of hydrous minerals such as serpentine, there are very few of the general type. Solid solutions are excepted, but the limitations of use are very restricted and interpretations are sometimes ambiguous. General methods for determining temperature and pressure conditions during the crystallization of minerals would have considerable scientific and economic value. It is not the purpose of this paper to discuss the various methods of geological thermometry and barometry, but to present one general method, applicable to all minerals, and to describe what progress has been made in the methods of measurement. The general method, in brief, is a study of the stress conditions in and around various types of foreign inclusions which are trapped in minerals during growth. The method depends upon the fact that any homogeneous gas, liquid, or solid will in general have coefficients of thermal expansion and volume compressibility different from those of any given mineral. Therefore stress must develop in and around all types of inclusions in minerals if the temperature or the pressure, or both, are changed from the conditions which prevailed during the growth of the mineral. The methods of measurement consist of determinations of the temperature-pressure conditions of fit of the inclusions in the host mineral. The types of inclusions in minerals are: I—gas, or liquid plus vapor, when observed at room temperature, due to crystallization under pneumatolytic conditions; 2—liquid, or liquid plus vapor when observed at room temperature, due to crystallization under hydrothermal conditions; 3—glassy solid, or devitrified glass, due to crystallization under magmatic or high temperature metamorphic conditions, in a siliceous liquid; and 4—crystals, due to overgrowth of other minerals crystallizing simultaneously or of other minerals which crystallized previously. A survey of the literature shows that much valuable earlier work on inclusions, especially that car- ried out in England in the last century, has dropped out of current knowledge. The following is a brief summary of the significant contributions to the problem up to the present day. Davy in 1822 asserted that fluid inclusions in minerals consist of an aqueous solution of salts and a gas bubble, the whole being either at lower or higher pressure than atmospheric.&apos; At intervals from 1823 to 1862, Brewster contributed information concerning other types of inclusion consisting of 1—aqueous solution, a much more expansible liquid, and a gas; 2—aqueous solution, salt crystals, and a gas; and 3—the very expansible liquid and a gas. The very expansible liquid fills the gas space between 20" and 30°C. Compression strain effects were seen around inclusions in diamond, topaz, and other minerals.&apos;-" Sorby in 1858 and 1869 further advanced the study begun by Davy, stating that fluid inclusions represent a sample of the mother liquor of crystallization and that the degree of filling of aqueous inclusions at room temperature defines the temperature-pressure relations during formation of the host. The degree of filling may be measured by determining the minimum temperature of filling of the inclusion by the liquid phase. The very expansible liquid in some fluid inclusions is liquid carbon dioxide. The temperature at which salt crystals in fluid inclusions completely dissolve in the fluid is the minimum temperature of formation. Inclusions of glass or devitrified glass indicate crystallization from a melt. Inclusions of crystals in minerals are often centers of strain, which may be seen by optical effects or by radial tension cracks. Sorby realized that an analysis of stress-strain relations about inclusions could be used to provide precise data on the temperature of crystallization, but the matter was never pursued.&apos; , Hartley (1876, 1877),".&apos;" Hawes (1881)," Wright (1881),&apos;V ohnsen (1920),I3 and Holden (1925)"

Jan 1, 1953

By F. Gordon Smith

ALTHOUGH several geological indicators of the critical type are known, including quartz inversions and decomposition of hydrous minerals such as serpentine, there are very few of the general type. Solid solutions are excepted, but the limitations of use are very restricted and interpretations are sometimes ambiguous. General methods for determining temperature and pressure conditions during the crystallization of minerals would have considerable scientific and economic value. It is not the purpose of this paper to discuss the various methods of geological thermometry and barometry, but to present one general method, applicable to all minerals, and to describe what progress has been made in the methods of measurement. The general method, in brief, is a study of the stress conditions in and around various types of foreign inclusions which are trapped in minerals during growth. The method depends upon the fact that any homogeneous gas, liquid, or solid will in general have coefficients of thermal expansion and volume compressibility different from those of any given mineral. Therefore stress must develop in and around all types of inclusions in minerals if the temperature or the pressure, or both, are changed from the conditions which prevailed during the growth of the mineral. The methods of measurement consist of determinations of the temperature-pressure conditions of fit of the inclusions in the host mineral. The types of inclusions in minerals are: I—gas, or liquid plus vapor, when observed at room temperature, due to crystallization under pneumatolytic conditions; 2—liquid, or liquid plus vapor when observed at room temperature, due to crystallization under hydrothermal conditions; 3—glassy solid, or devitrified glass, due to crystallization under magmatic or high temperature metamorphic conditions, in a siliceous liquid; and 4—crystals, due to overgrowth of other minerals crystallizing simultaneously or of other minerals which crystallized previously. A survey of the literature shows that much valuable earlier work on inclusions, especially that car- ried out in England in the last century, has dropped out of current knowledge. The following is a brief summary of the significant contributions to the problem up to the present day. Davy in 1822 asserted that fluid inclusions in minerals consist of an aqueous solution of salts and a gas bubble, the whole being either at lower or higher pressure than atmospheric.&apos; At intervals from 1823 to 1862, Brewster contributed information concerning other types of inclusion consisting of 1—aqueous solution, a much more expansible liquid, and a gas; 2—aqueous solution, salt crystals, and a gas; and 3—the very expansible liquid and a gas. The very expansible liquid fills the gas space between 20" and 30°C. Compression strain effects were seen around inclusions in diamond, topaz, and other minerals.&apos;-" Sorby in 1858 and 1869 further advanced the study begun by Davy, stating that fluid inclusions represent a sample of the mother liquor of crystallization and that the degree of filling of aqueous inclusions at room temperature defines the temperature-pressure relations during formation of the host. The degree of filling may be measured by determining the minimum temperature of filling of the inclusion by the liquid phase. The very expansible liquid in some fluid inclusions is liquid carbon dioxide. The temperature at which salt crystals in fluid inclusions completely dissolve in the fluid is the minimum temperature of formation. Inclusions of glass or devitrified glass indicate crystallization from a melt. Inclusions of crystals in minerals are often centers of strain, which may be seen by optical effects or by radial tension cracks. Sorby realized that an analysis of stress-strain relations about inclusions could be used to provide precise data on the temperature of crystallization, but the matter was never pursued.&apos; , Hartley (1876, 1877),".&apos;" Hawes (1881)," Wright (1881),&apos;V ohnsen (1920),I3 and Holden (1925)"

Jan 1, 1953

By M. A. El-Aouar

Drilling muds are complex colloidal systems. In an effort to analyze the physical properties that affect their performance at high temperature, an investigation was made of the hydrothermal stability of the monoionic hydrogen, calcium and magnesium forms of montmorillonite and attapulgite. Samples were heated to 400°1F in a simple autoclave and analyzed for cation exchange capacity, c-axis X-ray diffraction spacing and pH factor us sodium hydroxide titration dependence. Out of the six monoionic clay minerals tested, calcium montmorillonite and magnesium attapulgite showed the least base exchange capacities and were relatively unaffected to heat. It is concluded that, of the six clays studied, these two clays were the most suitable forms for use in drilling fluids. INTRODUCTION The drilling of deep wells is dependent upon the drilling mud, and problems frequently arise concerning the effect of mud temperature and prolonged heating. Drilling muds are complex colloidal systems. In an effort to analyze their behavior in deep, hot boreholes, it was thought desirable in this investigation to study suspension of pure, monoionic montmorillonite and attapulgite clays before and after hydrothermal heating. Chemically, clays are hydrated aluminum silicates, such as montmorillonites and attapulgites. The structure of montmorillonite is composed of units made up of two silica tetrahedral sheets with a central alumina octahedral sheet. Attapulgite consists of double silica chains running parallel to the c-axis with the chains linked together at their longitudinal edges. The apexes of the tetrahedrons in successive chains point in opposite direc-tions.4,7,10 One of the most important characteristics of clay minerals is their property of sorbing certain anions and cations and retaining them in exchangeable state. The property of exchange capacity is measured in terms of milliequivalents (meq) per 100 gm of clay. The cation-exchange capacities are taken at pH 7 from titration curves which show the relation between pH and milliequivalents of added NaoH.9,11 The following are causes of cation exchange: (1) broken bonds around the edges of the silica alumina units would give rise to unsatisfied charges, which would be balanced by adsorbed cations; (2) substitution within the lattice structure of trivalent aluminum for quadrivalent silicon in the tetrahedral sheet and of ions of lower valence; and (3) the hydrogen of exposed hydroxyls may be replaced by a cation which would be exchangeable. The primary techniques employed to study the effects of hydrothermal treatment were titration curves and X-ray diffraction analyses. Titration curves were useful in characterizing the nature of the clay mineral surface and the stability of the various samples to attack by acid or alkali before and after hydrothermal treatment. X-ray diffraction analyses were used to determine changes in the nature and crystal structure of the minerals present and to determine the cations in exchange positions on the montmorillonite samples. EXPERMENTAL PROCEDURE PREPARATION OF SAMPLES Selected samples of Wyoming bentonite and attapulgite clay from mines along the Georgia and Florida border were purified by centrifugation to remove all particles larger than 2 microns. In preparing magnesium attapulgite, first a 1.00-normal solution of magnesium acetate was prepared. The pH of the solution was adjusted to 7.0 by adding magnesium oxide. One hundred grams of attapulgite clay were mixed in a blender for 30 minutes with the magnesium acetate solution. The attapulgite solution was left overnight, and then filtered in a Buchner funnel. The suspension was leached with further quantity of magnesium acetate solution and then washed with distilled water for 20 days, after which the filtrate was free from magnesium ions. The concentration of magnesium attapulgite in the clay slurry was adjusted and determined quantita-

By Wayne L. Worrell

Thermodynamic data for the stable carbides and oxides of chromium, molybdenum, and tungsten have been critically eualuuted and are used to determine the stable condensed phases at 1 atm total pressure in each metal-carbon-oxygen system. Pourbaix-Ellingham diagrams have been constructed and are used to estimate the minimum temperatures neces.see)*y to oblain these metals by corbolhermic redlrction under reduced pressures. USING recently reported data for the carbides and oxides of chromium, molybdenum, and tungsten, the stable phases in the Group VI metal-carbon-oxygen systems at specific temperatures and pressures can be determined. The thermodynamically stable regions in each metal-carbon-oxygen system can be completely pictured using a Pourbaix-Ellingham diagram in which the two coordinates are oxygen potential (-RT In Po2) and temperature. Worrell and chipman&apos; have described and constructed such diagrams for the metal-carbon-oxygen systems of vanadium, columbium, and tantalum. Downing2 has used similar diagrams to provide a description of the reactions in submerged arc furnaces. One of the examples which he cited was the Cr-C-O system; however, an improved diagram is presented in this paper using recently obtained data for the carbides of chromium. Although some of the oxides and carbides of the Group VI metals exist over slight ranges of composition,3 all solid phases are represented as stoi-chiometric compounds in the subsequent diagrams. Any error introduced into the calculated equilibria by neglecting the compositional variation of these phases is less than the uncertainty in the thermodynamic data. In a Pourbaix-Ellingham diagram for a three-component system, the phase rule requires that at equilibrium four condensed phases specify a point, three condensed phases determine a line, and two condensed phases define an area. To obtain a well-defined area at high oxygen potentials and teniperatures, the initial carbon to oxygen mole ratio is specified to be more than sufficient to reduce the oxide but not enough to convert all the metal to carbide. Thus, the lowest regions in the subse- quent Pourbaix-Ellingham diagrams represent areas in which the metal and its carbide are the stable condensed phases. Cr-C-O SYSTEM A) Thermodynamic Data. To simplify the subsequent calculations, all thermodynamic data are presented in the form of linear Gibbs-energy-of-formation* equations which were derived for the temperature range 505° to 2000°K and are reasonably accurate over the entire temperature ranges of the diagrams. The thermodynamic data for the compounds occurring in the Cr-C-O system are summarized in Table I. Several volatile chromium oxides have been reported,4 but their influence is negligible in the Cr-C-O system. Thermodynamic data for CO and CO2 were taken from the compilation of Coughlin.5 The only solid chromium oxide sufficiently stable to exist below 1800°K in the Cr-C-O system is Cr2Os.11 The Gibbs-energy-of-formation equation for Cr2O3 was obtained by revising the data tabulated in Coughlin,5 using the heat of formation of ah&apos; and the high-temperature thermal data for chromium.7 In 1953 Richardson8 re-evaluated the original work of Kelley and coworkers,9 who measured the thermodynamic properties of the chromium carbides. However, the high-temperature equilibrium results of Boericke9 yield values for ?G°f of Cr2O3 which are 2 kcal more positive than the accepted data. His results for the chromium carbides are based on the same questionable equilibria. Gleiser10 recently determined the Gibbs energy of formation of Cr3C2 and the equation in Table I was calculated from her data and the thermal data of Kelley and coworkers.9 Thermodynamic data for the Cr7C3 and Cr23C6 carbides were obtained by combining the Cr3C2 data

Jan 1, 1965

By M. Cross

INTRODUCTION The development of refining processes involving gas injection into liquid metals has seen the evolution of a variety of designs [I]. During the last few years or so the top, bottom and side blown processes have, essentially, been combined to produce a series of mixed gas blowing schemes (eg. the LBE [2,3] and LD-OB processes [4,5]). Although, a variety of mixed gas blowing processes have been carefully engineered and recorded operational experience now attests to their effectiveness, optimising their design and operational performance is another matter. Detailed process analysis of such operations is impeded by the complex nature of the physical and chemical processes involved: . fluid flow behaviour of the liquid metal and injected gases, . chemical reactions in the vessel, . mass and heat transfer within the fluid phases in the vessel, . influence of the geometry (eg. location of the tuyeres, etc.), . monitoring the operational performance. By and large, design and operating parameters evolve from a combination of laboratory scale experiments, pilot plant work and, of course, practical engineering experience. Conventionally, the gas-liquid system used to model the industrial process is air-water, for a number of valid reasons including – . visual demonstration plus geometric and (to some extent) dynamic similarity, . the fact that the system demonstrates, at least, qualitatively many of the complex features of gas injection - interaction between the gas and the liquid - the nature of the gas flow behaviour - qualitative indications of mass transfer rates Important experimental contributions to the understanding of gas liquid systems include the work of Brimacombe [6], Wraith [7], T.Robertson [8], D.Robertson [9], Guthrie [10] and Szekely [11] (these references provide starting points to identify the major published work in this area). Although air-water systems have been instrumental in providing a great deal of qualitative insight into the fluid flow behaviour of the processes, reliable quantitative predictions are more problematic. A number of workers (notably Brimacombe and colleagues [6, 12])have demonstrated that different gas-liquid systems can lead to wide variations in behaviour. Other limitations of physical modelling include the assumptions of isothermal flow, scale- up of turbulent flow behaviour and the inability to directly include the effects of chemical reactions. The above limitations in physical modelling, therefore, yield a natural role for mathematical modelling provided that it can be made comprehensive enough whilst remaining cost-effective. The major physical processes involved in gas injection systems may be summarised as . behaviour of the gas envelope structure in moving through the liquid metal . the global movement of the liquid metal . momentum interaction between the liquid metal and the injected gas . heat transfer and mass exchange via the chemical react ions. The behaviour of gas envelopes forming at a tuyere and then moving through the liquid has been investigated thoroughly by a large number of workers. Such models are generally based upon simple physic- al concepts derived from visual observation. An early model by Davidson and Schuler (DS), based upon a force balance between the liquid inertia and a combination of the bouyancy and momentum forces on the envelope, has proved to be a useful tool in predicting bubble volumes at detachment [13]. The DS model has been successfully modified in a number of ways [14-17] to cope with heat transfer and reactive transfer rates, etc. Furthermore, the basic ideas summarised above have also been confirmed on industrial converters by both Guthrie [I8] and Brimacombe [19]. Models have also been developed to predict the behaviour of the "chemical" component of some gas-liquid systems. However, one thing is clear; all gas-liquid reactive systems are governed by the degree to which effective mixing occurs. As such, the processes are dominated by the behaviour of the heat and mass exchange of the two phase fluid flow system. It is the pattern of movement of the turbulent liquid metal plus its interaction with the injected gas which must be predicted reliably, if sufficiently comprehensive comp-

Jan 1, 1984

By J. G. Eenink, R. A. Cunningham

One phase has been completed of a laboratory invesrigation of formations with relatively high permeability under conditions of overburden, formation and mud coltrmn pressures. The following statements are based on these tests within the limitations described. Drilling rate decreased when mud column pressure was greater than formation pressure. The decrease was primarily due to a layer of cuttings and mud particles held to the hole bottom by the difference in pressure. Much of the force of the bit was wasted in this layer and was unavailable for penetrating virgin formation. Adequate jet velocities helped clean the filter cake and chips away and resulted in increased drilling rate — the higher the jet velocity the faster the drilling rate. Drilling rate increased slightly when formation pres sure was greater than mud column pressure. The i11(reuse resulted from cleaning the hole bottom as formarion fluid flowed into the borehole. Overburden pressure had practically no effect on drilling rate. INTRODUCTION Some formations are difficult to drill in a wellbore hole but are easily drilled under atmospheric conditions in the laboratory. The reasons for this apparent difference in drillability are not fully understood. An increase in formation compressive strength causes a reduction in drilling rate. It has been shown by Griggs1, Handin2 and others3-1 that formation strength increases under conditions of high hydrostatic stress. Based on these results, a decrease in drilling rate would be expected with increased drilling fluid pressure in formations of low permeability having low formation fluid pressure. Such a decrease has been shown by Murray. et a1 and Eckel" This effect holds even when using nitrogcn under high pressure as the drilling fluid (unrc-ported tests from the Hughes laboratory). Poor bottom-hole cleaning may account for decreased drilling rates. Even though it may not be easily explained, evidence of poor cleaning can easily bc seen. One evidence is the increased drilling rate that accompanies decreased balling obtained with addition of oil in some muds used to drill shales.4,7,3 Many variables can affect drilling rates. These may include mud composition and properties, formation properties, bit types and operating conditions, formation fluid, mud column and overburden pressures. These must be studied a few at a time under carefully controlled conditions which simulate actual conditions as closely as possible. The purpose of this work was to investigate effect of overburden, formation and mud column pressures on drilling rates in permeable formations. These tests were suggested by conditions shown in Fig. I. Represented is a borehole drilled into a formation with relatively high permeability. In the vicinity of the hole bottom there exists a mud column pressure, P a formation fluid pressure, P and an overburden pressure, The mud has low water-loss properties and sufficient weight to more than overbalance the formation pressure. Laboratory investigation included a range with formation pressure greater than mud column pressure. NOMENCLATURE Pm mud column pressure = pressure of drilling fluid or mud at hole bottom Pf, formation pressure = pressure of fluid in intersticies of formation

By Raymond D. Sloan

During 1944, substantial gains were recorded in practically every phase of the petroleum industry in Oklahoma. With the spotlight of activity focused on other states during the more recent years, the production momentum that was achieved in Oklahoma through the discovery and development of such outstanding pools as Oklahoma City, Fitts. and major pools in the Seminole area during the late 1920&apos;s and early 1930&apos;s was practically spent by the end of 1943. However, with much of the industry&apos;s exploratory efforts greatly reduced in Illinois, Indiana and Kentucky, the strong demand for Mid-Continent crude and the optimism created by some discoveries in the state toward the end of 1943, Oklahoma has again become the recipient of much activity and interest, with its 1944 drilling activity, production and refinery runs to stills registering substantial gains over the preceding year. Development and Exploration Since the discovery of the West Edmond pool in April of 1943, with the substquent drilling and development boom, the scene of exploratory activity has shifted from the central. south and east portions of the state to the general trend of the . Granite Ridge, where leasillg and wild-catting were greatly activated within a wide radius of this high subsurface feature. Leasing activity, however, jumped the confines of the Granite Ridge trend and extended into all of the northwestern counties in the state. Seismograph activity at once began following and keeping pace with this shift in the leasing and exploration locale, by the concentration of a large number of seismograph crews in the area. At one time during the year,35 crews were concentrated in Oklahoma, the majority being west of the Granite Ridge and extending into the northwestern counties. To illustrate the increase in geophysical work done in the state during the year, records indicate that 377 crew months were worked in 1944, which is a 21 per cent increase over the 311.5 crew months activated in 1943. With the discovery in 1943 of several encouraging prospects and the finding of West Edmond as an incentive, drilling operations increased some 54.3 per cent, for a total of 1986 completions during the year. This is compared with a figure of 1287 for wells drilled in 1943, an increase of 699 wells. The industry, with a total of 1064 oil producers, enjoyed a success percentage of 53.6 per cent. Gas wells completed were 223, while 699 dry holes were drilled. In 1943, completions included 587 oil wells, 127 gas wells, and 573 dry holes. The number of oil producers completed during the Year was almost double the numher of wells in 1943. Of the total drilling wells completed in the state, 350, or 17.6 per cent, were wildcat operations—an increase of 5.1 per cent over the 194.3 total of 333. The wildcat success percentage was 13.4 per cent, 47 successful wells being completed for an average initial potential of 160 bbl. In the state, 23 gas wells were completed with an average initial production of

Jan 1, 1945

By H. Siebel, N. A. D. Parlee

MIZIKAR, Grace, and parlee&apos; studied rates of absorption of oxygen in stagnant liquid silver in a new type of open-ended gas-liquid metal diffusion cell, where the capillary diameters used were 1.5, 2, and 3 mm. They reported that the rates appeared to be controlled by the rate of diffusion of oxygen atoms in liquid silver and they used their data to arrive at the following equation as a reasonably good representation of the diffusivity (D) of oxygen in liquid silver: D = (14.7)10-4 exp (-7100 ± 3200)/RT cm2/sec Since this type of diffusion cell offers a promising method of measuring diffusivities of gas atoms in liquid metals, the senior author of this note and his students are carrying on a critical examination of the method, and are attempting to improve and extend the technique. The research reported briefly here was designed to examine more carefully the earlier conclusion1 that the absorption rate is diffusion-controlled, and to find out whether the diameter of the cell could be greatly increased and thus extend the method to systems where gas solubilities and diffusion rates are much smaller than in the case of oxygen in silver. For the purpose of these experiments, the open-ended capillary cell used by Mizikar, Grace, and Parlee was replaced by a small cylindrical crucible (1 to 1.7 cm in diam) containing pure silver (to a depth of 1.3 to 3 cm) resting in the bottom of a closed-ended vertical quartz tube, with an evacuated quartz filler tube just above the crucible to reduce dead volume. The method of measuring the rate of absorption was otherwise very similar to that of Mizikar, Grace, and Parlee and all the same precautions were taken. A much smaller and shorter tube furnace was used for heating, but simi- lar precautions about temperature gradients were observed. This technique was nothing more than the application of a Sieverts apparatus to the measurement of rates of absorption in small crucibles of stagnant liquid silver. In each run the volume of gas absorbed (V) was plotted against the square root of time (vt) and apparent diffusivities were calculated from the slope of the line, as described by Mizikar, Grace, and parlee.&apos; Care was taken to use only the data obtained during the early part of the total absorption period while the cell functioned as one of infinite depth. The results at 1000°C are summarized in Table I and Fig. 1. The following conclusions may be drawn: 1) Diffusion cells up to 1.7 cm in diam yield, within experimental error, the same values for the diffusivity of oxygen (D) as the capillary cells of 1.5, 2, and 3 mm in diam employed by Mizikar, Grace, and Parlee. Cell diameters as large as 1.7 cm are probably not desirable for the most precise work. Intermediate sizes, perhaps 0.5 cm, will be most practical. 2) All absorptions measured were consistent with a model where the rate is controlled by the diffusion of gas atoms in the metal. The fact that runs at absorption pressures ranging from 0.1 to 1 atm (at any one temperature) give

Jan 1, 1965

By F. Bonaccorso, G. Venturello, C. Antonione

A study of the effect of small amounts of interstitial impurities on recovery and re crystallization in high-purity iron (99.995 pct) has been undertaken. This paper gives results on the effect of carbon, introduced in small-dosed amounts @om 0.0005 to 0.0086 wt pct) by heating iron specimens of high purity m a static atmosphere of CO + ,. The materials prepared in this way, cold-worked 80 pct and subjected to a series of isochronal and isothermal annealings, were submitted to examinations by X-rays, micrographs, and hardness tests. It was observed that effect of carbon is remarkable in the sense of blocking recovery of mechanical properties up to the temperature at which re crystallization begins. On the contrary, carbon has a negligible effect on primary re crystallization temperature, when compared with the known effect of substitutional impurities. This is in agreement with the high mobility of interstitials. In effect, only a slight decrease, -20 pct, of the grain-boundary motion rate was noted, due to the interaction between the grain boundary and the carbon atoms. On the other hand, in the samples in which the carbon content is above the solubility limit at the temperature at which re crystallization occurs, a slight increase of nucleation frequency is noted due to the presence of precipitated carbides. ThE effect of purity on recovery and recrystalli-zation phenomena has been known for a long time; however, recently, new attention has been given to the problem since new methods for obtaining metals with extremely low impurity contents have become available. Most of this research work essentially concerns the effects of impurities which cause precipitates or give rise to substitutional solid solutions. The works of Bolling and winegard1 and Aust and utter&apos; on lead, and Vandermeer and orddon&apos; on aluminum, 01sen4 on nickel, and Abrahamson and Blakeney on iron should in particular be referred to. On the effect of this type of impurities a quantitative theory has been formulated by Detert and Lucke6 and has lately been discussed critically by Cahn7 and Gordon and vandermeer.&apos; Interstitial solid solutions in iron have not as yet been studied. As the study of the effect of interstitials is of great interest both from a theoretical and practical standpoint, it was deemed useful to examine the effects of carbon, nitrogen, and boron on recovery and primary recrystallization of iron. There already is some work by Chaudron et al.,1-" on the effect of interstitials on iron; their work, however, mainly concerns secondary recrystallization. The present paper refers in particular to the effect of carbon. EXPERIMENTAL PART Preparation of Materials. The pure iron used for this research was obtained from FeCIS recrystal-lized and purified by extraction with isopropylic ether. From the ferric chloride purified in this manner, hydroxide was precipitated with a solution of very pure ammonia, and, by calcination in pure sintered alumina crucibles, oxide was obtained. Reduction of the oxide to iron sponge was performed in sintered alumina tubes with very pure hydrogen at 650°C; at the end of the operation temperature was increased to 900°C. Specimens for the experiments were obtained by sintering the sponge at 1480°C in pure Hz after a pre sinter ing treatment at 900° C. It is important to note that the treatment at 1480°C in Ha produces a further purification from more volatile elements such as zinc, cadmium, arsenic, lead, and tin. Details on the preparation and characteristics of this type of very pure iron are given in a previous work.&apos;&apos; Only the complete analysis performed by neutron-activation methodsz3 is given here, Table I. Some of the specimens prepared in this way were carburized in the 0 region with very low amounts of carbon by treating them at 700°C in a static atmosphere of Ha containing a definite amount of CO. The set-up used is described in Fig. 1. A gas-tight quartz tube containing the specimen to be carburized and an internal friction control specimen, after being evacuated, was filled with Hz

Jan 1, 1963

By W. O. Philbrook, H. W. Hosking, N. B. Melcher

A simple and inexpensive mercury-tracer method has been develohed to study rates and Patterns of gas flow in blast furnaces. A Pulse of mercury is injected into the hot blast, and its arrival at the sampling Point is recorded by an electronic circuit that measures resonant absorbtion of ultraviolet light by mercury vapor in the gas aspirated .through the detecting. cell. A rapidly developing interest in blast-furnace technology in recent years has created a strong demand for precise knowledge of the rates and patterns of gas movement through the stock column. Classical fluid-flow theory cannot supply the desired information because the complexity of the blast-furnace situation precludes the accurate determination of many of the parameters required for solution of the equations. In fact some of the terms employed in fluid-flow treatments become almost meaningless under blast-furnace conditions. The void fraction of the stock column, for example, varies widely from place to place and is affected by the segregation of charge materials, attrition of lumps during descent, and channeling. Gas density varies with changes in local pressure, temperature, and composition, while viscosity varies with temperature and composition. Even if mean values for these functions were acceptable, they would be either impossible or extremely difficult to determine with a high degree of accuracy since the value of the mean for all gas properties would depend upon the flow distribution itself. Attempts to study gas movement by the pitot-tube determination of kinetic pressures have met with little acceptance. Reproducibility of results is extremely difficult to achieve in a situation where alignment of the tube orifice with respect to the flow direction in tortuous channels cannot be guaranteed. Because of these difficulties it would appear that progress can be achieved only with a direct method of measuring gas velocities and determining the shape of flow paths through the stock. The establishment of such a method is described in this paper. SELECTION OF METHOD Direct measurement suggests the use of a tracer by means of which the progress of a "parcel" of gas can be followed through the furnace. Several requirements must be met before any selected material can be considered a suitable tracer. First of all, it must be a substance not normally present in blastfurnace gas or dust; otherwise its presence in a gas sample will have little significance. It should be volatile at temperatures normally existing within the furnace, to ensure that it will travel with the gas, as a gas. It should possess special properties enabling its presence in the gas to be detected readily, and it should not undergo exchange reactions with any material normally present within the furnace, including refractories and coke as well as all components of the metal and slag. Radioactive materials naturally are considered because their radioactivity renders them susceptible to detection. However, of the many radioactive materials available, almost all are eliminated from final consideration by their unsuitability in terms of other properties. Radioactive silicon, carbon, calcium, phosphorus, and many others have been considered and discarded, either because they are not volatile over the whole range of significant temperatures, or because they could be expected to undergo exchange reactions with burden or refractory materials. Others which might be suitable in terms of these properties are eliminated because of a low activity, short half life, or the necessity for fragile counting apparatus. Of the many elements available in radioactive form, perhaps radon is the only one suitable for blast-furnace gas-flow studies. This has been used in several investigations1-4 with considerable success. The use of radiotracers has certain inherent disadvantages. One is the cost. An appreciable sum must be allowed for materials each time a test is made. Also, special methods of handling must be employed, and very often the progress of an investigation must be governed by the availability of isotopes and the delivery arrangements that can be made for them. An advantage is immediately gained if the tracer material can be a substance which can

Jan 1, 1960

By W. T. Weller

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

Jan 1, 1967

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

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

Jan 1, 1963

By A. W. Sommer, H. H. Kellogg

It is shown that SO3-O2 mixtures react with sphalerite at an appreciable rate ill the temperature range of 361° to 527°C to fornz ZnSO4. The rate of reaction follows a parabolle lax. Oxygen, or O2-SO2 mixtures, have a negligible effect on sphalcrite in the same temperatuee range. PRACTICAL roasting of sphalerite is usually performed at 800°C or higher. The calcine is composed of zinc oxide unless the roaster atmosphere contains relatively large amounts of SO2 and the temperature is held to 800oC or lower, in which case zinc sul-fate may also form. Ong, Wadsworth, and assell&apos; have made a careful kinetic study of sphalerite oxidation under these "normal conditions" and have concluded that the rate of the reaction is controlled by the decomposition of an activated complex adsorbed on the sphalerite surface. They found the rate of oxidation to be small at 700°C. Extrapolation of their data to 400°C would indicate an almost negligible rate at this temperature. In our work sphalerite was reacted with air, mixtures of SO2, O2, and N2, and mixtures of SO3, O2, and N2 in the temperature range 361o to 527OC. Negligible rates of oxidation were found, except for the gas mixtures containing SOs. With this latter gas, the oxidation of finely divided sphalerite was fairly rapid and zinc sulfate was the product. Based on the limited evidence available, the rate of sphalerite oxidation by SO3 is postulated to be controlled by gaseous diffusion through pores or cracks in the zinc-sulfate coating. This evidence for the direct reaction of sphalerite and SO, at low temperature may prove of importance to the understanding of zinc-sulfate formation in the dust-collecting equipment usually associated with zinc roasters. The roaster gas carries off fine dust, much of which may be unreacted sphalerite. The temperature in dust-collecting equipment drops from the roaster temperature (about 800°C) through the range of temperature we studied, to ambient temperature. Such dusts are known to contain far larger amounts of zinc sulfate than the primary calcine. It has been assumed in the past that this sulfate is formed by reaction of ZnO dust with SO3 (or SO, plus 0,) in the partly cooled gases. Our work shows that an alternative possibility exists—the direct reaction of SO3 with ZnS dust. EXPERIMENTAL Reaction rates for sphalerite at 350° to 600°C are relatively small so that it was necessary to use finely divided material with a relatively large surface area in order to obtain measureable amounts of reaction. In the preliminary experiments, pure mineral sphalerite, ground to pass a 325-mesh sieve, and chemically precipitated ZnS (in the form of an impalpable powder) were used. In the quantitative rate experiments the ground mineral sphalerite was processed in an Infrasizer to obtain a product with particle size between 9 and 18 µ. The analysis of this material was 66.7 pct Zn, 0.09 pct Fe, and 32.7 pct S (theoretical sphalerite is 67.1 pct Zn, 32.9 pct S). The sized powder was carefully mixed and divided into l-g samples. The apparatus used for the quantitative rate measurements is shown in Fig. 1. For the preliminary experiments the following change was made: The external catalyst furnace (B, C)* was not used. *Letters refer to Fig._______1 .___- Rather, when catalysis of the SO2 + O2 reaction was desired, the glass beads in basket M were replaced by vanadium-oxide pellets. Procedure—A l-g sample of ZnS was placed on a shallow stainless-steel tray, I, and spread evenly to a depth of about 1mm. The tray was placed in the reaction assembly and the assembly inserted in the cold furnace. The furnace tube was flushed with dry argon and then heated to the desired reaction temperature. The temperature was controlled to ±0.5oC. Gas mixtures (O2 + N2 or O2 + SO2) were prepared from tank gases by means of a mixing device identical to that employed by Darken and Gurry.2 The total flow rate of gas was 205 ml per min to point A of Fig. 1. The accuracy of the individual and

Jan 1, 1960