Search Documents

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By W. R. Hibbard

Evidence is presented which confirms previous findings that models of solution strengthening depending solely on lattice parameter changes are incomplete. Direct evidence for the Suzuki interaction of chemical effects was obtained from alloys with constant lattice parameter, for which solution strengthening was nearly independent of temperature from 300° to 20°K. SOLUTION hardening principles have been reviewed recently by Parker and Hazlett.' In addition, the interactions of dislocations and solute atoms were reviewed by Cottrell.' It is not proposed to review this information again here, but to refer only to those details which are pertinent to the developments which have occurred since these papers were published. Literature Survey To set the stage, it has been known for some time that solid solutions are stronger than the pure parent metals, and that the strength increases with increasing amounts of alloying element. One of the most extensive earlier investigations was that of Lacy and Gensamer,3 who reported that the solid solution strengthening of iron is a power function of the composition, that the relative strengthening is a function of the solubility limit, and that the rate of strain hardening is a function of the amount of solution hardening. Thus, the slope of the plastic portion of the stress-strain curve increases with increasing yield strength. An approach to the development of principles resulted from work by Gulyeav,4 who reported that the strengthening effect of solid solution alloys is periodic in character when plotted as a function of the atomic number of the solute. A more quantitative approach was to relate the amount of solution hardening to the lattice parameter of a solid solution. This was done quite early by Norbury,5 and later by Brick, Martin, and Angier,6 Frye and Hume-Rothery,7 and French and Hibbard. This relationship is only approximate in character. The next development was to relate the amount of solution hardening to lattice stress, as proposed by Frye and co-workers.9 They calculated what they called ionic overlap from the lattice parameter changes represented as a change in volume, using Bridgman's hydrostatic compression equation. The resulting correlation was somewhat improved as compared to the simple correlation with lattice parameter. The next step was carried out by Dorn, Pietrow-sky, and Tietz,10 who suggested an effect of valence differences in addition to the effect of lattice parameter. Their proposed parameter shows flow stress as a single valued function of the concentration and a relationship involving these two atomic variables. In addition, they contributed new data on the temperature dependence of solution hardening which indicated that it was a strong function of temperature for aluminum alloys. Allen, Schofield, and Tate11 carried the valency effect even further by reporting data indicating that with a given solvent element the stress-strain curve is a function of electron density, entirely neglecting the size factor, which varied by a factor of two. This evidence led to a quantum mechanical calculation by Cottrell, Hunter, and Nabarro,12 based on the formation of a dipole around an edge disloca-

Jan 1, 1959

By Sumner M. Anderson

Any discussion of outside support to Latin American mineral industries must concede at once the pre-eminent role of U. S. industry and business. American capital has developed the great copper resources of Chile and Peru, the lead-zinc of Peru and Mexico, the bauxite of Surinam and Jamaica, the petroleum and iron ore of Venezuela. American consumers have provided the market for their products. U. S. companies operating south of the border have trained Latin American geologists, engineers, chemists, and skilled labor—including metal, mineral, and fuel technicians—to new and greater competence. The industrial demands of World War II drastically altered world economics and concepts of U. S. foreign mineral requirements. In this country discovery of vital mineral resources fell far behind the alarming depletion rates, and there was new appreciation of the actual and potential importance of the world's less extensively explored areas. During the war and afterward, the task of supplying the nation's mineral requirements became greater than private industry could carry alone, and the Government assumed a position of continuing support. In Latin America, as elsewhere, this support applies in three categories: banking, procurement, and technical assistance. The first two are now generally accepted as necessary Government functions; the third has often been misunderstood and misinterpreted, although it has given rise to legitimate controversy. There have been major accomplishments in all three categories, as well as some errors of judgment. Since the beginning of World War 11, U. S. Government financing, procurement, and technical assistance have been applied on a large scale to a wide range of commodities and services throughout the free world. The present discussion—under the separate headings mentioned—will concern only the support related to a single industrial group in a single geographic area. During the War Federal banking was conducted on an emergency scale by temporary agencies created to meet the abnormal circumstances. Unfortunately, there are no readily available data on these activities that would permit evaluating loans and credits extended for mineral projects in Latin America. The only agency operating in this field since World War II is the Export-Import Bank, under its own authority and to some extent from credits authorized by the Office of Defense Mobilization. The Bank can offer more capital—on more favorable terms—than is available from most Latin American banks, and its staff of engineers evidently can appraise loan applications more favorably than commercial banks in this country are willing to do. The risks seem to be justified, as defaults have been well under 1 pct of total loans. From V-J day 1945 to the end of 1957, mineral loans to Latin America have been authorized for projects totaling $541.63 million in eight countries. In order of magnitude of credits, they have been made to Brazil, Peru, Chile, Mexico, Argentina, Venezuela, Guatemala, and Bolivia. In order of magnitude by industries. they are for development and expansion of steel mills, copper, manganese, nitrate plants, iron ore, zinc plants, sulfur plants, cement plants, tungsten, coal, and lead. Loan Distribution: By far the largest loans have been $307.9 million for iron and steel projects in Brazil, Chile, Mexico, Argentina, and Peru: $100.4 for developing the great low-grade copper deposits of southern Peru; and $55.3 million for developing manganese deposits in the Territory of Amapa, Brazil. Others range from $20.8 million for a zinc smelter in central Peru to $19 thousand for iron-mine trucking equipment in Chile. It has been noted that 59 pct of the loans have been to foreign companies, 14 pct to companies of mixed foreign and American capital, and 27 pct to American mining companies. About 40 pct of the total authorizations has been disbursed so far. In addition to these amounts, the Strategic Materials Division of the Economic Cooperation Admin-

Jan 1, 1959

By C. F. Weinaug, R. W. Farley, J. F. Wolfe

A rapid, accurate method for predicting the dew points of gas condensate systems and their subsequent normal and retrograde phase behavior with pressure decline has been developed. The method predicts dewpoint pressure, volume percent liquid, condensate content of produced gas, and the liquid and gas compositions at pressures below the dewpoint pressure. Use of a modified form of the BWR equation of state, which was established by regression calculations on experimental data from several gas condensates, provides depletion behavior predictions that are significantly improved over predictions using K-values based on convergence pressures and Standing's correlation for density. INTRODUCTION Prediction of complex hydrocarbon-mixture thermodynamics is fundamental to reservoir and gas process engineering. In particular, accurate descriptions of phase behavior are needed to plan production operations for gas-condensate reservoirs. Although analysis of reservoir behavior during producing and cycling operations can be made for simple reservoir models, extensions of these studies to more complex reservoir systems are still difficult, both with regard to numerical techniques and the phase behavior of the fluids. This report presents a mathematical method for handling the complex phase relationships encountered in retrograde condensates; it represents an advance in reservoir behavior technology in that it will aid the development of closed-form mathematical representations of total production systems. A typical phase diagram for a complex hydrocarbon mixture is shown in Fig. 1. When the reservoir temperature T, is higher than the critical temperature Tc, and lower than the critical-condensation temperature Tct' the mixture is a gas condensate. Decreasing the pressure on a condensate while maintaining nearly constant temperature will result in two phases forming at the dew point. Further decrease in pressure will cause additional liquids to be formed. This phenomenon is called retrograde condensation. As pressure continues to decrease, normal phase behavior regions are encountered (liquid vaporizes as the pressure decreases) until the low-pressure dew point is encountered. When a gas-condensate reservoir is produced, gas is withdrawn and the volume of fluids in the reservoir remains nearly constant as the pressure declines. Heavier components tend to accumulate in the liquid phase as a result of retrograde condensation, while lighter components tend to concentrate in the vapor phase. Thus, as the gas is produced the composition of the gas and the condensate

Jan 1, 1970

By K. Leutwyler

The key to realistic casing stress analysis in thermal recovery installations is accurate knowledge of the temperatures involved. Much information leading to prediction of heat losses from tubing string to wellbore has been published. This .same information permits calculation of casing temperature under certain conditions. This paper reviews the basis for most of the published work, the line source solution to the diflusivity equation and examines some of its application criteria in the light of unsteady-state conditions. An alternate finite differences method is presented for correlation, and experimental data pertaining to the heat transfer mechanism in the annulus is reported. In the light of the information presented in the paper, it appears that the assumption of a line source for the tubing-wellbore system is valid even for reasonably short periods of heat flow, .such as "steam soak" operations. The empirically modified equation presented earlier matches computer derived results of the finite differences approach lor periods as short as 50 hours. Previously predicted cas ing temperatures and the associated stress levels are thus verified based on these correlative computations and the reported experimental information. Fall out from the computer output resulting in some studies of the influence of insulating cement properties on casing temperature, and experimentally studied effects of nitrogen placed in the annulus at various pressures., are described. INTRODUCTION Generally speaking, predictions of average casing temperatures have been based on an idealized model involving a centralized tubing string at uniform constant temperature radiating energy towards the casing under steady-state conditions.'. ' Heat is transferred away from the casing to the formation by unsteady-state conduction across thermal barriers consisting of cement and the mud sheath. Heat balances written on the casing permit prediction of instantaneous casing temperatures. During early periods of injection, casing temperature varies appreciably as a function of time; however, for long injection times this variation decays as heat flow from the wellbore to the formation stabilizes, and it becomes virtually negligible as conditions surrounding the wellbore become quasi-steady-state (Fig. 1). Carslaw and Jaeger -escribe a solution to the diffusivity equation for a line source (Fig. 2). A constant supply of heat is introduced into the medium from the line placed into the origin of the x-y coordinate system and parallel to the z-axis. Heat spreads outward in the infinite cylindrical solid and the temperature increase on a line going through point (x, y) parallel to the z-axis at a distance r can be predicted by: The following criterion must be considered: the Carslaw and Jaeger equation is obtained on the basis of a large

Jan 1, 1967

By N. A. D. Parlee, I. D. Shah, W. C. Phelps

THIS note reports on a continuation of a program on the thermodynamics of Fe-S-X systems in the "steel-burning" temperature range begun by Bock, Barloga, and parlee.' Small coils of 0.039-in.-diameter pure iron (Ferrovac-E) wire were equilibrated with Fe-S-0 melts in ingot iron crucibles under argon at 1324°C in the assembly shown in Fig. 1, the lower part of which was heated in a large Globar furnace. Fig. 1 shows (R) the coiled iron wire sample, (M) the Fe-S-0 melt prepared from desired proportions of ferrous sulfide and ferric oxide, (K) the iron crucible, (I) a mullite guide tube, (J) an alumina crucible, (L) alundum sand, (N) crushed mullite, (Q) hollow cylinders of K30 brick, (C) a steel rod, (B), (El, (F), (G), (P), and (O), details of brass water-cooled head, and (H) a mullite tube 2-3/4 in. ID and 26 in. long. The rod (C) bearing the sample could be removed at will and be replaced by a tube for sampling the melt, or by a thermocouple for comparing the melt temperature with a control thermocouple not far from the 2-3/4-in. mullite tube. After equilibration the iron coils were treated with acid to remove the outer layer of metal and then they were analyzed for sulfur by a standard combustion method. The melt samples were analyzed for oxygen by a variation of Lander's2 hydrogen-reduction method for mattes, and for sulfur by the ASTM combustion-iodometric method NO. E 30-47. When the results of some forty-two runs were examined it was found that the sulfur and oxygen contents of the equilibrated melts showed good consistency, in that the plotted compositions fell along the smooth curve shown in Fig. 2. The points on this curve represent the experimental compositions of melts in equilibrium with ?-iron (Fe-S-0) solid solutions. These solid solutions are saturated with sulfur and oxygen with respect to the melt and therefore have the same sulfur and oxygen activities as the melt. The curve drawn through the points therefore delineates the ? liquidus. It falls in such a position as to roughly confirm the liquidus curves of Hilty and crafts3 at 1200" and 1300°C which are

Jan 1, 1965

By G. W. Thomas

Jan 1, 1967

By D. F. Hewett

CONTENTS PAGE Introduction 1 Geologic features of eastern Oregon 2 Geology of the Sumpter quadrangle 2 Rocks older than the ore deposits 3 Rocks younger than the ore deposits 8 Lodes of the Sumpter quadrangle 8 Areal distribution 8 Geologic distribution 9 Major features of the veins 9 Mineralogy and texture of the veins 13 Ore deposits zonally disposed around the Bald Mountain batholith . . . . 13 Summary of mineralogy and texture 32 Zonal relations 33 Ore deposits zonally disposed around Greenhorn Mountain batholith . . . 36 Zonal relations 37 Relation of textural features to zonal arrangement of veins 37 Summary of zonal relations 38 Comparison with similar districts 39 Outlook for the district 40 INTRODUCTION IN an earlier report on the Sumpter quadrangle1 it was briefly stated that the quartz lodes of that region might be separated into several groups or belts on the basis of the kind and abundance of the minerals of which they are composed, and that these belts were disposed as areal zones around the principal intrusive mass of the region. When the opportunity to revisit the region was offered in 1929, as a result of a cooperative agreement between the State Mining Board of Oregon and the United States Geological Survey, it seemed advisable to pay particular attention to this feature of the lodes, in the hope that if the idea of zonal relations was confirmed it might serve as a guide in exploitation.

Jan 1, 1931

By A. G. Rozner

TRANSVERSE rupture tests have been commonly used in mechanical investigations of brittle materials. The specimens are simple, easy to prepare, and loading presents no difficulty. Owing to the complexity of stress distribution in bending1 one can determine stresses rigorously only for the elastic range. The bend-test results, however, are at times converted into tensile-strength data, determined as maximum fiber stress. The assumption is made that for small plastic deformations the linear stress-strain relation applies. This problem was experimentally investigated by Bach3 who tested cast iron in bending and in tension and found that the maximum fiber strength obtained in bending exceeded the experimental tensile strength by a factor depending on the geometry of the tested piece. The factor varied from 1.57 for a rectangular cross section to 2.12 for a circular cross section of the tested beam. Since this work was carried out 60 years ago, the cast iron used could be expected to have been highly brittle, thus the extent of plastic deformation prior to fracture very small. We have therefore carried out a number of tests on magnesium-modified cast iron with an average elongation of -5 pct. It was hoped that the data obtained from the tests would serve as a guideline for comparing data obtained in tension and bending for brittle materials. Mechanical properties were measured at room temperature on an Instron tensile machine at a crosshead speed of 0.05 in. per min. The three-point bend test employed cylindrical bars 0.20 in. in diameter with a 1.0-in. span. For the tensile test, cylindrical tensile specimens with a gage length of 0.85 in. and a diameter of 0.160 in. were used. Results are summarized in Table I. From Table I one can see that fracture strength

Jan 1, 1965

By R. B. Booth, W. F. Linke

The continuous interest of the American Cyanamid Company in producing superior polymeric flocculants and dispersants for the mining industry has resulted in a broad, general study of the physical chemistry of mineral-polymer systems. The aim of this program has been to establish a sound theoretical basis on which practical problems could be analyzed, and which would lead to the development of new and improved flocculating and dispersing agents. This study has required the cooperative effort of physical, analytical, and polymer chemists, as well as engineers familiar with mining problems. This report summarizes the results of some of our investigations over the past two years. Further studies along these and other lines are in progress and will be reported in the near future. MEASURING FLOCCULATION IN THE LABORATORY The foremost problem in studying flocculation in the laboratory has been to establish reproducible methods of comparing and evaluating the effectiveness of various reagents. Previous workers have measured settling -rates 1,8 filtration rates2 , and sediment volumes3 . Direct observations of the floc size, or the clarity of the supernatant liquid, can also be used. AIthough these parameters are all interrelated, they do not always measure the same physical properties of the suspensions, and hence are not equally useful or reliable. Whatever the method of comparison selected. the testing procedure must be carefully defined and standardized. The results must be satisfactorily reproducible, and for practical purposes must accurately reflect the results obtained in mill use. In order to compare the data obtained in these various types of test, the physical nature of the flocculated system must be examined in some detail. In an unflocculated system, only the larger or denser particles (the "sands") settle at a useful rate. The sediment can be observed to build up from the bottom of the cylinder, while the fines remain in suspension for a long period of time. If the system is not disturbed, even the microscopic fines can be seen to classify themselves as they fall over a period of days. When a high molecular weight polymer is added to the system, the colloidal particles are tied together randomly in flocs, and the sands and fines settle together. The line of demarkation between solid and liquid phases moves downward from the top of the cylinder. The settling rate of the flocs depends on their size and on their effective density. For a given mineral sample, the larger the flocs, the faster the settling. Since the large flocs cannot pack as closely as small ones, they produce a sediment with more pore spaces, and hence of larger volume, which is easier to pump and to agitate. Although it is nearly impossible to measure individual floc sizes, it is relatively easy to determine settling rates and sediment volumes. Of these two, the settling rate is simpIer to measure, and hence generally more practical. Reproducible sediment volumes can only be obtained by allowing the solids to pack for a rather long period of time (e.g., hours to days), and although the volumes thus observed are directly related to the degree of flocculation (average floc size), it is often impractical to wait so long for the results.

Jan 1, 1961

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967

By A. S. Odeh

During the last few years, several authors1-6 have advanced mathematical solutions, both exact and approximate, to the pressure behavior of a well producing from a region bounded by a circular discontinuity. The most recent solutions were reported by Carter, and by Bixel and van Poollen. Varter solved the appropriate equations analytically. Bixel and van Poollen reported a numerical solution and provided correlation charts for the interpretation of flow tests; their method uses an overlay technique and is valid for wide ranges of mobility and storage capacity ratios. The mobility ratio is defined as the ratio between k/p of the two zones, and the capacity ratio is defined as the ratio between $c of the two zones. We have developed independently an analytic solution, programmed for the CDC 1604, that agreed with Carter's.'1 The same differential equation, assumptions, boundary, and initial conditions were used. We ran 30 drawdowns and 30 buildups of a unit capacity ratio and of varied mobility ratios. We also varied the radial distances to the discontinuity, and the radii of drainage. We examined the results and obtained correlations that were simple to apply and useful in flow test interpretation. The results and the correlations reported here in graphical form allow the determination of the mobility of the two zones and the radial distance to the discontinuity. They eliminate the need to fit field data to a curve as is required in the overlay technique.6 Thus, the advantages of the present correlations are that they require less effort and they provide a unique answer. They also provide guidelines for selecting the parts of the flow test plot to be used in the interpretations. Moreover, the results obtained from the programmed analytic solution define the time intervals where the approximate solution reported by Hurstl and discussed by Larkin2 are applicable. Normal buildup and drawdown plots resulted in an early straight line, a long transition zone and a late (second) straight line. However, the second straight line will not occur if the drainage radius is not at least 10 times the radial distance to the discontinuity. The mobility of the imer and outer zones can be calculated from the slopes of the first and second straight lines, respectively, by the normal methods of the flow test analysis. We found that the radial distance to the discontinuity may be determined from the abscissa value for the point of intersection of the extensions of the two straight lines. Results similar to some of these were reported by Adams et aL7 The approximate solution compared favorably with the exact one during the first and second straight-line periods. However, during the transition zone, deviation occurred between the two solutions. Discussion of Figures Figs. 1 through 7 summarize the results of this investigation. On Figs. 3 through 7 there appears the equation of each respective plot. In this paper the mobility ratio M is defined as klp2/kzpl, where subscripts 1 and 2 refer respectively to the inner and outer zones, and dimensionless time = t, = 6.32 X In the investigation, M ranged from 0.25 to 50;

Jan 1, 1970

By E. T. Turkdogan

Jan 1, 1962

[ ]

Jan 1, 1938

By H. H. Dreissen, F. J. Fontein

The operating principle of the classifying cyclone has been described often but some special aspects of the cyclone receive further clarification here. In addition, several uses of the classifiers are summarized and made clear by means of practical examples. In coal and ore washeries, the washing cyclone is widely used for the separation by specific gravity of sizes up to about 12 mm. This application has lately been extended to the separation by specific gravity of particles up to about 50 mm. This use makes special demands on the construction and wear resistance of the cyclone. The sieve bend, a development of the Dutch State Mines which is less well known than the hydrocy-clone, is dealt with at some length. The sieve bend is an apparatus designed for wet screening. Attention is given to the variables influencing the performance of the sieve bend. Hydrocyclones are now so widely applied that it would be difficult to find an industry not, directly or indirectly, using them or knowing how to use them since publications on this subject abound.1-3,5,9-11,13-15 However, because of the large number of variables influencing the operation of classifying cyclones, it is not entirely possible to calculate the results of this type of cyclone. Therefore, a number of lesser known variables will be discussed here along with some special applications of the classifying cyclone. FACTORS INFLUENCING CYCLONE CLASSIFICATION Shape of Particles: The shape of particles may cause a much coarser classification than the one expected, considering the sieve analysis and the cyclone dimensions chosen. This was clearly shown by feeding a mixture of spheres and disks in water into a cyclone. The spheres were 0.36 mm diam and the disks, 3 mm diam and 0.03 mm thick. These dimensions were chosen in order that the disks and spheres might have the same rates of fall under gravity. Both spheres and disks had a specific gravity of 1.05. The results obtained with these particles in a cyclone are summarized in Table I. Test No. 2 shows that the catch of spheres was

Jan 1, 1963

By C. C. Harris, A. Jowett, S. K. Ghosh

An equation derived from first-order kinetics and accounting for return from froth to pulp has correlated continuous flotation data. The concentration of floatable material in the pulp can be altered either by changing the feed concentration or the feed rate. The latter does not affect the flotation rate constant, the froth to pulp return constant, and the concentration of non-floatable material in both the pulp and concentrate. Both the return constant and the concentration of non-floatable material in the pulp are directly proportional to the feed concentration, while the concentration of non-floatable material in the concentrate is related to the feed concentration. The investigation of the froth flotation process by continuous test methods, previously used by Brown and smith,1 Jowett and Safvi2 and Jowett and West,3 is being continued. Recent tests, primarily concerned with flotation cell design, have yielded kinetic data suitable for testing an equation proposed by Arbiter and Harris4 which postulates that a flotation process is defined by two rate constants: 1) the rate constant for mass transfer from pulp to froth. Previous workers have regarded this as the flotation rate constant. 2) the rate constant for mass transfer from froth to pulp, which has previously been ignored. THEORY Under continuous flow conditions at the steady state in Fig. 1 the equation developed4 is Symbols are defined in Table I. The air-liquid ratio terms4 have been omitted, since they are substantially constant for a given set of conditions. Also, for a system operating at the steady state Vf and V, remain constant and Eq. 1 can be written: Cp/Cc=b'/a,+Qc/a' [2]

Jan 1, 1963

By E. F. Johnson, V. O. Naumann, D. P. Bossler

A method is presented for calculating individual gas and oil or water and oil relative permeabilities from data obtained during a gas drive or a waterflood experiment performed on a linear porous body. The method has been tested and found both rapid and reliable for normal-sized core samples. INTRODUCTION Individual oil and gas or oil and water relative permeabilities are required for a number of reservoir engineering applications. Chief among these is the evaluation of oil displacement under conditions where gravitational effects are significant, such as a water drive or crestal gas injection in a steeply dipping oil reservoir. Numerous proposed methods of obtaining relative permeability data on reservoir core samples have been too tedious and time consuming for practical use, or have yielded questionable and sometimes inconsistent results. A method has been developed by which the individual relative permeability curves can be calculated from data collected during a displacement test. The method is based on sound theoretical considerations. Using this method, with a properly designed experimental procedure, relative permeability curves can be obtained using core samples of normal size (i.e., 2 to 3 in. in length and 1 to 2 in. in diameter) within a few days after receipt of the core. In a recent publication D. A. Efros1 describes an approach to the calculation of individual relative permeabilities that is based on the same theoretical considerations. We believe the approach described in the present paper is more adaptable to practical application than the method implied by Efros. In addition, comparisons with independently determined relative permeabilities are furnished to substantiate the reliability of the new method. DERIVATION Previously the theory of Buckley and Leverett as extended by Welge2 as been used to calculate the ratio of relative permeabilities. In the derivation which follows, this theory is further extended to permit calculation of the individual relative permeabilities. The theory assumes two conditions which must be achieved before the method is applicable. They are that the flow velocity be high enough to achieve what has been termed stabilized displacement,3 and that the flow velocity is constant at all cross sections of the linear porous body. In stabilized displacement the flowing pressure gradient is high compared with the capillary pressure difference between the flowing phases. The high pressure drop insures that the portion of the core in which capillary effects predominate will be compressed to a negligibly small fraction of the total pore space. The assumption of constant flow velocity at all

By M. Mortada

Nonsteady-state flow of slightly compressible liquids in porous media of non-uniform properties has been the subject of a number of recent studies. Most of these studies considered one-dimensional flow in systems with cylindrical geometry and with the properties of the porous medium varying with radial distance only. A. Houpeurt' studied the problem in its general form by examining the solutions of the differential equation, In Eq. 1, the permeability k and the porosity + are arbitrary functions of radial distance. Houpeurt concluded that, except for special forms of the functions k and the problem is not well suited for mathematical treatment, but lends itself to high-speed digital computation. Carslaw and Jaeger' discussed the solution of Eq. 1 for the case in which the permeability varies as a power of the radial distance. P. Albert3 studied analytically the pseudo steady-state solutions in systems of finite radial extent and non-uniform properties. He also obtained some numerical solutions for nonsteady-state flow using a high-speed computer. William Hurst4 examined the interference pressure drop due to a point sink located in an infinite system consisting of two permeability regions in series. In this case, the permeability changed at a radial distance r from the point sink. INTERFERENCE BETWEEN OIL FIELDS IN A NON-UNIFORM AQUIFER The interference pressure drop for oil fields located in an extensive aquifer of uniform properties is discussed in a number of places in the literature. This work considers the interference pressure drop for oil fields located in a non-uniform extensive aquifer comprising two regions of different properties as shown in Fig. 1. Region I of the aquifer extends from which is the radius of Oil Field A, to ar,,, where a>1. It has permeability k,, porosity +, and effective compressibility c,. Oil Field B is located in Region II which extends from ar, to infinity and has different properties from those of Region I, but the same formation thickness h. Region II has permeability k2, porosity and effective compressibility c2. The choice of this type of aquifer should not imply that it is of frequent occurrence. It is chosen primarily to illustrate the effects of non-uniform aquifer properties on oilfield interference. It is of interest to point out that a non-uniform radial aquifer in which both the permeability and the porosity are inversely proportional to the radial distance can be treated in a manner similar to a linear aquifer of uniform properties. The situation where only the permeability varies as the inverse of the radial distance is also amenable to analytical treatment. Both aquifers behave quite differently from a radial aquifer with uniform properties so far as the water influx and the interference pressure drop are concerned. METHOD OF SOLUTION The expression for the interference pressure drop in Oil Field B due to a constant rate of water influx in Oil Field A is developed in the Appendix. This is accomplished by the simultaneous solution of two partial differential equations describing the pressure in Regions I and 11. The two differential equations are solved subject to the following conditions: (1) at r,, the pressure gradient is constant; (2) at ar the pressure and

By F. W. Schaller

TETRAGONALITY is an important aspect of the models of transformation, strengthening, and tempering of ferrous martensite. When the c/a

Jan 1, 1963

By Douglas D. Donald

The New Jersey Zinc Co. successfully holed through a 2 1/2-mile haulage tunnel connecting its new Ivanhoe shaft with the Van Mater Shaft at Austinville, Va. This 8x 10-ft cross-section tunnel was driven from both ends and met at a point approximately a mile from one shaft and a mile and a half from the other. In order to make certain that the faces met accurately on line and grade it was necessary to refine conventional mine surveying practices. The mines that are serviced by these shafts are located in the southwestern part of Virginia approximately midway between Roanoke and Bristol at the west foot of the Blue Ridge. While the presence of lead-zinc sulfide ores has been known in this area for more than 200 years, it was in 1902 that New Jersey Zinc Co. acquired the Austinville property, which has been in constant operation since that time. About ten years ago the company undertook an extensive diamond drilling program at Ivanhoe several miles to the southwest. This program was successful in locating sufficient ore to justify the purchase of a number of the properties in this area and plans were made to develop the new orebodies. The discovery of this ore as well as substantial addi- tional reserves at Austinville resulted in a decision by the company to increase its daily milling capacity to 3000 tons, of which about 1000 tons would be supplied from the Ivanhoe mine. After careful consideration of all possible methods of transporting this volume of ore from Ivanhoe to the mill, it became evident that diesel locomotive haulage through a 2½-mile tunnel between the two mines would best suit prevailing conditions. Since the lowest operating level, known as the 1100 level, of each mine would be at approximately the same elevation, it was obvious that such a tunnel should be driven on this level. Not only would this tunnel facilitate feeding Ivanhoe ore into the underground crusher and skip loading installation at Austinville, but it would also open up a known orebody 5000 ft southwest of the Van Mater shaft. In addition it would provide a valuable diamond drilling base from which a large and a potentially favorable but untested block of ground might be explored. In order to hasten the completion of the tunnel, it was felt that it should be driven from both ends. ENGINEERING CONTROL PROCEDURES It was vital that the tunnel junction be made accurately on line and grade; thus it was felt that engineering control of an unusually high degree of reliability was required to prevent an accumulation

Jan 1, 1959