Search Documents

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 R. Karl Zipf

A method is presented for creating realistic numerical models for practical coal mine ground control. The method includes procedures to collect the necessary mechanical input parameters from a geologic core log, procedures to set up a model, and procedures to interpret calculation results. The input parameters come from a detailed geologic core log and extensive point load tests of estimate rock layer strength. A suite of material property input parameters is proposed that allows the user to go from core log to numerical model inputs. Rock bolt anchorage properties are also linked to the material properties of each geologic layer in the model. Following this procedure leads to very realistic calculations of the rock failure process and rock support system behavior. These calculations in turn enable realistic comparison of the effectiveness of alternative rock support systems.

Jan 5, 2007

By R. K. Zipf

A method is presented for creating realistic numerical models for practical coal mine ground control. The method includes procedures to collect the necessary mechanical input parameters from a geologic core log, procedures to set up a model and procedures to interpret calculation results. The input parameters come from a detailed geologic core log and extensive point load tests of estimate rock layer strength. A suite of material property input parameters is proposed which allow the user to go from core log to numerical model inputs. Rock bolt anchorage properties are also linked to the material properties of each geologic layer in the model. Following this procedure leads to very realistic calculations of the rock failure process and rock support system behavior. These calculations in turn enable realistic comparison of the effectiveness of alternative rock support systems.

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 A. P. E. Dirige

Extensive field study has been performed to assess the capabilities of thin sprayon liners (TSLs) and conventional spray-on support systems for preventing rock and support material damage that often results due to mining-induced rockbursts. Support performance was studied using field scale explosive detonation trials to simulate dynamic failure effects that are known to develop during typical rockburst events. Multiple seismic and high speed photographic monitoring techniques were used to provide detailed information concerning rock motion, surface fracturing, ejected fragment motion, and support liner survivability characteristics. Results of this study have demonstrated that TSL?s and variant layer combinations may be as effective as or better than conventional support materials for mitigating rockburst or like damage in highly stressed mine environments. In this study, field results have been verified using numerical modeling procedures to better understand the support behaviour of TSLs when subjected to highly stressed mine environments and mining-induced rockbursts, a factor that is essential when designing any rock support system. A series of numerical modeling assessments were conducted using FLAC3D modeling of a half circular TSL-lined tunnel, influenced by anisotropic stresses and simulated mining-induced rockbursts, to investigate the support potential of TSLs under prototype conditions. The effect of a mine rockburst was simulated using the three-dimensional dynamic analysis option of FLAC3D. Results of TSL-lined tunnels in highly stressed and mine rockburst conditions indicated exceptional effectiveness of TSLs for suppressing rock deterioration resulting from stresses and rockbursts. The rock support potential of TSLs on tunnel surfaces was compared with that of thin shotcrete linings, and results indicated that thin spray-on lining products currently available may be equally as effective in support as shotcrete materials.

May 1, 2009

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 Jimin Su, Enguang Zhu, Kangning Zhang, Yang Li

"The effect of interburden on overlying strata movement during the mining of a lower group coal seam is still a challenge in China. Based on the multiple coal seams in western China, numerical simulation and field tests are applied to study the interburden and overburden movement during the undermining process. This study focuses on the impact of the thickness and structure of the interburden, the proportion of hard rock, and the ratio of the overburden thickness to interburden thickness (OB/IB). In this paper, four typical underground coal mines are simulated for lower seam longwall mining with different interburden parameters. The results indicate that the ratio of the overburden thickness to the interburden thickness is not the dominant actor for multiple seam longwall mining, especially in the mines in western China,which feature shallow depth. The thick and hard interburden could strongly affect the lower coal seam mining process. INTRODUCTION The project background is based on four typical underground coal mines in western China.The information on the overburden and interburden is shown in Table 1. Different overburden to interburden thickness (OB/IB) (Ellenberger, 2009) values of the four coal mines are also studied: 6.2 in Bulianta coal mine, 7.2 in Kaida coal mine, 8.2 in Shigetai coal mine, and 9.4 in Halagou coal mine. The proportions of hard rock in the interburden of the four coal mines are also different: 0% in Kaida coal mine, 70% in Shigetai coal mine, 80% in Bulianta coal mine, and 100% in Halagou coal mine. To be specific: the thickness of interburden are listed as follows: 7m in Halagou coal mine, 6m in Kaida coal mine, 10m in Shigetai coal mine, 28m in Bulianta coal mine. The depth of the upper coal seam are: 66m in Halagou coal mine, 54m in Kaida coal mine, 79m in Shigetai coal mine, 182m in Bulianta coal mine. NUMERICAL SIMULATION UDEC 4.0 is applied in this paper to analyse the overlying strata movement and the interburden failure.(Liu, 2007) This is based on mining in the lower coal seam under the mined-out upper coal seam (Majdi, 2012, Morsy, 2006)."

Jan 1, 2017

By Andrew Theobald, Peter Dare-Bryan, Trevor Byers

The use of advanced modelling techniques and the timing flexibility of electronic initiation systems allows for unprecedented control over blast outcomes as shown in several complex blasting problems experienced by the open pit operation Xstrata Ernest Henry Mine (EHM) in north-west Queensland, Australia.

Jan 1, 2010

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 R. Pakalnis, D. Tenney, B. Lang

"This paper describes three case histories whereby numerical methods were employed as a tool in the identification of potential instability. Simple two dimensional (PCBEM) and pseudo three-dimensional (BESOL, MINTAB) codes were employed to identify zones of induced stress which were subsequently related to extensometer, stress meter and observational data as the stope was excavated. The following case histories will be discussed. MacLellan Mine - sill pillar stability; Ruttan Mine-Hanging Wall Stability; and Detour Lake Mine-narrow vein mining.The MacLellan case history showed how modelling and instrumentation enabled mining to be conducted safely and under controlled conditions. MINTAB and PCBEM were employed to predict sill pillar failure. The sill pillar failure was associated, with an unloading of stress which was recorded by instrumentation (stress meter) and a decrease in the over-all rock mass strength. The yielded pillar, however, maintained its over-all configuration enabling access through the pillar.The Ruttan case history is an example of hangingwall and footwall design whereby the existence of adverse structure coupled with a zone of relaxation results in dilution. BESOL and PCBEM were employed to model a test stope and to verify the observations and measurements through instrumentation . Empirical methods of dilution estimates are employed to determine the level of walls lough with respect to the rock mass rating, stope dimensions and rate of excavation.The Detour Lake case history in an example of open stoping a narrow steeply dipping deposit. Stope dimensions are 135m in strike length, 47 m in vertical height with an average width (F/W to H/W) of3.6 m. Stress meters coupled with visual observations are related to modelled estimates of stress and dilution."

Jan 1, 1991

By Ashok Kumar, Rajendra Singh, Petr Waclawik

High-grade coal is locked in the shaft protective pillars at 900 m below surface in the Czech part of Upper Silesian Coal Basin. To optimise recovery of coal from the larger size shaft protective pillar, room and pillar (called bord and pillar in India) is trialled to further develop these on smaller pillars. This development followed a constant gallery width of 5.2 m and 3.5 m height. The monitored pillar size in the first developed panel varied from 860 m2 to 1225 m2 of area. The stability of these pillars (i.e., natural supports) at this depth of cover is observed to be a vital issue to control ground movement and safety of the surface structures. A simple calculation of safety factor revealed that application of an empirical method for estimating size and shape of the pillar is inappropriate due to complex geo-mining conditions of the site. Accordingly, geo-technical instruments are installed in the panel for monitoring the stability of the diamond shaped pillar at different stages of the development. Extensive monitoring of the pillar behaviour by these instruments provided an important set of data for such optimum recovery. An established process of simulation-based pillar strength estimation is often used by the mining industry, but calibration is generally done with the help of empirical formulation. However, for the deep-seated pillar, the strength suggested by the existing empirical formulations (except the CMRI formula) is found to be much different from the instrument observations. Using a study conducted at Czech hard coal mines, an attempt is made to develop a strain-softening-based numerical modelling approach to assess pillar performance at deeper cover. The recorded value and nature of the pillar spalling in the field along with the CMRI formula are used to calibrate the developed simulation approach. An assessment of the stability of the developed pillars using simulated models did not indicate any collapse of the natural supports. This paper briefly presents difficulties in pillar design at great depth of cover and field and simulation results of pillar strength estimation in Czech hard coal mines.

Jan 1, 2018

By J. W. Latilla, J. J. van Wijk

In late 1997 a request was received to assist with determining the powered support unit (shield) load capacities for the planned 2 seam short wall. This project had a number of aspects which made it attractive to utilize numerical modelling in conjunction with other more conventional design tools in the design of the shields. Amongst these were: The 2 seam panels were to be laid out about 25 m below the previously short walled 4 seam, the 4 seam had in turn been over-mined by the 5 seam long wall, the proposed face height of the 2 seam short wall was 5.5 m while the 4 seam had been mined at 4 m. The top seam in the area (5 seam) had been mined at between 1.5 and 2.2 m high, the face conditions with about 95 m of broken material applying dead weight loading to the 2 to 4 seam inter-burden septum were hard to predict and cyclical loading due to the presence of strong sandstone bands in the 2 to 4 seam inter-burden. The basic thrust of the study was to take the known, and successful, situation (4 seam short wall) and carry out comparative modelling to enable prediction of likely 2 seam conditions. Some two dimensional aspects studied were: 5 m goaf hanging up behind shields, stresses in beams to determine likely first goaf roof span, stresses in cantilever beams to asses cyclical loading due to sandstone bands in the interburden being loaded by the overlying goaf and stresses on chain and crush pillars adjacent to the goaf edge. Three dimensional aspects that were investigated included: Load distribution along face for close goafing and 5 m goaf hang-up and load in a 1 m wide rib modelled in the position of the powered support units. From the modelling it was possible to determine the optimal shield leg capacity of the new 2 seam equipment. These outputs were used when drawing up the technical specifications for manufacture of the new equipment. The face was successfully commissioned in August 2002. Fin 1997, requête fut faite pour aider à la détermination des capacités de charge pour l’unité de soutènement mécanisée (bouclier) pour la couche 2 courte taille en projet. Certains aspects de ce projet le rendaient intéressant pour utiliser la modélisation numérique en même temps que d’autres méthodes conventionnelles pour la conception des boucliers. Parmi ces aspects se trouvaient : Les panneaux de la couche 2 devaient être placés environ 25 m au dessous de la couche 4 courte taille ; la couche 4 avait été remplacée par la couche 5 longue taille ; la hauteur de taille proposée pour la couche 2 courte taille était 5.5 m, tandis que la couche 4 avait été exploitée à 4 m. La couche supérieure de la région (couche 5) avait été exploitée entre 1.5 et 2.2 m de hauteur ; il était difficile de prédire les conditions de taille, avec environ 95m de matériaux cassés ajoutant une charge morte à la plaque de séparation entre les couches 2 et 4 et un chargement cyclique causé par des bandes de grès dans l’entre-fardeau 2-4. L’idée principale de l’étude fut de prendre la situation connue et réussie (couche 4 courte taille) et de faire une modélisation comparative pour prédire les conditions probables de la couche 2. Certains des aspects à double dimension qui furent étudiés étaient : Un accrochage de remblai de 5m derrière les boucliers ; contraintes dans les poutres pour déterminer le premier intervalle probable entre remblai et toit ; contraintes dans les poutres en porte-à-faux pour évaluer le chargement cyclique causé par les bandes de grès dans l’entre-fardeau qui est accablé par le remblai et les contraintes sur les piliers par chaîne et d’écrasement adjacents au bord du remblai. Aspects à trois dimensions qui furent étudiés : Distribution de la charge le long de la taille pour les remblais proches et l’accrochage de remblai de 5 m et la charge dans une nervure de 1 m de large modelée dans la position des unités de soutènement mécanisées. La modélisation a permis de déterminer la capacité de portée optimale du bouclier du dispositif nouveau pour la couche 2. Ces résultats furent utilisés pour dresser les spécifications techniques pour la fabrication de l’équipement nouveau. La taille fut commissionnée avec succès en août

Jan 1, 2003

By M. W. Grabinsky, M. S. Diederichs, E. Hoek

Provided that they are used correctly, numerical models can be of considerable assistance to mining engineers in designing underground excavations. Guidance is given to the engineer or geologist who is not a specialist in numerical analysis but who is interested in using numerical modelling tools or in employing someone to analyse the stresses, displacements and failure in the rock surrounding underground excavations. A description of the theoretical background to a number of these numerical models and discussion of the advantages and disadvantages of each are followed by guidelines on the selection of appropriate models for each stage in the development of a mine, from exploration to the completion of mining

Jan 1, 1991

By Kuruppu M

Geotechnical design plays a key role in mine design. 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, (ii) empirical, based on empirical design formulae; and (iii) 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 is provided. 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 mining structures. The authors hope that the design examples will demonstrate the use of the numerical modelling and analysis applications in excavation design.

Jan 1, 1995

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 O. K. Mahabadi

The scope of this study is to simulate the behaviour of a layered rock sample under standard laboratory Brazilian Disc test using an innovative combined finite-discrete element method (FEM/DEM) research code. The influence of layering and the direction of loading on the samples behaviour are studied through various layering and loading configurations. A parametric study including the effect of loading rate on the tensile strength of samples has been carried out. This paper demonstrates the suitability of the numerical approach to explicitly model rock deformation and failure.

May 1, 2009

By Lok M. Sharma, Deep K. Khatri, Bradfield. Brett

"A shored enclosure was required to access a control structure upstream below the dam of a deep lake. The shoring involves the interaction of the structure, soil, and water and the loads imposed include the hydrostatic forces of the water, soil pressures as well as dynamic forces due to currents and waves. Deep excavation support systems have traditionally been designed using classical methods that involve direct computations of earth pressures. However, traditional analysis does not take into account soil-structure interaction. Numerical modelling offers better means to analyze soil-structure interaction. This paper presents a 2D finite element (2D-FE) analysis of a deep excavation support structure that was required to remove and replace a control structure at the upstream toe of a dam. The sheet-piling shoring was designed to withstand the lateral pressures from 14 m deep free water and 5.5 m deep excavation below the free water to reach the required excavation level to replace the gate structure. Elastic-plastic analyses were performed to predict the stresses and deformations of the complicated shoring structure and excavation bottom. The soil parameters, procedures used to develop the 2D-FE model, and the results are discussed in this paper. The shoring system is currently under construction.INTRODUCTIONA deep shoring system enclosure is required to keep out water and soil and to provide a support for unbalanced loadings due to dewatering and/or, excavation in order to permit construction of the permanent structures in the dry. A braced shoring system consists of a single wall sheet piling which is driven into the ground to form a “box” around the excavation site and is braced on the inside. After bracing the sheet piling, seal concrete is placed into the bottom to prevent water from seeping into the open excavation from underneath the sheet piling, and the water is pumped out to provide the dry work area. The braced shoring involves the interaction of the structure, soil, and water and the loads imposed include the hydrostatic forces of the water, soil pressures as well as dynamic forces due to currents and waves. Deep excavation support systems have traditionally been designed using classical methods that involve direct computations of earth pressures, (active, passive, or apparent pressures). There have been some uncertainties regarding the actual distribution of earth pressures due to influence of the support wall flexibility, as these traditional methods do not take into account soil-structure interaction. Numerical modelling offers better means to analyze soilstructure interaction."

Jan 1, 2017

By Medhurst TP, Trueman R

The aim of the paper was to model numerically some aspects of the extraction of coal pillars by split and fendering. A two-dimensional elasto-plastic plane strain model with strain softening and interface mod- elling was used. The model was able to simulate the improvement of conditions experienced in lifting from the fender, when the roof adjacent to the fender has caved. In contrast, when a stook was left, providing support to a span of uncaved roof, the model pre- dicted the fender would suffer considerably more de- formation and strain.

Jan 1, 1992