Search Documents

By J. W. Leomrd, B. A. Donahue

Results of research are presented examining the extent to which the analytical characteristics of the relatively distinct coal bands from a variety of coal seams can be related to each other. This paper dis-cusses an approach for developing a practical coal petrographic quality control program based on conventional analyses, most of which is part of the standard A.S.T.M. procedure. The work is a final follow-up to Part I of this series which was prepared by the authors as an approach to applying conventional coal petrography to single coal seams. Recent published work by the authors entitled Petrography for Coal Mining and Coal Prepration: Part I dealt with interrelating various chemical and physical properties of coal1 measured using conventional analyses made on distinct petrographic bands taken from single coal seams. Since most coal production facilities process coal from a single coal seam or from very closely related coal seams in the same area,2 emphasis on interrelating the properties within a single seam appeared appropriate. The distinct petrographic bands were analyzed on the assumption that such differentiated data would be more representative of a heterogeneous coal seam than the single analytical value commonly used to characterize each property of an entire seam.34 Effort was directed at demonstrating the extent to which the interrelated chemical and physical properties 5,6 could be developed into nomographs or petrographic standardization graphs. Thus, one analysis, determined on a series of petrographic fractions separated from a single sample, was used to estimate numerous other properties in each fraction by referring to the previously established petrographic standardization graphs. This conventional approach to coal petrography was undertaken as a suggested feasible means by which a few coal analyses could be employed to develope a more penetrating knowledge of the properties of coal from any given seam in order to monitor more extensively its performance at the point of utilization. Such procedures can support development of the type of in- formation commonly sought through automated testing and through the use of computers.7 The broad knowledge which can be developed through these procedures is intended for application in the generation of an analytical profile or broad characterization of coal. These estimates were not intended as replacements for individual coal analytical tests. In this expanded second part of the research program, distinct petrographic bands from nine coal seams in the Central Appalachian Region were physically and chemically analyzed to elucidate the extent to which this concept of conventional petrography could be broadened for application to numerous coal seams. In presenting this second phase of work, the relationships developed are presented individually and not in a connecting series of nomographs or petrographic standardization graphs as in the previous work, thus leaving open the combinations of possibilities to individual interpretation and application. MATERIAL AND EXPERIMENTAL WORK Nine coal seams, representing a wide range of rank, from the Central Appalachian coal fields were used in this study. The distinct petrographic bands from the Kittanning, Pond Creek, Jawbone, Tiller, Poca-hontas No. 3, No. 2 Gas, Eagle, Winifrede, and Pittsburgh coal seams were separated by carefully removing a portion of each band at the face of the seam. The following were determined: ash, sulfur, free swelling index, heating value, bulk specific gravity, volatile matter, Hardgrove Grindability Index, and Gieseler Plastometer measurements. Determinations, where procedures were available, were carried out using ASTM standard procedures.8 Bulk specific gravity was determined using a kerosene volume displacement procedure modified from a method applied by Headlee and McClelland of the West Virginia Geological survey,9 Those bands with a bulk specific gravity greater than 1.60, which is generally above the practical specific gravity cleaning range of bituminous coals, were excluded and no analyses were performed. Much of the initial organization of this second phase of work was developed through the extensive use of a computerized statistical monitoring program (see Ref. 4). However, in order to achieve the closest possible interpretation of results, the final organization of data proceeded mainly from exhaustive trial

Jan 1, 1968

By S. Stopek, E. D. Hinkley, R. H. Rediker

Epitaxial heterojunctions have been prepared by melting the lower -melting-point semiconductor of the interface between two dijferent semiconductors. when the temperature is reduced, the melted material recrystallizes, having alloyed into the higher -melting-point semiconditctor. The electrical and elcctro-optical properties of such single-crystal heterojunctions between GaAs and Gash and between p-type InAs and n-type Gash are the subject of this paper. The forward current varies as exp (AV), where A is substantially independent of temperatuve. For Gds-Gash heterojunctions at temperatures above, 370°K, if the current-voltage relationship were to he expressed as exp (qV/nkT), then n would he less than unity. The injection luminescence associated with forward current is, for the most part, characteristic. of the lower bandgap semicondutctor. These results can he explained by carrier injection into the lower bandgap semiconductor by tunneling through a barrier at the interface. The photovoltaic effect measured for incident photons having energies in the range between the bandgaps of the two semiconductors is much smaller than that produced by higher-energy photon The smallness of this between -the-gap photovoltaic response can he explained by the low probability for penetration of the barrier by the carriers produced in the smaller -bandgap semiconductor. THE technique of interface alloying has been used to produce single-crystal junctions between dissimilar semiconductors.' Oriented wafers are placed on a carbon heater strip, Fig. 1. so that semiconductor S1, which has the lower melting point, is supported by semiconductor S2. Electrical current passed through the heater strip produces a temperature gradient such that S2 is at a higher temperature than S1. As the temperature is raised the lower face of S1 begins to melt. Before the entire wafer can melt, however, the heater-strip current is turned off and, as illustrated in Fig. l. the melted portion recrystallizes, having alloyed into S2. Junctions have been fabricated by the above procedure between GaAs and germanium, between GaAs and GaSb, and between InAs and GaSb. Mroczkowski, Lavine. and Gatos have described the metallurgical and chemical aspects of the GaAs-Ge junction.2 The transition from GaAs to germanium is not monotonic and a portion of the recrystallized region consists of the GaAs-Ge eutectic. Since gallium is an acceptor and arsenic is a donor in germanium, since germanium dopes GaAs, and since the electrical properties of the GaAs-Ge eutectic have not been investigated, any interpretation of the electrical characteristics in terms of simple heterojunction theory would be incorrect. That the rectification of GaAs-Ge heterojunctions is not a property of the impurity doping of the GaAs or the germanium, but is most probably due to the impurity distribution in the recrystallized region, is clear from the fact that forward conduction occurred for all the GaAs-Ge interface-alloy junctions (whether they be n-n,n-p. p-n. or p-p) when the germanium was biased positively. The electrical characteristics of these GaAs-Ge heterojunctions will not be discussed further in this paper. Electron-beam microprobe analysis of GaAs-GaSb heterojunctions showed that the transition from arsenic to antimony atoms was without structure. and that the transition occurred within a 2 to 3 region.' In this paper we will describe the electrical properties of the GaAs-GaSb heterojunctions as well as electrical and electro-optical properties of the InAs-GaSb heterojunctions. A band model for the junction will be proposed which can explain these properties. ELECTRICAL PROPERTIES OF GaAs-GaSb HETEROJUNCTIONS After interface alloying. in preparation for the electrical and electro-optical experiments, ohmic contacts were made by conventional means. For example. Kovar tabs clad with tin were alloyed to n-GaAs or similar tabs clad with Au-Zn were alloyed to p-GaAs. The units were mounted and then etched. All four combinations of conductivity types

Jan 1, 1965

By R. M. S. B. Horta, W. T. Roberts, D. V. Wilson

Evaluation of results obtained by the Harris method for inverse pole figures is discussed. Two existing analyses and a new approach are compared. In the most frequently used analysis, different reflections are accorded equal weight. A second method involves area-weighting of reflections, mainly to take account of nonuniformity of Pole distribution. In the new method, reflections are weighted according to their multiplicity factors. Intensity ratios obtained from a sample of steel sheet by the three methods are compared. INVERSE pole figures offer a convenient method of depicting the proportions of material with various orientations referred to a unique specimen axis. The technique was first applied to fiber textures, but it is now frequently used as a description of texture in sheet metal, the unique axis in this case usually being the sheet normal. The Harris method' for the determination of inverse pole figures has the great merit of simplicity. It suffers certain intrinsic limitations, but provided that these are borne in mind, the results allow useful conclusions to be drawn regarding the effects of texture components on certain mechanical properties. Much thought has been given to alternative methods of texture description, aimed at eliminating the ambiguity in the simpler methods.2"6 Refinements in the method of description are obtained at the price of increasing complexity, and it seems clear that a simple method, such as that of Harris, will continue to find application. The original derivation by Harris1 has been corrected by Morris7 and later by Mueller, Chernock, and Beck.2 Since the latter reference is more accessible, it is referred to more frequently than the work of Morris. Most users of the inverse pole figure method present their results on the basis of the treatment by Mueller et al., and it is the purpose of the present communication to comment on this treatment and to suggest a simple modification to the method of presentation which enhances the meaningfulness of the results. It will be convenient first to summarize the Mueller analysis as follows: The intensity Ihkl of a hkl reflection from a textured specimen is: Ihkl=CI0AL\Fhkl\2NhklPhkl [1] where C is a constant for a given sample I, is the intensity of the incident beam A is an absorption factor. With usual diffrac-tometer geometry and a flat specimen, A is inversely proportional to the density of the specimen. L is the Lorenz polarization factor Fhkl is the structure factor for the hkl reflection Nhkl is the multiplicity factor phkl is the fraction of crystals in the polycrystal-line aggregate with any particular (nkl) plane parallel to the surface. The corresponding equation for a specimen with randomly oriented grains is lR,hkl=CRl°ARL\fnkl\2Nnklk [2] the subscripts R indicating the factors that are different for the random sample, which will, in general, have a density different from that of the textured one. The intensity ratio is hkl , c A Phkl [3] IR,hkl CR AR PR The constants are eliminated by summing over all the measured reflections ^iRMi °R AR PR ljPkkl and from Eqs. [3] and [4] T IR,hkl The only unknown in Eq. [5] is Yj phkl and before this can be evaluated, it is necessary to adopt some form of normalizing procedure. Mueller et al. assume that when a large number, n, of reflections is used, the average value of p as defined by Eq. [6] below will be unity ? = £j**l3l [6] and the final result of their analysis follows: hkl Pkkl = i I&Jf~ [7] ±y* Ihkl ^ in,hki The normalizing procedure used by Mueller et al. involves an unrealistic assumption since it gives equal weight to each of the hkl reflections. Other normalizing procedures can now be examined with a view to selecting a more appropriate one. For an infinite number of reflections, the function p can be averaged over a continuous range of orientations covering the complete solid angle 4p by the integral relation: p =1 jp(a)dQ = l [8] This is the only exact procedure, but in practice a

Jan 1, 1970

By J. C. Purcupile, R. L. Morris

Present methods of roof support require that first a hole must be drilled, and then a toggle belt or grouted resin bolt be set. Automating this system is very difficult. Therefore, we decided to develop a system which would insert pins into the roof in one step. We think that this type of system can be automated much easier. It would then speed up the roof bolting operation and eliminate the need for a man to work under unsupported roof, making roof bolting a much safer job. There are two ways that a pin can be inserted in a mine roof without drilling a hole. It can either be pushed in, using a very high, constant axial force; or it can be hammered in like a nail. The former method involves very high compressive stresses in the pin, and requires complicated machinery to prevent buckling. Since hammering the pin does not require high axial forces buckling is not as likely. Also, experiments have shown that there are no problems associated with getting the pin to go in straight. We chose to develop the second method. While investigating methods for hammering a pin into a mine roof, we learned that piles were being driven 10 times faster than by conventional methods with a method called "orboresonance." We decided to turn the pile driver upside down, and drive pins up into the mine roof. Sonic Pin Setting Machine The machine consists of a steel bar supported at the nodal points of its first mode of transverse vibration (Fig. 1). A mechanical oscillator, driven by a hydraulic motor, is attached to one end of the bar. This provides the forcing function for the vibrations in the bar. The hydraulic power is supplied by a portable unit. When the forcing function is near the natural frequency of the bar, large amplitudes of vibration occur and require little power input. An anvil is attached to the bar at a point of maximum deflection, providing a means for striking the pin. Since the bar is being driven near its natural frequency, it takes little energy to accelerate the bar back and forth between each hit. This means that most of the power being put into the bar is going into driving the pin. This is unlike conventional pile drivers which use most of their power to accelerate the hammer back and forth between hits. The resonant frequency of the bar we used is about 250 Hz, which is in the audible range, hence the term "sonic pin setting machine." The system benefits are: (1) elimination of the two-step drilling and pin-setting roof bolt operation, (2) increased speed, (3) more intimate contact between the formation and bolts, (4) less expensive roof bolts, (5) higher integrity, (6) potential for automatic operation. Modifications During the months from Oct. to Dec. 12, 1974, under U.S. Bureau of Mines grant No. 00144104, "Projects on Coal Extraction," six members of the senior class at Carnegie Mellon University designed and built the experimental machine. They drove three pins into the Safety Research coal mine at Bruceton, Pa., and determined that the project was worth continuing. In January 1975 we began making modifications and repairs in preparation for a more extensive program of testing at Bruceton. The major items were: 1) Design and build a collapsible pin guide to keep the pin parallel with the line of travel of the lift table. 2) Move the anvil to the center of the bar. 3) Refinish the nodal holes and install new node pins and bushings. 4) Realign the bar in its frame. 5) Install a relief valve on the lift unit to control the upward force. The upward force that we exert on the pin is relatively small, so there is no tendency for it to buckle. This means that all we need to do is hold the pin straight during its initial penetration. This allows the bushing to travel up and down in a straight line. At the beginning of penetration, the bushing

Jan 1, 1977

By D. D. Smith, D. P. Squier, E. L. Dougherty

A system of differential equations describing the temperature behavior of fluid injected at constant surface temperature in a well is derived and .solved analytically. A formula for the fluid temperature at any time and depth is given, as well us a special formula valid for very large times. These formulas are used to calculate temperatures for several typical cases. The results indicate that, initially, the temperature of the water entering the formation is considerably lower than the injection temperature. This condition lasts for only a short period— less than three days for most cases of practical interest. Following this highly transient period, during which the temperature of the fluid entering the formation builds up to about 50 to 75 per cent of the injection temperature. the system enters a quasi-steady state in which the temperature changes are very slow. After severl years, the bottom-hole temperature will still be 50" to 100°F lower than the injection temperature, hilt the heat losses may he tolerable. INTRODUCTION Predicting the behavior of a hot-water flood requires that the temperature of the water entering the injection interval be estimated. This report describes the development and solution of a system of equations which describes the temperature behavior of the injected water in the wellbore with certain simplifying assumptions. The only previous means known to the authors for describing such a process is that of Moss and White.' Their results appear to be close to those obtained by our method in the practical cases which were compared; this agreement is largely due to the fact that in our method temperature soon approaches a quasi-steady state, as was assumed in their method throughout. However, our model covers all times, is continuous (whereas the Moss-White model depends on breaking the depth into discrete intervals) and. we feel. more closely describes the physical problem. FORMULATION OF THE PROBLEM PHYSICAL SYSTEM AND ASSUMPTIONS The injection procedure consists of pumping water at a fixed surface temperature T., down an infinitely long cylindrical well or tubing of inner radius Any material exterior to the water column such as mud, casing, or cement is regarded as part of the formation. The general behavior of the system may be described qualitatively as follows. When the hot water is first introduced into the system, the temperature difference between the formation and the water is large, resulting in a high rate of heat transfer. As a result, the temperature adjacent to the wellbore rises very quickly. Because the segment of the formation adjacent to the wellbore largely controls the heat transfer rate, the heat transfer rate will become relatively constant when this portion has reached a temperature close to that of the water opposite it. The temperature of the water and formation then increase very slowly with time. The length of the initial highly transient period and the temperature of the water at its conclusion will be functions of depth, injection rate, injection-string radius, surface injection temperature and the physical properties associated with the water-formation system. The following additional assumptions were made. 1. There is no heat transfer by radiation in the system. 2. There is no heat transfer by conduction in the vertical direction in either the injection stream or the formation. 3. The linear volumetric and mass flow rate of the water is constant throughout the injection stream. 4. No horizontal temperature gradient exists in the injection stream. 5. The product of density and heat capacity is constant for both the water and the formation, and the formation thermal conductivity is constant. 6. Initially, both the water in the wellbore and the reservoir are at a temperature given by the (constant) ambient surface temperature plus the product of depth and geothermal gradient (assumed constant). At large distances for the wellbore (r m), the formation will remain at this temperature. 7. The water temperature and the formation temperature at r — r,, are equal for all depths D. DERIVATION OF EQUATIONS The differential equation satisfied by the fluid temperature T,(D, t), which is obtained by writing a heat balance on a cylindrical differential of volume dV of the injection string between the depths D and D i dD, is

By H. Dykstra

Calculated wellbore pressures were obtained for parameters of radilcs ratio and permeability. In all cases bur two, after-production was allowed to occur for one day. The calculated pressure build-up data were compared with actual pressure build-up data from a condensate well. The paper discusses the comparison, gives reasons for, inability to match calculated (2nd actual data, and presents conc1usions derived from the comparison. INTRODUCTION Considerable information has been published on the decline in wellbore pressure for oil and gas fields resulting from production of fluids. Considerable information also has been published on pressure build-up for oil fields, but relatively little has been published for gas and condensate fields. Perrine' summarized the equations derived for pressure build-up in oil wells and showed how the different methods could be applied. On the other hand, very little theoretical information has been published on pressure build-up in gas or gas-condensate wells.2-4 Tracy,' by pointing out the similarity between the equation describing oil wellbore pressure decline and the equation describing gas wellbore pressure decline, showed how methods of pressure build-up analysis could be applied to a gas well. In analyzing pressure build-up data for gas or gas-condensate wells, it would be very desirable to compare actual pressure build-up data with calculated pressure build-up data. Such a comparison would lead to a better feeling for the quantitative picture of gas-well build-up curves and would help to establish a degree of confidence in a given build-up curve. This paper discusses the comparison and presents conclusions made from it. Calculated wellbore pressures are given for parameters of radius ratio and permeability. In all cases but two, after-production was allowed to occur for one day. The work was done in an attempt to evaluate the reserves of a gas-condensate field. As will be discussed later, it was not possible to make an evaluation. METHOD OF CALCULATION The method used for solving the depletion and pressure build-up behavior for a gas-condensate well was that developed by West, Garvin and Sheldon.5 or the IBM 704 computer program used, it was assumed that liquid condensing out of solution would have only a minor effect on the flow behavior. Therefore, the field was treated as having single-phase flow, with the gas saturation S, remaining constant at unity minus the connate-water saturation Sw. As will be discussed later, this may not have been a good assumption. It also was assumed that oil production could be considered as equivalent gas production by using a conversion factor based on an analysis of the liquid produced. An outer boundary condition of no-flow was assumed. BASIC DATA The data required for the study were reservoir properties, fluid properties and production data. Reservoir properties included the following: thickness h, 179 ft; porosity, 0.194; connate-water saturation, 0.43; gas saturation, 0.57; permeabilities k, 1.0, 0.5, 0.25 and 0.15 md; radius ratios re/rw, 800, 1,000 and 2,000; wellbore radius rw, 5 in.; initial reservoir pressure P, 6,529 psia; and reservoir temperature, 272°F. The selection of permeabilities was based on core-analysis data showing an average air permeability of about 2 md. With a relatively high connate-water saturation, the average effective gas permeability would be about one-half, or less, of the air permeability. The selection of radius ratios was based partly on the observed rapid decline in wellbore pressure during production and partly on the fact that the well was believed to be located in a relatively small fault block. Reservoir fluid properties were obtained from a reservoir fluid study on a recombined sample of gas and condensate. Gas formation-volume factors were calculated from pressure-volume measurements made at reservoir temperature. Gas viscosities were calculated by the method of Carr, et al,G from an analysis of the re-combined sample. These data are shown as a function of pressure in Table 1. The production schedule from the gas-condensate

By K. E. Beu

IN the measurement of retained austenite concentrations in steels using the integrated intensity method,1 Averbach has pointed out" that the absorption factor A(0) for a flat sample making a glancing angle 4 with the incident X-ray beam can be combined with his constant factor R to obtain our constant factor" G provided that "1—There is no preferred orientation in the sample, and 2—The geometric requirements [for the sample, film, and X-ray beam] have been met precisely."' If A(0) can be combined with R to give G according to the equation G = R . A(0) [14] then the possibility for making non-compensating errors in the austenite determination has been eliminated. This has been described previously in a technical note." Because of the brevity of the previous note,3 it was impossible to emphasize the fact that the two conditions on preferred orientation and geometry can be easily met experimentally, contrary to Aver-bach's statement that: ". . . the necessary conditions must be tested experimentally for each determination, and this is done most easily by observing whether the apparent absorption has the form of Eq. 2 [the theoretical equation for A(0)]."" he way in which these two conditions can be met experimentally will be discussed briefly to help clear up this point. In addition to these two conditions, other factors such as sample shape, homogeneity, and grain size which also affect A will be included in this discussion. A (0) depends on sample shape. The theoretical function for A was derived originally for a flat surface.' In general we have found that a flat sample surface is readily obtainable." If such a surface can * The effect of surfaces which are not flat on the measurement of retained austenite will be discussed later. be obtained, it has the following advantages: 1—it is easily reproducible from sample to sample, 2—it is the form required for metallurgical examination —this type of examination being frequently desirable for this work, and 3—it is an efficient shape for diffraction purposes. Perhaps the most important of these features is that a flat sample surface is easily reproducible; hence, this requirement on the repro-ducibility of A from sample to sample is met for all such samples. Fig. 1—Schematic diagram of the quartz crystal monochromator diffraction unit. The centerline of the main beam is at 6' to the target face. The tangent to the crystal face is at 16.8' to the moin beam. The centerline of the monochromatic beam is at 33.6 to the main beam. The sample surface can be rotated in its own plane. The angle of the sample surface and the monochromatic beam can be adjusted by rotating about the vertical axis, The film holder can also be rotated about B so that the film can be exposed over the desired angular range. For Fe K, X = 1.932A. 1011 planes of quartz have d = 3.35A. A(0) depends on homogeneity and grain size. If the sample is badly segregated or the grain size is large, the effect of micro-absorption and primary extinction1. ' must be considered. It has been shown, however, that for most plain carbon or low alloy hardened steels, neither micro-absorption nor primary extinction effects are present.' A (0) depends on camera geometry. For a given angle 0, A remains theoretically constant from a geometrical viewpoint only if the following factors are kept constant: 1—the angle of inclination .+ of the sample to the X-ray beam, and 2—the centering of the sample with respect to the film cylinder. These are mechanical problems which can be solved readily if the facilities of a good machine shop are available. The arrangement used to insure that the angle 4 remains constant and the centering of the sample is reproducible is indicated schematically in Fig. 1. The sample is clamped against a thin plate with a hole in it by means of a spring-loaded pres-

Jan 1, 1954

By H. R. Crawford, P. B. Crawford, A. F. Tragesser

A method has been developed to calculate wellbore temperatures during mud circulation and the actual cementing operation to aid in the design of cement slurries. The method agrees within 10F with previously measured values. The calculation technique provides temperatures, as functions of time, at varying depths in both the casing and annulus. The technique also provides this information if a relatively cool cement slurry is pumped into the well immediately following circulation of hot mud. Circulating bottom hole temperatures of brine and a bentonite mud were measured. INTRODUCTION As wells are drilled deeper, greater demands are being made on all phases of the industry, and new technology has been developed to provide satisfactory well completions. However, little or no work has been conducted on accurately determining bottom-hole, static and circulating temperatures.. In designing a cement slurry, such factors as density, fluid loss control, viscosity, deterioration from ternperature, compressive strength and pumping time must be considered. Individual well conditions often make it necessary to include still other factors. Pumping time is a primary consideration and, as wells are drilled deeper, encountering higher bottom-hole temperatures, this property becomes even more important. Cement slurries must be designed with sufficient pumping time to provide safe placement in the well; however, the slurry cannot be overly retarded as this will prevent the development of satisfactory compressive strength. The pumping time of a specific cement is currently obtained by subjecting the cement to simulated conditions of temperature and pressure. A reasonably accurate bottom-hole pressure may be obtained by considering hydrostatic heads of fluids, friction pressure and wellhead pressures. However, accurately determining bottom-hole temperatures is much more difficult. Bottom-hole static temperatures are estimated by considering several sources of information, including logging temperatures, published temperature gradient maps and field experience. This information is usually questionable due to disagreement of data from the various sources. Temperature gradient maps were constructed based on temperatures recorded many years ago while running bottom-hole pressure tests. These thermal gradients then represent an average of well conditions and cannot always apply to a specific well. Also, logging temperatures may be affected by the time since fluid was last circulated, rate of penetration, circulating rate and many other factors. Therefore, even though logging temperatures are available, the question still exists as to the correction factor that should be applied to obtain an accurate static temperature. After obtaining static bottom-hole temperature, it is then necessary to relate this to circulating temperatures actually encountered by the cement slurry. This is accomplished by selecting a test schedule from the API RP-10B corresponding to the estimated well conditions.' The API-recommended practice for testing oilwell cement provides testing schedules for various well depths and conditions. These schedules are intended to simulate down-hole conditions during cementing. They provide a rate at which both temperature and pressure are increased until the estimated circulating conditions are reached. These testing schedules represent circulating temperatures for an average well and, although there is flexibility in choosing the test schedule that most accurately simulates the temperature of an individual well, it still is not possible to consider all the well conditions that will affect the bottom-hole temperature. Many factors affect cement temperatures; for example, the length of time a well has remained static prior to running casing and cementing, the circulation time, the temperature of fluids used in cementing, fluid density and flow properties of fluids. The pumping time for a typical retarded cement could vary from 2 to 4 hours with a 10F change in testing temperature. Variations in pumping time are the most critical in highly retarded cements used in deep, hot wells; yet, predicting bottom-hole circulating temperatures is more difficult in these wells. This work was conducted to develop a means of calculating circulating temperatures as a function of well depth, casing and hole size, pumping rate and time, fluid and reservoir physical properties and thermal status of the well. PREVIOUS WORK In 1941 Farris reported on a study to develop information leading to a more practical laboratory evaluation of oilfield cementing mixtures and performance.' It was then recognized that the pressure factor was being neglected,

By John D. Price

COAL blending, in the preparation of coal before coke making, is so commonly practiced as to be almost universal. But the reasons underlying this practice, the benefits resulting from it, and the materials used in blending vary widely. This paper will outline the various phases of the subject and present information which may be correlated with work that has been done elsewhere. It will deal entirely with work done on the high-volatile coking coals of the western part of the United States, special emphasis being given to the coals of Colorado and Utah. A surveyL of the 86 coke plants in active operation in the United States during 1949 indicates that only 9 plants, or 10.5 pct of the total, charged one single rank of coal into their ovens, while the remaining 89.5 pcl made use of blending in some form. This report indicated that of these total plants 5 used straight high-volatile coal, 4 used straight medium-volatile coal, 47 used blends of high and low-vola-tile coals, 25 used blends of high, medium, and low-volatile coals, 2 used blends of high and medium-volatile coals, 3 used blends of medium and low-volatile coals. The fact that certain plants operated on a single kind of coal should not be interpreted to mean that no blending was practiced there, for invariably such plants secure their coal from more than one source and in the interest of uniformity do blend the coals as received. The general term coal blending covers two fields, the first of which is the mechanical mixing of a number of coals to secure uniformity. Often it is found necessary to secure coal for coke production from a number of different mines; these coals, though of the same general type or rank, may differ in their chemical composition or in the physical qualities they impart to coke made from them. Again, it is not unusual to find that coal from different sections of the same mine may show variations in quality. Under such conditions it may be necessary, in the interest of a uniform final product, to introduce a system of blending bins, a bedding yard, or other mechanical methods of securing a uniform mixture. Unfortunately this form of blending has received very little attention up to the present time; it has not received the consideration its value merits. The second type of blending, while also for the purpose of coke improvement, deals more particularly with the use of a blending agent differing in character from the base coal: it is this form of blending that will be discussed here. To consider only the western coals, for blending may be found necessary for other reasons with other coals, blending has been practiced experimentally or commercially under the following conditions: 1—When a single coal or mixture of coals of the same rank does not produce a satisfactory coke. For example, a high-volatile coal when used alone is likely to contract when coked so that a comparatively weak coke is formed. Or, if of very low rank, the coal may be deficient in the necessary bitumens required for good coke production. 2—When a product of some special quality is required, for example, when a plant ordinarily producing blast furnace coke must operate at slow coking rate to produce a high-grade foundry coke. Under this condition the reduced daily production of all products which accompanies slow coking time may be undesirable, and the use of some blending agent to increase the size of coke made at faster coking rates may be necessary. 3—When greater yield of coke or its coproducts is needed. Depending upon economic values of the products it may be found desirable to increase the yield of one or the other. 4—When supply of a particular coal must be used, either to protect reserves of high quality coking coal or to utilize a surplus or inferior product not otherwise usable. Many materials have been used for blending purposes, the exact agent to be used depending both upon the condition to be corrected and the nature of the base coal. No universal blending agent that can

Jan 1, 1954

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By Alfonso W. Ovalle

INTRODUCTION The block-caving mining method is one of the most economical ways existing today to extract ore from nature. It is undoubtedly the least costly of the underground systems and moreover it compares favorably to open-pit mining. In addition to its low cost characteristic, the block-caving mining method permits a great extraction volume as can be exemplified by El Teniente's output capacity of over 22 million tons of ore per year. These characteristics of block-caving give assurance that, even though it is a very old mining system, it will be a good answer to the many future large massive low ore grade deposits to come and a1 so to rep1 ace open-pi t mining as such operations get deeper with higher stripping ratios. These are the reasons that induce further study of block-caving. In this paper block-caving at El Teniente is analyzed with due consideration to the major problems that this mine will encounter. Emphasis is given to the factors influencing the mining of soft ore which accounts for most of the ore to be extracted in the next ten years. A separate paper is devoted to a major modification introduced to block-caving hard primary ore at El Teniente by means of LHD equipment extraction. DESCRIPTION OF THE PRESENT MINING Ore Production and Flow Over 62.000 tons per day of ore are extracted 359 days per year at the El Teniente Mine. The Teniente 8 train haulage level at 1.983 m takes 29.000 tons daily to the Col6n Concentrator. Two trains make the 8 km haul with 16 cars per train. Car capacity is 90 tons. The Teniente 5 train haulage level 298 m above takes 33.000 t per day to the old Sewell Concentrator. Ore is produced from three production areas in two different levels. The so called North Mine produces 42.000 tons daily from the Teniente 4 level 65 m above the Teniente 5 haulage level. The Central Mine produces 1.000 tons from the first hard rock block caved at the Te- niente 4 level. The combined 43.000 tons per day production is sent through inclined 2 m dia- meter raises to the Teniente 5 haul age level , where 33.000 tons are taken 2.5 km to the Sewell Concentrator and 9.000 tons are hauled 0.5 km to a 5 m diameter vertical ore pass which sends ore to Teniente 8. The haulage is by 7 trains, 15 gabble-bottom cars per train. Car capacity is 20 tons. The so called South Mine produces 19.000 tons per day from the Teniente 1 level 645 m above the Teniente 8 haulage level. Ore is transferred by 2 m diameter inclined raises to the Ten 1 Retram level 22 m below, where 5 trains with 20 Sandfor Day cars per train, 8 tons capacity per car, haul the ore 1.5 km to two vertical 5 m diameter ore passes which connect to Ten 8. Starting in 1982 the Central Mine will in- crease production to 6.000 tons per day. An in- dependent loop haulage system will start operat- ing at the Teniente 5 level to feed a new under- ground 1.372 mm (54 in.) gyratory crusher which will send minus 152 mm ore (9 in.) to Teniente 8. In two or three years time the North Mine will only feed the 33.000 tons to the Sewel 1 Concentrator with its old Teniente 5 haulage system. The Col6n Concentrator will be fed from the Central and South Mines. Figure 1 shows the transfer systems and ore flows. The Ore Body In plan view at the Teniente 4 level the ore body has a kidney shape 1.800 m long and between 400 to 700 m wide. The north-eastern section of the ore body was washed away by the Teniente River. At the Sub 5 Level, 125 m below Teniente 4, the mineralization completely surrounds the circular 1.000 m diameter steep-sided barren breccia pipe (Braden pipe) . The ore body rock is 73% andesite, 14% diori te, 13% dacite and 6% mineralized breccia. Most of this ore is secondary enrichment sulfides which apart from the dacite is strongly fractured being chal cocite, cove1ine , bornite and

Jan 1, 1981

By V. C. Williams

Some of the physical, chemical, and electrical processes for conversion of saline water to potable or industrial water are economically surveyed from an engineering viewpoint. Since all these processes require energy for drive and equipment for containment, the correlative economic factors are developed which indicate directive influences in the choice of particular regional processes. The supply of natural waters and its distance also affect decision. Any one process will probably not prove dominant in the field because auxiliary considerations such as the saline water source; types and continuing availability of fuel; electric power use or recovery; area economic status and advancement; and the political pressures of population, group demands, and land use tend equivocally to obscure capital and operation cost decisions. Basic engineering considerations, data, and economic factors are presented to assist in the direction of these decisions. An exploding world population, increasing industrialization, advancing standards of living, and the desire of less-privileged nations for betterment focus attention sharply on a major problem: water. *19 Up to now, in retrospect, people have had it relatively easy in the handling of this problem. All the better dams in the most advantageous sites, the better aquifers, the shortest aqueducts have been built. In another phase of the problem, concern is evident that wastes cannot indefinitely be disposed of merely by keeping them dilute and discharging them promiscuously. 7-9 And, perhaps, as past civilizations have done,l5 water, watersheds, streams, and irrigation may have been mismanaged or, at the least, not adequately studied.3,5,36,37 In this last is perhaps the core of the problem. As Gross states, "Ignorance and too often, indifference are contributing factors. Archaeology and theology both furnish ample testimony to the existence of rich lands where deserts now stand; it was man who ravaged his land. Unless education is a companion to water development, development might as well be forgotten. But without water, there is no beginning."13 The U.S. is showing increasing concern about its water for predictions are that by 1980 the daily withdrawals will be 494 billion gal, a figure nearly equal to the dependable supply.Is This is based on a conservative projected population of 230 million. The major categories of withdrawals are: To make available this per capita average of 2150 gal per day will require an expenditure of $219 billion over the next 20 years. The U.S. is not alone in this concern. The United Nations shows as arid zones of the world: all of Africa north of the equator and south of the 20's parallel; all of the Arabian peninsula; all of the middle east and Iran, Iraq, Pakistan, Afghanistan, northern and central India; a great band about 1000 miles wide along the 40'~ parallel from the Caspian Sea east across Russia through China to the Pacific Ocean; all of Australia except the coastal plain; the Caribbean Islands; the western nations of South America; and the western third of the United States and of Mexico. With one quarter of the earth's 57,500,000 sq miles of land thus suffering from lack of good water, increasing attention goes to the treatment of brackish and sea waters. The U.S. has been a leader in this field4,12, 16123,24 through its Office of Saline Water in the Dept. of Interior because even now some of its cities and regions are short of potable water. 11j'7,M Industrial water is also of vital concern as a result of ever higher industrialization1,14122 Other nations, among them JaPan, Israel,13188 Germany, Union of South Africa, Australia, Netherlands, France, Yugoslavia, Russia, and groups such as the Organization for European Economic Cooperation (OEEC)' are also diligent. The objective is low cost water, which means that both technology and economics have prominent roles in saline water conversion processes. TECHNOLOGY: SALINE WATER CONVERSION A number of reviews of methods have been made, principally by staff members of the Office of Saline Water (U.S. Dept. of Interior). Jenkins,31'32 Gillam,34p

Jan 1, 1962

By W. L. Phillips

Stress-train curves were obtained for single crystals of silver, Ag-5 pct Zn, Ag-10 pct Zn, and Ag-20 pct Zn tested in tension and shear at 78°, 195°, and 297°K. At room temperature the critical resolced shear stress gC increased, the length of' the easy-glide region increased, and the rate of' work hardening dwving easy glide decreased with increasixg zinc concentration. The change in the ratio of uc at room temperature to that at lower temperatures was significantly greater for the alloys than for pure silver. It was found that an increment in stress was necessary to continue slip when the slip direction was rotated 60 deg. The magnitude of this increment increased with strain for all alloys, increased with zinc concentration for a given strain, and for a given strain increased with decreasing -temperature. DESPITE its practical importance in improving the mechanical properties, alloying is not fully understood. Except for copper alloys few sets of systematic data are available. Von Goler and Sachs' studied the deformation of Cu-Zn alloys of increasing zinc content and found that, for dilute alloys, the critical resolved shear stress increases linearly with concentration. The range of easy glide was found to increase with increasing zinc content. Schmid and seliger,2 Sachs and Weerts,3 and Osswald4 have shown that with Mg-A1, Au-Ag, and Cu-Ni crystals, respectively, the critical resolved shear stress also varies linearly with concentration. More recently, Linde and his coworkers have investigated the variation of the critical shear stress of copper alloyed with tin , antimony, indium, germanium, silicon, nickel, and gold. They found that the slope of the critical resolved shear stress is related to the change of lattice parameter with composition, and also to the difference in Goldschmidt's atomic diameter between solvent and solute atoms. Garstone, Honey-combe, and creetham6 have shown that similar relationships can be found for small additions of silver, gold, and germanium to pure copper. They found that, with increasing silver or gold concentration, the critical shear stress for glide is increased by alloying, and so is the range of easy glide, which reaches as much as 60 pct for 0.50 pct Ag alloy and 0.62 pct Au alloy, as compared to 6 pct for pure copper of similar initial orientation. They also found that the alloying additions had little effect on the rate of hardening during easy glide, the slope scarcely changing with increasing alloy content. General secondary slip was detected only when the crystals began to harden rapidly. Although the slip appeared to be very fine in the early stages of deformation, coarser slip bands were formed towards the end of the extensive easy-glide range. The present investigation describes the deformation characteristics of single crystals of Ag-Zn containing different concentrations of zinc. Tension and shear testing were used for this study. EXPERIMENTAL PROCEDURE The method of growing the single crystals, sample preparation, and method of testing have been described in detail previously.' EXPERIMENTAL RESULTS A) Tension-Room Temperature. The initial orientations and stress-strain curves of single crystals of silver, Ag-10.0 pct Zn, and Ag-20.0 pct Zn are shown in Fig. 1. It is evident that there is considerable change in the stress-strain characteristics as a function of zinc concentration. The effects of zinc concentration on the critical resolved shear stress for both CU-zn8 and Ag-Zn alloys are summarized in Fig. 2. At all concentrations the resolved shear stress of the Cu-Zn alloys is higher than that of the Ag-Zn alloys. The resolved shear stress increases parabolically as a function of composition for both alloy systems. The length of easy-glide region increased with increasing zinc concentration, Fig. 3b). As the length increased the slope (do/de) decreased slightly, Fig. 3(b). Metallographic investigations demonstrated two significant effects of increasing zinc concentration. First, the amount of clustering increased, compare Figs. 4(a) and (b). The slip lines changed from uniform in pure silver to clustered in the Ag-20 pct Zn and Ag-30 pct Zn alloys. Second, the amount of cross slip decreased as the amount of clustering decreased.

Jan 1, 1968

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 C. S. Matthews, M. J. F. Rosenbaum

The purpose of thir work wax to lcarn it~lzut infori~lation could he obtained from various typs of pilot water floods and to attempt to find the optunum pilot patter11, for a revervoir which had previously been depleted by a solution gas drive. The study was made in the laboratory with mathemetical methods a dynamic analog and a potentiotnetric analog. Results werp tested against the field llistorics of a nrrnlber of pilot water floods. At a reasonable valrre of currzulative injection, the total production rate for the one-injector five-spot should reach about 6.5 per cent of injection rate, and for a four-injector five-spot, about 9 per cent. Accurate estimates of ultimate recovery cannot be made on the basis of such snzall prorluction rates. However, with a pilot composed of nine ir1jector.s and 16 producers the production rate is approximately 50 pcr cent of injection rate at a reasonable value of camulative injection. Sonle inforn~ation for extended performance predictions might he obtained from such a large pilot. These conclusions were drawn on the basis of results obtained for unit mobility ratio, and a sturly using tlre potentiometric analog was made of the effect of other mobility ratios to determine the range of applicability of these predictions. For the four-injector, five-spot pilot with the ratio of production to injection rate (before water breakthrough) is about twice that for with it is about two-thirds; and with M0= 10, it is about one-third For high mobility ratios, it was found that the production rate increased considerably as water-cut increased. These result can be used to modify, qualitatively, the inter.pretntions based on curves for the unit rnobilit\. ratio CaSeS. It was found that the maximum ratio of production rate to injection rate obseriled in field pilot floods was of rhe scime order as that prerdicted by these methods. The time required to reach thisr maximum did not generally agree with the time predicted for a homogeti~orir reservoir. The differcrlce between predicted and observed time of response gives an indication of the permeability profile and of the condition of the producin,g wells. Pilot water floods of the pattern type are generally carried out in reservoirs which have been depleted by solution gas drive and are at low pressure. Under these conditions, oil and water can be considered incompressible. It is assumed that, as the water is injected, an oil bank forms ahead of it and that there is a distinct interface between the water zone (or bank) and the oil zone (or bank) and between the oil zone and the region ahead of the oil zone. It is further assumed that only gas is mobile in the unflooded (gas) region, only oil is mobile in the oil bank and only water is mobile in the water bank. The saturations and the mobilities associated with each zone are assumed uniform. We idealize our reservoir to be homogeneous, horizontal and of constant thickness. Effects of gravity within the producing layer are assumed negligible. If the actual time-dcpendent flow problem is approximated by a acries of steady-state problems. the potential and stream function in the oil bank and water hank satisfy Laplace&apos;s equation in two dimensions. We can therefore use a poteiitiometric analog of this system. Potentiometric models have yielded uscful results in this laboratory&apos; and clsewhere in the study of a variety of secondary recovery problems. For the case where M = I, we generally prefer to use theoretical mcthods as well as a simpler dynamic analog. Except where otherwise noted, the ratio, side of five-spot/wellbore radius. is taken to be 3,600. a figure which corrcsponds to a normal-size wellbore in a 10-acre well spacing. THEORETICAL EVALUATION OF VARIOUS PILOT PATTERNS, Mw0 = 1 <&apos;he theoretical models which we used to examine the performance of various pilots are shown in Fig. 1. Image theory was used to determine the ratio of production rate to injection rate as a function of the volumc of the flood. The ratio of production rate to injection rate was chosen because this is an easily measurable quantity which is characteristic of a pilot

By C. R. Knopp

Gas-oil relative permeability ratio is an important relationship in oil reservoir predictive calculations. A correlation has been developed from 107 gas-flood k/k tests on Venezuelan core samples. The correlating parameter is based on restored-state water saturation tests and&apos; is applicable to both consolidated and poorly consolidated sandstone reservoirs. Data of the correlation show that there are no distinguishah1e differences between the mass-data groupings for the two c1assifications A procedure is recommended for running .sufficient relative. permeability analyses to compute a geometric mean of the sample group. The geometric mean is more representative of the total core, and probably the entire reservoir. For example, while only one in four of the k,,/k,., test curves agreed closely with the resultant correlation of this report, the geometric mean curves of the 16 suites (three samples or more). showed good agreetment ill three cases out of four. INTRODUCTION The gas-oil relative permeability ratio is an important, fundamental relationship in most oil reservoir predictive calculations. Predictive calculations are made to estimate future reservoir production characteristics and ultimate oil recovery. The k1,/k2, relationship is specifically needed to relate the surface gas-oil ratio to the reservoir oil and gas saturation, and to calculate the relative movement of these phases within the reservoir whenever some of the more complex driving mechanisms are present. Laboratory k1/k2, tests are not generally run as a routine analysis. Consequently, k1/k2 data often are not available when needed because the cost of laboratory work could not be justified or the need for such data had not been properly anticipated. When laboratory k1/k2, data are available, they are often very difficult to interpret. For example, wide divergence is sometimes shown in a family of k1,/k1, tests representative of the producing horizon in a single well. With these considerations in mind, a study was made to determine if a relationship might exist between the k1,/k2, curve and some other simple laboratory test criteria. The most probable k1/k2, curve correlation for Venezuela described in this paper is the result of the investigation. The presented correlation defines the most probable gas-flood k,,/k,, curve through the medium of air-water capillary displacement and centrifuge water saturation tests. The laboratory procedures of these tests are. relatively simple, and inexpensive; test data should be. widely available- from routine analysis. DATA AVAILABLE, LABORATORY METHODS The report correlation utilized 107 gas-Hood k1/k2, tests run on sandstone cores of Venezuelan reservoirs. Table 1 is a general tabulation of data pertinent to the tests, while Table 2 summarizes the data. Thetests include 96 from Western Venezuela and 11 from Eastern Venezuela. Eighty-two- of the 107 test samples were sandstones that varied from poorly consolidated to-unconsolidated; 25 were consolidated. The average sample porosity was 26.7 per cent and the average permeability was 1,121 md; these values typify the better sandstone reservoirs of&apos; Venezuela. The Welge gas-flood technique,&apos; based on fundamental Buckley-Leverett frontal displacement theory, was introduced in about 1952 and is widely accepted in the industry. The laboratory procedure is relatively simple, rapid, and can be performed on small core samples. While there have been some minor variations in sample preparation and laboratory procedure in the tests used for the correlation, these tests can be generally summarized as follows. The core sample was first sol vent-extracted and dried. Connate-water saturation was restored by the oil-flushing or evaporation-blow down methods. At the beginning of gas flood the hydrocarbon pore volume was completeiy saturated with the test oil phase. Unsteady-state gas-oil displacement then began with the injection of nitrogen or helium. while the displaced oil and gas phases were incrementally metered at the out-flow face. From the test data, the k,,/k,, curve was calculated by the Welge method.&apos; The individual oil and gas relative permeabilities were also calculated." CORRELATING PROCEDURES In attempting to establish a basis of correlation, we found that broad mid-range sections of 105 of the 107 k,,/k,, test curves could be closely duplicated by a straight line. Only two curves did not show a degree of linearity in this region. Correlation-curve definition parameters were subsequently developed from this observation of consistent mid-range linearity. Possible correlating variables were limited to the physical properties measured on core samples that (1) were widely available as common test data and (2) could be easily and cheaply obtained through future laboratory work. The more obvious possibilities were porosity, permeability and

Jan 1, 1966

By R. Al-Hussainy, P. B. Crawford, H. J. Ramey

The effect of variations of pressure-dependent viscosity and gas law deviation factor on the flow of real gases through porous media has been considered. A rigorous gas flow equation was developed which is a second order, non-linear partial differential equation with variable coeficients. This equation was reduced by a change of variable to a form similar to the diffusivity equation, but with potential-dependent diffusivity. The change of variable can be used as a new pseudo-pressure for gas flow which replaces pressure or pressure-squared as currently applied to gas flow. Substitution of the real gas pseudo-pressure has a nrtmber of important consequences. First, second degree pressure gradient terms which have commonly been neglected under the assumption that the pressure gradient is small everywhere in the flow system, are rigorously handled. Omission of second degree terms leads to verious errors in estimated pressure distributions for tight formations. Second, flow equations in terms of the real gas pseudo-pressure do not contain viscosity or gas law deviation factors, and thus avoid the need for selection of an average pressure to evaluate physical properties. Third, the real gas pseudo-pressure can be determined numerically in term of pseudo-reduced pressures and temperatures from existing physical property correlations to provide generally useful information. The real gas pseudo-pressure was determined by numerical integration and is presented in both tabular and graphical form in this paper. Finally, production of real gas can be correlated in terms of the real gas pseudo-pressure and shown to be similar to liquid flow as described by diffusivity equation solutions. Applications of the real gas pseudo-pressnre to radial flow systems under transient, steady-state or approximate pseudo-steady-state injection or production have been considered. Superposition of the linearized real gas flow solutions to generate variable rate performance was investigated and found satisfactory. This provides justification for pressure build-up testing. It is believed that the concept of the real gas pseudo-pressure will lead to improved interpretation of results of current gas well testing procedures, both steady and unsteady-state in nature, and improved forecasting of gas production. INTRODUCTION In recent years a considerable effort has been directed to the theory of isothermal flow of gases through porous media. The present state of knowledge is far from being fully developed. The difficulty lies in the non-linearity of partial differential equations which describe both real and ideal gas flow. Solutions which are available are approximate analytical solutions, graphical solutions, analogue solutions and numerical solutions. The earliest attempt to solve this problem involved the method of successions of steady states proposed by Muskat.&apos; Approximate analytical solutions&apos; were obtained by linearizing the flow equation for ideal gas to yield a diffusivity-type equation. Such solutions, though widely used and easy to apply to engineering problems, are of limited value bemuse of idealized assumptions and restrictions imposed upon the flow equation. The validity of linearized equations and the conditions under which their solutions apply have not been fully discussed in the literature. Approximate solutions are those of Heatherington et al.. MacRobertsl and Janicek and Katz.&apos; A graphical solution of the linearized equation was given by Cornell and Katz. Also, by using the mean value of the time derivative in the flow equation, Rowan and Clegg&apos; gave several simple approximate solutions. All the solutions were obtained assuming small pressure gradients and constant gas properties. Variation of gas properties with pressure has been neglected because of analytic difficulties. even in approximate analytic solutions. Green and Wilts8 used an electrical network for simulating one-dimensional flow of an ideal gas. Numerical methods using finite difference equations and digital computing techniques have been used extensively for solving both ideal and real gas equations. Aronofsky and Jenkins"I " and Bruce et al.11 gave numerical solutions for linear and radial gas flow. Douglas et al." gave a solution for a square drainage area. Aronofsky 13 included the effect of slippage on ideal gas flow. The most important contribution to the theory of flow of ideal gases through porous media was the conclusion reached by Aronofsky and Jenkins" that solutions for the liquid flow case&apos;" could be used to generate approximate solutions for constant rate production of ideal gases. An equation describing the flow of real gases has been solved for special cases by a number of investigators using numerical methods. Aronofsky and Ferris 10 onsidered linear flow, while Aronofsky and Porter 17 considered radial gas flow. Gas properties were permitted to vary as linear functions of pressure. Recently, CarteP 18 proposed an empirical correlation by which gas well behavior can be estimated from solutions of the diffusivity equation using instantaneour values of pressure-dependent gas

Jan 1, 1967

By Jeff Steinhoff

INTRODUCTION The Henderson mine utilizes a highly mechanized, continuous, panel-caving, mining system to extract ore from a deep, massive, molybdenite deposit. The mine is located 80.5 km (50 miles) west of Denver, Colorado. The mine surface facilities are located 3,170 m (10,400 feet) above sea level in a steep valley on the eastern side of the Continental Divide. Milling facilities are 24 km (15 miles) west on the western side of the divide at an elevation of 2,804 m (9,200 feet) above sea level. The ore- body is located approximately 3,000 feet south of the valley under Red Mountain. Access to the ore for men and materials is through a 915 m (3,000-foot) deep, 8.5 m (28-foot) diameter, vertical, concrete-lined, service shaft. Access from the mill is through a 15.5 km (9.6-mile) rail haulage tunnel. The mine is ventilated through an additional intake shaft and two exhaust shafts. Mine production at this time is 27,255 mtpd (30,000 stpd). The mine ventilation system supplies 1,038 cubic meters per second (2.2 million cfm) through approximately 60 miles of drifting or 2.7 tons of air per ton of ore mined. There are 130 fans in the mine in fixed locations and in vent lines with 6,900 connected horsepower in the mine. MINING METHOD AND LAYOUT The orebody is divided vertically into two major zones. The upper zone is the 8100-level production area. The bottom zone is the 7700- level production area which is in the early development stage. The rail haulage level at 7500 feet is common to both production zones. Each mine production zone consists of five associated sublevels. The cave undercut level is 16.8 m (55 feet) above the production level. Two boundary cutoff levels are located 44.2 m (145 feet) and 62.5 m (205 feet) respectively above the production level. The fresh-air level is positioned 15.2 m (50 feet) below the production level, and the exhaust vent level is 19.8 m (65 feet) below the production level. Horizontally, each production zone is divided into three panels each, 224 m (800 feet) wide. These panels are caved from south to north. As the caving in one panel nears completion, caving in the adjacent panel is initiated. Development for the caving panels is continuous so that the sublevels above the production level and the production level itself have a combination of development drifting and production-related activities. UNDERGROUND VENTILATION NETWORK The ventilation system is zoned in the same manner as the orebody itself. One major split of 600 cubic meters per second (1,270,000 cfm) ventilates the 8100-level production zone; one split of 100 cubic meters per second (210,000 cfm) ventilates the development of the 7700- level production zone; and one split of 165 cubic meters per second (350,000 cfm) ventilates the 7500 rail haulage level. The haulage tunnel requires an additional 188 cubic meters per second (400,000 cfm) of air. Development-drift ventilation is accomplished by hanging 1.0 m (3.5-foot) diameter steel ducting in the drifts with 40-horsepower, 0.96 m (38-inch) diameter fans supplying 9.4 cubic meters per second (20.000 cfm). The normal maximum length for these systems is 300 m (1,000 feet). The 8100-level production-area ventilation system is especially suited to a high level of mucking activity confined in a small area. Approximately 93 per cent of the mine&apos;s total production is transferred from the drawpoints to ore passes in 10 production drifts. The active area in each drift is 300 m (1,000 feet). Twelve 5-cubic-yard LHD units with Cat turbo- charged 170-horsepower engines are assigned to the area. Under these conditions, more than one LHD is assigned to a particular production drift. Adequate ventilation is maintained by making an air change every 97 m (320 feet) along the production drifts. Fresh air is brought into the production drifts from the fresh-air level through 1.37 m (4.5-foot) diameter raises. Air travels south along the production drift to the ore pass where it is exhausted down the ore pass to the exhaust level. The ore pass is followed by another intake which is followed by an ore-pass exhaust. At the south end of the production area, a series of exhaust fans maintain a southerly air- flow through the production level.

Jan 1, 1981

By A. Steiner, K. L. Komarek

Activities of aluminum in solid Ni-A1 alloys have been determined between 20 and 60 at. pet Al and 1200" and 1400°K by an isopiestic method in which nickel specimens, heated in a temperature gradient, are equilibrated with aluminum vabor in a closed all-alumina system. The activity of aluminum shows a strong negative deviation from Raoult&apos;s law at low concentrations but increases by three orders of magnitude within the ß(NiAl) phase. The partial molar enthalpy and entropy of mixing are negative. Using Wagner and Schottky&apos;s theory of ordered compounds, a degree of disorder of 4 x 10 -4 for NiAl and 1.25 X 10-2 for FeAl has been calculated THE Ni-A1 system has been studied by a great number of investigators, and the results, as far as the phase diagram is concerned, have been compiled by Hansen.1 The phase boundaries from 0 to 50 at. pet Ni are well-established. At higher nickel contents the boundaries are still in dispute and an additional phase, A12Ni3, has been reported.&apos; The phase diagram is dominated by a very stable high-melting compound, NiA1, with a relatively wide range of homogeneity. Heats of formation of solid alloys have been determined calorimetrically by Oelsen and Middel3 from 20 to 95 at. pet Ni and by Kubaschewski4 from 25 to 80 at. pet Ni. According to the most recent compilation5 no other thermodynamic investigations have been reported for the Ni-A1 system. Due to the corrosive nature and the low vapor pressure of aluminum, a method has been employed for determining activities of aluminum which was previously developed for the Fe-A1 system.= Nickel specimens, heated in a closed evacuated alumina system in a temperature gradient, were equilibrated with aluminum vapor from a source within the system kept at constant temperature. After complete equilibration the specimens were analyzed and activities calculated from the known vapor pressure of aluminum. APPARATUS AND EXPERIMENTAL PROCEDURE Materials. The nickel specimens were made from wafers of electrolytic nickel (International Nickel Corp.) of 99.99 pet purity which were rolled to a 0.001-in.-thick foil by Driver-Harris Co. and to a 0.005-in.-thick sheet in our laboratory. The aluminum (Aluminum Corp. of America) had a purity of 99.99+ pct. The alumina tubes and crucibles were made of impervious recrystallized alumina with an alumina content of 99.7 pet (Triangle RR, Mor-ganite Inc.). Experimental Procedure. Annular specimens were punched from the sheet, the punching burrs removed, and the specimens degreased in carbon tetrachloride and acetone and weighed on a micro-balance to within an accuracy of ±0.01 mg. The specimens were positioned with alumina spacers along an alumina tube, and the positions measured. Aluminum metal was machined into cylindrical shape, and placed into an alumina crucible. The tube with the specimens was then inserted into a hole drilled into the aluminum metal. An alumina tube with its closed end at the top was slipped over the specimens so that its lower end fitted snugly into the alumina crucible. The assembled reaction tube was inserted into a mullite tube with a water-cooled brass head which had an opening for a quartz thermocouple protection tube and a metal-to-glass connection to a conventional vacuum system. A Pt-Pt 10 pet Rh thermocouple could be raised and lowered in the quartz tube which was placed along the outside of the alumina reaction tube. The mullite tube was heated by two separately controlled resistance-tube furnaces so that in the experimental temperature range an over-all temperature gradient of approximately 150o to 250°C could be imposed on the reaction tube. The position of the mullite tube was adjusted so that the surface of the aluminum metal was always at the temperature minimum. The reaction tube was thoroughly evacuated before and during slowly heating the assembly up to the melting point of aluminum. A pressure of less than 2 µ (Hg) was maintained during an experiment. Once the aluminum had melted, it isolated the contents of the alumina tube from the surroundings. Several times during an experiment the temperature gradient was carefully measured. An experiment lasted from 3 to 6 weeks and it was terminated by air cooling the evacuated mullite tube. For further details of the experimental procedure the paper on the Fe-A1 system6 should be consulted.

Jan 1, 1964

By J. L. Smith, A. L. Albee

This article presents a series of guidelines by which the geologist can evaluate the likelihood of surface faulting and its probable extent at any given site in Southern California and Nevada. The information is intended primarily for geologists concerned with establishing design criteria for proposed nuclear plants. The geologic problems involved in the location of a nuclear power plant are fundamentally no different from those for other types of installations. They fall into four main categories: foundation stability, landslides and slope stability, shaking due to earthquakes, and surface faulting. For problems in the first three categories the foundation engineer, the geologist, and the seismologist can provide criteria for plant location and design, and these problems can generally be economically handled by appropriate design measures for a project of the magnitude of a nuclear plant. However, problems in the last category — surface faulting — are more difficult to handle and require a unique evaluation. Accordingly, this paper will deal primarily with the problem of establishing design criteria for surface faulting, particularly as it affects nuclear facilities. A nuclear reactor is a power source that for greater safety is contained in a heavy, air-tight structure, just as gas, oil, water, and other power sources must be contained. Surface faulting is significant in that it may reduce the integrity of the containment by affecting critical exterior piping or by breaching of the containment. A similar significance exists relative to dams or tanks for the storage of water, gas, or oil, except that in these latter examples the breaching of the container automatically releases the fluids to do their damage. This is not necessarily the case with the rupture of a reactor containment structure because the function of the containment is totally protective, i.e., it is necessary only in the event that radioactive products are released from the reactor, and there are many other safeguards to pre- vent this release even if the containment is ruptured. At the present time, nuclear power plants must be located near large sources of water for cooling the steam generated. The construction of an industrial facility of any kind on the coast line is esthetically distasteful to most people since, unfortunately, there is not enough coast to fulfill all the needs and all the desires of all the people. In most cases where industrial facilities encroach on the lives of citizens, there is no mechanism other than zoning laws by which a person can effectively protest. But in the case of a nuclear facility, the public hearing required by the Atomic Energy Commission provides a forum for dissent, as in recent case histories, and the question of safety provides objectors with a weapon for fighting the construction of the plant. The nature of a public hearing for a nuclear plant is such that the prospective owner and operator must prove that there is no undue hazard, whereas the objectors need only demonstrate that there is a reasonable doubt. It is in this situation that geology becomes the Achilles Heel of nuclear power plant location. For his investigation, the geologist has natural exposures of rock at the ground surface and a limited number of trenches and drill holes to give him a fairly complete picture of the distribution of the various rock types. From the surficial data, he must infer three more dimensions — depth below the surface, past time, and future time. Unlike many problems faced by engineers, the geologist has only this one set of data from which to reach a conclusion — he is unable to reproduce the natural sequence of events in order to obtain another set of data for comparison. Hence, by the very deductive nature of a geologic conclusion, it is difficult to prove a geologic conclusion beyond a reasonable doubt even to other geologists — and perhaps one should say especially to other geologists because the experience and background of a geologist will strongly influence his conclusion, and no two geologists have exactly the same experience and training. The engineer and the public official would like the geologist to conclude that faulting cannot occur at a given site or to assign a numerical value to the probability of its occurrence — but no responsible geologist can do either of these things. Since government officials and others must make decisions that affect the public safety, it would seem that the geologic profession must attempt to establish criteria and

Jan 1, 1968