Search Documents

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

By G. H. F. Gardner

The boundary conditions F = 1.0 and 4 = 1.0 (100 per cent) are dictated by the definition of F in the equation 1/F = aFp. These conditions have physical significance as the "limiting conditions" for the observed behavior in systems of dispersed particles. The other boundary conditions, cited by Gardner in his Eq. 3 as applying to this system, have no significance in the systems we studied. Maxwell's equation (Eq. 2) cited in Gardner's comments was developed for and limited to a regular, rather than random, dispersion of spheres which are sufficiently removed from each other so that they do not interfere with current flow about their neighbors. These restrictions also apply to Rayleigh's analysis. The title of Rayleigh's paper cited by Gardner is "On the Influence of Obstacles Arrayed in Rectangular Order upon the Properties of a Medium". He is not discussing randomly dispersed spheres. It might be inferred from Gardner's comments that Rayleigh specifically considered the differential equation (Gardner's Eq. 3). This is not the case. Mathematically, Gardner's Eq. 3 is exactly true for the system to which it applies, however, as both Rayleigh and Maxwell note that their equation for F in terms of F is not for randomly dispersed spheres. At the time we wrote our paper, we were aware of the work of Wyllie and Gregory as cited by Gardner. We cited their work in Ref. 5 of our paper. Our comment on page 287 of the paper about the validity of an m value of 1.3 for spheres with a porosity of 38 per cent remains as our answer to Gardner's question on this point. Our data, at 38 per cent porosity, show that the condition a = 1 and m = 1.3 holds for spheres near 40 per cent porosity.

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By Gregory Yates

INTRODUCTION Project Description The Silver Line is a Bus Rapid Transit (BRT) Line and is Boston?s newest transit service. Silver Line Phase I extends along Washington Street entirely on the surface in priority bus lanes from Dudley Station in Roxbury through the South End and Chinatown neighborhoods to downtown Boston. The Phase I opened in July 2002. Silver Line Phase II consists of an underground transit tunnel from South Station to the South Boston Waterfront with underground stations at Courthouse and World Trade Center to Silver Line Way east of D Street. Phase II continues through the Ted Williams Tunnel to all the terminals at Logan Airport, and on the surface to the Boston Marine Industrial Park and City Point. Phase II opened in December 2004. Silver Line Phase III is a one mile tunnel connection to the existing Phase II tunnel in the vicinity of South Station and will eventually portal on Tremont Street just north of Marginal Road. After leaving the portal, the buses will use dedicated bus only lanes on the surface roadway network to connect to the Phase I service. Figure 1 shows the alignment of the Silver Line upon completion of the Phase III project.

By R. J. Wygal, J. Naar

An equation for three-phase (water, oil, gas) imbibition oil permeability is developed, assuming the water to be the dominant wetting fluid. Oil isoperms are obtained for consolidated sandstones characterized by 1/pc2 = CS*. The evolution of an oil-gas system imbibing water from Swi to(Sw,imb )max is shown to proceed along a line of constant oil saturation with increasing oil permeability and decreasing gas saturation. When the gas satuation cannot be reduced further, the system evolves along a line of constant Sg with decreasing oil saturation and permeability. The initial gas saturation is shown to reduce markedly the effect of complete wetting by either oil or water on flow performance. INTRODUCTION Imbibition oil isoperms are required for performance prediction when a well is producing water, oil and gas. This situation occurs in multiphase displacements such as underground combustion, steam injection and the water flooding of highly depleted reservoirs. In a recent paper,1 a model was presented for the prediction of two-phase imbibition characteristics. This paper extends the imbibition model to the case of three phases by assuming that the water is the dominant wetting fluid. The following results were obtained from the model: (1) an analytical expression of oil isoperms; (2) oil isoperrns as functions of reduced water, oil and gas saturations, valid for all sandstones having a capillary pressure curve which can be approximated by 1/PC = CS*; and (3) evaluation of the three-phase flow performance as dictated by complete wetting by either oil or water. The agreement between predicted and observed oil recovery in the presence of a gas phase, reported in Ref. 1, is a partial support for the present development. However, experimental data are not available at this time to check fully the model predictions. Perhaps this paper will stimulate the collection of such data. The imbibition model of a porous medium has been described previously, and the reader is referred to the paper of Naar and Henderson1 for details. In brief, the model is formed by the random interconnection of straight capillaries, with a provision for the blocking of the non-wetting phase by the invading wetting fluid. The following expressions were derived from the model. 1. The amount of oil blocked by the invading water is

By R. N. Hoffman

A new technique for determining capillary pressure curves has been developed and tested. The technique differs from previously reported centrifuge techniques in that the centrifuge is slowly accelerated from zero to the maximum desired speed rather than being held constant at particular, progressively higher speeds. An important advantage of this technique over other methods of determining capillary pressure curves is the short time required to obtain the desired amount of data over the chosen pressure range. For example, to obtain 30 data points between 1.2 and 104 psig with a 1.55-in. long, 3/4-in. diameter core using a brine-air system, 6.6 hours were required with this technique. An equally important development of this paper is an analytic method for the conversion of the data from the centrifuge experiment to capillary pressure curve data. Previously there has been only an approximate conversion available. Although the capillary pressure curves determined by this technique appear to be as accurate as those determined by other techniques, the accuracy could be improved if certain variables, not treated in this experiment, were investigated. Among these are the dynamic distortion of the centrifuge equipment and imperfect initial saturation of the cores. INTRODUCTION Pirson 1 defines capillary pressure in a porous medium as "the differential pressure that exists between two fluid phases at their interfaces when they are distributed under static equilibrium within a porous material". Capillary pressure in rocks is known to be a function of fluid saturation, among other things, and a capillary pressure curve is defined for the purposes of this paper as a plot of the capillary pressure-wetting phase saturation relationship for a particular rock sample. Several methods are used for determining capillary pressure curves for small rock cores. Prominent among these are the semi-permeable barrier, mercury injection, and, to a lesser extent, centrifuge methods. The semi-permeable barrier method is currently the most popular. It features simplicity in both execution and the mathematical conversion of the experimental data into a capillary pressure curve. The main disadvantages of the semi-permeable barrier method are the time required — as long as two months — to obtain several points for the curve, and a fairly low maximum pressure before breakthrough of the non-wetting phase into the barrier occurs, for example, about 32 psig for a brine-air system. It is for these reasons that other methods such as the centrifuge method have been introduced. High accelerations and the absence of a barrier result in quicker attainment of saturation equilibrium at a given pressure. However, the centrifuge method involves much more expensive equipment and more difficult procedures and calculations than does the barrier method. The purpose of this investigation has been to improve the equipment and procedures of the centrifuge method and to develop an analytic method for the conversion of the experimental data into a capillary pressure curve. Hassler and Brunner did the original work in the determination of capillary pressure using a centrifuge. In their work the centrifuge speed was increased in a step-wise manner, each speed being held constant until saturation equilibrium was reached in each core. Saturation equilibrium was indicated when the volume of liquid collected in the graduated pipette of the core holder remained constant. According to Hassler and Brunner, equilibrium was reached in "a few minutes to one-half hour or more". In the centrifuge method, as opposed to the barrier method, the fluid saturation of the core is not a constant throughout the length of the core, but varies with the radius of centrifugation. Also, the capillary pressure cannot be read directly but must be calculated from a knowledge of the centrifuge speed and other parameters. For these calculations, Hassler and Brunner derived the following equations: