Search Documents

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

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

By M. Shevchenko

Integrated experimental and modelling research program on the phase equilibria and development of thermodynamic databases of the lead and copper metallurgical gas/ slag-matte-metal-solids systems (PbO-“Cu2O”-FeO-Fe2O3- ZnO-CaO-SiO2) is being undertaken to support improvements in the pyrometallurgical smelting and recycling processes. This is the first systematic investigation of phase equilibria of slag systems in equilibrium with Pb metal, providing information for systems in which copper coexists in slag with lead or zinc as major components. The experimental studies involve high-temperature equilibration of synthetic samples, rapid quenching, and measurement of the compositions of equilibrium phases using electron probe X-ray microanalysis (EPMA). FactSage-based thermodynamic database development is integrated with experimental research. Initial thermodynamic assessments are used to identify priority compositions and conditions for experiments, which are planned to provide specific data to assist thermodynamic optimisation. Significant improvement in the accuracy of the phase equilibria description is achieved. Continuous database improvement and extension to include new elements are targeted. Example of industrial application is given.

By J. W. Walker, G. J. Heuer, E. W. Hough

An improved apparatus is described for the deterrni-nation of interfacial tension* by the pendant drop method in the pressure range from 15 to 15000 psi and the temperature range from 100º td 280°F. The chief characteristics of this apparatus are reliability and ease of operation. Temperature control was accomplished by means of an electrically heated jacket which could he completely removed while cleaning the apptrratus. The cell was designed so that the camera could he placed about I in, from the window without overheating the carnera lens. Separate "0" rings employed to seal each window and corresponding window backup plate gave very reliable window operation at 10,000 psi. An extra pistorr-cylinder combination was added so that a high compression ratio could he attained when gas was employed as one phase in the pendant drop cell. Other essential parts of the system were similar to an apparatus already described.' The interfacial tension for carbon dioxide and water varies from about 69 dynes/cm at 15 psi to about 19 dynes/cm at 1,250 psi for 100°F. For pressure.r between t,500 and 3,500 psi, the interfacial tension is between 32 and 24 dynes/cm at this temperature. The interfacial tension at 280°F varies from about 50 dynes/cm at 200 psi to about 17 dynes /cm at 10,000 psi. The maximum surface excess for carbon dioxide, assuming the surface excess of water to be zero, occurs at about 1,500 psi jor the 100°F isotherni and Itas a value of about -79 X 1014 mol/sq cm. At 160°F this maximum IS shifted to about 2,000 psi and has a value of about 10 X 1014 111ol/.sq cm. Other maxima, us well as several minima, are present at these and other temperatures. Judging from the data on the solubility of water in carbon dioxide at elevtrted pressures2,3, it appears that the principal maxinza natned are in regions of increasing solubility of the water in the carbon dioxide rich phase. INTRODUCTION AND APPARATUS A pendant drop interfacial tension apparatus which would accommodate an opaque liquid in a gas phase was constructed several ycars ago' for industrial purposes. This apparatus was used as a starting point for the present development. The objective of the present development was to enhance reliability and convenience of operation. Other articles'" have described pendant drop apparatus which could not accommodate a gaseous phase adjacent to the windows. A schematic diagram of the apparatus is shown in Fig. 1. The water used followed the path; A, intro-

By Haijian Fan

Currently, deep foundation designs for the service limit state are deterministic in nature. The settlements and lateral deflections of deep foundations are calculated by means of numerical techniques whereby the foundation-soil interactions are characterized by discrete load transfer curves such as t-z or p-y curves. The noted deficiency of the deterministic approach is that the uncertainties arising from various sources, such as soil properties and model errors are not considered. To address this deficiency, a performance-based reliability design methodology is developed using Monte Carlo statistical methods. The performance based design criteria are defined in terms of the specified allowable displacement. The input to the numerical solution techniques, including the soil strength parameters, unit weight of the soil are relevant to the generation of the load transfer curves are treated as random field. The well established random field generator algorithm, local averaging subdivision, is used to generate random samples of input, where the mean values, variances, the probability distribution, and correlation structures are required statistical parameters for these random variables. The objective of the performance based reliability design is to determine the optimized foundation dimensions so that the probability of failure is less than the target probability of failure. The probability of failure by the Monte Carlo approach is simply the ratio of the number of unsatisfactory performance to the sample size. A numerical example for laterally loaded drilled shaft is given to illustrate the application of the methodology. The incorporation of the correlation length is shown to be critical for properly accounting for the site soil spatial variability.

When most people think of coal mining, they envision men in soot-covered clothes and hard hats digging at deposits in deep underground caverns. These mental images commonly include coal cars rumbling on rails as they emerge from the earth, filled to the brim with rocky ore. The world of modern, longwall coal mining is a far different place. Yes, it?s underground, the miners still wear hard hats and the environment is filled with soot. However, today?s mining operations are carefully monitored and the entire process is highly automated. Five kilometers (3 miles) underground, a sophisticated, high-powered mining machine moves steadily through a 2-m (6-ft) high tunnel at an average speed of 3.2 km/h (2 mph), carving out coal, then conveying it to a crusher and carrying it to the surface. Instead of digging out the ore, a longwall mining system slices it from the exposed tunnel face, using giant rotating shears. The shorn coal tumbles onto a conveyer below the shears and is conveyed to the stageloader, where it is crushed and dumped onto the main conveyer belt that eventually hauls it out of the mine. Integral to the longwall shear and conveyor is a moving metal roof that protects it from the thousand plus feet of ground above. The ceiling panels automatically advance above the shears, powered by hydraulic roof supports that stand on the other side of the tunnel, moving in a slow, sideways lock-step with the shears. This complex and massive machine makes longwall mining the most ef?cient and productive method for excavating coal. In fact, it produces 2.7 kt to 22.6 kt (3,000 to 25,000 st) of ore per shift, versus 544 to 907 st (600 to 1,000 st) for the conventional room-and-pillar mining. Using that method, miners develop one section of a mine at a time after carving out load-bearing columns of coal to support the roof above their heads. Longwall mining also is considered safer than the room-and-pillar method, largely because the process is mostly automated, requiring fewer workers. Still, longwall mining has its drawbacks: The equipment is expensive and the mining conditions are harsh. An entire longwall system can cost between $50 million and $100 million, while the electrical controls used to operate it can cost an additional $1 million to $2 million. These investments grind away all day, every day, in hot tunnels filled with coal dust.

By G. D. Boozer, W. H. Somerton

RESEARCH CORP., LA HABRA, CALIF. Thermal diffusivities of some typical sedimentary rocks have been measured by a rapid unsteady-state technique. Thermal data including diffusivity and conductivity are presented for the temperature range of 200" to 1,800°F. Diffusivities calculated from steady-state conductivity rneasurernents compared favorably with unsteady-state data at the lower temperatures. Unsteady-state data are less reliable at higher temperatures; reported values may be as much as 20 per cent high at temperatures in excess of 1,500°F. This is caused partly by the simplified method used in calculating diffusivity from experimental data and partly by the many thermal reactions which occur at the higher temperatures. The marked decrease in thermal diffusivity and conductivity at elevated temperatures indicates the importance of considering the temperature dependence of thermal properties in subsurface heat-transfer calculntions. Knowledge of the effects of liquid saturation, overburden and pore pressures, and heats of reaction of reservoir materials will be required before complete subsurface thermal analyses may be made. INTRODUCTION Heat-transfer calculations are becoming increasingly important in the analysis of in situ combustion and underground nuclear explosions. Thermal data for natural rock materials at elevated temperatures are generally unavailable. Calculations at high temperatures may be grossly inaccurate if temperature-dependent thermal properties of rocks are not considered. It is known that the heat capacities of rocks increase with temperature. Heat capacities may be calculated from chemical or mineralogica1 analyses of the rock system. Thermal bulk expansion of rocks is believed to be small, although few data are available. Thermal conductivities of most natural rock materials decrease with increased temperature. Thermal diffusivity, defined as the thermal conductivity divided by the product of the heat capacity and the bulk density of the material, is used in unsteadv-state heat-transfer calculations. It is apparent that thermal diffusivity decreases markedly with increased temperature. A new method, which gives a nearly continuous record of thermal diffusivity within the temperature range of 200° to 1,800°F and requires less than 1.5 hours for the run,' was used to measure the diffusivity of several sedimentary rocks. Thermal conductivities were calculated from measured bulk densities, heat capacities calculated from chemical analysis and the measured diffusivity values. Thermal conductivities were measured at lower temperatures by a steady-state technique1 to test the agreement of data obtained by the two methods.

By H. D. Hodges, J. K. Godbey

In order to better understand the fracturing process, bottom-hole pressures were measured during a number of typical fracturing operations. A recently developed system was used that allows simultaneous surface recording of both the bottom-hole and wellhead pressures on the same chart. The results from six fracruring treatments are summarized on the basis of the pressure data obtained. Al-though no complete analysis is attempted, the value of accurate pressure measurements is emphasized. Important characteristics of the bottom-hole pressure record do not appear at the wellhead because of the damping effect of the fluid-filled column. In four of the six treatments described, the formations apparently fractured during the initial surge of pressure with only crude oil in the well. The properties of the fluids used during the treatments are given and the fluid friction losses are obtained directly from the pressure records. This technique is also shown to be adequate for determining when various fluids, used during the process, enter the formation. INTRODUCTION Hydraulic fracturing for the purpose of increasing well productivity is now accepted in many areas as a regular completion and workover practice. Numerous articles have appeared in the literature discussing the various techniques and theories of hydraulic fracturing'. In general, three basic types of formation fractures are recognized today. These are the horizontal fracture, the vertical fracture, and fractures along natural planes of weakness in the formation'. Any one or all three of these fracture types may be present in a fracturing operation. However, with only the wellhead pressure record as a guide, it is difficult at best to determine if the formation actually fractured, and is almost impossible to determine the type of fracture induced. These difficulties arise in part because the wellhead pressure record, especially when fracturing through tubing, does not accurately reflect the pressure variations occurring at the formation. Several factors contribute to this effect and preclude the possibility of using the wellhead pressure as a basis for accurately calculating the bottom-hole pressure. These factors are: 1. The compressibilities of the fluids which damp the pressure variations. 2. The changes in the densities of the fluids or apparent densities of the sand-laden fluids. 3. The flowing friction of the various fluids and mixtures, which is dependent on the flow rates and the condition of the tubing, casing, or wellbore. 4. The non-Newtonian characteristics of a sand-oil mixture and its dependence upon the fluid properties, the concentration of sand, and the mesh size used. 5. The unknown and variable temperatures throughout the fluid column. Because of these reasons it was determined that in order to obtain a more accurate knowledge of the nature of fracturing, the bottom-hole pressure must be measured along with the pressure at the surface during a fracturing treatment. Even with accurate pressure data, a reliable estimate of the nature of fracturing is still dependent upon knowledge of the tectonic conditions. However, the hydraulic pressure on the formation is basic to any approach to a complete analysis. In order to accomplish this objective a system was developed to record the wellhead and bottom-hole pressures simultaneously at the surface. By recording both pressures on a dual pen strip-chart recorder, it was possible to greatly expand the time scale so that rapid pressure variations would be faithfully recorded. By such simultaneous recording, time discrepancies inherent in separate records are eliminated, thus overcoming one of the most difficult problems associated with bottom-hole recording systems. This paper illustrates the results obtained by using this system during six typical fracturing operations. All of these tests were taken in wells that were treated through tubing. By a direct comparison of the wellhead and bottom-hole pressures, the importance of obtaining complete pressure information during a fracturing treatment is emphasized. THE INSTRUMENTATION AND PROCEDURES The bottom-hole pressure measuring instrument consisted of a pressure-sensing element, a telemetering section, and a lead-filled weight or sinker bar. The pressure-sensing element used was an isoelastic Amerada pressure-gauge element. By using an isoelastic element, no temperature compensation was necessary in the tests described, since the temperature was believed to be well below the maximum temperature limit of 270°F. The rotary output shaft of this helical Bourdon tube element was coupled to a precision miniature potentiometer. The rotation of the pressure-gauge shaft thus changed the resistance presented by the potentiometer

By A. G. Weber, R. L. Ledbetter, A. L. Hicks

A procedure which makes use of both analog and digital computers has been developed for predicting the pressure-production behavior. of water-drive reservoirs. The electric analyzer is used for matching the past pressure and production history of the reservoir and inferring the physical characteristies of the aquifer. When the history has been matched, a .set of base curves called influence functionce" are derived and the future behavior of the reservoir is calculated from these curves. INTRODUCTION Electric analog methods are particularly well-suited for complex water-drive problems. Like all other methods of analyzing water-drive reservoirs, however, analog methods are subject to certain limitations. One of the more obvious and serious of these results from inaccurate or insufficient knowledge about physical characteristics of the reservoir and aquifer. However, it is not insurmountable because the unknown characteristics can usually be inferred from the reservoir's previous behavior. Simons and Spain have shown that even with scant reservoir data and limited production history, electric analyzer studies can provide useful and reliable results. Aside from this limitation on analyzer predictions, the main disadvantage has been lack of a technique for keeping predictions abreast of latest production forecasts. Pressure behavior of a reservoir depends on the rate at which fluid is withdrawn; it is therefore necessary to specify the future withdrawal schedule before the future pressure behavior can be calculated. Conventionally, the electric analyzer predicts pressure behavior for a specific production schedule. If this schedule is abandoned in favor of some new one, as is frequently the case, either the study must be rerun or results adjusted in some manner. The usual procedure is to predict pressures for several different forecasts covering a wide range of production rates. Pressures for specific rates other than those studied can then be approximated by interpolating between the analyzer predicted curves. However, reliability of this type of prediction is doubtful. This paper describes a technique for predicting future behavior of water-drive reservoirs that retains all advantages of electric analyzer analysis and eliminates the main disadvantage. That is, the analyzer is used conventionally to infer unknown reservoir parameters from the pressure history, but not to predict the future behavior directly. Instead, a derived set of base curves, called influence l'unctions, can be used to calculate future behavior for any production schedule.

By Cliff McCartney, Michael A. Trevits

The National Institute for Occupational Safety and Health (NIOSH) is conducting a program of research to evaluate and improve remote mine seal construction technology for mine fire control and suppression. The main focus of this work is to develop reliable technology that will remotely seal a mine opening from floor-to-roof and rib-to-rib. NIOSH, in partnership with Strata Mine Services Inc.1, is evaluating the capability of a two-component phenolic foam-based product, called Rocsil®1, that is intended for cavity filling, air and gas sealing, and consolidation of highly fractured strata. Rocsil® can be rapidly applied, has a quick cure time, a high expansion rate, good compressive resistance and excellent anti-static properties. When the resin and catalyst are mixed, an immediate foaming reaction occurs followed by rapid expansion up to 30 times the original volume. By virtue of its fire-resistance properties, Rocsil® foam has application for sealing underground mine areas undergoing heating from spontaneous combustion and it could be useful for remote mine seal construction to aid in mine fire control and suppression. This paper presents an overview of remote mine sealing technology, introduces the use of Rocsil® foam for remote mine seal installations, offers a case study of the use of this rigid foam material to remotely construct a temporary in-mine seal and discusses technology improvements. Disclaimer: The findings and conclusions in this report are those of the authors and do not necessarily represent the views of NIOSH.

By R. L. Dalton, L. A. Rapoport, C. W. Carpenter

Procedures have been developed to study the performance characteristics of unconfined pilot water floods using scaled laboratory models. The effects of operating conditions on pilot performance for variozrs well patterns have been investigated. The limitations of the laboratory model studies are outlined. Relationships between the oil recoveries of pilots and large-scale floods have been established. The applications of model studies to the prediction of large-scale waterflood performance and pilot flood design are discussed. INTRODUCTION Pilot water floods are undertaken to evaluate operational procedures and to gain advance information about the performance of an extensive water flood. In general, however, the performance of a pilot flood is not directly indicative of what can be expected from a large-scale pattern water flood. This difference stems from the fact that the over-all flow configurations are not the same for the two systems. A water flood normally consists of an array of identical well patterns. In such an array, the perimeters of the well patterns are axes of symmetry and act as impermeable boundaries. Thus, an extensive, pattern water flood can be visualized as a repetition of "confined" floods. In contrast, a pilot water flood involving only one or a few well patterns is "unconfined". Here the well patterns are not balanced by other flood units; hence, the perimeter of the pilot area does not act as an effective boundary. Accordingly, only a portion of the fluids injected (and also of those originally contained within the pilot area) is captured by the producers of the pilot. The rest escapes into the surrounding reservoir. This usually causes the amounts of oil and water produced in the pilot to be different from those recovered by a producer in an extensive water flood. Moreover, the recovery observed in the pilot can be expected to be greatly influenced by operating conditions, i.e., by the bottom-hole pressures at the injection and producing wells. The above considerations show that full utilization of pilot test results and reliable prediction of large- scale water floods require a knowledge of the relationship between the recovery performance of confined and unconfined well patterns. To date, only limited basic information is available on the behavior of pilot water floods. Craig, et al,' have investigated the effects of initial gas saturation, oil-to-water mobility ratio, and well damage on the oil recovered from a single, "normal" five-spot pilot pattern (Fig. 1). Rosenbaum and Matthews2 studied the effect of initial gas saturation and mobility ratio on the ratio of production rate to injection rate for various groupings of five-spot patterns. Paulse113 investigated the effect of mobility ratio and distance from the pilot to the edge of the reservoir on the areal sweep in a single "inverted" five-spot (Fig. 1). However, none of these studies considered the effect of operating conditions on pilot flood behavior. Moreover, hardly any of the work to date dealt explicitly with the oil-recovery performance of the pilot floods. The purpose of this paper is to introduce the concept of interrelation between pilot flood behavior and operating conditions, to illustrate this interrelation by means of experimental data and to compare the recovery performances of dif-

By R. F. Hinds

A pressurizing system was designed and built to apply a radial pressure of 5.000 psi to rock samples. Samples of the Bradford, Weir and Kirkwood sandstones were subjected to radial pressures parallel and perpendicular to the bedding in the system and the changes in conductivity, porosity and permeability were deter-mined. Asymptotic decreases in conductivity, porosity and Klinkenberg permeability were observed over the 0- to 3,500-psi range. From these data, the increase of formation factor over the same range was calculated. The cementation exponent expressed in the Archie equation was found to increase with pressu,re. The rate of change varied from formation to formation. Only small differences were observed for samples taken parallel to the bedding as comnared to those taken perpendicular to it except for the Bradford sand. Above 3,500 psi the data are not good enough to tell if the properties continue to fall-off asymptotically or in a different manner. A possible explanation for the probable litnit of the asymptotic rate around 3,500 to 4,500 psi is that this represents the maximum pressure which the formations have unrdergone during their geological life as a result of burial. INTRODUCTION For many years measurements have been made in the laboratory of the properties of oil-bearing sandstones with little consideration being given to whether these properties had been changed when the samples were moved from their original position in the ground to the laboratory. Relatively little has been reported in the literature on the effect of pressure on the properties of sandstones. The changes in permeability for eight typical consolidated oil-bearing sandstone samples with overburden pressure (0- to 15,000-psi range) was presented by Fatt and Davis' in graphical form. In 1953 Fatt' also reported the relationship between porosity and overburden pressure for four typical consolidated oil-bearing sandstone samples (0- to 5,000-psi range). Fatt later described the changes of formation factor and porosity of 20 typical oil-bearing sandstone samples with overburden pressure changes (0- to 5,000-psi range) giving the values of the coefficients and exponents in the Archie' equation, and the more general Winsauer equation', where C, k and m are empirical constants, F is the formation factor, and + is the porosity. During the experiments of this investigation the core was contained in a neoprene rubber sleeve as shown in Fig. 1. It was held between two steel end members which were mounted in a steel cylinder. The unassembled pressure cylinder with its component parts is shown in Fig. 2. The core was insulated from the two steel end members by a Teflon sleeve with a small hole in the

By R. V. Ramani, J. M. Mutmansky, P. C. Luna, V. A. Marple, Honglin He, B. A. Olson

The development of a sampling instrument for the assessment of an underground coal miner’s exposure to airborne respirable coal mine dust on a continuous and real-time basis has been long identified as essential for timely initiation of control actions. A respirable dust dosimeter (RDD) which operates on the principle of increasing pressure across a filter media with increasing dust mass on the filter has been developed by the National Institute for Occupational Safety and Health. Comparisons of the performance of the RDD with an U.S. approved gravimetric respirable dust sampler (GRD) have shown that the RDD may meet the requirements for measuring miners’ exposures on a cumulative, near real-time basis such as mid-shift and end of shift exposures. Recently, the performances of the RDD and GRD were compared by conducting side-by-side area sampling experiments in a controlled experimental environment in the laboratory and in underground coal mines under normal operating conditions. In this paper, the results from the in-mine study are presented and discussed. The experimental design for the in-mine study included two mines in different coal seams in different geographic coal regions. To ensure a wide range of concentrations, the experiments were conducted in a continuous miner section in Mine 1, and a longwall face in Mine 2. Airborne dust sampling was conducted for four different shifts at each mine. In each shift, two samples of different exposure durations were collected from two different locations. Each sample consisted of the readings of three dosimeters and the exposed filters from the three paired GRD samplers. As a result, a total of 96 comparative samples were obtained, 48 for each mine. The data were analyzed for the precision of the RDD and GRD measurements and for the relationships between the differential pressure of the RDD and the mass of dust on the paired GRD filter. The analyses indicate that, on the average, the RDD and the GRD track the airborne respirable dust concentrations well, exhibiting both a high correlation between the two measurements and a reliable predictive relationship of one from the other. Details of the experimental design, data analyses procedures, results, conclusions and recommendations of the underground experiments are presented in this paper.

By Richard P. Rothwell

In January last I published in the Engineering and Mining Journal a table giving, with a considerable degree of accuracy, the production of anthracite coal for the year 1874. At that time it mas impossible to give any correct estimate of the production of bituminous coal, for the returns I had been able to secure from such distant and scattered districts were not always forwarded punctually, and the work of .obtaining anything like a true estimate of the coal mined in districts from which I was unable to obtain official returns was both laborious and tedious. As will be seen from the detailed reports furnished herewith, I have succeeded in obtaining tolerably accurate returns of 20,093,704 tons; and upon careful revision by numerous well-informed local authorities, my first estimate of the amount of coal not officially returned to me has grown very considerably; indeed, as I get more information on the subject I find that our bituminous coal production greatly exceeds the amount we have been accustomed to hear it placed at, aggregating as it does, 26,465,704 tons. Mr. Andrew Roy, State Inspector of Mines for Ohio, has reported the production of Ohio at 4,168,445 tons, while his first estimates were but 3,820,000 tons; and similarly, in other cases, the more full the returns direct from the mines the more incorrect do I find all previous reports. I am able now to furnish, not an accurate statement, indeed, but by far the most full and reliable one that has yet been published of our coal production. The amount of labor involved in obtaining this admittedly imperfect report has been simply enormous. In the following tabular statement it will be perceived that Pennsylvania produces nearly one-half of all the bituminous coal mined

By J. Litsenberger, T. E. Marshall, P. R. Coyle, L. J. Prosser, A. T. Iannacchione

For the first time in an underground stone mine, the relationship between roof movement and microseismic emissions was examined in conjunction with two distinct roof fall areas. As roof monitoring increases in acceptance and monitoring technology advances, the goal of providing reliable roof fall detection systems to enhance the safety of underground mine workers moves closer to reality. Instrumental to reaching that goal is the ability to interpret accurately and completely roof movement and microseismic emissions, which can serve as precursors to roof falls. This paper examines the capabilities of convergence and microseismic monitoring systems to better understand roof rock failure mechanics and to anticipate roof falls. An understanding of how these techniques are used and how they interact with each other is critical in developing the most effective ground control strategy for our nation’s mines.

By Yong Pan, Purushotham Tukkaraja

Block caving is a preferred underground mining technique due to its high production rate and low operation cost. Before its final step into a mature caving system using gravity to break rock, drilling and blasting are required to form drawbells to connect the production level and undercut level and induce the initial caving propagation in undercut level. A large amount of toxic fume generated by blasting activities poses a big threat to people working in a confined environment and ventilation is assumed to be the most economical and efficient way to deal with this issue. Past research shows that it’s common to use traditional methods like using a formula or using numerical method (Computational Fluid Dynamic (CFD)) simulation to predict re-entry time after each blasting, but it fails to consider the percentage of toxic fume entrapped in muckpile, which could reach up to 70%. In this paper, CFD was used to investigate the blasting toxic fume characteristics under different entrapping percentages and under changing muck pile properties like broken rock size and porosity. It showed that CFD was a very useful method to study blasting fume in an underground mine and it took much longer time to dilute toxic fume into a safe level with higher entrapping percentage. The dilution time increased with a decrease of broken rock size and porosity.

By B. A. Eaton, K. E. Brown, C. R. Knowles, D. E. Andrews, I. H. Silberberg

This paper presents the resitlts of an investigation of two-phase, gm-liquid flow in horizontal pipelines. Experimental data were taken in three field-size, horizontal pipelines, two of which were constructed for this purpose. The data were obtained using water, distillate and crude oil separately as the liquid phase, and natural gas as the second phase. Experimental pressure-length traverse, liquid holdup and flow-pattern data were collected for each set of flow rates. These data were used to develop three correlations that are presented herein. The liquid-holdup values correlated with various flow parameters without regard to the existing flow pattern. The same was true for the energy-loss factors. A new flow-pattern map is presented that appears to be quite reliable, but not required for the pressure-loss calculations. The liquid-holdup correlation and the energy-loss factor correlation are used in conjunction with a two-phase flow power balance, developed during this study, to predict the pressure losses that occur during gas-liquid flow in horizontal pipelines. A recommended calculational procedure is given, as well as a statistical analysis of the results. This procedure lends itself to computer application, since several small pressure decrements are needed to calculate a pressure-length traverse. The correlations are shown graphically, but may be curve fitted with existing curve-fitring computer programs. INTRODUCTION Due to the frequent occurrence of gas-liquid flow in pipelines and the desire to accurately calculate the pressure losses that occur in these lines, two-phase flow is of considerable interest to the petroleum, chemical and nuclear industries. In the petroleum industry, gas-liquid mixtures have been transported over relatively long distances in a common line due to the advent of centralized gathering and separation systems. Long two-phase flowlines are usually accompanied by large pressure drops which influence the design of the system. Gas-lift installations are designed on the basis of known tubing pressures at the wellheads. The horizontal flowline connecting the wellhead and the separator system must be correctly sized in order to minimize the horizontal flowline pressure losses and the wellhead tubing pressure. Practically all oilwell production design involves horizontal two-phase flow in pipeknes. All of the flow processes of oil and gas production must be studied simultaneously to insure good well design. Since the beginning of offshore oilfield development, long horizontal flowlines have been constructed. Because pressure losses greatly influence the performance of producing wells, a method is desired that can be used to predict such pressure losses and select optimum flowline size. Several types of gas-liquid flow exist, and many of these are discussed by Gouse The study of pressure gradients, fluid distributions and flow patterns that occur in horizontaI multiphase flow is made difficult by the great number of variables involved. The various flow regimes give rise to changing velocities of the fluid particles in all directions. These instabilities of the interface between the gas and liquid prohibit the determination of actual vector velocities of fluid particles in each phase. Also, it is practically impossible to arrive at correct sets of boundary conditions. Therefore, most investigators have concluded that a solution to the problem by the classical fluid dynamics approach, whereby the Navier-Stokes equationsM are formulated and solved, is far too complex. Other methods must be utilized to develop general correlations that will predict the behavior of gas-liquid horizontal flow systems. Multiphase flow studies have sought to develop a technique with which the pressure drop can be calculated. Pressure losses in two-phase, gas-liquid flow are quite different from those encountered in single-phase flow; in most cases an interface exists and the gas slips past the liquid. The interface may be smooth or have varying degrees of roughness, depending on the flow pattern. Therefore, a transfer of energy from the gaseous phase to the liquid phase may take place while energy is lost from the system through the wetting phase at the pipe wall. Such an energy transfer may be either in the form of heat exchange or of acceleration. Since each phase must flow through a smaller area than if it flowed alone, amazingly high pressure losses occur when compared to single-phase flow. Most investigators of horizontal two-phase flow phenomena have chosen to separate their experimental data into several groups of observed flow patterns or regimes.

By D. R. Wilburn

This summary of international mineral exploration activities for the year 2013 draws upon information from industry sources, published literature and U.S. Geological Survey (USGS) specialists. The summary provides data on exploration budgets by region and mineral commodity, identifies significant mineral discoveries and areas of mineral exploration, discusses government programs affecting the mineral exploration industry and presents analyses of exploration activities performed by the mineral industry. Three sources of information are reported and analyzed in this annual review of international exploration for the year 2013: 1) budgetary statistics expressed in U.S. nominal dollars provided by SNL Metals & Mining (SNL) of Halifax, NS, Canada (formerly Metals Economics Group (MEG); 2) regional and site-specific exploration activities that took place in 2013 as compiled by the USGS; and 3) regional events that affected exploration activities, including economic, social and political conditions, that were derived from published sources and unpublished discussions with USGS and industry specialists. Commodity and regional compilations are presented in this summary. Because multiple sources are used to develop commodity and regional compilations, reported statistics may vary depending on the source and type of data that are being reported. The SNL data summarize planned company budgets for worldwide exploration activities in 2013 for 20 nonfuel mineral commodities based on surveys returned by companies primarily focused on precious metals (gold, platinum-group metals and silver), base metals (cobalt, copper, lead, molybdenum, nickel, tin and zinc), bulk commodities (coal, iron ore, phosphate and potash), and specialty commodities (diamond, lithium, niobium, rare-earth elements, tantalum and uranium). Information on uranium exploration activities was included in the MEG overview for the first time in 2007. MEG included data on the commodities lithium, niobium, phosphate, potash, rare-earth elements and tantalum for the first time in 2010 because of the increased topical significance of these commodities. Since 1999, companies with exploration budgets of $100,000 and greater were included in the MEG and SNL compilations. SNL budget estimates exclude aluminum, coal, and oil and gas. Budget estimates for iron ore have been compiled since 2011 but reported separately. The post-1999 budget surveys by MEG and SNL cover at least 90 percent of world?s reported nonferrous, nonfuel mineral exploration budgets; the 2013 survey is reported by SNL to cover an estimated 95 percent of these budgets. The 5 percent that is not covered is accounted for by companies that chose not to participate in the SNL survey, private companies that do not publish their budget data, and government-funded exploration activities.

By Stéphane Bourg

The SCRREEN project gathers more than 50 European initiatives, associations, organisations or independent experts working on Critical Raw Materials (CRM) into a long-lasting network including stakeholders, public authorities, and civil society representatives. SCRREEN contributes to improve the CRM strategy in Europe by (i) mapping primary and secondary resources and substitutes of CRMs, (ii) estimating the expected demand of various CRMs in the future and identifying major trends, (iii) providing policy and technology recommendations for actions improving the production and the potential substitution of CRM, (iv) addressing specific waste electrical and electronic equipment (WEEE) and other End-of-Life (EoL) product issues related to their mapping and treatment standardization, and (vi) identifying the knowledge gathered over previous years and making the data accessible to those beyond the project. As a first step, SCRREEN worked at providing stakeholders with a comprehensive analysis of the current use of the CRM, the mapping of potential primary and secondary resources, and identification of the key technologies used in their production. Substitution profiles, analysis of European business and policy issues pertinent to CRMs, and an identification of the standards, policies, and regulatory frameworks concerning these CRM was also completed. In a second step, to be conducted between mid-2018 and mid-2019, SCRREEN will provide guidance for improving the CRM market in Europe by identifying opportunities from research, development, and innovation (R&D&I) that could reduce the supply and/or economic risks around the most relevant CRM for Europe.

By Joseph E. Chilton, Robert J. Timko, Charles D. Taylor

No regular monitoring of methane is required in areas of gassy coal mines outby the mining face. Methane ignitions that have occurred in mine outby areas indicate the need to provide better protection to workers. Handheld methane monitors are now used by some miners to make periodic measurements of methane at the working face. The IYONI II gas detector which is incorporated into a miner’s cap lamp and worn on a miner’s helmet can continuously provide an alarm signal whenever methane levels exceed a set level. Tests were conducted to evaluate the performance characteristics of this methane detector by measuring response times with methane gas supplied through a calibration fixture or adaptor. Other response time tests were performed with the detector in an environmental test box. Performance was also evaluated in a full scale test gallery where face methane emission and underground ventilation were simulated. Procedures for calibration by response time measurement of the IYONI II detector have been developed. In limited testing, the IYONI II detector was found to reliably detect the presence of 1 percent by volume methane.

By David Coleman

This paper presents a case history involving the installation of drilled in piles with rock sockets for the Bicycle and Pedestrian Path project located west of the Henry Hudson Parkway within the Riverside Park between 83rd Street and 91st Street at the Hudson River in New York, NY. The paper specifically addresses access and construction procedures for a marine installation when working in an urban and environmentally sensitive setting.