Search Documents

By Harry M. Mikami, A. Gene Sidler

Chemistry of the evolution of materials in contact with copper converter tuyeres is delineated by means of analyses of periodic punch rod samples taken during a converter cycle. Lining samples from known peripheral and end locations from three converters were investigated chemically and petrographically. One set of eight specimens representative of all lining parts was removed for study before completion of the campaign. Purely surface removal by slag action is minor. Slag filling of interconnecting cracks extending to the hot face results in a "stoping'' action which shapes the hot face contour and is a significant factor. Major wear mechanism is believed to be copper metal slabbbing caused by metal first penetrating joins and invading fractures or impregnating brick structure laterally from the joints. Importance of lining design and joint fit are emphasized. New laboratory tests are indicated. In the copper converter, molten matte, which is essentially a solution of iron and copper sulfides, is blown with air to produce blister copper. During the first stage of the process the FeS component is oxidized to SOz and FeO, the latter combining with a siliceous flux to form a predominantly ferrous silicate slag. Successive increments of matte and flux are charged, blown, and skimmed of slag until sufficient CuzS ("white metal") has accumulated. This white metal is then blown to copper metal in the final or blister-forming stage by the oxidation of sulfur to SOz and the reduction of combined copper to metal. The reactions of both stages are exothermic so provide the required heat. Temperature is controlled by cold scrap copper, flue dusts, smelter reverts, and so forth. Considerable vari- ation in practices occur depending on the type of converter, grades of the matte and flu, and other factors. Lathe and Hodnett's' survey of world converter practices gives much information on these variations and their effects on the refractories. Principles of converting chemistry are presented by Newton and wilson,' Ruddle,' and schuhmann. Most copper converting occurs within the temperature range 2100" to 2370°F although individual smelters report1 ranges from as low as 1700" to 2240°F to a high of 2200" to 2400°F. Temperature intervals reported from beginning to end of blow varied from 45" to 540°F. Average temperature range was 2030" to 2260°F and the average interval 250°F During the matte-blowing stage varying amounts of magnetite (FeO.Fe,O,) are formed in the slag as suspended crystals and also the liquid portion of the slag will contain, in addition to ferrous silicate, some Fe20s which adds to the potential magnetite formation if the temperature falls. This is a component of the so-called "dissolved magnetite" often referred to in the literature. Although the slag may contain minor amounts of CaO, A1209, MgO, and alkalies, the principal components are Fe0-Fez03-SiO,. The phase equilibria relationships of this system have been elucidated by Muan~ and Schuhmann, et a1.' These data and estimates of converting atmospheres by schuhmann4 are useful tools in studying slag behavior. It would appear that in many cases reported magnetite contents of slags are apparently higher than actually is present at operating temperature. This is probably because much of it is formed during slag solidification. The quaternary systems which have Mg-0,1,0,' and addition to iron oxides and silica are also of interest in considering the reactions of the slags with basic refractories. The formation of magnetite has an important interrelationship throughout the pyrometallurgy of copper (Mossman,' Gronningsater and Drummond," Ruddle,' Elwood and Henderson") but its presence in the converter particularly relates to operating efficiency and lining life in certain instances. Archi-

Jan 1, 1963

By I. H. Silberberg, R. L. Reed, O. K. Kimbler

interfacial films have frequently been observed at interfaces between certain crude oils and water. Several investigators have postulated that the presence of these films should influence the efficiency of oil recovery in water drive or waterflood operations. They may also influence the stability of emulsions which are sometimes a problem in petroleum production, and may be a factor in the fonation of paraffin deposits in oil well tubing and flow lines. This paper presents a technique with which a modified Langmuir film balance may be used to study the compressibility and collapse pressure of these natural interfacial films. Experimental data are presented for several nude oil-water systems. Data developed are used to infer the phase state of the film as a function of such variables as rate of reduction of interfacial area, ionic composition of the subtrate and pH of the subtrate. A film of known physical characteristics is shown to have a sinnificant effect on oil recovery from an unconsolidated sand pack. Possible applications of these results to petroleum production are discussed. INTRODUCTION The use of water to displace petroleum from reservoir rocks is of major importance both as a primary and a secondary recovery process. As water invades the rock, oil is completely displaced from some pores and left as a discontinuous phase in other pores. The manner in which water moves from pore to pore is strongly influenced by capillary forces. In view of the complexity of reservoir fluid systems, there can be little doubt that complicated interactions take place at both the liquid-solid and oil-water interfaces. One of the more interesting, and least understood, of the phenomena which take place at the oil-water interface is the formation of interfacial films. These films are believed to result from the adsorption of high molecular weight polar molecules at the interface.l.2 Presence of such molecules may cause a striking alteration in interfacial tension. When the oil-water interfacial area of certain crudes is rapidly reduced, a thin region (film) about the interface assumes the appearance of a solid membrane, and striations, wrinkles and gross distortions may occur. If such a film is solid, it should greatly alter the interfacial tension normally assumed to exist between the oil and water phases. If the membrane is continuous, a solid phase would separate the oil and water. Interfacial films between crude oil and water were observed in 1949 by Bartell and Niederhausers who commented upon the apparent rigidity of the films and their possible importance in the petroleum industry. Morrell and Egloff4 had earlier attributed the extreme stability of emulsions of sea water in fuel oil to very stable asphaltic films. Numerous investigators have observed rigid films in the course of crude oil-water interfacial tension determinations by the pendent drop method. Several investigators5,6,2 have separated inter-facially active materials from crudes and attempted to characterize them chemically. Reisberg and Doscher,2 using Ventura crude, showed the interfacial tension against water (as measured by the pendent drop method) to be affected by aging, contraction and expansion of the interface, and the pH of the water. These investigators attributed the adhesion of oil to a water-wetted surface and the distortions of flow paths in glass capillaries to the presence of rigid films. Dodd7 has studied the interfacial viscosity of adsorbed films and found them to be non-Newtonian in behavior. Craighead and Harvey8 reported a series of displacements in tubes packed with 60 mesh glass beads. They interpreted the results as indicating an effect of stearic acid films on waterflood recovery and imply that natural films may produce similar results.

Jan 1, 1967

By F. C. Bond

Since the Third Theory of Comminution was presented eight years ago (I) it has found increasing use in crushing and grinding problems. The practical utility of its wok index equation is quite generally acknowledged (2). However, its theoretical basis has been questioned in at least three technical articles (3) (4) ('). The purpose of this paper is to present experimental proof that it is scientifically correct. Particles under compressive stress are strained and deformed. They absorb strain energy, and when this locally exceeds the breaking strength, a crack tip forms. The surrounding strain energy flows to the crack tip, which rapidly extends and splits the rock, releasing the strain energy as heat. The initial energy flow causes additional crack tips in highly strained areas. If the compression is rapidly applied by impact, crack tips may form before the strain energy has reached equilibrium in the particle, thus decreasing the total work input required for breakage. The energy necessary to break is essentially the energy necessary to produce crack tips, since the energy necessary to extend the cracks to breakage is already present as strain energy in the deformed particles. After breakage nearly all of this energy appears as heat. The crack length cannot be measured directly. However, in particles of regular and similar shape the crack tip length is considered as equal to the crack depth, or crack extension necessary to break, so that the crack length equals the square root of one-half of the surface area. The Third Theory states that the useful work done in crushing and grinding is directly proportional to the total length of the new cracks formed. It can be confirmed by showing that a constant work input produces a constant length of new cracks when reducing the same material to different product sizes. This is done in the present paper on a wide variety of material. The constant work input was supplied by one revolution of the 12" x 12" laboratory ball mill used in making grindability tests by the Allis-Chalmers method (12) (13) The new crack lengths produced per mill revolution were measured from all available grindability test results at 28, 35, 48, 65, and 100 mesh on fifteen different ores, and were found to remain substantially constant for each ore at all mesh sizes. A new technique is used for the measurement of crack lengths. Size distribution plots of the mill feed and product are made by the Third Theory method (9) and the crack lengths are obtained from these plots by the procedure described in the present paper. The energy input required to produce a unit length is of fundamental importance in the size reduction of brittle solids. The crack length Cr is expressed in centimeters per cubic centimeter of solid material. It bears a definite relationship to the external surface area of the crushed or ground solid. A uniform particle shape must be assumed before the surface area and crack length can be evaluated. In this paper it is assumed that the relationship between the surface area and the particle volume of a particle d microns in diameter is the same as that of cube d microns on a side. The external surface areas of particles with a cubical breakage probably agree approximately with this rule, and correction factors can be applied when physical measurements of the surface areas are available for comparison. However, the assumption of equivalent cubes has been found satisfactory for most calculation purposes. Assuming equivalent cubes, one cubic centimeter of particles d microns in diameter will have a crack length Cr of v30.000/d centimeters, and a surface area of 60,000/d square centimeters. The specific crack length is thus equal to the square root of one-half the specific surface area. Where Sa is the surface area in square centimeters per gram and Sg is the specific gravity of the ground solid, then Cr = vSg . Sa/2 = 173.2/ vd (1)

Jan 1, 1961

By J. A. Herbst, D. W. Fuerstenau

This paper examines the zero order production phenomenon in the context of the size discretized batch grinding model. A restrictive interrelationship between the selection and breakage parameters of the model, which is mathematically sufficient to ensure zero order behavior, is delineated. Based on this interrelationship, a scheme is developed for predicting the values of the selection and breakage parameters for all size fractions from a minimum of experimental data. To test this scheme, dolomite was ground in a laboratory batch ball mill. Successful simulation of the comminution behavior of this system was achieved using the batch grinding model with the parameter values obtained from the proposed scheme. Historically, a preponderance of comminution research has centered on attempts to relate breakage energy to the performance of comminution machines. Concomitantly, grinding data were interpreted almost exclusively in terms of inherently empirical energy-size reduction relationships1-3 or "laws of comminution"4- 6 which were based on highly oversimplified descriptions of the fracture process. In some instances, these relationships provide a crude basis for the correlation of experimental data, but, invariably, this approach is inadequate for meaningful process simulation. The control and optimum design of comminution circuits require a mathematical model capable of depicting the size reduction behavior of every size fraction for grinding conditions of technological importance. Energy-size relationships do not provide this detailed information. For example, an energy-size analysis for batch ball-milled dolomite has recently been completed in the authors' laboratories.7 This analysis provided a satisfactory correlation between the hypothetical size modulus, which characterizes only the finer portion of the size distribution, and the input energy. However, this relation did not constitute an adequate description of the complete system. In general, one of two distinct viewpoints toward simulation has been adopted in recent comminution research. The first approach focuses on the fracture of single mineral specimens, with the essential aim of representing the over-all process in terms of the breakage characteristics of individual particles and the characteristics of the stress field which the particles experience within the particular size reduction device. As illustrated by Harris's review,' the actual incorporation of single-specimen fracture information into a description of the behavior of a multiparticle comminution system has seldom been attempted. The only model to accomplish this incorporation was the one developed by Schönert.8 The Schonert model treats single-passage grinding machines in terms of a distribution of effective loads acting on a particle, a distribution of required particle breakage energies, and the breakage product distribution. However, the use of single-particle fracture behavior to derive models of the size reduction process is at present limited to machines in which the particle residence time is approximately equal to the time required to apply stress. As a result, Schönert's analysis cannot be used to advantage to simulate the complex environment which prevails within a tumbling mill. For multiple-passage systems a second approach to the description of comminution processes has been more fruitful. This approach entails the formulation of a mathematical model which is phenomenological in nature in that it lumps together the entire spectrum of stress-application events which prevail in a system under a given set of operating conditions. The appropriately defined average of these individual events is then considered to characterize the over-all breakage properties of the device. Thus, to analyze the performance of a tumbling mill, the manner in which the particles of a particular size (or size fraction) are stressed need not be distinguished. Instead, a single parameter is assumed to represent the resistance of that size to fracture, given the average grinding environment which exists in the mill. The isolation of such a parameter and a related set of quantities, which constitute the breakage product size distribution for the average event in this size fraction, allows the formulation of physically meaningful descriptive equations capable of yielding precise and detailed information for simulation.

Jan 1, 1969

By L. N. Cosner, C. F. Austin, J. K. Pringle

Shock wave attenuation studies with three igneous rocks are reported: a spessartite elastic to low intensity stress waves, a diorite which is anelastic, and a porous scoria. The level of shock passing through a given thickness of rock has been determined by donor-receptor methods using high speed streak cameras. Test results include shock propagation velocities adjacent to the explosive, the dependence of test results upon test geometry, and the variation in shock wave attenuation in these three common igneous rock fabrics. The utilization of explosives in rock mediums involves a series of interrelated phenomena. These are the detonation process within the explosive, the transfer of energy from the explosion products to the surrounding rock medium in the form of sound or shock waves, and the transmission of the resulting waves through the rock medium to create damage in the rock mass. The subject of this paper is the investigation of the close-in phenomena that occur in the intensely shocked rock immediately adjacent to a detonating explosive charge. The specific aspects studied include propagation velocities, free surface velocities, pulse attenuations and the determination of the equations of state for the test rocks employed. In addition, a series of tests was conducted to indicate the dependence of the test results on the charge diameter utilized. The equations of state for rock materials have been investigated by a number of authors. Lombard, ',' Grine and Fowles,3,4 Zharkov, 5 Hughes and McQueen,6 Chabai,7,8 and Bass9 have published equation of state data on rock materials similar in structure or composition to those investigated in the study. Rinehart10 has prepared a voluminous tabulation of the available equation of state data from these publications and in addition has included data on many other rock and earth materials in his compilation. The reader is referred to the extensive tabulations and bibliography of Rinehart for classical studies involving the equations of state of metals and other solids and liquids. The problems of sympathetic detonation between drill holes in hard rock environments has received the attention of various authors as have the problems of explosive rock coupling close-in stress wave attenuation and explosive diameter effects. 11 to 18 In particular, the don or-receptor geometry of this report is very similar to the geometry of adjacent drill holes in rock. Most published pulse velocity and attenuation work in rock has been devoted to velocities at much lower amplitudes than those associated with close-in explosives phenomena and most attenuation studies involve transmission and attenuation across large spans of material. Excellent work in these areas has been published by the U. S. Bureau of Mines and by many individuals and groups engaged in geophysical research. SELECTION OF TEST MATERIALS A series of rock materials was desired for this test program that would be meaningful industrially, and that would also yield comparative results from a research viewpoint. Based on the availability of a series of related petro-fabrics, for which there is an abundance of physical and shock response data, three igneous rocks were selected as the test series. The applicability of shock wave studies in igneous rocks to industrial problems is readily demonstrated. Hardrock tunneling, quarrying and mining take place in igneous rocks throughout the world. From a research point of view, igneous rocks exhibit a wide variety of grain sizes and considerable diversity in composition, yet the common igneous rocks have a general similarity in their constituent mineral structures and in their grain-to-grain bonding mechanisms. The latter is based on the concept of rock formation through magma cooling and crystallization, though the degree of liquefaction in granitic plutonic masses is indeed a moot question. Three specific types of igneous rock were chosen for study. A coarse-grained granitic textured rock was desired that would serve as a representative of the widely abundant intrusive rocks of the world and

Jan 1, 1967

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By P. F. Gnirk, J. B. Cheatham

Single bit-tooth penetration experiments under static load were conducted on six rocks at confining pressures of O to 5,000 psi using sharp wedge-shaped teeth with included angles ranging from 30 to 120° In general, the force-displacement curves for all rocks exhibit an increasingly nonlinear and discontinuous behavior with decreasing confining pressure. The confining pressure at which a rock exhibits a macroscopic transition from predominantly ductile to predominantly brittle behavior during penetration vmies from about 500 to 1,000 psi for the limestones to greater than 5,000 psi for dolomite. The correlation between calculated values of force per unit penetration based on plasticity theory and experimental values is quite encouraging, even at confining pressures as low as 1,000 psi A qualitative correlation between volume of fragmented rock per unit energy input for a single bit-tooth and drilling rate for microbits appears to exist over a confining pressure range of 0 to 5,000 psi. INTRODUCTION Laboratory experiments utilizing a small-scale drilling apparatus have demonstrated that penetration rates are reduced considerably as a result of increasing the confining pressure from atmospheric to a few thousand psi?-3 This undesirable situation can, in general, be attributed to a combination of decreased efficiency of chip removal at the bottom of the borehole, increased rock-failure strength, and a possible change in the mechanism of chip generation and rock fragmentation with increasing confining pressure. To more fully understand the principles underlying the last circumstance, it is the purpose of this investigation to experimentally study the mechanism of single bit-tooth penetration into dry rock at low confining pressures and, in particular, to establish the confining pressure at which the penetration mechanism may undergo a brittle to ductile transition for various rock types commonly encountered in drilling. Confining pressure as used here refers to the differential pressure between the borehole fluid pressure and the formation-pore fluid pressure. EXPERIMENTAL PROCEDURE Using an experimental apparatus previously described: a single, sharp wedge-shaped tool was forced under a "statically" applied load into an effectively semi-in finite dry rock sample subjected to a prescribed confining pressure. To prevent the invasion of the con fining-pressure fluid into the pores of the rock sample during penetration, the exposed surface of the rock was jacketed with a layer of silicon putty.* Electrical instrumentation incorporated into the apparatus yielded a graphical plot of force on the tool as a function of penetration or displacement of the tool into the rock during an experiment. During the course of the experimentation the following conditions were maintained constant: (1) pore pressure — atmospheric (i.e., the rock was dry); (2) temperature — 75F; (3) rate of loading — essentially static (approximately 0.002 in./sec); (4) bit tooth — a sharp wedge-shaped tool loaded normal to the rock surface; (5) rock surface — smooth and flat; (6) drilling fluid — hydraulic oil; and (7) maximum depth of penetration — approximateiy 0.1 in. In addition, each experiment was performed on a different rock sample so the rock surface is free of a layer of cuttings and of any previous indentation craters. The influence of the corners of a borehole was neglected, since each rock sample was cemented inta a section of aLuminum tubing to simulate a semi-infinite body. Indentation experiments were made on a sandstone, two limestones, a marble, a dolomite and a schist with a 60' bit-tooth over a range of confining pressures from atmospheric to 5,000 psi and at a constant confining pressure of 1,000 psi for a variety of bit-teeth with angles ranging from 30 to

Jan 1, 1966

By Giovanni Barla, Stefan H. Boshkov

This paper examines the application of gelatin mixture models to the study of factors such as gravity effects, tectonic and residual stresses, anisotropy, physical nonlinearity, time-dependent behavior of the rock mass and connectedness of the geologic space. The technique of calibrating gelatin mixtures proposed by Richards and Mark in 1966 is applied to the particular problem considered. Similitude relationships between model and prototype are discussed with reference to the constitutive equations describing the physical behavior of the rock mass. The interaction between two openings in a linearly elastic medium, under conditions of plane strain, is used as an illustrative example with the view to describing the history of stress and deformation changes on the rock medium and lining around a slusher drift while undercuts are made in block caving. The prediction of the displacement, stress field and mode of fracture associated with a given rock structure, its initial stress field, and its physical properties is essential to rock mechanics and ground control. The basic approach is to assume the rock mass to behave as a continuous, homogeneous, and isotropic medium. The classical theory of elasticity is the simplest theory based on such a concept. Specification of the elastic constants, the initial state of stress in the ground, and the boundary conditions on stress and displacement allows one to predict stress and displacement within the body. Subsequently, provided that a criterion which governs rock fracture under various stress conditions has been established, a stability analysis of the rock structure can be performed. A natural extension of the present knowledge consists of removing some of the simplifying assumptions which characterize the analytical solutions based on the classical theory of elasticity. Complex geometry, the significance of tectonic and residual stresses, anisotropy, physical nonlinearity, time-dependent behavior of the rock mass and connectedness of the geologic space are some of the important factors which are neglected in the classical approach. A better insight into the real problem can be attained with either the finite element method of stress analysis or an experimental technique. The former has already been applied with success on several occasion 1-4 and is destined to render great service in the field of rock mechanics and ground control. The latter, which will be considered in this paper, has been used extensively in applications too numerous to be enumerated. Experimental technique involves either prototype testing, if the area to be studied is accessible, or model testing. Two- and three-dimensional photo-elastic analyses are very convenient means of obtaining the stress distribution in models of complex geometry. In general, such analyses are confined to studies of stress fields around openings which are sufficiently far removed from the surface boundary so that the stresses in their neighborhood are practically equal to those produced by an initial uniform stress field; i.e., so that gravitational effects are negligible. There are, however, many rock structures where these effects must be taken into account to attain a more realistic representation of the problem. Surface and near surface rock structures, systems of openings interacting with one another, and structural components loaded only by their own weight are some of the obvious examples. This paper was written to show how factors neglected in the usual analysis based on the classical theory of elasticity, can be taken into account when gelatin mixture models are used inconjunction with photoelastic methods. The interaction between two openings in a linearly elastic medium, under conditions of plane strain, is used as an illustrative example, with a view to describing the history of stress and deformation changes on the rock medium and lining around a slusher drift while undercuts are made in block caving. The birefringence of gelatin mixtures has been known for some time. A series of experiments were carried out by Rossi in 1910,5 but after that very

Jan 1, 1970

By L. W. Holm

This study shows that in the presence of foam, gas and liquid flow separately through porous media representative of reservoir rock. These results were obtained by using tracer techniques to measure the flow of the gas and liquid comprising the foam. Foam does not flow through the porous medium as a body even when the liquid and gas are combined outside the system and injected as foam. Instead, the liquid and gas forming the foam separate as the foam films break and then re-form in the porous system. Liquid moves through the porous medium via the film network of the bubbles and gas moves progressively through the system by breaking and re-forming bubbles throughout the length of the flow path. The flow rates of the gas and liquid are a function of the number and strength of the films in the porous medium. There is no free flow of gas; i.e., no continuous gas phase. On the basis of these results, foam can be expected to improve a water flood or gas drive by decreasing the permeability of the reseruoir rock to a displacing liquid or gas. This improves the mobility ratio and thus the conformance of the flood. INTRODUCTION Foam is formed when gas and a solution of a surface active agent are injected into a porous medium either simultaneously or intermittently. During the past few years, several papers have been published on the subject of foam flow in Dorous media. Foam has been used successfully in the removal of capillary water blocks from producing formations.''2 The use of foam in gas storage reservoirs to reduce gas leaks and to increase storage capacity has been considered in recent years. Foam has also been investigated as an oil displacing agent,3-5 and as an agent to improve the mobility ratio in a waterflood.6-10 However, the mechanism by which the gas and liquid phases comprising the foam flow through a porous medium has nor been described adequately. Normally, when two immiscible phases (gas and liquid) flow concurrently through a porous medium, each phase follows separate paths or channels. At given saturations of the two phases, a certain number of channels are available to each phase, and as saturations change, the number and configuration of the channels available for each phase also change. The effective permeability of each phase is a function of the saturation of that phase only, and the flow of each phase can be described by Darcy's law. When foam is present, the effective permeability of the porous medium to each phase is greatly reduced compaced with permeabilities measured in the absence of foam.6,11-13 Based upon the observed flow of surfactant solutions and gas in capillaries, it has been concluded that the gas and liquid may flow separately or they may flow combined as foam. At least four mechanisms have been postulated to explain how fluids flow with foam present: 1. A large portion of the gas is trapped in the porous medium and a small fraction flows as free gas, following Darcy's law.9 2. The foam structure moves as a body; the rate of gas flow is the same as the rate of liquid flow. 6,13 3. Gas flows as a discontinuous phase by breaking and re-forming films. Liquid flows as a free phase.12 4. A portion of the liquid and gas move as a foam body while excess surfactant solution moves as a free phase.14 It also has been suggested that different flow mechanisms exist for high quality (dry) foams made from dilute surfactant solutions and for foams made from more concentrated solutions. Studies conducted on the flow of foam through capillaries have shown that a plug-type flow occurs and that foam flows as a body.6,15 Fried also has described the flow of foam through porous media as

Jan 1, 1969

By H. S. Choi, I. Iwasaki

The effect of temperature as a parameter in ore flotation has not been systematically studied, although for some ores it has been known for many years that selectivity and grade of concentrates can be improved by conditioning or flotation at moderately or substantially elevated temperatures. In 1934 Coghill and Clemmer' stated that "elevated temperature seems to make hard water more tolerable" in soap flotation, and a cursory examination of the literature indicates that some earlier investigators were aware of beneficial results obtained by floating above room or mill-water temperatures. Tartaron,2 Mitchell et aL.,3 and particularly Hamilton et al.' have shown that high pulp temperatures give high recovery and selectivity in soap flotation of fluorite from a variety of gangue minerals. As a result of the work of these and of other investigators, high-temperature flotation has been successfully applied to industrial flotation of fluorite. Falconer" stated: "The use of high temperature conditioning (above 35°C and preferably above 60°C) in connection with the soap flotation of nonsulphide ores, as covered in B. Kalinowski's French Patent No. 847,-215, Dec. 7, 1938, is claimed to provide more intense activation, better separation of values from gangue, and reduction in quantity of reagents." The various investigations quoted above clearly demonstrate the improvements to be realized by increasing the temperature at which fluorite is floated, and it is surprising that little or no information has been published on systematic investigation of the same method as applied to other ores. The present article describes comparative results obtained at room temperature and at elevated temperatures by otherwise conventional flotation of iron ores, by single-mineral flotation in the Hallimond tube, and by contact-angle measurements. It is not the purpose of this article to evaluate the commercial feasibility of high-temperature flotation of iron ores, although preliminary calculations show that the cost should not be excessive and that the gains to be realized may well outweigh the extra cost of heating the pulp. LABORATORY FLOTATION TESTS Experimental: All flotation tests reported in this article that use a laboratory Fagergren cell were made on ore A (either containing the primary slimes or deslimed at 20µ as specified) following the preparatory procedures described elsewhere- or this particular ore. Ore A contained hematite, goethite (limonite), quartz, chert, and minor quantities of magnetite. Confirmatory tests made on other ores show that the same general principles apply to ores of similar mineralogy and structure. Conditioning was performed by adding the collector, as a soap, to the vigorously stirred and heated pulp. Conditioning time was 5 min, the pulp contained 60 pct solids, and the equivalent of 0.5 Ib of fatty acid was added per ton of original (unde-slimed) ore. The pH of the pulps ranged from 6.8 to 7.1 (25°C measurement) during conditioning. The pH was adjusted in the flotation cell by addition of either NaOH or H2SO4, the values reported being those of the pulp cooled to 25°C. During conditioning and flotation the pulp temperatures were controlled within l° to 2° by immersion of Variac-controlled Calrod elements in the conditioning beaker and in the cell, and by the use of water, heated to the appropriate temperature, for transferring pulp, adjusting pulp volume, and washing. Selectivity indexes have been calculated on the same basis as those reported in the article by Cooke, Iwasaki, and Choi." Results: Table I gives the results obtained by floating deslimed ore A with pure oleic acid at 25" and 50°C, and Fig. 1 shows the selectivity indexes plotted against pH of flotation. Inspection of Table I

Jan 1, 1961

By G. A. Stamboltzis, N. Arbiter, C. C. Harris

Fracture under low-velocity free-fall and double impact and under slow compression have been investigated. The pattern of breakage and the size distribution of resulting fragments of sand-cement and glass spheres have been determined. Photoelasticity methods were used to simulate the stress distributions in free-fall impact in order to explain the observed patterns of breakage. Oblique fracture planes, occurring only in free-fall impact, develop along lines which coincide with the trajectories of maximum compression as determined by mathematical analysis. Breakage efficiencies for different modes of fracture were compared for both types of spheres. For the same specimen and loading system, static loading and low-velocity dynamic loading induce geometrically similar stress fields resulting in reasonably similar fracture patterns and shapes of fragments. This work had its origin in a desire to understand three fundamental processes involved in autogenous comminution: namely; free-fall impact, double impact, and slow compression. To simplify the investigation, it is preferable to reduce the entire multi-stage fracture process to its most elementary form; that is, to study products resulting from a single-stage operation, such as the breakage of a specimen under single fracture conditions. The goal of this research is to study the energy utilization in the single fracture of brittle specimens and to relate the pattern of breakage and the resultant fragment size distribution with the nature of the material, the specimen size, the manner of load application, and the rate of loading. Published works on single fracture have mainly been concerned with tests on large (> 1-in.) irregular mineral pieces, especially coal and coke. These were mainly friability tests employing a single blow on closely sized pieces. More recently, small mineral specimens in the sieve range and glass spheres have been investigated.5-7 Kick8 has reported free-fall and double impact tests with large large cast-iron, cement, and clay spheres. The establishment of a minimum breaking height independent of size in the free-fall impact tests, and the direct proportionality between the minimum work required for fracture and the volume of the specimen in double impact tests were used by Kick as an experimental proof of his classical law of "Proportional Resistances." In the present work, spherical shapes were chosen because of their simple body geometry and consequent impact and stress field symmetry. Because the study involves several physical principles in connection with brittle-fracture, it may be of interest in fields where the strength of materials is of importance. APPARATUS The equipment for the free-fall impact testing is illustrated in Figs. 1 and 2. This consisted of a massive (17 x 17 x 3-in.) hard steel plate onto which specimens were dropped and a specimen release mechanism which could be set at any desired height up to 10 ft above the plate so that the impact velocity could be varied up to about 25 ft per sec. Double impact breakage is characterized by two points of loading situated at opposite poles. For this mode of breakage the apparatus was modified (Figs. 3, 4 and 5) so that a falling mass provided an impact of predetermined magnitude on a specimen resting on the plate. A spring-loaded device operated immediately after fracture to arrest the falling mass and to record its residual kinetic energy. Precautions were taken to avoid secondary breakage. Slow compression tests were performed on a conventional hydraulic testing machine, the specimen being held between two hardened parallel bearing blocks. In the case of glass spheres, carbide bits were used. Photoelasticity studies were performed on two-dimensional models held in a loading frame equipped with a force gauge (Fig. 6). The models were viewed in a polariscope and isochromatic fringe patterns were recorded photographically. SPECIMENS Sand-Cement Spheres: A series of sand-cement spheres were prepared by molding in spherical glass flasks.

Jan 1, 1970

By Robert Maes, Robert de Strycker

The Fe-Ni-As phase diagram has been established by the study of about a hundred alloys, by microscopic observation, and by thermal analysis, with arsenic contents up to 50 pct. The iron and nickel arsenides present extensive solid-solution fields, owing to the substitution of nickel by iron and vice versa; the extent of the solubility field of each compound has been determined with an accuracy of ±1 pct. In the investigated range of compositions, the solidification reactions were established, and the temperatures of the invariant reaclions detevmined with a preciston of iZ°C in the most favorable and ±5°C in the least favorable cases. The isothermal lines of the liquidus surface have also been drazum, with an accuracy estimated at i5°C. Reactions in the solid state, which take place for the formation or the decomposition of certain phases , were investigated in detail. ThE constitution diagram of the Fe-Ni-As system is fundamental for the understanding of the properties of nickel speiss; these by-products of the extraction of certain nonferrous metals indeed often contain the three elements iron, nickel, and arsenic as main constituents. The Fe-Ni-As system has already been the object of earlier investigations. In 1932, Guertler and Savels-berg' have presented some elements of the ternary phase diagram, for arsenic contents up to 55 pct (all percentages in this paper are given in weight percent), including the vertical section FezAs-Ni&s2. This investigation is, however, incomplete and certain anomalies suggested the necessity to verify the results published by these authors: the section Fe2As-Ni5AsZ, for example, is presented as quasi-binary, with a field of complete miscibility in the solid state, even though these compounds do not have the same structure. Recently, ~useck' established an isothermal section at 800°C in the Fe-Ni-As system by X-ray diffraction and microscopic examination of water-quenched alloys. These techniques are sometimes inaccurate for the determination of the fields where alloys are liquid at the investigated temperature, and they may lead to erroneous conclusions when a compound exists at the investigated temperature but is not stable at room temperature and cannot be maintained by quenching. LIMITING BINARY PHASE DIAGRAMS The Fe-As constitution diagram has been the object of several investigations which have been reviewed by Hansen and Anderko.3 For the Ni-As phase diagram, a recent study has been effected by Yund.4 In this system, Heyding and calvert5 have determined the existence of a compound of unidentified structure at arsenic contents slightly higher than those corresponding to Ni5As2 and at temperatures lower than about 200°C; by analogy with iron and cobalt arsenides, this compound could correspond to the formula Ni2As, as suggested by Kulle-rud; although this is not definitely established. In the region of arsenic contents from 35 to 55 pct, Fried-rich7 detected anomalies in the solidification reactions; an interpretation of these anomalies was given by Hansen ,8 which assumed that the solidification reactions in practice are not in equilibrium, but are meta-stable. The main features of the Fe-Ni constitution diagram are the existence of a complete miscibility field, at least at high temperatures, for the fcc phase (which will be designated by My in the remainder of this work) and of a limited solubility field for the bcc phase (designated by Ma). EXPERIMENTAL METHODS The fundamental technique used in this investigation was microscopic observation, which allowed the determination of the reactions occurring during solidification of the alloys, and possible reactions in the solid state.

Jan 1, 1968

By C. D. Thompson, D. C. Violetta, C. A. Czako

The mechanics of the continuous sintering process are briefly reviewed and the application of this process for the beneficiation of principal minerals used in making cement is shown. Pilot plant tests using pyrobeneficiation techniques were conducted with mixtures of limestone and shale in various ratios. Results are illustrated and show the percentage removal of carbon, carbonates, sulfur, alkalies, and volatiles as a function of the limestone-shale ratio. Tests with varying quantities of solid fwls, added to mixtures containing limestone and shale in a 4 to 1 ratio, were made to attain optimum elimination of the volatile constituents. Experimental results show that 90 pct of these constituents were removed during the calcination sintem'ng process. The advantages of using the beneficiated material for making cement are discussed with regard to prepam'ng, handling, and clinkering operations. A sintering process for treating raw materials or kiln waste products which contain high amounts of alkaline materials is described. Tests were performed with a high alkaline waste dust mixed in a 1 to 1 ratio with raw cement meal. The data from these tests and others using different halides for aiding alkali removal are evaluated by the alkali content of the cement clinker. Test results showed that the Dwight-Lloyd process could be applied with facility in beneficiating raw materials for cement making. Substantial quantities of deleterious volutile constituents were removed by calcining and sintering techniques. The continuous sintering process has been applied for a number of years as a means of beneficiating minerals for subsequent conversion to merchantable products. The process was developed and first used for nonferrous metallurgical applicatfons. However, it has risen to prominence for ferrous metallurgical applications where last year in the U.S. and Canada nearly 40 million tons of iron ore were treated by grate firing processes. Through the years of process development many investigators have realized the versatility of the sintering process and have directed experiments toward the use of it for applications with nonmetallic minerals. These have included calcining limestone, carbonizing coal, expanding shale, and burning cement raw meal. Today commercial sintering plants produce cement clinker from mixtures of limestone and shale, lightweight aggregate from shale or clay, and calcined products from phosphate shale. The mechanics of the downdraft sintering process have been presented by many authors. In general practice the materials to be agglomated are mixed as moistened fines with a solid fuel and this mixture is spread uniformly on the grates of the sintering machine. The charge is then momentarily directed under an ignition furnace. After ignition a continuous downdraft of air promotes combustion of the fuel and through heat transfer particles become partially fused into a coherent mass. This is cooled from continued draft induced after combustion. Many investigators have shown that distinct zones are formed during the process of sintering iron ores. At the midway point of the sintering process these are prominent and are, respectively, from the top of the bed to the bottom, comprised of a cooling zone of sintered material, a combustion zone of sinter formation, a calcination zone, a drying zone, and a zone of unreacted materials. Completion of sintering occurs near the discharge end of the machine where the firing zone approaches either the grates or grate dressing material called hearth layer. Here the peak draft temperature usually can be noted in the windboxes. The zonal pattern has been illustrated with the use of graphs which show thermal gradients as they progress through the bed of materials during sintering. Fig. 1 shows, in simplified form, thermal gradients existing when clinkering cement on a sintering machine. The thermal gradients illustrated are for periods of initiation, midway and completion of the clinkering process. The occurrence of these zones has importance for the beneficiation of cement raw materials. A research program was performed at the Dwight-

Jan 1, 1961

By W. A. Riggs

In a single year the total transportation cost equals nearly 30 pct of the value of mineral commodities, the largest single cost from the deposit to consumer. The magnitude of this economic factor calls for more complete understanding of cost and operational problems between producer and carrier than now exists. Details are given of the costs of providing transportation as well as freight rates for selling transportation. In 1958, transportation of mineral commodities in the United States required 600 billion ton-miles to move 1 3/4 billion tons of products of mines comprising solid fuels, crude oil, industrial minerals and rocks, and ores. Possibly 17 pct was handled by more than one mode of transport, since combinations of several modes of transport may, at times, achieve better service/cost results than any single mode. Total transportation cost of nearly $5 billion equaled roughly 30 pct of the value of the mineral commodities, the largest single cost from the mineral deposit to the consumer. The tremendous magnitude of this transportation undertaking demands much more complete understanding than now exists between mineral producers and carriers about the economic and operational problems of each other. A substantial portion of the daily work of the writer is spent in trying to act as a combination interpreter and catalyst between potential shippers and carriers, hoping to direct both to points of mutual understanding and advantage. For the purpose of this paper, it is necessary to ask mining people to accept some terminology and to try to appreciate the mental approach of the freight rate making people of the carriers. This paper will consider costs of providing transportation as well as freight rates for selling transportation—although these are fields of limited certainty. There are no accurate statistical breakdowns for all modes of transport and types of carriers indicating commodities, volume of traffic, distances, costs, and rates. Much information is available for railroads, and some for the regulated portions of air, motor, pipeline, and water transportation. Unregulated carriers, handling about half of the transportation of products of mines, are not required to report, and very little is known of their activities. Figures used herein are national in scope, unless otherwise indicated, and have generally been derived from Interstate Commerce Commission reports and files. Tonnage figures have been reduced to net tons of 2000 lb. The very generalized comparison of average mileage freight rate levels for various modes of transport provided by Fig. 1 affords advantageous perspective for the entire field of transportation costs. Dashed guide lines show cents per ton-mile, and dotted guide lines show percentage of first class railroad freight rate. While costs of other than railroad transportation more nearly follow constant cents per ton-mile functions, costs of railroad transportation show characteristic taper, or decrease in cents per-ton mile with increasing length of haul. This fundamental difference in cost is due largely to the fact that railroads and pipelines provide their own fully taxed roadway facilities while highway, waterway, and air carriers use government-provided tax-free roadway at no cost other than fuel, license, and certain excise taxes. Thus railroad and pipeline transportation might be said to carry heavy threshold costs in addition to those transportation costs which are directly variable with ton miles of transportation produced. The original basis of railroad freight rate structures lies in the classification of all articles of commerce into railroad class ralings ranging from 7 to 40 pct of first class (100 pct) rate level, dependent on value, hazard, loading characteristics, service requirements, and whether carload or less than carload. Classification tables and mileage-based (but arbitrarily grouped), rate scales are published as tariffs subject to approval by regulatory bodies. Used together, these determine class rates in cents per 100 lb between points in the U. S. Where local commercial or operating conditions render it advantageous for carriers to depart from

Jan 1, 1961

By H. L. Hartman

The estimation of the magnitude of pressure changes which occur in mine ventilation circuits is of primary importance to the mining engineer in making changes in an existing mine or in projecting the ventilation requirements of a new mine. It is a formidable task in many cases, because no adequate theory is available and empirical data are unreliable or even lacking. The common practice of increasing one's estimate of the mine static pressure by an arbitrary "safety factor" to compensate for the unknown is a poor substitute for accuracy. Generally little difficulty is encountered in calculating friction losses in mine airways. As long as the friction factor can be reliably estimated, the well-known Atkinson formula gives accurate results. However, it is in the determination of shock losses that the ventilation engineer encounters difficulty and is apt to resort to an arbitrary allowance. Shock losses may constitute as little as five or as much as twenty-five percent or more of the overall pressure drop in the mine, and such allowances are generally completely unsatisfactory. Since other shock losses have been covered by McElroy (1), this paper deals with a long-neglected source of shock and pressure change in ventilation circuits, divergence and convergence of airflows, here referred to by their more common mining terms, splits and junctions. PREVIOUS WORK There is practically no fundamental information available in ventilation literature on the splitting and junction of airflows, either in mining or mechanical engineering periodicals or handbooks. Fluid mechanics textbooks seemingly avoid the subject of divergence and convergence, except for brief empirical treatment. The only investigation reported for mine ventilation systems was by Weeks, et a1 (21), in 1933. Although applicable to the case of line brattices in coal mines, it cannot be applied to the general case of mine splitting and junctions involving change of direction with little or no area change. In dealing briefly with the subject of shock losses at splits and junctions, McElroy (1) suggests only that they be considered comparable to bends, with losses at junctions being increased 50 pct to allow for interference of merging streams. He recognizes the importance of quantity ratio but concludes that exact quantitative data are lacking to compute losses at splits and junctions. A very capable summary and analysis of the state of knowledge in industrial ventilation and fluid mechanics regarding losses due to divergence has recently been reported by Gilman (3). In considering all earlier work dealing with both air and water as fluids, he compares results on a common basis and obtains good empirical agreement. For divergent flow, the controlling variable in each branch is the quantity ratio, although the deflection angle is of secondary importance. Empirical equations are presented for straight and 90-deg branches with varying velocity ratios, assuming branches smaller than the main duct. These formulas are modifications of the Borda equation for abrupt contraction and expansion. The only information available on shock losses in convergent flow are approximate experimental data presented by Alden (4) and the Manual of Recommended Practice for Industrial Ventilation of the American Conference of Governmental Industrial Hygienists. Alden considers the effects of both quantity ratio and deflection angle. The data indicate an effect opposite to that observed in divergence, the shock loss varying directly with quantity ratio. SHOCK LOSS THEORY Because pressure losses due to shock have been found to bear a constant relation to the mean velocity of flow in a given conduit, they are frequently expressed as a dimensionless function of the velocity head, termed the shock loss factor. For airflow in ducts and mine airways, the shock loss may be represented by the formula (1): Hx - XHv, (1)

Jan 1, 1961

By P. E. Baker

An experimental study of heat flow in steam flooding was carried out with steam displacing water in a plane-radial fluid-flow model. Temperature distributions in the model reservoir, overburden, and substratum were measured at frequent intervals during flooding by means of a fixed array of thermocouples. Heat losses to overburden and substratum and the division of heat between the steam and hot-water zones in the reservoir were calculated from the temperature data by performing numerical integrations over the appropriate volumes. Temperature profiles and heat distribution results are presented. Some results were (1) the fraction of injected heat that is lost to overburden and substratum, when expressed as a function of time, did not depend on injection rate; (2) a significant portion of the injected heat was contained in the hot-water zone ahead of and underlying the steam front in the reservoir. The division of heat between the steam and hot-water zones did depend on injection rate; and (3) considerable gravity override was observed. The experimental results are compared with those of several published theoretical studies of heat transfer in hot-fluid injection; approximate agreement on heat losses is found with some of the most idealized theories. A valid comparison cannot be made in regard to steam us hot-water zone; basic interpretation of the hot-water zone, in terms of a vertical steam front, is clouded by gravity override in the experiments. INTRODUCTION The experimental study reported here primarily concerns two problems of heat transfer in steam flooding that have economic and engineering importance. These are (1) the fraction of injected heat lost to the overburden and substratum, and (2) the division of the remaining heat between a zone containing steam and a zone of hot water at lower-than-boiling temperature. The volume of the steamed region, which has a direct effect on oil recovery, can be calculated readily from these results. Several valuable theoretical studies have been published on this problem, but no experimental tests of the theories have previously appeared. Tests were needed because the effects of the necessary simplifications in the theories could not be determined otherwise. The basic, untested theoretical problem is the effect of "decoupling" heat and mass transfer, and the resulting inability to include the effects of heat transfer and condensation on fluid flow. The theoretical problem solved so far has been a heat transfer problem in which fluid is present as a heat carrier and flows in a stipulated manner, either linear or plane-radial. The simplifications that have sometimes been made in regard to thermal conductivity by, for example, Marx and Langenheim,l Lauwerier2 and Thomas3 apparently do not have an important effect on heat loss calculations, although the results of the different methods disagree to a certain extent, as shown by Spillette4 and by Ramey.5 A recent paper by the author6 compared the results of experiments in hot-water flooding with several theoretical results and showed a rather poor agreement. It will be seen that the experimental heat loss results for steam flooding, however, show acceptable agreement. These comparisons are discussed in a later section. EXPERIMENTAL APPARATUS AND PROCEDURE The system studied was one in which a sand pack, initially water saturated and at a uniform temperature, was flooded with steam at a constant rate and temperature. The temperature distribution in the "reservoir" and adjacent "strata" was measured with an array of thermocouples, and heat losses were determined from these measurements. The sand pack contained no oil because the only concern of the experiments was heat flow, and it was important to have the flow pattern as simple

Jan 1, 1970

By A. L. Benham, W. E. Dowden, W. J. Kunzman

The mechanism involved in creating a miscible fluid displacement of typical reservoir fluids through porous media by light hydrocarbon mixtures usually involves a number of mixing and separation stages in the contact area. This mechanism has been illustrated on ternary diagrams for miscible displacement of reservoir fluid by rich hydrocarbon gases. A relatively simple method has been devised for calculating the approximate conditions for a miscible displacement of reservoir fluid by rich, light hydrocarbon gases or by LPG mixtures. Miscibility is favored by increased pressures, decreased temperatures, light reservoir fluids and light hydrocarbon-displacing fluids rich in the heavier components. The conditions for miscible displacements of five typical reservoir fluids by various light hydrocarbon mixtures were determined for a range of temperatures from 70" to 260°F and pressures from 1,000 to 3,000 psia. A correlation of these conditions was obtained using the variables of temperature, pressure, C,+ molecular weight of the reservoir fluid, C,+ molecular weight of the displacing fluid, and mol per cent methane in the displacing fluid. This correlation has been applied to data appearing in the literature with very good success. A partial check of the correlation has been made using a step-wise procedure with a windowed PVT cell. INTRODUCTION The term, "miscible fluid displacement", may be defined as any oil-recovery displacement process where there is an absence of a phase boundary or interface between the displaced and displacing fluids. A water flood would not be a miscible fluid displacement, although cycling in a condensate reservoir would be. Another example is displacement of oil by gasoline. Not so obvious, however, are the conditions under which fluids (such as propane, mixtures of propane and methane, or other similar systems) will give a miscible fluid displacement of oil. There are four standard miscible-fluid-displacement processes being used for the recovery of oil from an "nderground reservoir—(1) high-pressure dry-gas miscible displacement, (2) rich-gas miscible displacement, (3) miscible slug using a slug of LPG followed by a lean dry gas under miscible conditions, and (4) gas driven-slug process using a slug of LPG followed by a lean dry gas under immiscible conditions. The first problem faced by the engineer designing a miscible displacement process is to determine the conditions for miscibility. A number of papers have included discussions of the process of obtaining miscibility in these processes.1-8 However, methods for calculating the conditions for miscibility have not been published. This paper includes a method for calculating the conditions for a miscible displacement of reservoir fluid by a rich gas (LPG diluted with methane or natural gas). In this paper, a rich-gas miscible displacement process is defined as an oil-recovery process in the reservoir. The cost of a rich-gas miscible displacement will be decreased by diluting the LPG with methane or natural gas as much as possible. Therefore, it is important to be able to predict the maximum dilution for miscibility. MECHANISMS The mechanisms for obtaining miscibility between a displacing mixture of light hydrocarbons and displaced reservoir fluid often involve a number of displacing and mixing stages. As stated previously, the cost of a rich-gas miscible displacement is decreased by diluting available LPG with natural gas. The problem is to predict the maximum dilution allowed while maintaining miscibility. To get a picture of the displacing and mixing stages, it is helpful to consider these complex hydrocarbon mixtures as pseu do-ternary systems. Graphical representation of a complex system as a pseu do-ternary system can be thermodynamically rigorous if enough variables are used to describe the character of each pseudo-component in one of the phases. However, even if the characterization of each pseudo-component is not complete, pseudo-ternary representation does serve to illustrate the mechanism for obtaining miscibility and can be used to make some general conclusions concerning the conditions for obtaining miscibility. A phase diagram for the ternary system, methane-

By P. H. Turnock, R. D. Pehlke

The effects of mutual interactions between elements on the solubility of nitrogen in liquid iron alloys containing chrormium, columbium, molybdenurn, nickel, or silicon have been determined. The equilibrium nitrogen solubilities in the liquid iron alloys were measured by the Sieverts method. The second-order effects caused by the presence of chrormium, columbium , or silicon in the melt were found to be significant. The solubility of nitrogen in liquid iron alloys containing several alloying elements has also been measured as a function of melt temperature and nitrogen pressure. The heat of solution of nitrogen in microcomponent iron alloys has been found to be a function of the logarithnz of the activity coefficient of nitrogen, irrvespectilje of the composition of the melt. The heat of solution of nitrogen in pure liquid iron was determined to be 1265 ± 450 cal per g-atom of N over the temperatxue range 1600° to 1800°. Sieverts ' law was obeyed for all melt compositions studied in the pressure interval 0.4 to 1.0 atm. The solubility of nitrogen in liquid pure iron and in many dilute binary liquid-iron alloys has been the subject of many investigations. Pehlke and Elliott' made a comprehensive investigation on this subject in which they summarized and compared previous researches. However, the practical application of these solubility studies requires prediction of the nitrogen solubility in the multicomponent alloys encountered in operating practice. The interaction parameter as defined by wagner2 has been used to predict the behavior of a solute in a complex alloy, and a modification of this approach has met with reasonable success in several instances.'1"5 The necessity for extending the nitrogen solubility data on binary alloys is caused by the lack of direct measurements of the solubility of nitrogen in molten alloy steels. A direct comparison between measured solubilities and those predicted for multicomponent alloys based on data in binary systems has been made in only a few instances, and most of these were reviewed by Langenberg,6 who developed a graphical technique for predicting nitrogen solubilities in multicomponent iron alloys. The solubility of nitrogen has also been studied over the entire Fe-Cr-Ni ternary system by Humbert and Elliott.? Both studies compared the solubilities computed from binary iron alloy data with those measured in ternary alloys in the liquid state at 1600°C. Most of the nitrogen solubility measurements referred to above were for iron alloys at 1600°C; thus, the ability to predict nitrogen solubilities in molten iron alloys has been restricted primarily to alloys at 1600°C. Recently, Nelsson8 devised an approximate method for calculating the change in nitrogen solubility in molten iron alloys as a function of temperature, but an assumption of regular solution behavior limits the applicability of the method. Chipman and Corrigan9,10 have developed an empirical correlation. based on the activity coefficient of nitrogen in the iron alloy at 1600°C, which facilitates the estimation of nitrogen solubilities at higher temperatures. This method also has limitations, due to the necessity of calculating the activity coefficient for nitrogen in a multicomponent iron alloy at 1600°C from the nitrogen solubility data for binary iron alloys. The shortcoming of such a straightforward approach to the problem of predicting nitrogen solubilities is the disregard of mutual interactions between elements in a complex iron melt. This investigation was conducted in an attempt to provide a more reliable method of predicting nitrogen solubilities in multicomponent iron alloys. The effects of the mutual interactions between chromium, nickel, silicon, molybdenum, or columbium in iron were studied, using the Sieverts' technique1,7 to measure the effects of alloy composition, temperature, and pressure on the solubility of nitrogen in ternary liquid iron alloys.

Jan 1, 1967

By B. C. Allen, W. D. Kingery

Surface tension and its temperature dependence have been determined for pure liquid Fe, Cu, Co, Ni, and Sn and for Fe-C, Co-C, and Ni-C alloys. The temperature coefficient of surface tension is negative for all these liquids. The surface tension and wetting behavior of tin and tin-titanium alloys on A120,, Si,N,, MoSi,, and Sic, as well as the behavior of titanium-containing nickel alloys on A1 20, have been determined over a wide temperature range. For all the solids investigated titanium additims markedly lower the contact angle. Some limits for the solid-surface energies have been calculated from liquid surface tension and wetting behavior. SURFACE tension and interface energy are important variables in a number of metallurgical and ceramic phenomena, and have been discussed at some length in connection with determining micro-structure of polyphase systems,' solid-state sintering,' thermal etching of grain boundaries, heats of solution and fine powders,'1 grain-growth phenomena,' the mechanical behavior of fine foils and fibers, nucleation, soldering and brazing,' and other phenomenz of academic or practical interest However, until recently few quantitative data have been available for high-temperature systems. Among available data are some indicating positive temperature coefficients for the surface tension of liquid copper,*,11,'2 cast iron,l3 and iron-carbon14 alloys. Although carbon was found to be not surface active in iron at 1570" some preliminary studies had indicated the possibility of a positive temperature coefficient for the eutectic composition at lower temperatures.16 One part of the present study has been largely concerned with determining the temperature coefficient of surface tension for several metals at elevated temperatures, In addition to surface-tension studies, the wetting behavior of alloys containing titanium has been investigated over a wide temperature range. Previous measurements indicated that titanium was selectively adsorbed from nickel at a liquid nickel-solid alumina interface,17 and the use of titanium hydride or titanium-alloy techniques for metal-ceramic soldering is well known. In the present study, previous data have been extended to new alloy systems, and measurements have been carried out with other oxidation-resistant ceramics. EXPERIMENTAL The sessile-drop method was employed to measure liquid surface tension and contact angle between liquid and solid, utilizing apparatus and techniques which have been described previously.'g The shape of a sessile drop depends on equilibrium between forces of surface tension tending to form a spherical surface of minimum area and forces of gravity tending to flatten the drop. Typical sessile drops are shown in Fig. 1. From precise measurements of maximum drop diameter and drop height the surface tension and drop volume were calculated. Four solid materials were employed as supporting plaques for the liquid-metal drops. Aluminum oxide (99.9 pct, J, T. Baker, reagent grade) was pressed and sintered at 1800°C. Molybdenum-disilicide powder (99.0 pct, Electro Metallurgical Co.) was pressed and sintered at 1650°C. Plaques were cut from silicon carbide (99.4 pct, hot-pressed at Alfred University) and silicon nitride (95.1 pct, with Al, Fe, 0 as major impurities, Norton Co.) rods. All plaques were polished, cleaned, and dried prior to use. A number of metals and alloys were employed. Chemical analyses of vacuum-melted ingots of iron, nickel, cobalt (National Research Corp.), commercial nickel-base alloys, electrolytic copper (International Smelting and Refining Co.), and tin (Mal-linckrodt Chemical Works) are given in Table I. Alloys were prepared in the laboratory by vacuum melting the base metal with either titanium hydride (Metal Hydrides, Xnc.) or spectroscopic-grade carbon (National Carbon Corp.). Nickel-titanium alloys contained between 0.001 and 0.005 pct oxygen and the ferromagnetic metal-carbon alloys contained 0.001 pct. Cylindrical metal specimens (1 to 2.5 g) were machined and bevelled at the bottom to insure a uniform advancing contact angle on melting. Before use the samples were acid-etched, cleaned with acetone, and rinsed in petrol ether.

Jan 1, 1960