Search Documents

By S. Finsterle, R. C. Trautz

Numerical modeling studies were performed to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is an important factor affecting the performance of a nuclear waste repository, such as the one proposed at Yucca Mountain, NV. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network, as well as on the capillarity of individual fractures intersecting the opening. A high-resolution numerical model of an unsaturated fracture network was developed, and the appropriateness of seepage predictions using a fracture-continuum model was examined. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2002

By S. Finsterle, R. C. Trautz

We performed numerical modeling studies to investigate water seepage into underground openings excavated in unsaturated fractured rock. Water seepage is a key factor affecting the performance of a nuclear waste repository such as the one proposed at Yucca Mountain, Nevada. The amount of water dripping into an underground opening depends on the connectivity and permeability of the fracture network as well as on the capillarity of individual fractures intersecting the opening. We developed a high-resolution numerical model of an unsaturated fracture network and examined the appropriateness of seepage predictions using a fracture-continuum model. Model calibrations to liquid-release tests performed at Yucca Mountain are also discussed.

Jan 1, 2000

By R. L. McNearny, K. A. Barker

Four numerical models written with the "PFC2D " program were compared against four similar physical models. The numerical models were constrained numerically to closely follow the caving behavior of the physical models. Results from the numerical analyses closely simulated both the granular and blocky behaviors of the physical models, including bridging and interlocking of the simulated orebody during draw. The drawdown patterns and the rate of draw developed during numerical modeling also closely simulated the physical models. The same assumptions and governing parameters used in analyzing the first four models were applied to an additional model. The fifth model employed a completely different, uneven, draw sequence than the other four models, allowing the material to flow in a more granular manner with a demonstrated instability of the undercut sections between each drawpoint opening. Drawdown patterns and behavior of the numerical material of the first four models progressed at uneven and inconsistent rates, similar to those seen in the physical models. By inspection, the overall shape and flow lines of both the numerical and physical models were close in the area removed and in the flow characteristics. The fifth model proved to be the most unstable of all the models tested, with poor correlation to the physical models. The uneven draw sequence, coupled with an additional frictionless wall condition and an added drawpoint, caused the material to flow in a very granular manner. Information generated by this study can be used as a first approach in determining the proper draw rate, drawpoint spacing and undercut thickness for a block-caving operation. Numerical models are often difficult to predict and interpret. However, this study should increase the confidence level for use of the PFC2D program in the future.

Jan 1, 2000

By S. Guo

In this paper, a model for the resin anchor bolt system is developed in which the smooth bar is modeled by using one-dimensional truss elements and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined together and then incorporated into the' existing computer code ROOFBOLT. After this development, all of the components in the resin anchor bolt system including smooth bar, bearing plate and resin annulus can be correctly modeled.

Jan 1, 1995

By J. C. Stankus, S. Guo, S. S. Peng

In this paper, a model is developed for a resin anchor bolt system in which the smooth bar is modeled by using one-dimensional truss elements, and the resin annulus is modeled by using axisymmetric elements. These two submodels are combined and then incorporated into an existing computer code called Roofbolt. After this development, all of the components in the resin anchor bolt system, including smooth bar, bearing plate and resin annulus, can be correctly modeled.

Jan 1, 1997

By Matthew E. Fowler, Jean Habimana, Roy F. Cook

The West Side CSO Tunnel, Pump Station & Pipeline Project is currently beginning construction along the downtown waterfront of Portland, Oregon. This4.3 meter diameter, 5,500 meter-long soft ground tunnel will be driven along the westbank of the Willamette River adjacent to buildings and utilities and several major pile-supported bridges. High groundwater pressure, extended lengths of gravel to cobble-size alluvium, and the need to control ground loss and surface settlements were key considerations in the selection of slurry tunnel boring machines for the project This paper presents an overview of the numerical modeling that was performed during design to estimate the strains and additional loads in the pile-supported structures and to assess the efficiency of proposed mitigation measures.

Jan 1, 2003

By D. G. Fredlund, J. A. Wilson

Conventional disposal of waste rock results in the construction of benches with fine and coarse layers dipping at the angle of repose. This interbedded structure influences flow pathways within the waste rock. A numerical modeling program was conducted to study a four-layer system consisting of alternating fine and coarse waste rock. The system was inclined from the vertical to forty-five degrees with heights varying from one to twenty meters. The results show the effects of gravity as preferential flow developed under unsaturated flow conditions.

Jan 1, 2000

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of the steel mills, which contains around 70% iron and other allied impurities. In the present context, a dynamic prototype counter current reactor (CCR) has been used for experimentation and simulation of particle bed by using discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/ + 270 (53 μm), and then it was oxidized in CCR to make single oxide Fe2O3 phase by maintaining required thermodynamic conditions such as pure oxygen at 1100 °C. Then this oxidized powder was subjected to reduction reaction in CCR having gas mixture ( H2:N2) ratio of 1:4 at 875 ± 5 °C for filling degree (fd) and reactor revolution (ω) ranging from 5.38 to 16.14% and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re ˂ 1) and natural heat convection (Gr/Re2 ˃˃ 1) has been maintained while undergoing solid–gas reaction. Besides particle bed behavior has been quantified by Froude number (Fr) for the optimum operating window. The simulation was validated by the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images. Prevailing heat transfer mechanism and reduction reaction mechanism has been established for solid–gas counter current conditions.

Nov 22, 2021

By S. S. SHINDE, P. C. JIKAR, N. B. DHOKEY

Mill scale is an oxide waste product of steel mills, containing around 70 percent iron and other allied impurities. In the present context, a dynamic prototype counter-current reactor (CCR) was used for the experimentation and simulation of a particle bed using the discrete element method (DEM). Mill scale was pulverized to ASTM mesh range 140 (106 μm)/+ 270 (53 μm) and oxidized in a CCR to make a single oxide Fe2O3 phase by maintaining the required thermodynamic conditions, such as pure oxygen at 1,100 °C. This oxidized powder was then subjected to a reduction reaction in the CCR, with gas mixture (H2:N2) ratio of 1:4 at 875±5 °C for filling degree, fd , and reactor revolution, ω, ranging from 5.38 to 16.14 percent and 2 to 5 rpm, respectively. A creeping fluid flow condition (Re < 1) and natural heat convection (Gr/Re2 >> 1) were maintained during the solid–gas reaction, and the particle-bed behavior was quantified by Froude number for the optimum operating window. The simulation was validated with the experimental results of reactor reaction product by X-ray diffraction (XRD) and scanning electron microscope (SEM) images.

By S. K. Sahoo, S. K. Sharma, U. K. Singh, G. S. P. Singh

This research focuses on the strata control issues in a Bord and Pillar depillaring working under the influence of dead load condition of voluminous and fragmented softcover of overburden dump material. It reports the modeling of a mine working through the FLAC 2D software under the Indian geo-mining conditions. The modeling results of an actual mine working under softcover were compared with that under the intact overburden condition to assess the influence of the soft overburden over the depillaring working. The modeling work involved progressive caving of the strata and the cyclic goaf filling following the modelsimulated occurrence of main fall and periodic caving. Mohr-Coulomb Failure criteria along with a simplified strain-softening (MCSS) model was used to study the failure and caving mechanism. Salamon’s compaction model was used to simulate compaction and resultant stress recovery in the goaf material. Comparison of the model findings was done in terms of load on supports, abutment stress, and convergence, apart from the stress redistribution and failure mechanism. The load transfer mechanism of the soft overburden was also studied for an overall understanding of the strata behavior in such workings.

Jun 16, 2020

By SK SAHOO, UK SINGH, SK SHARMA, GSP SINGH

Surface mining accounts for most coal production in Indian geo-mining conditions compared with underground methods, whose implementations are in deeper coal reserves and for lower stripping ratio in a congested operating space. The bord-and-pillar method is the preferred underground method in India for reasons of flexibility and investment cost. While previous researchers conducted field investigations on intact overburden conditions, this study focuses on the strata control issues in a bord-and-pillar depillaring working under the influence of the dead-load condition of voluminous and fragmented softcover of overburden dump material. The modeling results helped in ascertaining the criticality of strata behavior and rib-pillar stability in actual mining conditions when compared with the conventional intact overburden.

By Juan Andres Jarufe, Carlos Cifuentes

Numerical modelling has been extensively used to estimate plastic displacement in faults and rock mass surrounding mine excavations. Particularly, for fault slip seismicity, modelled plastic displacement is usually utilized to estimate seismic potential, however these methods assume that all the displacement caused during an specific simulated excavation step, is released as a single event at each step. This representation corresponds to a conservative estimation since rock mass displacement does not occur as one single event but as an exponential distribution of many small displacements and a few large ones. When displacement is expressed as seismic magnitude, this exponential distribution is addressed as the Gutenberg Richter (GR) relationship. This paper present a method to decompose the results from numerical models into a distribution of seismic events The reconstruction of the GR relationship from modelling displacement permits not only a more realistic approximation of the expected maximum magnitudes but also provides a tool to develop a probabilistic estimation of the maximum possible magnitudes.

Apr 26, 2022

By A. J. Basu

Geotechnical design plays a key role in mine design. I A proper geotechnical design provides a safe working environment and facilitates efficient mining operation. There are mainly three design approaches I) judgmental, based on observational and experiential facts, 2) empirical, based on empirical design formulae, and 3) analytical, incorporating numerical analysis. Most design applications are based on the design methods one and two as mentioned earlier. However, complexities arise in deep underground mining operations and deep open pit mines. These complexities and difficult geological environment introduce more uncertainties and simple empirical formulae would not be adequate as a design tool. These situations require numerical modelling and analysis of the design problems including interpretation of the results obtained from the analysis. There are several tools available for performing the numerical modelling and analysis. A brief review of the available tools was addressed earlier (Basu and Kuruppu, 1995). The authors selected ANSYS (a general purpose finite element analysis tool) for use in mining geomechanics and geotechnical design. This paper includes several examples demonstrating the applicability of this software in designing structures applicable to mining geomechnics. The authors hope that the design examples will demonstrate the use of the numerical modelling and analysis application in excavation design.

Jan 1, 1996

By W. Keilich

A methodology of subsidence prediction using the Distinct Element code UDEC has been developed as an alternative for subsidence modelling in the Southern Coalfield, New South Wales, Australia. The models have been validated by comparison with empirical results and observed caving behaviour. At this stage, modelling capability is limited to flat lying terrain. It is planned to apply the methodology to areas of high topographical relief to investigate the mechanics of valley closure.

Jan 1, 2006

By Andrzej Walentek

This paper presents the results of numerical modeling of the fractured zone around longwall gateroads. The numerical computations have been based on the results of underground measurements of the fractured zone around longwall gateroads, which were performed for various geological and mining conditions and measured with the use of the endoscopic method with a digital camera. Having results of the underground investigations, trials of duplicating them by the Finite Element Method using software PHASE and the Hoek-Brown failure criterion were performed. In order that the fractured zone predicted by the model-based computations closely simulate that measured underground, some selected parameters in the Hoek-Brown failure criterion were modified based on the earlier determined relations used for prediction of the roof-to-floor convergence in longwall gateroads. As a result, the model-based computations produce results very close to those obtained from the underground measurements. Determination of the fractured-rock zone around the mine roads is of great importance to mining practice, especially in planning and designing supports for longwall gateroads. The influence of the front abutment pressure causes the fractured rock zone around the roads to increase and results in an increase of support load. Moreover, proper determination of the fractured rock zone is one of basic elements in the process of planning roof bolting for mine workings.

Jan 1, 2009

Floor heave events in underground mining often happen suddenly and cause risk to workers and equipment. Reducing the risk of floor heave must be addressed in mining layout design. To improve that design one must understand the significant parameters that are correlated to the occurrence of floor heave. To accomplish that goal, researchers from the National Institute for Occupational Safety and Health (NIOSH) completed a numerical study using FLAC3D, taking advantage of its strain softening and ubiquitous joint capabilities to better simulate the floor heave and caving mechanisms involved at an actual mine site.

Jan 7, 2020

By S. P. Girrens

A coupled finite element geomechanical- hydrology code is currently under development for application to the problem of predicting groundwater disturbances associated with mine subsidence. The structural-fluid coupling is addressed by calculating the subsided mine geometry, with emphasis placed on determining the strata disturbance and locating damaged regions, for input into a hydrology code, which determines localized volume flow rates and aquifer fluctuations. Benefits from coupling will be best realized when field measurements, an additional aspect of the study concurrent with analytical investigations, indicating the relationship between increasing rock strain and increasing permeability are incorporated into hydraulic material descriptions. Hydrologic and struc- tural calculations are presented to demonstrate computational capabilities applicable to mine subsidence

Jan 1, 1982

By W. M. Ma, Daniel W. H. Su, K. Centofanti, Yi Luo, W. L. Zhong, Syd S. Peng

NTRODUCTION Surface subsidence monitoring is the most effective means for investigating the nature of, and developing the techniques for pre¬dicting the surface movement and deformation process. However, because of imperfection in survey equipment and field practices, the surface movement data collected from field survey are not error-free. Also due to lack of proper numerical techniques for computing surface deformation from the measured surface move¬ment data, a relatively large portion of the subsidence data col¬lected could not provide information other than surface movement. In a paper by Luo and Peng (1992), the sources and magnitude of errors in the measured surface movements using an electronic total station in subsidence survey are discussed. If the errors in the measured surface movement components are left untreated before calculation for deformation indices (i.e., slope, strain and curva¬ture), the errors will be carried over to the deformation indices and enlarged if improper numerical procedures are used in the calcu¬lation. It should be noticed that in protecting structures affected by ground subsidence, accurate surface deformation indices are much more important than surface movement. In order to improve the accuracy of the collected surface movement data and the calcu¬lated deformation indices, a number of numerical procedures are also proposed in the same paper. The proposed techniques include: a data smoothing technique for data pre-processing; and a number of numerical methods and/or procedures for calculating surface deformation indices from the measured surface movement data. DATA SMOOTHING Because of various factors, the measured surface movement data always contain errors, small or large. In order to better utilize the collected data, a data smoothing method is proposed to correct the errors in the measured movement data. The method uses the same concept as the influence function method (e.g., Knothe, 1957) on which many existing subsidence prediction methods are based. The data smoothing method makes necessary corrections on the measured movement for the point of interest according to not only the measured movement at this point but also those at the surrounding points. The contribution of a surrounding point in making the correction for the point of interest is inversely propor¬tional, in exponential form, to the distance between them. The new magnitude of a movement component at the point of interest (Zi*) is calculated by

Jan 1, 1992

By H. Huang, J. C. Wu, C. H. Xie, Y. Q. Xue, Z. H. Zhang

Two- and three-dimensional groundwater flow models for dewatering in an opencast mine are presented. The models focus on solving special problems using numerical simulations of the opencast mine. These special problems include the drying up of the aquifer, the movement of the model boundary and the prediction of the residual head. These models are used to describe the decline of groumdwater levels in the Yuanbaoshan Opencast Coal Mine and its neighboring area during the drainage period. The simulated values agree very well with the field data. There have been several successful examples of numerical simulation for predicting the mine dewatering rate in an opencast mine. However, The hydrogeological conditions for the Yuanbaoshan Opencast Coal Mine are rather complex, and the total average flow rate is very high at 4.676 x 105 m3/d (1.65 x 107 cu ft/day). In this paper, a two-dimensional groundwater flow model in a large area is developed to solve the problems mentioned above by simulating regional drawdowns of the groundwater level during the drainage period, and a three-dimensional flow model is developed to simulate the residual head near the open-pit bench.

Jan 1, 1998

By I. R. McGowan, J. Tilk, J. I. Hershman, R. L. Lewis

The White Pine Mine is a large underground copper mine located in the Upper Peninsula of Michigan. A numerical ground water flow model of the White Pine Mine was constructed to evaluate potential engineering controls that would reduce the rate of mine-water discharge after closure. Of the passive controls evaluated, the placement of bulkheads within the mine workings was predicted to offer the greatest benefit, reducing the rate of discharge by as much as 35%.

Jan 1, 1998