Search Documents

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By G. W. Doepel, W. P. Sibley

A three-dimensional, multilayer, mathematical model has been developed for predicting performance by mis-c.ible displacement. Areal and vertical coverage, total gas-oil ratio, propane-oil ratio and oil-recovery factor for the multilayer system are computed as functions of hydrocarbon pore volumes of solvent injected. The technique has been programed for use on the IBM-7090 digital computer and was used in calculating perfortnance for continuous miscible displacement of a five-spot, 50-layer .system. Results reported are based on the use of published laboratory data for a single, uniform five-spot layer. However, the technique is sufficiently generalized to calculate performance of any system of layers of any pattern configuration and mobility ratio, provided data are available to describe perfortilance of a single layer. The technique has also been adapted for predicting oil recovery resulting from various gas-driven slug sizes. The results presented for both continuous and slug-type injection show, that permeability variation has 0 significant effect on oil recovery but that solvent composition, within the range studied, has a lesser effect on oil recovery. INTRODUCTION Experience from many of the miscible-phase displacement projects currently in operation indicates that the injected solvent flows preferentially through the more permeable zones in the reservoir. This preferential solvent flow in the high-permeability streaks results in a reduction of vertical coverage (C,) and an attendant decrease in oil recovery. The multilayer technique presented in this paper uses areal coverage (C.,) data for a single layer from the literature and from laboratory measurements of displacement efficiency (herein referred to as microscopic coverage, C,,) to calculate areal coverage (C,) and vertical coverage (C,) for the multilayer system. Oil recovery for the multilayer system is obtained as the product of (Ca) (Cm) (Cv). It is assumed that C, in a single layer is 100 per cent; therefore, (C,) (Cm) is a measure of oil recovery from that layer. Data presented in the literature for a single layer generally deal only with areal coverage and/or microscopic coverage —but seldom specifically with oil recovery. Nevertheless, this technique is flexible in that either Ca and C,, of an individual layer or a direct measurement of oil recovery for the layer may be used. In developing this technique, the problem arose concerning a method of determining the relative amounts of solvent that will enter each layer during any increment of solvent injected. To solve the problem, a dimensionless equation was derived for expressing solvent injected into each layer as a fraction of total injection. A finite-difference method was applied in integrating the equation to determine oil recovery for the system of layers as a function of pore volumes injected. The technique subsequently was programed for solution on a digital computer. The method was applied to predict performance of a deep California reservoir under continuous miscible displacement and has been extended for calculating oil recovery as a function of slug size for gas-driven slug injection. The technique may be applied to any petroleum reservoir that can be adequately represented as a layered system, by considering permeability variation,1-3 permeability range and average permeability. Each layer is then assigned its appropriate permeability, porosity and hydrocarbon saturation. RESERVOIR VOLUME CONTACTED BY SOLVENT It has been recognized for some time that high microscopic coverage (C,) can be obtained by injecting a miscible agent, such as propane, into an oil reservoir. These high values of C,,, have been demonstrated in the laboratory and are normally between 80 and 100 per cent, depending on solvent and reservoir oil composition as well as on reservoir temperature and pressure. In miscible-phase displacement, because there is no interface between solvent and oil, effective permeability to solvent is the same as that to oil. (This assumes that oil saturation has not changed and that water is immobile.) Consequently, mobility ratio (M) simply reduces to that of viscosity ratio . Perrine4 and Habermann5 have shown that large

By R. A. Morse, C. R. Holmgren

The production of oil by water flooding can be substantially increased by the maintenance of free gas saturation in the reservoir during the flooding operation. This effect is accomplished by the alteration of oil relative permeability characteristics and the occupation by gas of pore space that would otherwise be filled with residual oil. The amount of reduction In residual oil can be calculated from appropriate water-oil relative permeability characteristics. This paper presents experimental data in support of the foregoing conclusions and an example of the calculations. The microscopic pore saturation concepts of the mechanism are discussed. A method of practical application to field floods is presented together with discussion of certain limitations. INTRODUCTION The presence of free gas has been reported by a number of investigators to significantly affect the oil recovery which can be obtained from sandstone flow systems by water flooding.1,2,3,4,5 The effect of gas, noted in every instance, has been to cause lower residual oil saturations than could be obtained by water flooding the same systems in the absence of free gas. The degree of improvement in recovery has been observed to vary widely, depending on the systems used and the conditions of the tests. The increased oil recovery obtained because of the presence of gas during a water flood has been variously attributed to changes in physical characteristics of the oil, selective plugging action of the gas, inclusion of oil mist in the free gas phase, and the additional sweeping or driving action of the free gas. All but the first of these suggestions imply changes in the displacement mechanism. The change in viscosity and inter-facial tension of the oil phase, within the pressure range used for all the experimental work, is certainly not sufficient to account for the differences in residual oil saturation noted dess there is a drastic change in the displacement process. One other effect which logically seems capable of causing differences in residual oil saturation of the magnitude noted in the experimental work is that of simple replacement. In a water-wet system containing oil, water, and gas, it is to be expected that the gas will exist inside the oil. This is the position of minimum free surface energy, since the gas-oil interfacial tension will be less than the gas-water interfacial tension. There is no apparent reason to expect that the existence of free gas within the oil phase should alter the saturation at which the non-wetting phase (now oil and gas) should become discontinuous and hence trapped so as to be unrecoverable by direct displacement by water. If this is the situation, then trapping of a certain percentage of gas saturation during water flood should result, at infinite water-oil ratio, in a like reduction of oil saturation below that attainable by flooding in the absence of free gas. It is visualized that the gas will exist as bubbles inside the discontinuous residual oil as illustrated in Fig. 1, with the size of the oil bubbles being substantially unchanged due to the presence of the gas. As a practical matter, it can be anticipated that the presence of a free gas saturation inside the oil phase will reduce the relative permeability to oil which will exist at any particular water saturation. This reduction will be caused by two factors — -the addition of the gas-oil interface, and the reduction of area available for oil flow in the pores containing gas. This reduction in oil permeability at any particular water saturation will result in water breakthrough at a lower water satu-

Jan 1, 1951

By R. A. Morse, C. R. Holmgren

The production of oil by water flooding can be substantially increased by the maintenance of free gas saturation in the reservoir during the flooding operation. This effect is accomplished by the alteration of oil relative permeability characteristics and the occupation by gas of pore space that would otherwise be filled with residual oil. The amount of reduction In residual oil can be calculated from appropriate water-oil relative permeability characteristics. This paper presents experimental data in support of the foregoing conclusions and an example of the calculations. The microscopic pore saturation concepts of the mechanism are discussed. A method of practical application to field floods is presented together with discussion of certain limitations. INTRODUCTION The presence of free gas has been reported by a number of investigators to significantly affect the oil recovery which can be obtained from sandstone flow systems by water flooding.1,2,3,4,5 The effect of gas, noted in every instance, has been to cause lower residual oil saturations than could be obtained by water flooding the same systems in the absence of free gas. The degree of improvement in recovery has been observed to vary widely, depending on the systems used and the conditions of the tests. The increased oil recovery obtained because of the presence of gas during a water flood has been variously attributed to changes in physical characteristics of the oil, selective plugging action of the gas, inclusion of oil mist in the free gas phase, and the additional sweeping or driving action of the free gas. All but the first of these suggestions imply changes in the displacement mechanism. The change in viscosity and inter-facial tension of the oil phase, within the pressure range used for all the experimental work, is certainly not sufficient to account for the differences in residual oil saturation noted dess there is a drastic change in the displacement process. One other effect which logically seems capable of causing differences in residual oil saturation of the magnitude noted in the experimental work is that of simple replacement. In a water-wet system containing oil, water, and gas, it is to be expected that the gas will exist inside the oil. This is the position of minimum free surface energy, since the gas-oil interfacial tension will be less than the gas-water interfacial tension. There is no apparent reason to expect that the existence of free gas within the oil phase should alter the saturation at which the non-wetting phase (now oil and gas) should become discontinuous and hence trapped so as to be unrecoverable by direct displacement by water. If this is the situation, then trapping of a certain percentage of gas saturation during water flood should result, at infinite water-oil ratio, in a like reduction of oil saturation below that attainable by flooding in the absence of free gas. It is visualized that the gas will exist as bubbles inside the discontinuous residual oil as illustrated in Fig. 1, with the size of the oil bubbles being substantially unchanged due to the presence of the gas. As a practical matter, it can be anticipated that the presence of a free gas saturation inside the oil phase will reduce the relative permeability to oil which will exist at any particular water saturation. This reduction will be caused by two factors — -the addition of the gas-oil interface, and the reduction of area available for oil flow in the pores containing gas. This reduction in oil permeability at any particular water saturation will result in water breakthrough at a lower water satu-

Jan 1, 1951

By Pol Duwez

The effect of the rate of cooling on titanium, zirconium, and thallium has been measured. For titanium-molybdenum alloys, it has been shown that for molybdenum concentration up to 8 pct the rate of cooling has no effect on the transformation temperature, but this temperature decreases progressively with increasing molybdenum content. HIGH strength to weight ratio and remarkable corrosion resistance have often been advanced as the main reasons for the present widespread interest in titanium. For the physical metallurgist, the existence of an allotropic transformation in titanium is another important factor. As in the case of iron, an allotropic transformation may be the fundamental property upon which the response of Ti-base alloys to heat treatment will depend. It is therefore important to study the mechanism of the allotropic transformation in pure titanium and the effect of alloying elements on this transformation. The existence of an allotropic transformation in titanium was reported in 1931 by Schulze.' Later, De Boer et al.' determined the transformation temperature from a discontinuity in the electrical resistivity at 882"42O°C. According to recent work," the transformation temperature is 882.5"&1°C. The crystal structure of the room temperature a phase is hexagonal close-packed and the high temperature fi phase is body-centered cubic.' It is well known that the high temperature /3 phase of pure titanium cannot be retained by quenching to room temperature. Several investigators'-' have shown that in alloys of titanium with metals which are soluble in the @ phase (for example, nickel, molybdenum, chromium, or manganese), the /3 structure can be retained, providing the concentration of the alloying element is greater than a certain critical value. Below this critical value, the alloys transform during quenching, at least partially, into a supersaturated solid solution having the a. structure. This behavior is very similar to that found in some iron binary alloys (Fe-Ni and Fe-Mh, for example) and is an indication that the transformation might occur by means of a shear mechanism analogous to that observed in the austenite-martensite reaction. The present investigation is mostly concerned with the effect of the rate of cooling on the /3 to a transformation. In the first portion of this paper, the results obtained on pure titanium are described. Two other pure metals, zirconium and thallium, were also investigated, because they both have an allotropic transformation like titanium, from an hexagonal close-packed to a body-centered cubic structure. For purposes of comparison, the effect of rate of cooling of the y to a transformation of iron was included in the experimental program. The second part of this paper deals with Ti-Mo alloys. The equilibrium phase boundaries involving the solid state reactions on the Ti-rich side of the system were determined and a study was made of the effect of the rate of cooling on the transformation of the high temperature /3 solid solution. Experimental Procedure Since the main object of the present investigation was the determination of the effect of the rate of cooling on a transformation temperature, the best technique available was that used by Greninger" in his study of the martensite arrest in carbon steels. This method consists of heating a very small sample of metal to which fine thermocouple wires are .attached, cooling the sample rapidly by means of a blast of helium, and recording the temperature vs. time curve on an oscillograph. The apparatus used in the present experiments was very similar to that described by Greninger, except that the heating coil was made of molybdenum wire (0.040 in. in diam) and the assembly was placed in a bell jar under vacuum of 10" mm Hg or less. The specimen con-

Jan 1, 1952

By M. van Swaay, C. E. Birchenall

A. S. lBrling (Johnson, Matthey & Co. Ltd., Laboratories)— Because of its initial emphasis upon the production of membranes by powder metallurgy and by normal casting and rolling techniques, this paper is of considerable interest to those concerned with the design and application of industrial diffusion units. As the subsequent performances of the three grades of palladium were not compared in detail, it can perhaps be concluded that specific permeabilities were very similar. A short table in which the absolute diffusion rates obtained by the Authors with their three grades of palladium were compared with those of other investigators18-25 would, however, have been of particular value in view of the very variable results obtained with this membrane material. The Authors state early in the paper that the permeability of the thoriated samples decreased in much the same way as that of the nonthoriated samples. Later, they emphasize that sintered palladium without thoria was very susceptible to poisoning and did not respond as well as the other materials to the decontamination treatment. These results, which suggest that sintered membranes are less effective than those produced by normal methods, have considerable industrial implications, and amplified experimental details would be very welcome. The results of the Authors differ in one important respect from those of other investigators. Thus, no "bleeding" was required to maintain a constant rate of diffusion, although Davis, Darling, and De Rosset found this to be necessary. Davisz3 found that "bleeding" was necessary even when the hydrogen had already been purified by diffusion through a palladium tube, and it would be of value to know the purity of the gas stream emerging from the purification train shown on Fig. 1 of the Authors' paper. The Authors show that diffusion rates varied according to the square root of the pressure difference providing that the palladium membrane was not contaminated in any way. This finding is in agreement with the results of those workers2,18,19,23,24 who diffused their hydrogen through the membrane into a partly or completely evacuated system.Those workers who did not find this square-root law dependence did not have a partial vacuum on the downstream-side of the membrane.'' The difference between the two experimental con- ditions appears to be significant, as the absolute diffusion rates per sq cm of diffusion area are lower when the square-root law dependence is not observed. While the effect could be attributable to poisoning of the membrane it is conceivable that a positive pressure on the downstream side of the membrane might influence the mechanism of diffusion. In view of the capricious behavior of palladium, factors of this type, which are likely to influence its industrial significance, should be carefully considered. M. van Swaay and C. E. Birchenall (authors' 9-eply) — Our Eq. 191 may be compared with several others for hydrogen diffusion in palladium: D = 5.95 x exp (-5770/1)' D = 3.55 x 10"3 exp (-5550/)' D = 4.3 x 10"3 exp (-5620/1')' Other values are given in the literature, but in many cases surface poisoning may have affected the results. No analysis was performed on the hydrogen gas stream after purification. It should be noted that the presence of liquid nitrogen traps near the heated membranes may have set up convection currents and provided the equivalent of "bleeding." The question of whether a partial vacuum existed on the downstream side of the membrane is not clear. Does it imply that there is something unique about the permeability measured near one atmosphere pressure? We do not believe that there is. When square root pressure dependence is not observed, diffusion rates are not measured. Is "absolute diffusion rate per sq. cm. of diffusion area" equivalent to permeability? The terminology is new to us. In any case we doubt very much that pressures of the order of several atmospheres can affect the mechanism of diffusion.

Jan 1, 1962

By Herbert H. Kellogg

Graphs representing the standard free-energy of formation as a function of temperature for 21 fluorides are presented, along with estimated values for the standard free-energy of formation of 20 atedadditional fluorides. A few of the many possible uses of these data in metallurgical calculations are discussed, including the fluorination of oxides, sulphides, metallurgicaland chlorides, and the reduction of metal fluorides. RECENT papers which have presented free-energy data for oxides, sulphides, chlorides, carbonates, sulphates, and silicates in the form of free energy vs. temperature plots,1-5 have found important uses in teaching and metallurgical development work. This method of presentation of thermo-dynamic data serves to make available in readily usable form thermodynamic constants which are scattered throughout the literature, and to show by graphical means the stability relations which exist between the compounds of a given type. It is certainly desirable to have these free-energy-temperature diagrams for as many types of compounds as are of present or potential use in metallurgical processes. Thus, even though there are but few uses of metallic fluorides in extractive metallurgy today, the potential uses of fluorine processes are of great importance and serve as the principal justification of this paper. Recent developments in the production of fluoro-carbons have resulted in vastly improved methods for the production and handling of fluorine and corrosive fluorides, with consequent reduction in costs.".' The extractive metallurgist should no longer overlook possible fluorine processes when development of a new process for a relatively high unit-price metal is being considered. In a previous paperv he author outlined in some detail the importance and use of free-energy data in metallurgical calculations and also discussed the properties of a standard free-energy vs. temperature diagram. The same method of presentation is adopted in this paper, and the reader is referred to the earlier work for explanations of the use of free-energy data. Thermodynamic Data This paper presents standard free-energy vs. temperature diagrams, Figs. 1 and 2, for those fluorides for which the major thermodynamic constants are known. In the calculations of the standard free-energy of formation (AGO), values of the enthalpy of formation at 298°K (?Ho,) were taken from "Selected Values of Thermodynamic Constants"8 most of the values of standard entropy at 298°K (S°298) were taken from "Selected Values"," but in a few cases the values given by Kelley9 ere preferred; data for the Phase changes came from "selected Values" and several of Kellley's publics- tions;0-12 and most of the specific-heat data were taken from Kelley.l0 The appendix gives a detailed list of the sources of the data. As explained in an earlier paper,6 he accuracy of some of the standard free-energy relationships is not high. An error of +1 kcal at 1000°C is expected in most of the curves of Figs. 1 and 2. The error may be as high as +3 kcal in some of the curves. Rather than attempt to estimate the reliability of each curve, the source of the data used is reported in the appendix, and the reader may compare these data with newer values as they become available. Figs. 1 and 2 fail to show values for the fluorides of many important metals. In most cases this is because the entropy of the metallic fluoride is unknown. For preliminary calculations it is often useful to have available the best estimate of an unknown datum, rather than no information at all. With this in mind, Table I lists estimated standard free-energy values for a number of important metallic fluorides not included in Figs. 1 and 2. The values reported in Table I were, calculated by Brewer et al.13 and are based on estimated entropy values. The estimated uncertainty of the values for each fluoride were calculated from the uncertainty values given by Brewer1' for ?H°, and the entropy function.

Jan 1, 1952

By C. H. Gorski

A method for preparing titanium tetrachloride is described which consists of reducing rutile with coke and chlorinating the reduced product at 200° to 500°C. The crude distillate is purified by treatment with copper powder and distillation. Recoveries, exceeding 90 pct of the titanium in the rutile, were obtained. IN recent years considerable interest in metallic titanium has been aroused because of its outstanding physical and chemical properties. The present commercial processes for preparing titanium are modifications of the Kroll1 process and consist essentially of the reduction of titanium tetrachloride with magnesium in an inert atmosphere to produce metallic titanium with magnesium chloride as a byproduct. At present,' commercial titanium tetrachloride costs $0.32 a pound, which makes the price of the included titanium $1.28 a pound. As the titanium tetrachloride must be purified further before processing into metal, its ultimate cost obviously will be higher. Consequently, one of the prime requisites for the production of lower cost titanium metal is a source of low-cost titanium tetrachloride. Recently, Knickerbocker, Gorski, Kenworthy, and Starliper3 escribed a method for preparing titanium tetrachloride from low-grade Arkansas rutile by smelting the rutile with pyrite and coke to form a titanium matte and then chlorinating the titanium matte to prepare a crude titanium tetrachloride. This process had several inherent advantages. Low-grade rutile concentrates were utilized, the chlorination temperature was relatively low (approximately 200°C), and the reaction was strongly exothermic. However, this process had several disadvantages, chief of them being the rather low recoveries of titanium tetrachloride (65 to 81 pct) and the difficulty in purifying the crude titanium tetrachloride obtained. Because of the high-titanium and low-iron con- tents of rutile as compared to ilmenite, rutile should be of greater potential value as a source of titanium tetrachloride than ilmenite. Rutile, however, is not chlorinated directly in commercial practice, because reaction temperatures above 700 °C are required, at which point chlorine attacks most of the common materials of construction. The process now in use consists of converting the titanium in rutile to the cyanonitride by high-temperature reduction in an electric furnace, and then chlorinating this compound. The disadvantages of this process are the difficulty of producing the cyanonitride, the use of high-grade rutile concentrates in the charge, and the presence of several percent of silicon tetrachloride in the finished product. It would be of considerable interest to develop a new process for the preparation of titanium tetrachloride from rutile by chlorination at 200" to 500°C. Process Development Several investigators'-' have studied the reduction of titanium dioxide with carbon in an effort to prepare titanium monoxide. Invariably, the titanium monoxide formed in this manner was not pure, as it contained titanium carbide and one or more of the lower oxides of titanium. As these lower oxides of titanium are more reactive chemically than rutile, and as it is known that titanium carbide may be chlorinated readily, a study was made of the chlorination of carbon-reduced rutile. A sample of rutile concentrates from Hot Spring County, Ark., containing 76.19 pct TiO,, 5.65 pct Fe2O3, 4.18 pct SiO2, and 2.88 pct Al2O3, was used for this investigation. The experimental procedure consisted of making an intimate mixture of —100-mesh rutile and petroleum coke, adding water and molasses as a binder to foi-m a paste, charging into a graphite crucible that was covered, except for a vent to allow the

Jan 1, 1952

By R. L. Smith, I. G. Greenfield

The binary system Ga-Sb has been investigated by thermal, X-ray, and metallo-graphic methods. 'The intermetallic compound GaSb melts at 705.9OC and forms a eutectic with antimony at 11.8 atomic pct Ga. This eutectic melts at 589.8OC. A eutectic appears to be formed at 29.8°C but is too close to the composition of pure gallium to be detected. IN the course of an investigation of the semiconducting properties of the intermetallic compound, GaSb, the phase equilibrium between gallium and antimony has been investigated. Recently, the phase diagram of similar systems, such as In-Sb,' In-AS,' and Bi-As3 have been determined. These diagrams have the same general appearance in that an intermetallic compound is formed at 50 atomic pct and a eutectic is formed on both sides of the intermetallic compound. Jenny,& who has been investigating the semiconducting properties of intermetallic compounds, has concluded from thermal data of several alloys that the phase diagram of the Ga-Sb system resembles those mentioned previously. This paper represents a more detailed study of the Ga-Sb equilibrium. Experimental Procedure The gallium that was used to prepare the alloys was obtained from A. D. Mackay Inc. and the Aluminum Co. of America; the latter reported the purity as 99.95 pct. The antimony was obtained from Bradley Mining Corp. who specified a purity of 99.997 pct. Since the intermetallic compound of GaSb was prepared primarily to be used for further studies of the semiconducting properties, additional purification by a standard zone-melting technique under an atmosphere of hydrogen was performed. Hall-effect measurements on this material indicated an impurity content of less than 3 ppm.5 The antimony also was purified by zone melting and used in conjunction with the compound GaSb to make the necessary alloys for the determinations of the antimony-rich side of the phase diagram. The antimony-rich alloys were chill cast in Vycor tubes which were evacuated to a pressure of less than one micron. Since these alloys expand on freezing, precautions were taken to prevent oxidation of the ingot in the event the tube fractured. To insure homogeneity, the brittle alloys were ground and then placed in the furnace for the thermal-arrest determinations. Since the cost of gallium is high, the experimental procedure was complicated somewhat because of the necessarily limited sizes of the melts. In order to economize, the gallium-rich side of the diagram was determined by using a 35 g ingot of pure gallium initially and adding increasing amounts of antimony in order to obtain additional alloys. The furnace used in the determination of the cooling curves is illustrated in Fig. 1. A 2¼ in. quartz tube, sealed at one end, was fitted with a brass flange which was attached with DeKhotinsky cement. The brass collar was water cooled to prevent the cement from softening. An "0" ring was used to provide a seal between the brass collar and the furnace cover. Nichrome ribbon was wound around a 2 ½ in. asbestos-covered copper core which was used to reduce the thermal gradient in the furnace. The insulation consisted of alternate layers of asbestos and aluminum foil surrounded by sil-o-cel refractory and

Jan 1, 1956

By T. M. Geffen, J. L. Sanderlin, F. F. Craig, D. W. Moore

Scaled reservoir models have been used to study the effect of gravity on oil recovery performance in frontal-drive operations; namely, water, gas, or solvent flooding. The difference in density between the reservoir oil and the injected fluid causes their segregation, resulting in a non-uniform advance of the fluid front. In the laboratory flow tests, which simulated both five-spot and linear injection operations in flat reservoirs, the viscous, capillary, and gravity forces present in these operations were scaled. Dyed fluids were used so that the gross movement of the injected fluid could be observed. The studies covered a range of injection rates, formation thicknesses, and rock and fluid properties normally encountered in field operations. The results of the model tests indicate that the volume of the reservoir contacted by the injected fluid at its breakthrough into the producing well is less than that expected based on information which neglects gravity effects. This difference can often be as much as 80 per cent by gas or water injection in uniform sand bodies. Preliminary flow tests on a non-uniform sand body indicate that the uniformity of the flood fronts may in some situations be influenced to a much greater degree by permeability variations within the rock body than by gravity effects. The magnitude of fluid segregation due to gravity is controlled by the average injection rate, rather than day to day or week to week variations. INTRODUCTION One of the important factors controlling the oil recovery from a frontal-drive operation, such as water, gas, or solvent injection, is the volumetric sweep efficiency. This factor is a measure of the gross portion of the reservoir that is contacted by the displacing fluid. The volumetric sweep efficiency is influenced by gravity effects, well arrangements, and variations in rock permeability within the reservoir. The gravity effects are due to the displacing fluid being of different density than the reservoir oil. This causes the displacing fluid to move preferentially toward either the top or bottom of the formation. Non-uniform advance of the flood front can be caused by the relative positions of the injection and production wells. Variations in rock permeability also result in an uneven advance of the displacing fluid. Gas zones, representing paths of low resistance, may result in non-uniform movement of the injected fluid. The more uneven the advance of the displacing fluid, the lower is the volumetric sweep efficiency and thus the lower the oil recovery efficiency. The effects of well arrangement and permeability variations, upon the gross movement of the injected fluid have been the subjects of extensive investiga-ti0ns.1,2,3,4,5,6 Gravity effects have been recognized as being present in frontal-drive operations'.' but there has been little quantitative measure of their magnitude. Very likely the reason gravity effects have not been studied is that this problem does not lend itself to simple experimental or mathematical analysis. The development of reservoir modeling techniques has made possible a quantitative determination of the effect of gravity on fluid segregation and thus upon the oil recovery performance of frontal-drive operations. This paper presents a progress report on a laboratory investigation of gravity effects in fully liquid-saturated, horizontal, uniform as well as non-uniform systems. These results were obtained from scaled experiments on both linear and five-spot systems within the range of conditions normally encountered in secondary recovery operations. This study is concerned with the gross movement of the injected fluid to the time it first reaches the production wellbore. Therefore, it should be emphasized that the results of this paper can not be used to determine the effect of gravity segregation on oil recovery efficiency at abandonment conditions. MODEL SCALING The science of scaling models is relatively new as applied as petroleum production problems, Many of the scaling techniques have been adaptions of those used in basic studies of heat and fluid flow. There are

Jan 1, 1958

By J. D. Doherty

That America's great coal deposits eventually will be our principal source of liquid as well as solid fuels is generally accepted. Moreover, the day when synthetic oil from coal will begin to supplement our petroleum supplies does not appear to be far off. The American petroleum industry has not spared effort or expenditure in its attempt to keep pace with the demand for oil, and recent production has exceeded expectations. However, the rapidly rising demand for liquid fuels plus the increasing difficulty of finding and developing the new oil needed even to maintain present production rates make it highly unlikely that oil from domestic wells can fully meet future oil requirements. Therefore, it appears that one or more of the following courses must be followed: 1. Curtail oil consumption. 2. Incr ase imports. 3. Develop new sources such as synthetic liquid fuels. All three courses may well be needed, as it is doubtful that any one of the three alternatives alone can satisfactorily solve the problem indicated by the prevailing trend in demand for oil. Curtailment of oil consumption is not a constructive policy, except where oil is being wasted or coal can be readily substituted as in stationary steam plants. Reduced oil supplies for most uses will curb our economic advancement and, in addition to inconveniencing consumers, will have a pronounced adverse effect on American industry. The oil-burner-manufacturing industry already has been dealt a hard blow. Shipments in the first 5 months of 1948 decreased 70 pct from those in the same period in 1947.1 The effect of a similar cut in the automobile industry can be imagined. A severe prolonged shortage of gasoline could seriously disturb our present economic well-being. It might be pointed out also that even the prospect of a shortage or lack of assurance for the future can discourage buyers of equipment, as occurred in the case of oil burners. An increase in our imports of oil seems to be a necessary expedient for the immediate future. Its advantages are that it will require a minimum of manpower and steel in a time of critical shortages and will conserve domestic reserves. Any major dependence upon foreign sources for our oil supplies should be avoided, however. Other nations are increasing their oil consumption more rapidly percentagewise than the United States and eventually will compete with us for the not unlimited world supplies. In peacetime, foreign countries might restrict ship- ments of oil to this country or impose exhorbitant prices. In the event of war, there is serious danger that such sources would be cut off, and American-developed oil fields overseas might even be used against us. Although synthetic fuels cannot relieve any oil shortage that might develop within the next year or two, an immediate start toward a substantial synthetics industry to supplement petroleum not only could help relieve shortages several years hence but also would help establish the much needed "know-how" for large-scale production. By assuring ample future supplies, such a start would permit sound planning by industries dependent upon oil and maintain or restore the confidence of prospective customers for oil-con-suming equipment. If it is agreed that domestic production of crude oil cannot be expected to meet future domestic demands, that oil rationing in general is undesirable. and that too great dependence on foreign oil should be avoided—and the weight of the evidence supports such conclusions—it is inescapable that commercial synthetic fuel production should be developed to supplement domestic petroleum supplies, which, of course, will furnish the bulk of our requirements for yean to come. Sources of Synthetic Oil The three principal sources of synthetic oil are natural gas, oil shale, and coal. A commercial plant using natural gas is under construction, and expansion of this industry is expected.

Jan 1, 1950

By Russell G. Haworth

Three companies, United States Potash Company, Potash Company of America, and International Minerals and Chemical Corporation, are now operating potash mines and reGneries in the Carlsbad, New Mexico, area. The three mines are located approximately twenty miles east of Car1sbad. The deposits have been found ideal for mechanized mining methods. This summary has been compiled from information made available by the three operating companies. Prospecting Since outcrops of potash are nonexistent, prospecting is accomplished by drilling and coring of the salt section to establish the grade, thickness, and outline of the potash beds contained in the salt deposits. Information regarding the nature of the formations in the cover of 400 to 600 ft thick over the salt and potash beds is also obtained from drill cuttings in the same operation. By proper interpretation of assembled data the best location for shafts may be determined. The overburden varies in character but consists in general of interbedded shales, sand, gravels, limestones, and anhydrite beds. Water, in unpleasant quantities, especially in the sands and limestones, often completes the geologic section. Some of the shaftas in the area have been sunk without undue difficulties, others have been delayed by flows of water. Most of the shafts are concreted from the surface to the salt section. Each of the three mines in the Carlsbad area has one shaft for hoisting ore and one for hoisting men and supplies. Mining Methods Information acquired from initial drilling and prospecting affected the choice of mining methods. Since the tabular, slightly rolling deposits were similar to coal seams, the room-and-pillar system used in coal mining was adopted by all the companies. Pillars have not been removed in any of the mines on a large scale. One company has conducted experimental robbing of pillars in one restricted area. Where it is known that overlying strata contain water in large amounts, it is not considered practical to risk flooding of the mines by removing pillars without fill. Some thought is now being directed toward development of economical methods of filling for pillar recovery. The problem is not identical at the three mines and local conditions will govern the methods employed for eventual pillar recovery. The important point is that, with large reserves, it has not been necessary to resort to mining of pillars for adequate production. The ore deposits are laterally extensive but are irregular in shape and lie at a depth of 700 to 1100 ft below the surface. The major operations of all mines are on the same bed strati-graphically, except for the langbeinite level at International Minerals and Chemical Corporation. The deposits are apparently not connected and in many parts of the deposits being mined, "salt horses" or barren zones occur. From 400 to 700 ft of salt with minor beds of polyhalite and anhydrite lie immediately above the ore horizon. The ore differs from the salt only in the potassium chloride content and color due to small amounts of iron associated with the sylvite, giving the ore a somewhat mottled red and white cast. Thickness of the ore varies from 5 to 14 ft. Rooms and break-throughs in the mines vary from 24 to 40 ft in width, the overlying salt forming an excellent roof and providing safe mining conditions. Almost no timber is used in any of the mines for roof support. Drilling is normally the first step in mining the face of a room. Electric motor driven auger drills are used in the area. Jack hammers were used to some extent when the mines were first opened but the ore is easily drilled if bits are sharp and the coal type

Jan 1, 1950

By Albert S. Trube

The purpose of this paper is to clarify the definition of compressibility and to present a uniform basis upon which instantaneous compressibilities of liquids and gases can be compared. The equations gaverning the instantaneous compressibilities of imperfect gases are derived and the concept of pseudo-reduced compressibility is introduced. Part of the data presented by Brown, Katz et a1 on compressibility factors for natural gases has been rearranged. A graph of pseudo-reduced compressibility vs pseudo-reduced pressure for various pseudo-reduced temperatures is presented. The need for additional work in relating the compressibilities of liquids and gases is discussed. This information should be of value to reservoir engineers in making non-steady state performance calculations in gas reservoirs. It should be of further use irz pointing the direction for additional research in the nature of liquid and gas compressibilities. INTRODUCTION With the increasing use of steady and non-steady state well and reservoir data, there is a corresponding increase in the importance of the various factors entering into such calculations. Increasing emphasis is being placed on the necessity for obtaining reasonably accurate estimates of the physical properties of the reservoir fluids well in advance of the more accurate laboratory data. One such factor is the isothermal coefficient of expansion of the media which are transmitting and attenuating the non-steady state pressure waves. The average isothermal coefficient of expansion, or "compressibility" is a complex function controlled by the physical properties of the formation and the fluids contained therein. The isothermal expansion coefficients for reservoir gases are usually quite variable, in many cases being highly-pressure sensitive. The coefficients for reservoir liquids tend to be pressure sensitive, but not nearly so much as reservoir gases. The coefficients for solids, usually expressed in terms of a "modulus of elasticity" are relatively insensitive to pressure variations within their elastic limits. For this reason, and also because many previous applications have been limited to rel- atively small pressure ranges, there has been a tendency to ignore the variable nature of isothermal expansion coefficients and treat them as constants. Also, the term "compressibility" by which these coefficients are generally designated is commonly confused with a similar term, z, used to define the deviation of an imperfect natural gas from the perfect gas laws. A clear distinction should be made at the outset between the term "compressibility", which is an isothermal coefficient of expansion of a substance, and the term "compressibility factor", z, which refers to the deviation of a gas from the perfect gas laws. Although the scope of this paper is limited to the compressibility of single phase natural gases, it is definitely related to the problem of accurately estimating the compressibilities of single phase hydrocarbon reservoir liquids, which will form the basis of a future presentation. BASIC PRINCIPLES The coefficient of isothermal compressibility of a substance, c, is usually determined from pressure-volume or pressure-length -measurements depending upon whether the substance is single phase gas, liquid, or solid. A convenient method for making such estimates for a finite change in pressure and volume at constant temperature is to use the well known equation V1-V2/V1 (p2 - p1) .....(1) Eq. 1 is negative because the volume of a confined substance decreases as the pressure is increased. In this case V1 > V2 and p2 > p1. This equation is useful in approximating the compressibilities of single phase gases and liquids undergoing small pressure changes. It is evident, however, that this equation is almost identical with the determination of Young's modulus of elasticity for solids. If the assumption is made that change in length is directly proportional to change in volume, as would very nearly be the case for a steel rod in tension within its elastic limit, then E5=-L1 (p2 - p1)/L1 - L2 .......(2) in which E. is the isothermal expansion coefficient, or Young's modulus of elasticity, for a solid. And further, for this special case L1 (p2 - p1)/L1 - L2 .......(3)

Jan 1, 1958

By R. E. Smallman, I. Milne

The deformation behavior of single crystals of Nb-Mo alloys has been investigated with particular reference to the influence of composition, orientation, and temperature. Strong solid-solution hardening was observed reaching a maximum at the equiatomic cotrlposition and can be attributed to the difference in atomic size between niobium and molybdenutrz. Changes in the form of stress-strain curve, as shown by a high work-hardening rate and restricted elongation to fracture, were observed at a composition of Nb-85 pct Mo and are attributed to the presence of MozC DreciDitate. Conjugate slip was only extensive in dilute alloy samples; at the 50/50 composition deformation rnainly occurred by primary slip, and the onset of conjugate slip gave rise to failure by cleavage on (100). The variation of yield stress of Nb-50 pet Mo with orientation was consistent with slip on (011)(111) slip systems. The temperature deperndence of the yield stress between -196" and 250°C was similar to that of pure bcc metals, but at a much higher stress level; no evidence for twinning %as found. IN recent years the deformation behavior of various pure metals in groups VA and VIA has received considerable attention, but surprisingly little work has been carried out on binary alloys made by mixing metals from the two groups. Such an investigation would be of interest since single crystals of metals of group VA have been shown to deform characteristically with a multistage deformation curve1"3 while a parabolic type of deformation curve has been reported for most of the group VIA metals.4'5 It has been suggested by Law ley and Gaigher~ that the difficulty encountered in obtaining multistage deformation curves for molybdenum in group VIA was possibly because of the presence of a microprecipitate of MozC which they observed even at carbon contents as low as 11 ppm. Recently a multistage deformation curve has been reported for molybdenum ," although the stages are not so definitive as those for group VA metals. The binary alloys of the particular refractory metals which have been investigated in single-crystal form include Ta-w,' Ta- Mo,' and Nb- Na." While a large amount of hardening was observed for alloys of the Ta-W and Ta-Mo systems, associated with room-temperature brittleness for alloys approaching the equiatomic composition, Ta-Nb remained ductile over the complete composition range with little or no solution hardening. Other systems have been investigated by hardness measurements on polycrystalline material and a discussion of the hardening of these alloys has been presented by ~udman." The purpose of the present investigation was to examine the deformation behavior of Nb-Mo alloys in detail, with particular reference to alloy composition and single-crystal orientation. In this way it was hoped to shed some light upon the restricted ductility of these alloy specimens. 1) EXPERIMENTAL PROCEDURE The starting materials were obtained in the form of beam-melted niobium rod and sintered molybdenum rod of suitable dimensions. Since niobium and molybdenum form a complete solid-solution series at all temperatures, alloy single crystals were produced by melting the two constituents together in an electron bombardment furnace (EBM). To produce specimens free from segregation a molten zone was passed over the length of each rod six times in alternate directions at a speed of 10 in. per hr. Typical specimens were analyzed for interstitial impurities by gas analysis and for metallic impurities by spectrographic analysis. The results of this analysis are shown in Table I. Many of the tensile specimens were also analyzed (after testing) by scanning the gage length in an electron beam microanalyzer, from which it was found possible to predict the approximate composition of a specimen from the original proportions of each element in the EBM. The tensile specimens were made with a gage length of 0.5 in. and diameter of 0.075 in., using a Servomet Spark machine. By careful machining on the finest range for the final i hr of this technique, surface cracks could be reduced to the level where they were easily removed by electropolishing in a solution of nitric and hydrofluoric acids. The specimens were strained at a rate of 10 4 sec-' using friction grips designed to prevent accidental straining and maintain a good alignment before straining. The orientations of the individual specimens tested are shown in Fig. 1 and the corresponding compositions listed in Table I1 together with collated experimental data. 2)RESULTS a) General Deformation Behavior. The effect of composition on the room-temperature deformation curves of similarly oriented specimens is shown in Fig. 2. The yield stresses of the pure constituents, while not the lowest reported to date, were at least comparable with existing data. Although the solution hardening was large for alloys at either end of the phase diagram, and comparable with the Ta-W solution-hardening data of Ferris et a1.,8 the low work-hardening rate characteristic of niobium was sustained until a composition of Nb-85 pct MO had been reached. Associated with the peak yield stress ob-

Jan 1, 1969

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By C. E. Taylor, L. V. Fegan

With on-stream electrical-resistance measurement as a basis, Bethlehem Steel developed a system for the completely automatic control of moisture in the feed to pelletizing. This system was later successfully adapted to the sintering process, with automatic control of moisture serving to maintain the air permeability required for effective burnthrough. Control of moisture content is an important factor in the agglomeration processes by which iron ores are converted into either pellets or sinter for use in blast furnace burdens. In both pelletizing and sintering the iron ore feed must be mixed with the proper amount of water so that capillary attraction and surface tension forces can be utilized to cause the ore particles to coalesce and form larger agglomerates. The success or failure of the entire pelletizing or sintering process depends on maintaining the proper moisture. Therefore, continuous automatic control of moisture is one of the conditions for the automation and optimum use of both pelletizing and sintering processes. In pelletizing, variation in moisture can cause ball size to change and is a common cause of cycling which results in a fluctuating feed to induration. Excessive moisture can cause lumps to form at the cutter bar and contributes to the formation of buildup on the green pellet screens. Green pellets with excess moisture deform and reduce bed permeability. In sintering, accurate control of moisture is also important because moisture must be within the limits that give the bed the required permeability for good burnthrough. In response to these problems, a moisture meter, based on the principle of electrical resistance measurement, was developed. This moisture meter has been a key factor in the development of systems for providing uniform feed of concentrates in most of Bethlehem's pelletizing and sintering plants. The automatic moisture control system for pelletizing was successfully plant-tested and resulted in the following: (1) Control of moisture in the green pellets to +0.2%. (2) Up to 50% reduction of bentonite consumption. (3) Control of unscreened green pellets so that 85% of the furnace product at Cornwall was within the size range of -5/8 +3/8 in. (4) Improved quality of product with a tumble test of 94 and a "Q" index of 90 as the result of feeding a consistently uniform feed to a shaft furnace. In sintering, where sinter-mix temperature change and variable feed could cause errors in moisture measurement, the system was designed to maintain the correct moisture under virtually all plant operating conditions. It was successfully plant-tested and is now in use in most of Bethlehem's sintering plants. The system achieves: (1) Control of moisture within the desired limits for optimum production. (2) Continuous measurement and control of moisture even with intermittent feeding of mix, so that attention by the operator is seldom required. The two automatic control systems provide smoother startup, minimize upsets and contribute to the consistent control of process variables in pelletizing and sintering. MOISTURE-CONTROL SYSTEM FOR PELLETIZING Fig. 1 shows the general arrangement of the system for automatic control of moisture in pelletizing.1 In this system a bin level controller, using a single probe which regulates filter speed, maintains a uniform low concentrate level in the bin.' The moisture meter, comprised of an analyzer and control components, continuously measures the moisture in the bed and adds water to adjust to the desired level. The meter is placed in an open-loop (feed-forward) control, i.e., water addition occurs downstream from the probes. Filter vacuum must be set so as to produce a filtercake moisture which is below the optimum for pelletizing. Bentonite is added after the

Jan 1, 1969

By B. M. Larsen, T. E. Brower, J. W. Bain

A mass of circumstantial evidence is presented to indicate that the main source of alloy losses in open-hearth tapping is oxidation by air, with the steel apparently reacting with an amount of oxygen equivalent to about 30 times its own volume of air. The effect is erratic from heat to heat, depending largely on turbulence and distance of free fall of the stream of liquid metal. THIS paper reports another phase in a general study of oxidation in the open-hearth, various aspects of which have already been discussed in previous papers.1,2,3 It is based on experiments conducted on commercial, basic open-hearth furnaces at various times when opportunity offered, beginning about 1938, in an attempt to determine the mechanism and measure the amount of air oxidation during the flow of the tapping stream through the air into a ladle. This effect is closely connected, not only with alloy losses and control of chemical analysis, but possibly with various questions regarding inclusion content and steel quality and with a true picture of the deoxidation of steel in general. Thus it is believed that discussion of a large part of the results to date may be of assistance to others who are studying problems related to deoxidation or to steel cleanliness and quality. The literature reveals that strangely little attention has been devoted specifically to the effect of pouring steel through open air. Bardenheuer and Henke4 in 1939, in a paper devoted largely to overall losses of manganese in the basic open-hearth, reported that for a number of heats tapped into a tilted ladle held just beneath the runner so that the maximum distance of drop of the stream through the air was held to less than about 30 in., there was in most cases almost no loss of manganese from furnace to ladle. Hultgren,5 in 1945, showed that in high-carbon basic electric steel deoxidized in the furnace and containing no large inclusions as it left the furnace, large silicate inclusions, often rich in manganese oxide, were present in samples dipped from the pool in the ladle; also that the rise in content of large inclusions (above 0.01 mm dim) was greater when the metal stream was made to spray or flow in a turbulent manner. Our interest in this phase of oxidation in the open-hearth was first aroused by a study of alloy efficiencies on ladle additions. The frequency curves of manganese efficiency (net recovery of metal added in ladle) given in fig. 1 are of interest in this connection. Lowest, and also most variable (47 to 95 pct), recovery was in low-carbon heats having no coal additions to the ladle, progressing upward to a maximum average and least variable (82 to 97 pct) recovery in high-carbon heats killed with aluminum and silicon. Both coal additions and higher carbon contents in the tap stream, as well as additions of silicon and aluminum, help to lower the amount of manganese loss. But the most interesting point is the amount of manganese oxidized in many heats tapped at 0.50 to 1.06 pct carbon. It has been shown in a previous paper3 that the amount of oxygen in the tap stream in such heats is very small, usually below 0.01 pct, too small to oxidize any appreciable amount of manganese. Also, the fact that manganese should be oxidized at all in a metal solution containing such more stable oxide-forming elements as aluminum and silicon is an indication of some rather intense oxidizing effect such as could be caused by reaction with the air during ladle filling. Oxygen Balances from Tap Stream to Teeming Stream: Data in the form of alloy efficiencies, such as those shown in fig. 1, are not a proper measure of such a general oxidation effect. Having available

Jan 1, 1951

By L. L. Handy

The most common method of identifying hydrocarbon-bearing strata in a well that penetrates many different formations involves measurement and interpretation of the electrical properties of the formations as determined by electrical logs. Even though this method is used extensively, and even though in a great many instances it is capable of indicating presence of oil or gas, situations arise for which it is extremely difficult, if not impossible, to deduce the presence of hydrocarbons. These situations may involve the following. 1. A thin formation, bounded by highly resistive formations, in which it is impossible to obtain the actual resistivity of the uninvaded zone with existing logging devices. 2. A formation in which invasion has been so extensive that a value for the uninvaded zone resistivity cannot be obtained. 3. A very shaly formation in which the resistivity index, I, is lower than that usually associated with productive formations. 4. Laminated formations comprised of thin productive sands separated by thin shale streaks in which the individual sand and shale streaks are too thin to permit measurement of uninvaded-zone resistivity with existing logging devices. 5. Productive formations in which the water saturation is high. To extend the utility of electric log interpretation to identification of hydrocarbons in all types of formations, there is strong incentive to find a method not subject to these limitations. Some time ago, in connection with research on the wettability of reservoir rock, an investigation was conducted in which the resistivities of cores were measured shortly after they were removed from a core barrel, and again after they had been extracted and restored to their original oil and brine saturation.' The resistivities after extraction were generally lower. Other tests made on the cores indicated that they were more nearly water wet after they were extracted; thus, it was assumed that the observed changes in resistivities were due to a change in wettability of the cores. Other experiments'," have shown that resistivities of rock samples are sharply dependent on wettability. These experiments have shown that oil-wet samples are more resistive than water-wet samples. To obtain an understanding of how the wetting properties of the surfaces of core material affect electrical resistivity, a series of experiments was conducted. Two groups of core samples were prepared for testing. One group contained brine, but no residual oil. The other group was saturated with brine, flooded with oil to a low water saturation, then flooded with brine to a final residual oil saturation. Resistivity measurements were made on each group. Both groups were then flooded with the original brine to which a chemical had been added that renders sand and clay surfaces preferentially oil wet, a so-called reverse-wetting agent. very little change in resistivity was observed in cores containing only water. The group containing residual oil, however, showed resistivity increases of 100 to 200 per cent. These experiments showed that the resistivity of a core containing oil could be altered by changing wet-tability of the core. Moreover, the possibility was introduced that reverse-wetting agents might be employed as the basis for a logging method for identification of oil-bearing strata. Since behavior of a porous rock containing gas and water might be expected to be similar to that of a rock containing oil and water, such a method should also be applicable to identification of gas-bearing zones. In principle the wettability of the invaded zone could be reversed without altering conductivity of the interstitial water or the hydrocarbon saturation therein. Those strata showing significantly increased invaded-zone resistivities would, therefore, contain hydrocarbons; those with no significant change would be filled only with water. Addition of a reverse-wetting agent to a hydro carbon -bearing zone which is, by nature, already preferentially oil wet would not result in an enhancement of its resistivity. It is generally believed, however, that most hydrocarbon-bearing strata are preferentially water wet.