Search Documents

By E. C. Barfield, D. G. McCarty

Two-dimensional analyses offer considerable promise in providing the basic information required to effect more precise control of petroleum reservoir performance. This paper describes a method for conducting some engineering analyses of this type using a high-qreed digital computer. The general approach is to (I) develop a representative numerical model of the reservoir and (2) employ a suitable numerical technique to solve the basic equation of flow and to develop the required engineering information for the particular case represented by the model. This technique has been used to calculate pattern performance in connection with several field projects involving water injection into oil reservoirs. This type of analysis involves sequentially (I) calculating the pres-' sure distribution in the reservoir and (2) tracking the progress of the interface between the displaced and displacing fluids in a step-wise manner to provide a depletion history of the operation. The particular analysis presented in this discussion is subject to the restriction that the mobility of the displaced and displacing fluids be equal and assumes that the fluids are incompressible and that gravity and capillary pressure do not affect the shape of the flood pattern. Complete two-dimensional .flexibility is maintained with regard to definition of reservoir rock and fluid characteristics, placement of physical restrains and boundary conditions, investigation of flow characteristics in the reservoir, and movement of the displaced and displacing fluid interface. The results of these studies indicate that the highspeed digital computer is well suited for conducting reservoir performance studies in two-space dimensions and that dependable numerical techniques are available for making such analyses. INTRODUCTION One of the problems facing petroleum engineers is that of providing more precise engineering control over reservoir operations at all stages during the depletion cycle. The most direct approach to this problem is to simulate the reservoir with some type of physical or mathematical model and use the information developed from the behavior of the prototype to predict the performance of the actual reservoir. This principle is, of course, not new, and a variety of simulation techniques have been developed to provide information on the behavior of oil and gas reservoirs. One straightforward refinement of present reservoir analysis techniques is the progress from one-dimensional reservoir models to those where the reservoir can be represented in two- and three-space dimensions. Multidimensional reservoir analyses are desired in order to show not only what is happening in the reservoir but in addition to show where in the reservoir the various phenomena of interest are occurring. Thus far, two-dimensional analyses have been confined predominantly to potentiometric and electrolytic model studies. However, when large numbers of wells are involved and where large-scale models are required, such problems become very tedious and instrumentation problems become acute." It is possible through use of numerical methods and general purpose digital computers to avoid many of the instrumentation difficulties and to relax some of the restrictions of previous techniques. The purpose of this paper is to discuss some recent experiences with the use of high-speed digital computers for predicting pattern performance in fluid injection operations and to describe the numerical procedures which have been used for making such analyses. THE NUMERICAL MODEL The problem of representing the reservoir with an appropriate numerical model is composed of at least two important considerations. The first requires that an adequate description of the reservoir mechanism be formulated in the form of the basic differential equation describing the process. The second consideration is definition of the physical characteristics of the reservoir. The latter requires that representative basic data be obtained, examined, and interpreted to provide sufficient detailed definition of the pertinent physical properties throughout the system. The Differential Equation One of the basic assumptions made for the purpose of performing the two-dimensional studies being dis-

By T. Freeman

In the past ten years, the Geotechnical Baseline Report (GBR) has arguably become the key document for tunnel construction. This report not only allocates much of the risk involved with the work, it serves as the basis for bid preparation and is used extensively in resolving disputes during construction. This paper discusses some important issues related to the GBR and presents suggestions for improving these vital reports.

By Michael A. Piepenburg

The East Side Access (ESA) Project will extend Long Island Rail Road (LIRR) into midtown Manhattan via 11 km of tunnels, a new terminal station and associated support facilities. In addition to the hard-rock tunneling currently underway in Manhattan, soft-ground tunnels in Queens will traverse beneath Amtrak?s Sunnyside Yard and Harold Interlocking, the busiest rail interlocking in North America. This paper details ESA?s coordinated approach to interfacing with both Amtrak and LIRR with regard to the planned tunnels, shafts and related construction beneath or adjacent to the operating railroads. This interface includes integration of railroad staff, extensive presentations, reviews and agreements with the railroads regarding settlement mitigation, real-time automated geo-instrumentation monitoring, ground improvement, and geotechnical/environmental controls.

Relatively few longwall mines in the U.S. operate under multiple seam conditions where the two seams are less than 200 ft apart. This paper describes the experience of six that do. These operations are located in PA, WV and UT, and include examples of both undermining and overmining. Some operate above or beneath their own workings, others are in historic mining districts and must contend with abandoned mines that are decades old. The lessons that these mines have learned cover a broad range of topics, including: • Whether to stack gateroads or place them under old gob areas; • How to size pillars and select artificial support to cross longwall stop lines; • How to use yield pillars to minimize multiple seam stresses and coal bump potential, and; • When to anticipate the creation of pathways for gas, water, or oxygen between current and abandoned gobs.

By Michael Schmitz, Eugene Ahrenholtz, William A. Heitbrink, Jeanne A. Zimmer, Andrew B. Cecala, Ernest S. Moyer, Earle H. Andrews, Chris C. Coppock, John A. Organiscak

Many different types of surface mining equipment use enclosed cabs to protect equipment operators from health and safety hazards. The overburden removal and mining process can be extremely dusty and can cause excessive dust exposure. To study this issue, a cooperative research effort was established between the National Institute for Occupational Safety and Health, U.S. Silica Co., Clean Air Filter Co., and Red Dot Corp. in an effort to lower respirable dust levels in an enclosed cab on an older surface drill at a silica sand operation. Throughout this research effort, a number of modifications were incorporated into the drill's filtration and pressurization system, as well as in other areas, to improve its design and performance. An average cab efficiency of 93.4% was determined with gravimetric sampling instruments when comparing the outside with the inside cab dust levels on the final design. Although this study considered just one operation, the goal was to identify cost-effective improvements that could be implemented on all types of enclosed cabs to lower respirable dust concentrations. Two critical components for an effective enclosed cab system are having a properly designed, installed, and maintained filtration and pressurization system, along with a method for maintaining structural cab integrity, which allows the cab to be positively pressurized. Another important component is maintaining cab cleanliness. Although this research was originally directed toward the mining industry, it is also applicable to agricultural or construction equipment.

By W. S. Gatley, F. C. Appl

This paper presents a combined analytical and experimental study of chisel penetration vs time during chisel impact on rock, a problem of fundamental importance in improving the performance of roller-cone bits or percussion drilling toots. For a given force-time relationship between chisel and rock, the problem of determining the penetration (displacement) vs time of the chisel is formidable. This is so because the rock is a nonlinear system with distributed mass and distributed damping (friction, dissipation of energy due to rupture, etc.). Since the literature does not contain adaptable solutions, the rock behavior to impact was simulated approximately by an "equivalent" lumped system, that is, an "equivalent" mass, spring, dashpot system. With this assumption, an analytical solution was found for chisel penetration vs time due to a sinusoidal load between chisel and rock. From this solution were found curves, in terms of dimension-less variables, for the maximum depth of penetration vs the frequency of the sinusoidal loading and for the energy transfer vs frequency. The results of this analysis were used to predict the penetration rate of rotary rock bits us rotary speed. The curve indicated that an optimum speed exists. To verify this analysis, an experimental apparatus was constructed and used to apply a sinusoidal pulse to a chisel penetrating a rock specimen under atmospheric conditions. Strain gauges were mounted on the chisel shank and a velocity transducer was mounted between the chisel and the rock surface. The velocity was integrated electrically and picked up simultaneously with the strain gauge signal on an oscilloscope. Permanent records were made photographically to provide simultaneous records of force us time and penetration us time. In comparing the experimental results for limestone and dolomite with the theoretical results, good agreement was found in the frequency range of the experiments. Unfortunately, the inertia effect (peak penetration) indicated by the theory occurs at a frequency much higher than could be obtained experimentally with the apparatus constructed. A "rate-of-loading" effect is indicated theoretically, but has not yet been verified experimentally. INTRODUCTION The process of drilling with percussion tools or rotary rock bits is basically related to the transient response of rock to surface impact. Each time a bit tooth contacts the rock, high stresses are developed which result in penetration and rock removal. As the tooth moves on, stresses are relieved and a new cycle begins as the next tooth contacts the rock. Thus the drilling process, which consists of an endless succession of these cycles, can be studied in terms of a single cycle. It is apparent, therefore, that the study of single-chisel impact on rock is fundamentally important in improving the performance of roller-cone and percussive-type drills. Previous studies in this area have been conducted by Simon1 and Hartman2 by means of drop tests. In these tests a chisel was attached to a weight and allowed to fall, due to the force of gravity, so that the chisel was driven into a rock specimen upon impact. Strain gauges were attached to the chisel shank and the resulting force-vs-time curves were recorded photographically from an oscilloscope screen. The depth of penetration and crater dimensions were also measured. These tests have provided much valuable information but, as mentioned by the investigators, have not provided complete information on the effect of "rate of loading'! This is partly due to the fact that the chisel motion during drop tests is not a controlled motion which can be varied in form and frequency. Therefore, it seemed that additional information could be obtained by studying chisel impact under conditions where both the motion and the frequency of loading could be controlled. This paper presents a combined analytical and exper-

By R. Unger, L. Mallett

The use of virtual reality (VR) technologies to train miners is an evolving area in research and development. In July 2006, representatives of 14 international organizations met to discuss current VR-related efforts and explore ways to employ VR in the future. This paper gives an overview of the various approaches presented at the meeting and discusses directions for future implementation. An example of a training package that teaches mine map-reading skills is presented. The meeting was the first step to building collaborations for furthering VR applications for mining. Planning has begun for a 2008 meeting open to broader participation. Disclaimer: The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Institute for Occupational Safety and Health.

By T. Seldenrath

In the development of a fluid -operated hammer drill' for accelerated penetration of hard rock formations in oil wells, a research investigation was conducted to evaluate the percussion effects obtained with different design characteristics and to determine the possibilities of percussion in the low frequency range. Tests were conducted on granite blocks and comparable impact forces were measured with a load cell under selected test conditions. These full-scale laboratory tests provided an evaluation of the effectiveness of percussion which was further supported by actual downhole field performance. While much pertinent data have been presented by other investigators the method of evaluation described in this article resulted in good correlation between laboratory and field performance and may be applicable in other low frequency percussion developments. EXPERIMENTAL METHOD For these tests the bit was held stationary and the test block was rotated and forced upward against the bit by a hydraulically-lifted rotary table. "Weight on bit" was calculated from pressure readings on the fluid system of the lifting cylinders—corrected for tare and friction. Percussion was obtained with several tubular hammers of different weights up to 320 lb. These were lifted mechanically and allowed to drop by gravity from various heights onto a slidably-mounted bit sub or "floating anvil" for calibration of the load cell response with known kinetic energy of the hammers. During evaluation of percussion effects, the hammers were hydraulically operated to provide a range of percussion frequencies up to 1,020 blows per minute delivered to the anvil. To obtain comparable test readings, the load cell of the strain gauge type, was attached to the bottom of the anvil or bit sub in place of the bit. The load cell rested on a typical granite test block mounted on the drilling table. The latter was hydraulically supported as in the actual drilling operation. The effect of support resistance was the same for all tools when calibrating. Impact forces so determined werethusdirectly comparable with one another for the purposes of this investigation. When a roller bit was substituted for the load cell for actual drilling, the impact force would probably be reduced due to the resiliency of the bit body, legs. cones and teeth and the penetration of the teeth into the rock. However, since the bits and rock specimens were uniform in these tests, an evaluation of actual tool performance, in terms of parameters previously determined with the load cell at the selected standard conditions of calibration, proved satisfactory. Throughout these tests, 8¾-in. standard tricone, hard rock toothed roller bits and carbide-insert roller bits were used. Rotary speeds were 50 and 100 rpm; static weight on bit was applied up to 40,000 Ib, and percussion frequencies up to 1,020 blows per minute were available. During this investigation it was found that toothed bits on which the teeth had developed flat ends about 3/16 in. wide did not change appreciably for the duration of a test. Hence, all toothed bits used in this study were in this condition. Conventional reference curves obtained in these tests show that increasing weight on hit by 2:1 yielded an increased penetration rate of approximately 2.7:1, while an increase in rotary speed of 2:1 netted an increased penetration rate of roughly 1.7:l. Results are thus of the same order as those reported for hard rock drilling by other investigators." CALIBRATION Under the selected standard test conditions previously described, with the lifting table hydraulically supported, and the anvil interposed between hammer and load cell, the hammers were dropped from known heights. Thus, with the kinetic energy of the hammer known, the corresponding maximum impact force shown by the cell was determined. (It may be of passing interest to note here that special drop tests with the hammers showed that no appreciable reduction in impact force resulted from introduction of the anvil or hit sub.) From these calibrations, it was possible to determine the theoretical kinetic energy available in the hammer from the maximum impact force shown by the load cell while operat-

By Linda L. Chasko, Ronald S. Conti, John D. Cool

Mine operators often rely on mine rescue teams to save lives during an underground emergency such as an underground fire, explosion, roof fall or water innundation. It is extremely important that team members are provided with adequate exploration equipment and that training simulations are conducted in a realistic manner. A series of mine rescue training exercises was developed, conducted, and evaluated by the National Institute for Occupational Safety and Health (NIOSH) in cooperation with the Pennsylvania Department of Environmental Protection, Bureau of Deep Mine Safety and several mining companies. The exercises were conducted at NIOSH’s Lake Lynn Laboratory and operating mines during 1995 to 1999 and resulted in improved technology and training for mine rescue teams. For example, existing technologies were identified to help rescue teams during exploration. These included various chemical light shapes and laser pointers to identify team members and mark crosscuts and mine materials, strobe lights for mapping out escapeways and marking crosscuts, and thermal imaging systems to see through smoke. A hands-free communication system showed potential for enhanced communications between team members and the fresh air base. Anew team lifeline that allows for flexibility of movement between team members was developed. A positive-pressure inflatable escape device (IED) was used to isolate the “hazardous” environment from fresh air and allow rescue team members to traverse through. An inflatable feed-tube partition that can rapidly block large openings such as underground passageways and simultaneously provide a feed-tube for high-expansion foam generators was also deployed for several simulations. During the simulations, mine rescue team members donned self-contained breathing apparatuses (SCBA) and traversed more than 400-m of mine passageways filled with nontoxic smoke. Visibility ranged from 0.3 to 0.9-m. In one area of the mine, a kerosene heater was used to simulate a fire and reduce the oxygen concentration. In another area, a 1.5 pct methane zone was established. The team members searched for "victims,” administered first aid, mapped the passageways, built temporary stoppings, and ventilated the smoke-filled entries. Advanced training included the extinguishment of a conveyor belt fire in a surface fire gallery. Two teams, each carrying charged waterlines, were escorted into the fire gallery through both doors. These teams controlled the rollback smoke and extinguished 7 to 9-m of burning conveyor belt. Team members indicated that the exercises were extremely beneficial due to the realism of thetraining simulations.

By Thomas M. Barczak, Joe Bower, Jinsheng Chen

Pumpable roof supports provide an alternative longwall tailgate roof support and have grown in usage during the past few years. Heintzmann Corporation has been installing pumpable roof supports at the RAG Resources Emerald Mine in western PA since 1998, where they have provided effective roof control in their longwall tailgates. Despite the success of these supports in this application, questions remain regarding critical design issues for optimizing the use of this support technology. The support loading profile is characterized by a high initial stiffness with peak loading occurring at less than one inch of convergence, followed by significant load shedding with a post failure capacity comparable to that of wood cribbing. Therefore, the key to optimizing the support utilization is to provide sufficient load density to prevent convergence from occurring beyond the peak loading capacity. This requires an understanding of how the supports interact with the ground conditions, hence measurement of the ground reaction curve. In order to obtain this information, pumpable roof supports were instrumented to measure support loading. Roof deformation and roof-to-floor convergence measurements in the vicinity of the instrumented supports were also made. The experimental parameters for the installation were the support spacing and water-to-solids ratio of the grout, which controls the grout strength and ultimately the maximum capacity of the support. The study clearly shows that at this mine a 24-inch diameter support is fully capable of providing adequate ground control under depths of cover of 750 ft as only 50 pct of the available support capacity was utilized outby the longwall face. It was also shown that the 2.00 to 1 water-solids ratio, despite providing a slightly softer support response, is sufficient for maintaining the same degree of roof control provided by the traditional 1.75 to 1 grout mix. The 10-ft spacing of the supports did not cause any ground control problems outby the face. However, inby the face the performance of the support is degraded once the peak capacity of the supports is exceeded. It appears that the large load shedding behavior, which is characteristic of this support following peak loading, allows the immediate roof to separate. When this happens, the support is unable to regain control of the roof in time to prevent failure of the immediate roof beam.

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

This paper presents the results of the data obtained in the first stage of a long-range study at high pressures of the system, vapor-hydrate-water rich liquid-hydrocarbon rich liquid. The data presented are for the three-phase systems in which no hydrocarbon liquid exists. Tests were performed on 10 gases at pressures from 1,000 to 10,000 psia. One of these was substantially pure methane, and the remainder were binary mixtures of methane with ethane, propane, iso-butane and normal butane. Several conclusions may be drawn from the data. 1. Contrary to previous extrapolations, the hydrocarbon mixtures tested form straight lines in the range of 6,000 to 10,000 psia which are parallel to the curves for pure methane, when the log of pressure is plotted vs hydrate formation temperature. 2. The hydrate formation temperature may be predicted accurately at pressures from 6,000 to 10,000 psia by using a modified form of the Clapeyron equation. The total hydrate curve may be predicted by using the vapor-solid equilibrium constants of Carson and Katz' to 4,000 psia and joining the two segments with a smooth continuous curve between 4,000 and 6,000 psia. 3. The use of gas specific gravity as a parameter in hydrate correlations is unsatisfactory at elevated pressures. 4. The hydrate crystal lattice is pressure sensitive at elevated pressures. INTRODUCTION Prior to 1950 many studies had been made of the hydrate forming conditions for typical natural gases to pressures of 4,000 psia.""'"'"" Most of these attempted to correlate the log of system pressure vs hydrate formation temperature, with gas specific gravity as a parameter. One of the more promising correlations was made by Katz, et al, which utilized vapor-solid equilibrium constants. The only published data above 4,000 psia are those of Kobayashi and Katz7 for pure methane to a pressure of 11,240 psia. In the intervening years, most published charts for the high-pressure range have represented nothing more than extrapolations of the low-pressure data, with the methane line serving as a general guide. The reliability of these charts has become increasingly doubtful (and critical) in our present technology as we handle more high-pressure systems. The portion of our high-pressure hydrate research program reported here was designed to: (1) investigate the reliability of existing charts; (2) obtain actual data on gas mixtures to 10,000 psia; and (3.) develop a simple hydrate correlation that was more reliable than those which simply used specific gravity as a parameter. Binary mixtures of methane and ethane, propane normal butane, or iso-butane were injected into a high-pressure visual cell containing an excess of distilled water. Hydrates were formed and then melted to observe the decomposition temperature of the hydrates at pressures from 1,000 to 10,000 psia. EQUIPMENT The equipment consisted of a Jerguson 10,000-lb high-pressure visual cell, a 10,000-1b high-pressure blind cell and a Ruska 25,000-1b pressure mercury pump. The visual cell was placed in a constant-temperature water bath controlled by a refrigeration unit and an electric filament heater. A Beckman GC-2 gas chromatograph was used in analyzing the gas mixtures after each run was completed. EXPERIMENTAL PROCEDURE After evacuating the gas system, the heavier hydrocarbon was injected into the high-pressure mixing cell to that pressure necessary to give the desired composition. This cell then was pressured to 1,100 to 1,200 psia by methane from a high-pressure cylinder. The mixing cell holding the gas contained a steel flapper plate and was shaken intermittently over a period of 15 minutes. After mixing, the valve to the high-pressure visual cell containing excess distilled water was opened, and the gas mixture was allowed to flow into the cell. The temperature in the water bath was lowered 10" to 15'F below the estimated hydrate decomposition point. As a first check, the temperature was increased at a rate of 1°F every six minutes to find the approximate point of decomposition. It was again lowered 1.5° to 5°F to form hydrates. The temperature was raised to within l° of the estimated decomposition point and then increased 0.2F every 10 to 15 minutes until the hydrates decomposed. This procedure was repeated at various pressures to obtain 7 to 13 points for each mixture between 1,000 and 10,000 psia. After completion of the hydrate decomposition tests, the gas mixture composition was analyzed with a calibrated gas chromatograph. These gas analyses have an estimated error of ± .1 per cent.

By John R. Eckel

A Iaboratory drilling rig has been devised and placed in operation which permits the application of hydrostatic, terrastatic, and formation pore pressures to a rock sample for drilling under controlled conditions. Calibrated rests of this equipment in-dicate qualitative agreerment with field results on the effect of weight on bit, rotary speed, and rate of circ~~lation on drilling rate. Tests under pressure, using water and air as the drilling fluids, indicate that when drilling saturated limestone samples, a pressure differential between hydrostatic and formation pressure is the only pressure which affects drilling rate. A differential in either direction causes a reduction in rate, but the same diflerential from the wellbore into the formation causes the greater reduction. Also, in lime.stone, air favorably affected drilling rote only when compared to water at high hydrostatic pressurre, all other conditions remaining constant. Large increases in drilling rate with air could be obtained in the laboratory only while drilling a shale, which offered a limited permeability to nitrogen and which could be drilled with a differential from the formation into the wellbore. In general, the magnitude of the changes due to the pressure effect were small compared to changes in drilling rate obtainable by changing either bit weight or rotary speed. INTRODUCTION In the course of Humble's drilling research program, the desirability of investigating several factors which lend themselves to laboratory rather than field study became obvious. Principal reasons for this are the expense involved in field studies and the difficulty of exercising the required degree of control. Of primary importance are investigations of the effect of fluid pressure on rock rate, the effect of pressure on rock drillability, and a study of the fundamental mechanics of rock failure and chip generation. Several years ago some work in the laboratory on the effect of mud properties on drilling rate was re- ported to the API', but the completion of this work and the investigation of the other variables affecting drilling rate required equipment somewhat more elaborate than that used on these earlier tests. This equipment has been assembled to form Humble's Drilling Engineering laboratory. To date only the effect of pressure on rock drillability while using water and air has been investigated, and this work is not yet complete. Facilities are available to investigate the effect of fluid properties on rate of penetration, and information is being assembled on a current basis on the mechanics of drilling. The sequence of presentation in this progress report is somewhat different from that in which the data were obtained since the order of ex-

By Harry Asche

The last 12 months has seen the construction of the first Earth Pressure Balance (EPB) TBM excavated tunnel in New Zealand on the Project Hobson sewer storage tunnel for Watercare Services Limited, with the concurrent excavation of the second EPB tunnel on the Rosedale Outfall project for North Shore City Council. Together these tunnels involve a total length of 6km and follow an 8 year pause in TBM activity in Auckland. TBM performance on these two projects is thus of significant interest as there are four major tunneling projects totaling over NZ$8B being planned for Auckland, in similar ground conditions. Auckland?s bedrock is a weak sandstone/siltstone (East Coast Bays Formation), subject to infilled paleovalleys. The geology in the region has been significantly influenced by volcanic activity; 49 volcanoes are present within 20km of the city centre. The paper describes the design and construction challenges for the projects, including TBM performance, ?sticky spoil? management, and risks associated with basalt lava flows and two paleovalleys. Watercare?s Project Hobson also saw the successful implementation of a Geotechnical Baseline Report?one of the first uses of this type of document in Australasia. Rosedale Outfall is designed for high internal water pressures with varying hydraulic flows, with the risk of exfiltration explicitly addressed.

By E. B. Brauer, B. F. Pennington, E. W. Hough, G. L. Stegemeier

Interfacial tension divided by the difference in density between the liquid and the vapor phases was determined experimentally by the pendant drop method on several isotherms in the two phase region below the critical point for the methane-normal decane system. The density difference data of Sage and Lacey was used in the calculation of interfacial tension. Both interfacial tension and and interfacial tension divided by density difference were found to vanish at the critical point. Interfacial tensions of less than one dyne/centimeter were observed as far as 1,000 pounds per square inch below the critical pressure. EXPERIMENTAL PROCEDURE The interfacial tension divided by the density difference for the methane-normal decane system was determined at the 100°, 130°, 160' and 190°F isotherms, from pressures of about 1,000 psi to the critical pressure, which is more than 5,000 psi for these isotherms. Particular emphasis was placed upon the investigation at pressures slightly below the critical pressure where the interfacial tension is less than 0.5 dyne/cm. Volumetric properties in the two-phase region, including the critical pressures and temperatures, were taken from the work of Sage and Lacey.1 The experimental pendant-drop technique used for the determination of interfacial tensions at high pressures incorporated the ideas of Michaels and Hauser,2 Hough, et a1,3 Walker4 and Heuer.5 In addition, the technique for determination of extremely small interfacial tension by Jennings6 was utilized in the region near the critical points. A detailed description of the apparatus is given in a dissertation by one of the authors.7 Cleaning operations on the stainless-steel sample system included successive washings with chromic acid, tap water and, finally, distilled water. Subsequent cleanings were performed with re-distilled normal pentane, which had an extremely low residue upon evaporation. Specific composition requirements necessitated a fairly precise sample introduction — although, fora two-phase, two-component system, the composition of each phase is completely determined if pressure and temperature are controlled. The normal decane was delivered into the evacuated sample system as a liquid from a burette. The methane was then introduced into the system from a calibrated isothermal container, so that pressure differentials could be used to determine the amount introduced. High pressures were obtained by compressing the sample with a mercury injection pump until the critical pressure was reached for the particular isotherm being studied. Experimental data were then obtained for specific pressures by first decreasing the pressure slightly so that two phases would appear, and then photographing a drop at that pressure. Subsequent photographs were made at increments throughout the pressure range. Calculation of interfacial tension divided by density difference was made from measurements

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the

By Takuji Yamamoto

High-speed TBM excavation requires rapid ground assessment ahead of excavation, particularly in complex grounds. The Tunnel Reflector Tracing technology (TRT) developed by authors, uses seismic reflections to produce volumetric images of ground anomalies. To further minimize interruption in the TBM excavation the authors modified surveying technique switching from cable (wired) to radio (wireless) data acquisition system. This paper presents the outline of the Wireless TRT system, and compares the results for the wired and the wireless surveys. The comparison shows that the wireless system offers high precision data, is more efficient and economical and significantly less intrusive in tunneling operations.

Underground coal mining continues to evolve in the U.S., and more reserves are being mined under deeper cover, with worse roof, or with interactions from previous workings. At the same time, the mining community is responding to higher safety standards and intense competitive pressures. The need for effective ground control design has never been greater. Ground control safety issues that have been addressed by recent the National Institute for Occupational Safety and Health (NIOSH) research include: Improving roof support performance; Maintaining safe tailgate escapeways from longwalls; Optimizing pillar design for retreat mining; Controlling multiple seam interactions; Predicting roof conditions during extended cuts, and; Preventing massive pillar collapses. As funding from both government and the private sector has diminished, the emphasis in research has focused on providing the mining community with practical techniques for improved ground control design. Many projects have successfully employed empirical methods that emphasize the statistical analysis of case histories from underground mines. Other projects have employed numerical models and large-scale laboratory testing of roof support elements. Using these data, NIOSH has developed an entire toolbox of computer programs that have been effectively transferred to the mining community.

By R. J. Blackwell