Search Documents

By T. K. Redden

Protective coatings based on the formation of a surface coating of nickel aluminide (NiAl) were applied to the nickel-base superalloys IN 100, SEL 15, and U-700. Coated specimens were exposed to an oxidizing environment at temperatures between 1600 and 2200 F for times up to 1000 hr. The oxidation resistance and stability of the coating were evaluated by weight gain measurements, metallographic examination, and X-ray diffraction study of surface oxides and coating. The composition of the coating and diffusion zone was determined by electron microprobe traverse of samples before and after high-temperature exposure. Intermediate phases formed in the coating and diffusion zone were identified by X-ray diffraction in situ and after electrolytic extraction. The outer coating was found to consist of the inter-metallic compound, NiAl, while the diffusion zone contained MC, M23C6 or M6C carbides, and a phase in a matrix of NiAl + Nidl. Oxidation resulted in formation of an A1203 n'ch scale containing some Tz02. Depletion of aluminum during oxidation resulted in degradation of the outer coating to Ni3Al and the nickel alloy matrix. Diffusion of aluminum into the base metal was found to be slight and did not influence coating life significantly. The o formed in the diffusion zone during coating decomposed during elevated-temperature exposure to form stable carbide phases characteristic of the base metal. Diffusion zone phase changes were found to have no effect on the life of the aluminide coating in the oxidizing envzron?nent. THE oxidation resistance of many high-strength nickel-base superalloys is inadequate for extended exposure at temperatures above about 1600°F. In addition, some applications for these materials require that they be exposed to environments containing sulfur compounds and sodium salts which can cause surface attalk known as sulfidation or hot corrosion. In order to provide the necessary corrosion resistance to the high-strength alloys, protective coatings based on an aluminizing process have been developed. These processes, usually based on a pack cementation technique, result in the formation of a NiAl-rich outer coating layer either during the coating process or by a subsequent diffusion treatment. The performance of the aluminide coatings is affected by interactions between the coating layer and the base metal both during the coating process and during subsequent exposure at elevated temperatures. Knowledge of these interactions is required to guide the development of coatings capable of longer life and improved reliability. Goward et al.' recently reported the metallurgical factors which influence coating per- formance on MAR-M200. The present work is concerned with correlating the interactions and performance of coating compositions on several representative materials. EXPERIMENTAL PROCEDURES Materials. Three cast nickel-base superalloys which are used for turbine buckets in air-breathing engines were studied: IN 100, U-700, and SEL 15. Their chemical compositions are given in Table I. The alloys were vacuum-induction-melted and cast to slabs approximately 0.3 in. thick from which rectangular specimens 0.25 by 0.5 by 1 in. were machined. Coating Procedures. The machined specimens were coated by CODEP processes which were developed at the author's laboratory. These are based on pack cementation in various media to deposit either aluminum or aluminum in combination with titanium. The coating process which deposits only aluminum is designated CODEP-C, while the CODEP-D process deposits titanium in combination with aluminum. The CODEP-D process was applied only to IN 100. Both CODEP processes are applied at 1950" or 2000°F for 4 hr without need for a subsequent diffusion treatment. An outer coating about 1 mil in thickness is produced by these processes. Test Procedures. Coated specimens were exposed to static oxidation for periods ranging from 24 to 1000 hr at temperatures of 1600" to 2200°F. Terminal weight gain measurements and visual examination were used to evaluate oxidation resistance. including oxide spalling and coating failure. Both as-coated and exposed specimens of each alloy were studied by metallographic examination, electron microprobe analysis (EMA), and X-ray diffraction analysis either of the exposed surfaces or of phases extracted from the coating and diffusion zone. RESULTS As-Coated Condition. The microstructures of as-coated conditions were generally similar, irrespective of base materials or the particular coating process. They are typified by IN 100 coated by CODEP-D as shown in Fig. 1. The predominately single-phase outer layer, area A, Fig. 1, was identified by X-ray diffraction as the intermetallic compound NiA1. The NiAl zone extended inward to the original base metal interface. The diffusion zone, area B, Fig. 1, included carbide phases, a lamellar phase oriented perpendicular to the base metal surface, and a matrix phase consisting of a mixture of NiAl and Ni3Al. The phases in the diffusion zone were electrolytically extracted using a 10 pct HCl in methanol solution at approximately 1.3 amp per sq cm. The extracted phases were found to be M6C, MC, or M=C6 carbides and o as shown in Table I1 for each of the alloys. The d spacings from a typical diffraction pattern are

Jan 1, 1969

By H. Hooykaas

A few years ago Marsal published a method of predicting oil recovery as a function of time for an edge-water-drive reservoir with several rows of wells.' The method is based solely on oil- and water-production data obtained during a short period for the first row of wells and the gross productions (past and future) of all wells. This paper presents a streamlined version of the original publication and experimental evidence showing the reliability of this version. For this purpose scaled experirnents were performed in a large model (750 X 500 X 23 cm, 24.4 ft X 16.2 ft X 8.9 in.) in which edge-water-drive reservoirs with three rows of production wells were simulated. From the results of the experiments it is concluded that the method gives a fairly accurate prediction of the future production behavior of the reservoir and enables us to select the most profitable production plan. Results obtained with the method were always slightly conservative. INTRODUCTION The calculation method developed by Marsal aims at forecasting future reservoir behavior.' It can be applied to floods in which a number of down-dip wells are rapidly going to water. Essentially it is a method for predicting the future behavior of the up-structure wells on the basis of the observed behavior of the down-structure wells. As such, the method is especially useful for predicting behavior of edge-water drives toward several rows of wells. The requirement to be met for justifying the approach is identity of production mechanism in all regions of the reservoir. The attraction of the method then lies in the fact that it applies even when we have no complete description of the prevailing production mechanism. The original publication was in the German language and for this reason, perhaps, did not attract widespread attention. In our opinion it deserves to do so because it does offer in a simple way results which are in fair agreement with observed data. The aim of the present paper is to bring the method to the attention of a wider audience. The method in its most simplified form is reviewed and experimental data in support of it are presented. PRINCIPLES To show Marsal's principle. we will first write the continuity equation for very simple conditions, linear flow (in the positive x-direction) in an incompressible porous body [cross-section A(x), porosity +(x)], saturated with incompressible fluids (oil saturation S). When the oil in this formation is displaced by water, injected at x = 0 (inlet face), we get where q(t) is the total gross displacement rate and f the fractional oil flow. The equation can be solved if we make the assumptions: (1 ) that the fractional oil flow (f) is a univalent function of saturation (S); and (2) that the function is the same throughout the porous medium. Then: f = f(S) Solving this equation with the method of characteristics, we obtain

Jan 1, 1966

By Samuel G. Lasky

In a recent discussion of the mineral resource position of the United States, published in the September 1946 issue of Mining and Metallurgy, I expressed my conviction that ore deposits of major magnitude remain to be found in this country and that they can be found without undue cost or risk. I stated with equally strong conviction, however, that if the search for such deposits is to be notably successful there must be, among other things, a reorientation of philosophy and mental habits on the part of many who would be connected with the search. I would now like to enlarge upon that statement. For such a purpose. we may consider that the search for concealed deposits involves four steps: (I) the selection of areas for further study; (2) the geologic study of those areas in order to verify their general favorableness and in order to choose specific patches worthy of being tested; (3) the decision whether exploration of the patches chosen should actually be undertaken; and (4) the exploration itself. The first two steps and the last rest upon the judgment and skill of geologists. The third step—the decision whether exploration should be undertaken—rests upon management. In my opinion, the philosophies or principles that have guided these steps in the past are now inadequate. Philosophy of Appraisal If concealed deposits are to be looked for, how are the places of their possible concealment to be recognized? Every known mineralized area still in the prospect stage has already been repeatedly examined and condemned. As for those deposits concealed under later rocks, there may be not even the clue of a near-by outcrop for engineer or geologist to consider. flow then, are we to know where to look? Geologic Probability In the early days of our mining industry, ore deposits were found mainly by "following the ore." Later, as such opportunities decreased, and as geologic concepts became accepted, deposits were found more and more through use of the notion of geologic probability. Many instances can be cited wherein a courageous operator ignored the implications of the unsuccessful efforts of predecessors and, solely by "following the ore," converted a discredited prospect into a good mine. Similar experience doubtless will occur in the future, but "following the ore" will no longer in itself find new deposits in the amount needed. Nor can enough new ore be expected from further use of the notion of geologic probability, for, as mentioned above, our known prospects have already been ex-

Jan 1, 1949

By Samuel G. Lasky

In a recent discussion of the mineral resource position of the United States, published in the September 1946 issue of Mining and Metallurgy, I expressed my conviction that ore deposits of major magnitude remain to be found in this country and that they can be found without undue cost or risk. I stated with equally strong conviction, however, that if the search for such deposits is to be notably successful there must be, among other things, a reorientation of philosophy and mental habits on the part of many who would be connected with the search. I would now like to enlarge upon that statement. For such a purpose. we may consider that the search for concealed deposits involves four steps: (I) the selection of areas for further study; (2) the geologic study of those areas in order to verify their general favorableness and in order to choose specific patches worthy of being tested; (3) the decision whether exploration of the patches chosen should actually be undertaken; and (4) the exploration itself. The first two steps and the last rest upon the judgment and skill of geologists. The third step—the decision whether exploration should be undertaken—rests upon management. In my opinion, the philosophies or principles that have guided these steps in the past are now inadequate. Philosophy of Appraisal If concealed deposits are to be looked for, how are the places of their possible concealment to be recognized? Every known mineralized area still in the prospect stage has already been repeatedly examined and condemned. As for those deposits concealed under later rocks, there may be not even the clue of a near-by outcrop for engineer or geologist to consider. flow then, are we to know where to look? Geologic Probability In the early days of our mining industry, ore deposits were found mainly by "following the ore." Later, as such opportunities decreased, and as geologic concepts became accepted, deposits were found more and more through use of the notion of geologic probability. Many instances can be cited wherein a courageous operator ignored the implications of the unsuccessful efforts of predecessors and, solely by "following the ore," converted a discredited prospect into a good mine. Similar experience doubtless will occur in the future, but "following the ore" will no longer in itself find new deposits in the amount needed. Nor can enough new ore be expected from further use of the notion of geologic probability, for, as mentioned above, our known prospects have already been ex-

Jan 1, 1949

By John P. Dolan, Charles A. Einarsen, GiIman A. Hill

his paper discusses how the following formation characteristics can be determined mathematically from drill-stem test pressure charts: true formation pressures, effective permeability of the entire section tested, well productivities, wellbore damage, and possible detection of barriers (faults, pinchouts, etc.). This paper also presents a practical method for im-mediate determination of eflective permeability and well-bore damage from successful double shut-in pressure tests. A list of recommendations for improving the reliability of drill-stem rest pressures is also presented. INTRODUCTION A drill-stem test is a temporary completion designed to sample the formation fluid and to establish the possibility of commercial production. Early pressure recording devices were used merely to verify proper operation of the testing tool. Until recently the accuracy of the pressure gauges has been insufficient for any reliable quantitative use of the recorded pressures. In view of the need for more reliable formation evaluation and as a result of the recent interest in exploration work employing the concept of hydrodynamic en-trapment,1,2,3,4 better pressure recording gauges are now in use. These devices can record pressure within 1 per cent above 1,000 psig and can detect differential pressures as low as 1/2 psig. In addition to formation pressure, several other reservoir characteristics can be determined from DST charts; namely, well productivity, formation permeability, wellbore damage, and the possible existence of bar- riers (faults, pinchouts, changes in permeability, etc.). This paper presents a practical method to interpret DST pressure charts for formation pressures and many other reservoir properties, a method that has been developed in analyzing approximately 4,000 DST charts during the last five years. The techniques used are a composite of published articles on drill-stem testing5,6, together with well-known pressure build-up analysis methods. THEORY It has been shown7,8 that the following equation may be used for analysis of pressure build-up curves: Po - Pw = 162.6 qu/kh log ( to + ?/?) .....(1) When this equation is applied to the curves obtained in drill-stem testing, the assumptions and boundary con-ditions are more nearly realized than in conventional flow and build-up tests on producing wells. Zak and Griffen9 have recently discussed in detail the use of this equation in analyzing DST charts. One of the problems with DST curves is the lack of reservoir data for precise analysis. Therefore, it is necessary to develop empirical rules and field methods for analyzing DST charts in quantity. For this reason, the empirical methods presented in this paper have been developed, and their derivation is found in Appendices A and B. METHOD USED FOR ACCURATE READING OF DRILL-STEM TEST CHARTS In order to apply the pressure build-up theory to DST charts, it is necessary to obtain a digitized expression for the pressure and time data recorded by the pressure gauge. These data may be either provided by the service company or interpolated by projecting a

Jan 1, 1958

Jan 1, 1935

By I. Fatt, D. N. Saraf

A method is described for measuring two- and three-phase relative permeabilities in sandstones or sand packs using a nuclear magnetic resonance (NMR) technique to determine fluid saturations. Two- and three-phase relative permeabilities have been determined on Boise sandstone using the NMR technique of saturation measurement. Three-phase relative permeability to water was found to depend only on the water saturation, whereas three-phase permeability to oil depended on both the water and oil saturations. Relative permeability to gas in three-phase flow was found to depend only on the total liquid saturation. INTRODUCTION Three-phase relative permeabilities are extremely useful in calculating field performance for reservoirs being produced by simultaneous water and gas drives. Three-phase relative permeability data are also needed for analyzing solution gas-drive reservoirs which are partially depleted and are being produced by water drive. Some thermal recovery processes involve three-phase flow which require three-phase relative permeability data for predicting reservoir behavior. Unfortunately, three-phase relative permeability measurements have rarely been made. Also, because of the scarcity of three-phase data, it has not been possible to date to relate other measured rock characteristics to the relative permeabilities with a great certainty. Leverett and Lewis,l Reid2 and Snell3 have reported three-phase relative permeability data on unconsolidated sands. Leverett and Lewis used ring electrodes spaced along the length of the sand sample to measure the resistivity of the sample which was assumed to be monotonically related to brine saturation. Gas saturation was determined from pressure-volume measurements. Oil saturation was obtained by material balance on the cell containing the sand sample. This method is involved and time consuming. Another difficulty arises from the fact that the resistivity of the sand is a function not only of saturation of brine but also of the distribution and saturation history of the brine in the pore spaces. Reid2 used a gamma ray absorption technique for measuring liquid saturation. This method has the disadvantage that total liquid saturation rather than oil or brine saturation is all that can be measured and still another method is required to determine the saturations of individual components. snel13 used a neutron bombardment method which also required a separate determination of the individual component saturations. Caudle et al.4 measured three-phase relative permeability on consolidated sandstones using vacuum distillation for determining fluid saturations. Distillation after each reading makes this technique very lengthy and time consuming. Corey et a1. 5 and Naar and Wygal 6 measured three-phase relative permeability on sandstones by the capillary-pressure method. Semipermeable diaphragm assemblies were used at each end of the core specimen to keep the water phase in the core. Gravimetric methods were used to determine fluid saturations. Sarem7 recently reported an unsteady-state method for measuring three-phase relative permeability on sandstones. This method is an extension of Welge's method8 for measuring two-phase relative permeability. Although Sarem's method is simple and comparatively fast, the assumptions involved may oversimplify the problem. Sarem's assumption, that in all rocks relative permeability to each fluid will depend only on the saturation of that fluid, seems to be rather unrealistic. Neglecting capillary effects at the end of the core is also a weak assumption. Donaldson and dean9 measured three-phase relative permeability using a method similar to Sarem's.

By J. P. Heller, H. B. Bradley, A. S. Odeh

A potentiometric model technique is presented for determining the areal sweep efficiency of a five-spot well pattern, at and beyond breakthrough. A sharp interface between displaced and displacing fluids is assumed. Although the prototype svstem is referred to as a waterflood operation, obvious changes in the notation and computations will adapt the results to other displacement processes. The results of the study include the following. 1. Areal sweep efficiencies for a five-spot well pattern, at and beyond breakthrough, for mobility ratios (displacing to displaced fluid) of 4:1, 2:1, 1:1 and 1:4. 2. Extension of potentiometric analysis to the investigation of the areal sweep efficiency beyond breakthrough in a five-spot pattern by the application of conformal mapping and conductive-solid models. 3. The use of layers of conductive fabric in representing mobility ratio changes in potentiometric models. 4. The development of a probe mechanism for probing conductive solids. The results obtained by conductive-cloth models agree with earlier areal sweep efficiencies at breakthrough obtained by Aronofsky and Ramey 1 on the potentiometric analyzer using electrolytic-tank models. Results beyond breakthrough differ from those obtained by the X-Ray Shadowgraph technique. Data from this study show that, for mobility ratios greater than one, water cut rises rapidly as fluid is produced after breakthrough. However, for mobility ratios smaller than one, a large increase in area swept resulted with only a small increase in water cut. INTRODUCTION In calculating reservoir performance of waterflood-ing operations and other fluid-injection programs, it is necesssary to estimate the areal sweep efficiency * before and after injected-fluid breakthrough into production wells. Influence of mobility ratio on oil-production history, before and after breakthrough for a five-spot well pattern, has been studied by X-Ray Shadowgraph techniques 2-5 and by gelatin models.6 In none of these investigations was the transition zone controlled experimentally. Steady-state anaiog techniques assume that a vertical and discrete interface exists between the displacing and displaced phase. Thus, when using a potentiometric analysis, it is assumed that no transition zone exists. To determine the areal sweep efficiency beyond breakthrough for a five-spot well pattern using the potentiometric analyzer, four principal problems had to be resolved. These were the following. 1. The laborious method of representing changes in mobility ratio in electrolytic tanks by means of contoured wooden blocks. 2. Vapcrization of the electrolyte solution during the study causing concentration changes and, thereby, variable conductivity. 3. The scale-limitation problem 7 that is accentuated by abrupt eiectrolyte depth changes. 4. A wellbore geometry problem. EXPERIhlENTAL EQUIPMENT AND PROCEDURES CONDUCTIVE-SOLID MODEL The theory, equipment and procedural techniques of the electrolytic-tank potentiometric model have been treaced excellently by Lee8 and by Muskat.9 By the use of a conductive solid to represent the fluid conductivity of the porous media, the need for contouring wooden blocks to represent mobility ratio was eliminated. It was found that a carbon-black coated fabric, Uskon D-16,** could be employed as a fairly uniform conductive media. This material, when stacked in appropriate layers and pressed together, can be used to represent changes in mo-

By J. M. Tobin, L. H. Cadoff, L. M. Adelsberg

The reaction between liquid zirconium and graphite at temperatures above 2000 °C has been investigated. The reaction products were found to be carbon-saturated zirconium metal and ZrC which formed between the graphite and the metal. Parabolic growth behavior was observed for the ZrC Phase at all temperatures of this investigation. The parabolic growth constant at temperatures between 2000° and 2860°C was measured to be 1.83 exp 84,300/RT) sq cm per sec. The reaclion mechanism was proposed to be the rapid carbon saturation of the liquid metal and the formation of ZrC at the metal-carbide interface with diffusion of carbon through the ZrC, the "rate - determining step" of the reaction. The concentration-independent diffusion coefficient of carbon in ZrC, (DZrC), was expressed as 0.95 exp? 78,700/RT) sq cm per sec. This value mas calculated using the temperature-invariant ZrC phase fields proposed in the literature. The [ZrQiq)] — [ZYQiq) + ZrC] phase boundary over the temperature range 2000° to 2800°C was determined and the ZrC + C eutectic temperature was found to be 2890° ± 50°C. ThE Group IV and VB transition-metal refractory carbides are of interest because of their high melting points, high temperature strength properties, and relative inertness in certain corrosive environments. The present-day understanding of these materials, however, is limited by the general unavailability of accurate and reliable physical and chemical property data. This is due primarily to the difficulties associated with the preparation of suitable, high-density, high-purity carbide samples, and the achievement and control of uniform high-temperature environments. Accurate measurement of temperature is also an important factor limiting the reliability of reported data. The Zr-C reaction was selected for investigation because of the high melting point (>34000C) and favorable nuclear properties of the reaction product, ZrC. The direct reaction method afforded an opportunity to obtain kinetic data on the fully dense carbide. In this paper, layer-growth techniques were used to estimate the diffusivity of carbon in ZrC and to investigate the phase equilibria in the Zr-C system at temperatures above 2000°C. EXPERIMENTAL PROCEDURES AND DATA Crystal bar zirconium (99.9 pct Zr), purchased from the Nuclear Materials and Engineering Corp., and ZTA graphite (99.9+ pct) were used in this study. The analyses of the materials are listed in Table I. A schematic of the high-temperature carburization apparatus is shown in Fig. 1. The sample was a zirconium metal charge in a 1/2 -in.-ID graphite crucible which was capped with a tightly fitting graphite stern. The crucible and stem were designed to closely approximate black-body conditions. The graphite crucible was packed in lampblack (outgassed 1 hr at 2100°C) which provided insulation and thermal stability to the system. The inert atmosphere was maintained by a constant flow of high-purity argon or helium gas. The nozzle-diffuser section on the top flange was sufficient to prevent back-diffusion of air into the system. Chemical analysis of the carburized samples revealed oxygen and nitrogen concentrations of less than 44 and 20 ppm, respectively. The sample was heated by induction with a Westing-house 5 kw-450 kc power source. Temperature was measured with a Milletron two-color pyrometer which was sighted into the crucible by reflection from an overhead front surface mirror. Optical losses due to glass absorption and light reflection at air-glass interfaces were minimized by the employment of the dynamic-gas seal. It was found that no correction was required for the front surface reflector. The pyrometer was calibrated periodically against a U.S. Bureau of Standards tungsten ribbon secondary standard. Temperatures inside the crucible were also calibrated against the Zr-ZrC eutectic (1850°C)1 and the Nb-Nb2C eutectic (2335°C)2,3 temperatures and agreed to within ±10°C of these values. The temperature was controlled by manually adjusting the power input; variations of

Jan 1, 1967

By D. H. Yardley

ENGINEERS and geologists find the magnetic compass unreliable in many mineral bearing areas because of local magnetic attraction. This is the case in the vicinity of many contact deposits, and particularly true in the Lake Superior district where iron formation is wide spread and causes wide fluctuations of the magnetic needle. One of the most useful instruments for geologic traversing or running grid and claim lines under these conditions is the dial compass which does not depend upon magnetism. Difficulties encountered in the field because of local attraction can be overcome through transit observation of Polaris or the sun, or use of the solar attachment. However, a field party often lacks the equipment or time for work of such accuracy and it is here that the dial compass finds its use.

Jan 3, 1954

By Jim Douglas, H. H. Rachford, D. W. Peaceman

The problem of unsteady-state gas flow through porous media has been solved numerically only for the case of linear or radially symmetric reservoirs. A recently introduced numerical method for solving the unsteady-state heat flow equation in two dimensions is applied to the calculation of the depletion of a square region containing a perfect gas. Solutions are presented in graphical form for various values of dimensionless parameters. The solutions are compared with published .solutions for radial reservoirs. I NTR ODUCTION The problem of unsteady-state flow of gas through porous media gives rise to a second-order non-linear partial differential equation for which no analytical solution has been found. Numerical approximations to solutions of the gas flow problem have been obtained by the stepwise solution of an associated difference equation1,2. However, the methods so far developed have required that the reservoir be either linear or radially symmetric. This restriction in shape has been necessary so that only two independent variables be considered, namely, one distance variable and time. In order to deal with reservoirs having more realistic shapes, it is necessary to develop numerical procedures for the solution of the gas flow problem involving two distance variables. A numerical procedure, denoted as the alternating-direction implicit method, for the solution of the heat-flow problem in two dimensions has recently been introduced" By the use of this procedure, approximate solutions have been obtained for heat-flow problems in a square8, and in regions having various non-rectangular boundariesa. Because of the similarity between the equation for heat conduction in solids and the equation for gas flow in porous media, it is reasonable to expect that the alternating-direction implicit method should also be useful for solving the gas flow problem in various two-dimensional regions. Solutions have been calculated for the simplest two-dimensional region, a square reservoir with a single well in the center. While reservoirs of such simple geometry seldom, if ever, exist, the solution of this problem is of some practical importance because it is also the solution of another problem, that of determining the depletion history about an individual well in an infinite uniform reservoir containing wells spaced throughout on a square lattice and producing equally from each well. Finally, it is of interest to compare solutions for a square reservoir with those for a circular reservoir and to determine the effect of the shape of the boundary for that particular case. METHOD OF CALCULATION Basic Differential Equation By combining the equation of continuity, the perfect gas law,

Jan 1, 1956

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.

By B. E. Eakin, A. L. Lee

Atmospheric pressure viscosity values For light hydrocarbons were calculated by a simple mixing rule based on the Sutherland constants for pure components. This method was further extended so that the atmospheric pressure viscosity of a mixture may be calculated based on the molecular weight information only. INTRODUCTION A proven method for estimating the gas-phase viscosity of a mixture is the one by Herning and zipperer.l This mixing rule has the following form: One method for calculating the pure component viscosity i is the Sutherland equation,2 which has the following form: The Sutherland constants B and S are determined from experimental data. A set of Sutherland constants for light hydrocarbons based on the data of several investigators3-5 is presented in Table 1. APPLICATION The Herning-Zipperer equation and the Sutherland equation were used to obtain an equation for calculating the gas-phase viscosity of light hydrocarbons. This equation has the form, Eq. 3 was used successfully in the study of methane-n-butane mixtures.6 Further analysis indicated that the Sutherland constants are a function of molecular weight.7 The gas-phase viscosity may be calculated based on molecular weight information only. The following equation was developed: The atmospheric pressure viscosity of ethane, propane, n-butaneand four methane-propane mixtures calculated on the basis of Eq. 6 are compared with the data of other investigators in Tables 2, 3, 4 and 5, respectively. A plot of atmospheric pressure viscosity vs molecular weight based on Eq. 6 is presented in Fig. 1.

Jan 1, 1965

By R. A. Greenkorn, R. C. Smith

Hele-Shaw cells are used to model creeping flow through porous media (where Darcy&apos;s law is valid). The effects of inertia on flow about obstructions in a Hele-Shaw cell can be calculated by a perturbation method if one can determine a solution to Laplace&apos;s equation. Results of a computer solution for flow about circular, square and elliptical obstructions are presented. These results show that or a modified Reynolds number of less than I, the inertia terms are small; and for values of less than 3, the average streamline predicts the ideal flow. Therefore, the analogy might be used for studying flow in porous media up to a modified Reynolds number of at least 3. INTRODUCTION The nature of fluid flow in porous media is of interest in the fields of soil mechanics, ground water flow, petroleum production, filtration and flow in packed beds. Because it is very difficult to study the phenomenological behavior of flow in porous media, homologs and analogs are used to study flow characteristics. A Hele-Shaw model, made of two closely spaced plates - usually glass — is often used as an analogy to two-dimensional flow in porous media. Hele-Shaw2 showed experimentally that the streamline configuration for creeping flow around an obstacle located between two closely spaced parallel plates is the same as for two-dimensional, ideal flow about the same obstacle. Stokes7 verified these observations mathematically. The usual equation of motion for flow in porous media is Darcy&apos;s law. The form of the mathematical statement of Darcy&apos;s law is identical, within a multiplicative constant, to the expression for the average velocity over the plate gap in the plane of a Hele-Shaw model. These models may be used to describe flow in both homogeneous and heterogeneous porous media. In the mathematical proof of the Hele-Shaw analogy it is assumed that the convective terms in the Navier-Stokes equations are negligible and that the equations of motion degenerate to Laplace&apos;s equation, with pressure the dependent variable. Whenever a Hele-Shaw model is used as an analogy to flow in porous media, .the validity of this assumption is in question. Riegels5 showed that if convection is not neglected, the velocity distribution around a cylindrical obstruction in the flow field depends on a Reynolds number, the plate spacing, and a dimension characteristic of the obstacle. Riegels&apos; solution, a perturbation solution, uses the boundary condition that the flow rate into the obstacle averaged over the plate gap at any point on the obstacle is zero. The method requires that a solution to Poisson&apos;s equation for the perturbation pressure be found. Riegels evaluated this solution for the case of the cylindrical obstruction. His method may be simplified by eliminating the need for solving Poisson&apos;s equation for the perturbation pressure. Instead, an analytic expression for the perturbation pressure gradient is obtained (valid for arbitrary shapes) and used to eliminate pressure from the equations for the perturbation velocities. The results show, for symmetrical shapes, that if NRe < 1, the convective acceleration terms are small, and that the average velocities represent ideal flow up to at least NRe = 3, where: L is a characteristic dimension of the obstacle perpendicular to flow, b is the plate spacing, u is viscosity, v, is velocity of approach and p is density. THEORY* Incompressible fluid flow through a Hele-Shaw model and a porous medium are analogous because the equations for the average velocity components in two-dimensional flow are identical if:

Jan 1, 1970

By V. J. Berry, J. D. Pendergrass

Well pressure transient tests provide a means for directly obtaining information about formation pressure and reservoir flow capacity. Such tests have also been proposed for determining presence and location of faults or other reservoir closures and for measuring oil in place. For mathematical convenience, most theoretical studies have considered the reservoirs to be homogeneous. Definitive information is not yet available to show whether the actual presence of nonuniformities will make pressure transient behavior different from that of a uniform reservoir. The conclusions reached from actual transient tests are questionable, therefore, insofar as they rely on the original assumption of homogeneity. One type of nonuniformity commonly assumed to exist is that of stratification. In most reservoirs the strata are thought to be in vertical communication. Equations for the transient flow of a single-phase, compressible fluid in a one-well, bounded, circular reservoir have been solved for several situations involving crossflow between multiple strata of various thicknesses and permeabilities. The results show that except for the very early flow period, which usually is too short to be analyzed, the transient performance observed at the well is substantially identical with that of a homogeneous reservoir having the same dimensions and having the same steady-state flow capacity. Thus, stratification does not adversely affect interpretation of well transient tests. This conclusion holds for all commonly encountered combinations of reservoir thickness and external radius. Deviations are observed for unusually thick reservoirs whose outer radii are relatively small. The results of these studies also show that the presence and the amount of stratification cannot be simply diagnosed from reservoir pressure transient data when there is crossflow between strata. INTRODUCTION The last decade has brought wide acceptance of the transient well pressure test for determining reservoir parameters. Following the original work of Hurst and van Everdingen,1 the mathematical theory was thoroughly explored. Numerous authors have suggested how to determine static reservoir pressure, permeability-thickness product, original oil in place and reservoir limits for different reservoir geometry. The same mathematical techniques have been used to predict the transient performance of a reservoir over a long period of time. Most of the theoretical work has been for homogeneous, isotropic systems. Some results have also been presented for a homogeneous, anisotropic reservoir.2 Petroleum reservoirs are not homogeneous. The deposition process seems to favor creation of a stratified formation. This concept is sufficiently well accepted so that the most natural extension of the transient flow theory beyond the homogeneous case is to a stratified formation. Results for a stratified reservoir with no vertical communication between layers can be obtained from the results for a homogeneous reservoir. Lefkovitz, et a1,3 have given a thorough treatment of the two-layer case. Two recent papers4,5 have treated the case of a two-layered reservoir with vertical communication or crossflow, between layers. Russell and Prats4 find that, after a relatively short time, the two-layered reservoir with crossflow exhibits a simple exponential pressure decline. From this time forward, the behavior is not distinguishable from the behavior of a homogeneous reservoir having the same steady-state flow capacity. The results of Katz and Tek 5 are equivalent. Russell and Prats4 also speculate that a multi-layered reservoir with crossflow will behave as a homogeneous system after long enough production time "... providing the contrast in kh between layers is not too great". They also suggest that, at intermediate times, "... the relative positions of the layers with respect to each other will have a great influence on the production behavior and on the time at which the previously mentioned large-time approximation might be valid". Katz and Tek5 remark upon the mathematical difficulty of treating a reservoir having many layers or strata. If relative position of the layers is also

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval&apos;s heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy&apos;s Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved

By H. B. Bell, P. C. Lidster

A study has been made of the use of ThO2-Y2O3 solid electrolytes to determine activity of chromium in Fe-Cr and Ni-Cr alloys in the temperature range 1300° to 1700°K. This method has been shown to give results which agree with the best vapor pressure measurements. SEVERAL investigations1"5 have been made of chromium activities in Fe-Cr alloys. Most of these have been made by vapor pressure measurements using knudsen cell techniques while one has been made using gas equilibration with H2/H2O gas mixtures. In the present investigation electromotive force measurements have been made on this system using a thoria yttria solid electrolyte. This system was -chosen to establish the technique of measurement at low oxygen potentials, before extending it to Fe-Ni-Cr alloys. The cell used was Pt Cr • Cr2O3 ilThO2 • Y2O3 llalloy • Cr2O3 IPt Since the oxygen pressure of the chromium-chromic oxide mixture lay in the range 10"19 to 1027 atm in the temperature range investigated, thoria yttria electrolytes were used. It was soon found during the investigation that considerable care was required to prepare satisfactory electrolytes and a study was made to find the best method of preparing pellets of TI1O2-Y2O3 solid solution of as near theoretical density as possible. Various compositions in the range 2 to 15 mole pct Y2O3 were studied but the experimental measurements were made with electrolytes containing 4 to 5 pct Y2O3. The method of preparation finally adopted was to use thorium oxalate and yttria as raw materials. The thorium oxalate was decomposed at 900°C for 6 hr to give a reactive powder. This powder was then mixed with the appropriate amount of yttria and ground in an agate mortar for about 4 hr. The powder was moistened with absolute alcohol as a lubricant and pressed in a metal die at a pressure of 7 tons per sq in. The pellets were then heated at 1100°C for 24 hr followed by 24 hr at 1250°C. After this firing cycle the material was reground and repelleted and again fired at 1100° and 1250°C followed by heating for 7 days at 1600°C. It was found in later experiments that a further firing of 4 hr at 2000°C improved the density further and gave pellets of only 2 pct porosity. After firing, the pellets were polished on emery paper and finally with the 2 µ diamond paste. The polishing was followed by a final firing of 24 hr at 1500°C. Investigation showed that this final treatment minimized reaction between the electrolytes and chromium oxide. It was believed that this was because of a lowering of the surface energy of the electrolyte surface. Two types of investigation were carried out to check that the pellets were satisfactory oxygen ion elec-

Jan 1, 1970

By Oliver Bowles

WHILE vast quantities of limestone and dolomite are used in metallurgy, the estimated production in 1926 being 23,860,000 tons, there are many problems connected with their use which have not received adequate study. The literature of metallurgy is notably lacking in comprehensive discussions of fluxing or furnace stone. Approximately 900 lb. of limestone is used for every long ton of pig iron produced in the blast furnace, but this important constituent of the charge receives little attention compared with the intensive study that has been applied to ores and fuels, the other important constituents of the charge. In an effort to elucidate some of the doubtful points, the writer has made a wide review of metallurgical literature, and has had the disappointing experience of finding not more than a paragraph or two on metallurgical stone for each two or three thousand pages of literature reviewed. Furthermore, the limited information supplied is in some instances contradictory and even absurd. For example, a correspondent in a reputable metallurgical magazine classifies magnesia and silica together as impurities in flux "that must be melted and run off as slag." From the furnace operator&apos;s viewpoint it is desirable that as complete knowledge as possible be obtained concerning the effects of chemical composition and physical character of the stone on its action in the furnace and on the quality of the metal produced. A better knowledge of these features is advantageous also to the quarry operator in that it enables him to work more intelligently in producing stone that will best satisfy the metallurgists&apos; requirements. The following pages comprise a brief compilation of data now available. Much of the information is well established and probably acceptable to most metallurgists. On some points there is room for debate, and the chief purpose in presenting a paper at this time is to arouse discussion which will throw more light on controversial points and thus assist the author in preparing a more comprehensive and accurate dis¬cussion of metallurgical limestone problems.

Jan 1, 1928

By Earl C. Henry

Magnetic iron ore was mined in Essex County, N. Y., during the American Revolution; Benedict Arnold is said to have mined ore near Port Henry to secure iron for chains and spikes for the Lake Champlain fleet. Details regarding the operations prior to 1804 are not easily obtained but work on a commercial scale appears to have been started about 1838. The Witherbee, Sherman & Co. was established in 1849. The property of the district is held in fee almost entirely by the mining companies. The Mineville iron-ore district, as well as the entire Adirondack region, belongs to the pre-Cambrian complex. The readily identified rocks consist of gabbro, anorthosite, syenite, and granite, and are classed as igneous, eruptive, and intrusive. In addition, there is a great variety of gneisses, rocks with a variable mineral composition but all more or less distinctly foliated. The Grenville series are the oldest rocks in the district; this series consists of graphitic white limestones, graphitic sandstones, and quartzites and gneisses with or without graphite. The orebodies occur in the gneissic rocks and are comformable to the structure of the rocks. The typical occurrence of the ore in this district is in lenticular bodies with a parallel pitch to the southwest; this structure is best exemplified in the New Bed mine. The Harmony and Old Bed orebodies are faulted by two sets of diabase dikes, one set having a southwest strike and the other set an east-west strike. The ores occurring in the acid gneissic rocks are non-titaniferous. The ores occurring in the gabbros and anorthosites are high in titanium; these deposits are irregular, have no structural relation with one another, and are unquestionably of igneous origin. The origin of the gneissic rocks and the ores therein contained is not settled; both the sedimentary theory of ore concentration and the igneous theory of magmatic segregation have been advanced. The iron mineral in all the orebodies is crystalline magnetite; some very perfect and large specimens have been obtained. The gangue is

Jan 1, 1925

By Owen R. Rice

Changes in coke hardness affect the working of the blast furnace, for soft coke is an obstacle to proper furnace operation. Soft coke is due to a low hydrogen-oxygen ratio in the coal charged; increasing this ratio increases the hardness of the coke. The best means of determining the hardness of coke is the combined tumbler and ball-mill. What is Good Coke? It is not always appreciated that the correct answer may be such a reliable forecast to the physical condition of the blast furnace, and incidentally its entire operation, that among the ordinary variations in other raw materials and also in operating conditions, it assumes a position of primary significance. Inasmuch as coke forms two-thirds of the total bulk of material charged, this assertion will not seem exaggerated, and in the subsequent discussion its applicability will be presented. Coke is the ideal blast-furnace fuel because it is strong and hard enough to resist the crushing and grinding effects of the furnace burden, and its porosity allows the ready passage of gas upward through the charge, and rapid combustion before the tuyeres. Coke deprived of its strength and with a cell structure physically frail and chemically susceptible to solution by carbon dioxide has its chief asset eliminated. When the melting zones are approached the fluid cinder trickling toward the hearth will envelope the small pieces of coke and form a sticky mass that will retard the free passage of gas. Blast pressure will then rise, driving will slow up and fragments of coke will lie inert before the tuyeres, rendered practically incombustible by their slaggy envelop. It is believed also that coke messes at cast are attributable to quantities of such coke floating unburned in the hearth. It has been correctly claimed that coke rendered small outside the oven possesses greater combustibility than that so developed in the

Jan 1, 1922