Search Documents

By G. W. Tracy, R. D. Carter

A method based on an approach by Hurst' is developed for calculating water influx behavior. Using this method, superposition calculations may be eliminated. The principal diflerence between this method and that of Hurst is that, over finite intervals of time, constant oil-productton rates are assumed by Hurst whereas constant water influx rates are assumed in the presently described method. By adopting the assumption of constant water influx rate for finite time periods, a method is developed for calculating water influx behavior which lends itself to convenient combination with the material balance equation in the Schilthuis form. Above the bubble point, an equation is given for the explicit step-wise calculation of pressure history for a prescribed oil-production history. Use of this method of calculating water influx behavior should result in a considerable saving of time as compared to superposition calculations when the calculations must be made without using a digital computer. The connection is shown between the method described herein and the method developed by van Ever-dingen and Hursr2 (which requires superposition cal-culations). Also, results obtained using the two methods are compared.

By I. Fatt

A recently developed system for measuring electrical resistivity of liquids without use of electrodes offers some interesting possibilities in electric logging technology. The equipment as supplied by the manufacturer is satisfactory for continuous mud logging on a drilling rig or for measuring mud or filtrate resistivity in the laboratory. A simple modification of the commercially available instrument makes it suitable for measuring resistivity of core samples in the laboratory. The continuous measurement of mud resistivity on a drilling rig is a convenient means for detecting mixing of formation water with the drilling mud. Such information is useful to the geologist, the mud engineer and the logging engineer. However, continuous mud resistivity logging by conventional electrode-type resistivity cells is beset with difficulties. The mud, sand and rock chips abrade the electrodes, thereby changing the cell constant and eventually destroying the cell. Also, additives and crude oil in the mud may poison the electrodes by coating them with a nonconductive material. An electrode-type resistivity cell. therefore, may give erroneous readings under certain conditions. Electric logging companies circumvent the electrode poisoning problem by using a four-electrode resistivity cell for measurement of mud resistivity. In this cell, change in electrode area does not change the cell constant. However, the four-electrode cell is difficult to adapt for continuous reading and does not solve completely the problem of electrode abrasion by the sand and cuttings in the mud. The measurement of electric logging parameters on core samples in the laboratory encounters some of the same problems discussed in connection with mud logging. Ideally, the electrical resistivity of a core sample should be measured by placing platinum black electrodes in direct contact with the plane ends of a cylindrical or rectangular sample. Platinum black electrodes however, are much too fragile and easily abraded to be brought in contact with a rock sample. Also, oil or other constituents in the fluid contained in the sample will poison platinum black. In practice, gold-plated brass electrodes, in an AC bridge circuit operating at about 1,000 cps, are used for routine core analysis. For more precise work in research studies, a four-electrode scheme is used.',' Preparation of the samples for the four-electrode method is much too involved for routine core analysis. An apparatus for measuring resistivity of liquids without use of electrodes was described by Guthrie and Boys3 in 1879. They suspended a beaker, containing the electrolyte, by a torsion wire and rotated a set of permanent bar magnets around the vessel. The eddy currents induced in the electrolyte reacted against the rotating magnetic field to develop a torque, which was measured as a deflection of the torsion wire. In 1879 this method could not be made precise or convenient because of the lack of strong permanent magnets. The writer described a very greatly improved apparatus similar to that of Guthrie and Boys, but it was not suitable for continuous measurements or core samples.' Many electrodeless resistivity devices using radio frequency current are described in the literature.5, 6 These generally are suitable only for noting the end-point in a chemical titration. They do not measure resistivity, instead measuring a complex quantity which includes the dielectric constant and the magnetic permeability. The first description of the apparatus to be discussed in this paper was given by Relis.7 Improvements and modifications are described by Fielden,s Gupta and Hills,> and Eichholz and Bettens.10 DESCRIPTION OF APPARATUS The apparatus used in this study is based on the principle that the solution under measurement can form a loop coupling two transformer coils, as shown in Fig. 1. For a fixed AC voltage applied across Coil A, the voltage appearing across Coil B is a function of the resistance of the liquid-filled loop. The details of the voltage generating and measuring circuits are given in Refs. 7, 8, 9 and 10. A block diagram of the equipment is given in Fig. 7. Special features worth mentioning are the operating frequency of 18,000 cps and the automatic temperature compensation which results in the given resistivity readings being automatically correlated to 25°C. The liquid loop supplied by the manufacturer, shown in Figs. 1 and 2, was modified for use in core analysis (Fig. 3). The core sample under test is substituted for a section of the original loop. As shown in Fig. 3, the unit accepts only plastic-mounted cylindrical core specimens. A Hassler-type sleeve easily can be designed for the unit if unmounted cores are to be measured. EXPERIMENTAL PROCEDURE MUD LOGGING A simulated mud line was set up in the laboratory.

By Michael S. Wieland

This report investigates whether it is tractable and thereupon reasonable to remove the booster influences when ranking the toxic fumes of non-cap-sensitive mining explosives. The hazard potential is represented by relative fume toxicity RFT, the resultant influence computed from a formulated rule with a restricted set (tally) of fume components and multiplying constants. The RFT result is compared to the rule criterion that represents the worst case tolerable toxicity stipulated by regulations or otherwise. The typical nonstoichiometric booster composition can render notable unwanted influence on the RFT ranking results. Historically, the remedy was to shoot the booster separately and reduce the trial concentrations by the relevant subtractions. Unreal negative concentration results were noted on rare occasions, revealing the faulty nature of the reduction procedure. With thermodynamic reaction chemistry codes TDRC, on the other hand, the booster ingredients can be retained in or removed from the reactant tabulation, yielding reduction factors taken as the concentration ratios for the two circumstances. The underlying presumption regarding reduction factors is that the ratios of component concentrations with and without the booster ingredients when resolved theoretically ought to equal those rendered experimentally, when tractable. Upon rearranging the presumed relationship, the unadjusted fume concentrations divided by theoretical reduction factors yield the wanted readjusted fume concentrations, thereby removing the booster influence without generating unreal negative concentrations. The work principle from thermodynamics restricted by the constraint of zero-net interaction loss is utilized in conjunction with the TDRC to resolve the wanted (Z-state) fume spectrum. This circumvents the recognized difficulty with traditional hydrodynamic detonation theory that is related to the nonideal reaction characteristics of typical mining explosives. The tacit characteristic of the Z-state resolved from the non-trajectory work principle is that it can not be remedied or tampered with, unlike the traditional trajectory techniques where the reaction products can be readjusted by changing formula representations, nullifying specific reactions or otherwise manipulating the thermodynamic trajectory process. Though the work principle theory is an unfinished formulation, the Z-state results prove worthy enough for rendering trends and comparisons and reaching otherwise useful conclusions. Further refinements should yield more reliable fume spectra and indicate explosive composition readjustments that would reduce toxic effects within underground mines and the wider world environment.

By F. M. Perkins, R. H. Jamison

Publications concerning scaling laws for laboratory flow models of oil reservoirs indicate that the relative permeability and capillary pressure relations must be the same functions of saturation in the model and its prototype. In this paper, the relative permeabilities and saturations are redefined in a way which permits one to have diflerent relative permeability and capillary pressure relations in the model and prototype. The development of the new scaling criteria is demonstrated by deriving scaling laws for one simple flow problem. INTRODUCTION The successful application of information gathered through bench-scale model studies to the prediction of multi-phase fluid flow behavior in petroleum reservoirs requires careful consideration of the scaling criteria. The requirements which must be met for a small-scale model to duplicate the behavior of a much larger reservoir prototype have been discussed by Leverett, et al,' Rapoport,2 and Geertsma, et al.3 According to these authors, the relative permeabilities to water, oil and gas and the dimensionless capillary pressures must be the must have the same values in the model and prototype. Furthermore, they point out that the viscosity ratios must have the same values in the model and prototype. In addition to these requirements, other scaling criteria are given which, when considered, usually dictate a much higher permeability in the model than exists in the reservoir. Since the more permeable model sands may have connate-water and residual oil saturations quite different from those in the prototype, the scaling requirements pertaining to the dependence of relative permeabilities and dimension l ess capillary pressures on saturation usually are not met. Also, if the model fluids are selected to give equal viscosity ratios in the model and prototype, one may find greatly different mobility ratios* in the two cases and, therefore, may expect grossly different flow patterns. For this reason, Craig, et al,' correlated the results of their model studies with mo- bility ratio rather than with viscosity ratio. As yet, however, no satisfactory method has been advanced which permits use of anything other than the scaling requirements outlined in the previously mentioned papers. The purpose of this paper is to present a method of developing a set of scaling criteria which permits different relationships between saturation and relative permeability or capillary pressure in the model and prototype. Rather than attempt to report scaling criteria for a variety of practical situations, this paper demonstrates by use of one simple example how the modified scaling criteria are derived.

By M. J. Sapko

The Pittsburgh Research Laboratory (PRL) of the National Institute for Occupational Safety and Health (NIOSH) has conducted full-scale explosion experiments, to evaluate the strength characteristics of various seal designs, used for safely isolating worked-out areas in underground coal mines. Large-scale explosion tests conducted within the multiple entry section of PRL's Lake Lynn Experimental Mine (LLEM) is currently the only accepted test method used by the Mine Safety and Health Administration (MSHA), for deeming a seal design suitable for use in American mines. These explosion tests are labour intensive, expensive to conduct, and can interfere with other critical underground safety and health research programs conducted by NIOSH. The PRL has developed an alternative seal evaluation method, based on a hydrostatic pressure loading concept that can facilitate the in situ testing of seals in an operating mine. Two chambers within LLEM and one within PRL's Safety Research Coal Mine (SRCM) were used for hydrostatic pressure loading various seal designs. The results from these chamber tests compare favourably with those from the large-scale explosion tests in the multiple entries. In addition to testing seal designs at the required 20 psi static pressure level, the chamber test approach also allows for the determination of the seal's ultimate design strength. Six-scaling relationships for predicting the strength of seal designs as a function of entry size are also presented.

By Daniel T. Eyde

The U.S. Geological Survey (USGS) estimates world production of natural zeolites to be 2.7 Mt (3 million st) in 2013. China ranked ?rst with 2 Mt (2.2 million st), followed by the Republic of Korea at 230 kt (253,000 st) and Turkey at 50 kt (55,000 st). Japan?s production is reported to have increased as a result of zeolites, reportedly low grade local mordenites, being used to treat the low-level radioactive contaminants in the areas surrounding the damaged reactors at the Fukushima Daiichi power plant. China and other countries? high-volume, low-value construction applications, which include dimension stone, pozzolans and lightweight aggregate, are also included in production. This means the production data does not accurately indicate the quantities used in high-value added specialty applications. Consumption in the United States, according to USGS estimates, increased slightly, to an estimated 71.3 kt (78,400 st). Domestic uses for natural zeolites were, in decreasing order by tonnage, animal feed, pet litter, cement (primarily down-hole cement applications by the drilling industry), water purification, odor control, waste water cleanup, fungicide or pesticide carrier, gas absorbent, fertilizer carrier, oil absorbent, desiccant, catalyst and aquaculture. The six leading uses accounted for more than 70 percent of the domestic natural zeolite sales tonnage.

By Dennis Dolinar, Russell Frith, David Oyler, Christopher Mark

Unusual circumstances may require that a longwall retreat into or through a previously driven room. The operation can be completed successfully, but there have been a number of spectacular failures which exposed miners to serious roof fall hazards. To help determine what factors contribute to such failures. the National Institute for Occupational Safety and Health (NIOSH) compiled a comprehensive international database of 130 case histories. The cases include five failures where major rock falls occurred in front of the shields, and six even more serious failures involving major overburden weighting. This suggests two room failure mechanisms The first is a roof fall type failure caused by loading of the immediate roof at the face as the fender narrows. The second is a weighting type failure caused by the inability of the roof to bridge the recovery room and face area, and affecting rock well above the immediate roof The data indicate that the roof fall type of failure is less likely when intensive roof reinforcement (bolts, cables, and trusses) is employed together with higher-capacity shields The overburden weighting failures, in contrast, occurred when the roof was weak and little standing support was used

By Todd M. Ruff

Researchers at the National Institute for Occupational Safety and Health (NIOSH), Spokane Research Laboratory, are investigating technology and methods to monitor the blind areas around large haulage equipment used in surface mines. On average, six fatalities a year can be attributed to equipment colliding with other vehicles or pedestrian workers, or backing over the edge of a dump point. Several off-the-shelf technologies exist for monitoring areas that an equipment operator cannot see from the cab, including radar, video cameras, and radio signal detection systems. NIOSH researchers have tested many of these systems on large dump trucks both at their test facility and at surface mines. This article discusses the technologies available for this application and their effectiveness in detecting objects near an off-highway dump truck. Also, guidance is presented for mine safety personnel who are planning to evaluate collision warning systems or cameras on their haulage equipment. Available Technology Over the past few years, the number and types of systems available for monitoring the blind areas around heavy equipment have increased. Systems marketed for the automotive and recreational vehicle industries have been around for a while. Small video cameras can be found

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By J. B. Agan, R. J. Fernandes

This paper utilizes the layered-system approach, modified to include areal sweep efficiency, to determine the miscible-slug size required and to predict the perfort7zance of a miscible-slug process in a highly stratified reservoir having a permeability variation of 0.76. The slug size in each layer was determined from theoretical calculations and laboratory data which showed that, to maintain miscibility in the reservoir, the slug size had to be 8.2 per cent of the floodable hydrocarbon pore volume. Economic considerations, however, such as the value of the miscible product injected, the cost of injection and the cost of reprocessing the miscible product, limited the amount of miscible product which could be injected into the lower-permeability layers. This resulted in subjecting only 60 per cent of the reservoir to a slug size sufficient to maintain miscibility. Utilizing the relationship of producing gas-oil ratio and areal sweep efficiency as a function of pore volume injected (as obtained from data published by Dyes, Caudle and Erickson2), the theoretical performance was determined. This prediction is compared with actual performance to date. Actual performance reflects an increasing percentage of tniscible product in the produced fluid to a particular level, which remains static for a short period and then increases again to a higher level. This performance suggests that the reservoir has performed as a layered system. It is concluded that the layered-system approach can be used to predict performance in a highly stratified reservoir and that economic consideration is a limiting factor in determining the amount of the reservoir which can be miscibly swept. INTRODUCTION The Paloma zone in the Paloma oil field, Kern County, Calif., is comprised of several sand members separated by siltstone, shale and/or cherty sections. Rapid facies changes within the zone indicate that the sediments were deposited in a "near shore" environment. As a result, the Paloma sands have a high permeability variation, and the Paloma zone is comprised of several sand reservoirs. Since one of the sand reservoirs containing crude oil was not receiving the benefits of gas injection, consideration was given to recompleting certain wells to isolate the sand in order to perform a miscible-phase injection project and, thereby, to increase its ultimate oil recovery. Before proceeding with the project, a review of the literature was made to determine the miscible-slug size that would be required to maintain miscibility throughout the displacement process. Although the review showed that considerable information was available for homogeneous systems, insufficient information was available to determine slug-size requirements and to predict performances for a highly stratified reservoir. This paper has been prepared to describe the method used to determine the miscible-slug-size requirement and to predict the performance of the miscible-phase injection project initiated in the Paloma oil field. Although the Paloma project has only been in operation since Aug., 1958, sufficient production history has been obtained so that a comparison can be made with the predicted performance. DESCRIPTION OF RESERVOIR The reservoir in which the miscible-phase project was initiated in the Paloma field is shown on Figs. 1 and 2. Reservoir properties are as follows: average sand thickness, 43 ft; porosity, 19.1 per cent; water saturation, 47.2 per cent: oil viscosity, 0.328 cp; reservoir temperature, 250°F; reservoir pressure, 3,500 psi: gas saturation, 21.2 per cent at time of miscible product injection; and total net hydrocarbon pore volume, 19,993,000 reservoir bbl. LIMITATIONS The method described herein utilizes a layer-system approach, where injection into and production from each layer is assumed to be proportional to the layer's fractional capacity. The effects of increased fluid mobility on injection and production rates, in layers which have had breakthrough, are not reflected in this method. There has been considerable controversy over the slug size necessary to maintain miscibility throughout the displacement process. In this method, the slug-size

By W. E. Webb, R. H. Church, W. G. Hopkins

The development of systems for location of miners who have become entrapped following a mine disaster has been one of the major research goals in the Bureau of Mines. This study details the preliminary development of such a system. Current technologies for locating trapped miners from the surface are either simple direction, finders or time difference pulse systems that have accuracy Limitations. Developing technologies using full vector field measurements at several receiving points are currently being studied. The necessity of measuring field strength with great accuracy can be a drawback in these systems. A system that operates with continuous wave radio signals, whose phase differences can be measured accurately even in weak signal conditions, greatly improves the location accuracy. A microcomputer would be used, in conjunction with the system, to gather and manipulate the data to automatically pinpoint the location of the trapped miner. Comparison of this system with existing technology indicates a favorable accuracy/complexity ratio for the phase difference of arrival technique.

By J. M. Gaines, C. H. Herty

IRon oxide (FeO) plays an extremely important part in the manufacture of steel. In the open-hearth furnace and the Bessemer converter it is the chemically predominant compound and controls to a large extent both the operation of the furnaces and the composition and quality of their products. Almost all of the refining action is carried out through the medium of FeO, which oxidizes carbon, manganese, silicon and phosphorus. At the end of any open-hearth or Bessemer heat the steel contains dissolved iron oxide, which must be wholly or partly eliminated. The type of deoxidation product depends on the amount of dissolved iron oxide and on the deoxidizer used. After the heat has been deoxidized and while it is being teemed, iron oxide diffuses from the ladle slag into the metal, causing losses of deoxidizer and increase in iron oxide content of the steel. It is generally believed that non-metallic inclusions are responsible for a great many failures of steel in service and for decreased yield of steel during rolling. The basis of most non-metallic inclusions is dissolved iron oxide. This investigation was undertaken in order to determine the amount of iron oxide that will dissolve in iron under various conditions. Steel-making slags are complex and it has been thought best to determine solubility under simple slags as a first step. The study is being continued under slags whose compositions approach those of operating slags Previous Investigations The solubility of oxygen in iron has been determined by a number of investigators (see Table 1) but their effort's have been directed toward the

By Robert H. Peters, George R. Bockosh, Barbara Fotta

This paper was commissioned by the Japan Technical Cooperation Center for Coal Resources Development. It is to be published in a proceedings following a series of presentations to Japanese coal officials in February 1997. The paper focuses primarily on summarizing the past 20 years of U.S. research on promising techniques for ensuring that underground coal miners work safely. These techniques can be grouped into three major categories: management, workplace design, and training. Each technique is described, and the results of efforts to apply the technique at underground coal mines are summarized. Recent statistics on accidents and productivity in U.S. underground coal mining operations are presented before the overview of research findings is given. Because Japanese coal mines are using the longwall method of coal extraction, separate statistics are given for U.S. longwall versus room-and-pillar coal mining operations.

By Robert R. Stein, Roy C. Bartholomae, William W. Aljoe

This Bureau of Mines report describes in detail the design and operation of a test apparatus for measuring the sound power levels of drills used by the mining industry. The two major components of the test apparatus are a computer-controlled automated drill test fixture (ADTF) and a large (45,000-ft3) reverberation chamber that houses the ADTF. Design specifications and performance capabilities of the ADTF and the reverberation room are given. Initial test results for three types of drills--a pneumatic percussion drill, a hydraulic percussion drill, and a pneumatic rotary drill--are given to illustrate the types of experiments that can be conducted with the test appartus.

By Timothy P. Smirnoff

Sustaining the region?s mobility needs, achieving projected economic growth, and maintain Midtown Manhattan as a center of regional, national and global importance has led to continued investment in commuter rail and transit infrastructure in New York City and resulted in design and construction of large underground caverns. The current plans for stations in Manhattan include New York Penn Station Expansion (NYPSE) caverns for Trans Hudson Express (THE) Tunnel Project for New Jersey Transit (NJ Transit) and the Port Authority of New York and New Jersey (PANYNJ), Metropolitan Transportation Authority Capital Construction Long Island Railroad (MTACC/LIRR) twin caverns for East Side Access (ESA) Project, and station caverns of the No. 7 Line Subway Extension Project and Second Avenue Subway Project for Metropolitan Transportation Authority Capital Construction/New York City Transit (MTACC/NYCT). This paper provides a general overview of these large caverns now in design or construction in Manhattan, the area?s geologic setting, the geotechnical factors considered and approaches commonly applied for cavern design.

By Robert E. Mclvor

Classification system efficiency (CSE) is a simple metric that can be used to measure and increase ball milling circuit efficiency. It is defined as the percentage of the ball mill energy being used to grind target (oversize) particles. This paper describes CSE analysis conducted for a variety of operations with related improvements in circuit efficiency ranging from 3 to 30 percent. The achievable CSE for any given ball milling circuit is described, providing the means to determine the extent to which every plant can benefit from increasing grinding circuit CSE.

By Anthony Iannacchione, Christopher Mark

Residents of Maxwell Hill, a suburb of Beckley, West Virginia, reported from January until March of 1988 a series of earth tremors. The tremors rattled dishes, and the noises associated with tremors were sometimes loud enough to wake people at night. These events were puzzling because Beckley is in an area of low seismic activity. Maxwell Hill is located above abandoned coal mines, though no direct subsidence damage to surface structures has been documented. A research program was initiated to determine whether subsurface ground movements associated with the abandoned mines could be responsible for the seismic activity experienced at the surface. A hypothesis is proposed which states that seismic disturbances can be caused by abandoned mine failure when two conditions are met. First, some instability must exist at the mine level. The most likely cause of instability is the failure of pillars of coal that were left for support at the time of mining. Second, subsurface movements can result in significant seismic activity only when the rock above the mine contains at least one very strong and brittle layer that is prone to fail violently through rupture or fault slip movement. Analysis of the data indicates that both hypothesized conditions for mine-related seismicity appear to have been satisfied in the case of Maxwell Hill. It has therefore been concluded that the earth tremors experienced by Maxwell Hill residents probably are related to the abandoned mines. Analysis of the potential magnitude of these events indicates that they are unlikely to cause significant surface danlage. The probability that future events will disturb the residence of Maxwell Hill should be quite low based on the identified mechanism.

By J. A. Zimrner, R. J. Tirnko

The roles of main mine fans in underground mines are to induce airflow and continuously re- move hazardous gases and dust. While most larger mines use multiple fans to accomplish these tasks, many smaller mines employ only a single fan. This paper concentrates on those mines having only one fan to provide ventilation needs. If this fan should cease to function, it is likely that methane concentrations will in- crease in some areas of the mine. As a safety precaution, Federal Regulations require that personnel must be- gin evacuating the mine within fifteen minutes after fan stoppage. Powered haulage can be used to transport workers to the mine portal as long as travel is through areas where hazardous levels of methane are not expected. To determine where methane will accumulate along haulageways during a fan stoppage, air samples must be taken at specific areas and analyzed for hazardous gases. During this study, procedures were developed for measuring methane levels along the haulegeways that are used to exit a mine following stoppage of the main mine fan. Methane levels were measured in four different mines at selected underground locations during four planned fan stoppages. In these mines the measurements taken only at the specified locations showed that, if workers had used electrically powered haulage to exit the mine following a fan stoppage, the methane accumulated along haulageways during fan stoppages would not have significantly increased the risk of an ignition.

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

A system of partial differential equations describing miscible displacement of fluids in porous media is derived The system takes into account the influence of gravity, the spatial distribution of permeability, diffusion, and fluid viscosities and densities. A numerical procedure for approximating solutions to the differential systems has been tested for a horizontal two-dimensional geometry. In end-to-end displacements of oil with less-viscous solvents, the numerical solutions exhibited fingering qualitatively similar to that observed in laboratory models. Small random spatial variations in permeability about the mean value are sufficient to initiate fingering. Quantitative comparisons of computed results with laboratory data show good agreement. INTRODUCTION Miscible flooding of oil by solvent is receiving increasing consideration for field use.1-4 Unfortunately, such floods potentially present severe problems in loss of recovery through the by-passing of oil by expensive solvent. Consequently, their economic evaluation requires sound techniques for predicting recovery. The purposes of this work are to present a finite-difference method for calculating the multidimensional displacement of oil by solvent and to investigate the validity of the method by comparing results of calculations with data from displacements in laboratory models.5 The formulation of the method to simulate the model experiments treats a case of limited scope in the description of solvent flooding in the reservoir. The model experiments were carried out with fluids that were assumed to form an incompressible, ideal, two-component system with constant diffusivity. Establishing the validity of the method even in its present form provides a major step toward the goal of quantitatively evaluating individual solvent flooding projects. First, it demonstrates the fea- sibility of calculating the course of displacements which are dominated by an inherent macroscopic instability, i.e., viscous fingering. In addition, the method provides almost the only practical means of examining the effect of the size and extent of reservoir inhomogeneities on the development and propagation of the fingers, taking into account the important influences of diffusion and gravitational segregation. It might be inferred that such studies cannot yield results of practical value because in displacements dominated by macroscopic instability the pattern of finger development and the resulting performance should depend critically on small variations in the geometric distribution of reservoir inhomogeneities. The present method offers the capability inherent in computational techniques of predicting performance reproducibly with any arbitrary distribution of reservoir properties and, as such, provides a means of evaluating the sensitivity of behavior to uncertainties in reservoir definition and, thus, assessing the reliability of prediction of performance. THE PHYSICAL PROBLEM DATA USED FOR QUANTITATIVE TESTS OF CALCULATIONS In this work, a calculation for treating the miscible-displacement process is tested by comparing calculated results with experimental data. The experiments chosen for comparison have been described in detail,5 and the quantitative features are summarized in a later section. Briefly, oil was flooded by solvent of equal density from a thin rectangular channel in Lucite packed with uniform Ottawa sand. Under these conditions, two-dimensional geometry is considered adequate to represent the process. Inasmuch as the present work is concerned with testing of the method by comparison with specific two-dimensional data, the system of equations presented in the succeeding section will be oriented to defining the physical system specifically. THE DIFFERENTIAL EQUATIONS Suppose x, y to be a Cartesian co-ordinate system and define h(x, y) to be the height of a point above a horizontal reference plane. It will be assumed

By E. F. Johnson, F. H. Brinkman, H. J. Welge, A. L. Hicks

A method is presented for predicting the character of gas or water displacement in a radial system, which can be either horizontal or inclined. The latter case would comprise cone-shaped or dome-shaped symmetry. The method conriders the one dimension of radial distance, taken parallel to the bedding planes. Approximate allowance can be made for segregation of the two fluids perpendicular to the bedding planes. A complete saturation history can be obtained in about a day's time with the use of a desk calculator. Results of an example calculation agree well with the performance of simulated gas drive in a large model representing a section of a cone. This pie-shaped, sand-packed model was 8-ft long and 4-in. thick. The method of solution is useful in handling other partial differential equations of first order and first degree. INTRODUCTION A new analytical method has been developed for predicting the performance of gas or water drive in a reservoir having radial or cone-shaped symmetry (see Fig. 1). In addition, this new development offers a general method of handling partial differential equations of first order and first degree. This work extends those methods developed previously'- "or predicting the displacement of oil by gas or water in linear systems. The method of calculating oil recovery in a radial system involves two major parts. The first part is a calculation of the distribution of saturations at various times of depletion. The second part is the calculation of oil recovery and producing gas-oil or water-oil ratios through volumetric balance on the saturation distributions obtained in the first part. This second part differs only slightly in principle from similar calculations made with linear systems, for instance. The basic assumptions used in the calculations are as follows. 1. One-dimensional radial flow. An approach is explained later for approximating a second vertical dimension. (With- drawals of oil, and of the driving fluid after breakthrough, are assumed uniform around the periphery of the dome.) 2. Constant-volume fluid flow. Thus, in a gas drive there is no change in gas solubility, or in oil or gas volume due to viscous and gravity pressure drops. 3. Capillary forces are negligible. 4. Constant rate of injection. A program has been written for carrying out the calculations needed on a digital computer. However, a prediction can also be arrived at using a desk calculator; this latter method of work requires about a day's time for each prediction made. CALCULATION OF SATURATION DISTRIBUTION For definiteness, we describe the development of our method for the case of gas drive from the crest of the structure. The reader will note, however, that the method is just as applicable to an advancing radial water drive. The equation describing gas flow over an infinitesimal