Numerical modelling of the longwall mining and the stress state in Svea Nord Coal Mine

This thesis presents numerical and analytical investigation of the geomechanics underlying longwall mining. It was tried out to study the disturbances induced by longwall mining in nearby rocks and their influence on the stability of the gates, pillars and main tunnels of longwall mines. The thesis consists of two major parts: numerical and analytical investigations. The study site is the Svea Nord coalmine, Svalbard, Norway. A novel algorithm was proposed for numerical simulation of the longwall mining process. In the proposed algorithm progressive cave-in and fracturing of the roof strata, consolidation of the cave-in materials and stress changes are simulated in detail. In order to outline the caved-in roof rocks a criterion based on maximum principal strain (in tension) was used. The critical tensile strain of roof cave-in was determined through back-calculation of the surface subsidence above a longwall panel at the mine. The results of the simulations were then used to analyse stress changes induced by longwall mining and the stability of gates. The simulations revealed that the stability of the gates and the loading to the rock bolts are closely related to the width of the chain pillars. With slender pillars, shear displacements along weak interlayers and bedding planes result in heavy loading to the rock bolts. Therefore, the locations of weakness zones should be taken into account in rock bolt design. The developed algorithm was implemented to study the loading and stability of the barrier pillar of the mine. The barrier pillars protect the main tunnels and border area of the mine from disturbances induced by longwall mining in the panels. The simulations show that the stresses in the barrier pillars fluctuate up and down during mining because of periodic cave-in events behind the longwall face. A failure zone of about 12 m exists in the wall of the barrier pillars. A large portion of the barrier pillar is still intact and is, thus, capable of protecting the border area. The results of the detailed simulations of longwall mining via the developed algorithm were, also, implemented in a large-scale numerical model. The model consists of all of the longwall panels and the border area of the mine. It is intended that the coal in the border area on the other side of the longwall panels will be mined after completion of the longwall mining. There is concern about how the longwall mining affects the stress state in the border area and how stress changes would affect future mining in the border area. A failure zone of about 20 m developed in the wall of the main tunnels on the side of the border area after all the longwall panels were mined out. The stress state in the remaining portion of the border area remains unchanged. Therefore, it will be possible to mine the border area in the future. In order to investigate the roof strata cave-in mechanism in detail a discontinuous numerical simulation of roof cave-in process was conducted by UDEC code. The block size in the roof strata and the mechanical parameters of the discontinuities were obtained through back-calculations. The back-calculations were conducted with a statistical method, Design of Experiment (DOE). Numerical simulations revealed that jointed voussoir beams formed in the roof strata before the first cave-in. Beam bending results in stress fluctuations in the roof strata. The maximum deflection of a roof stratum at the study site before the first cave-in is about 70% of the stratum thickness. The simulations and field measurements show no periodic weighting on the longwall shields in this mine. Numerical sensitivity analyses show, however, that periodic weighting may occur in strong roof strata. Roof strata with a high Young’s modulus and large joint spacing are not suitable for longwall mining. The maximum sustainable deflection of the roof strata before cave-in depends upon the horizontal in-situ stress state. It slightly increases with the in-situ horizontal stress in the stratum beams, but the horizontal stress would increase the possibility of rock-crushing in deflected roof beams. The implemented numerical method would be useful in assessment of the cavability of the roof strata and in selection of longwall shields with adequate load capacity. As shown through discontinuous numerical simulations, the roof strata above the underground opening constructed in the stratified rocks form voussoir beams. The stability of those beams is the major concern in the study of the gate stability and roof cave-in assessment in the longwall panels. Two different analytical methods were developed for cases with and without the in-situ horizontal stress acting along the beams. In the analytical model for the beams without horizontal stress a bilinear shape was assumed for the compression arch generated within the voussoir beams. The stability of the compression arch is governed by the energy method. The model requires an iterative procedure for convergence, and an algorithm was proposed for it. The analytical method was verified with numerical simulations by means of a discrete element code, UDEC. For the beams subjected to in-situ horizontal stress, the classic beam theory was employed to drive the analytical solution for it. The superposition method was used to obtain bending/deflection equations of the beam. The validity of both the assumptions and the developed method were, also, investigated by numerical simulations. The developed analytical method revealed that high Young’s modulus of a beam rock increases the stability of the beams against buckling but it causes higher stress within the compression arch which increases the probability of crushing failures in the beam abutments and midspan. In-situ horizontal stress along beams increases their stability against buckling and abutment sliding failure, but it raises the possibility of crushing failure at the abutments and the midspan.PhD i geologi og bergteknikkPhD in Geology and Mineral Resources Engineerin

Coupled hydro-mechanical stability analysis of Åknes rock slope

Stability of Åknes rock slide was assessed by coupled hydro-mechanical numerical analysis. The analysis was carried out by considering explicitly the rock joints network and the sliding planes. Numerical analysis shows that with draining all the groundwater from the sliding surface the factor safety will improve ca. 0.5 and 4 % in the east and west flank of the slope, respectively. Fast infilling of the backscarp due to flooding can destabilise the west flank: decreasing the factor of safety to 0.89. It was shown that in hypothetically highest groundwater level inside the rock mass, the west flank has higher tendency of becoming unstable than the east flank.Norges vassdrags- og energidirektorat (NVE

A Simplified Approach to Estimate Anchoring Capacity of Blocky Rock Mass with Pressure Arch Theory

In this paper, a simplified method for predicting rock mass resistance against tensile forces from rock anchors (anchor) is presented. The interaction of anchors and rock mass was investigated using three-dimensional discontinuous numerical modelling. Several patterns of rock discontinuities were assumed in the numerical modelling while a single anchor is embedded in it. The numerical results show that the existence of a discontinuity set sub-parallel to the anchor significantly improves the rock mass resistance against the tensile force from the anchors. This phenomenon is due to the rock block interlocking at the sub-parallel discontinuity set. Rock block interlocking generates a zone of stress concentration inside the rock mass which has an arch shape (i.e., a pressure arch), resisting against the anchor's tensile force. The load-bearing capacity of the pressure arch plays a significant role in the resistance of the rock mass against the forces from the rock anchor. A voussoir beam analogy was utilised to study the load-bearing capacity of the pressure arch. A simplified analytical approach was developed to assess the load-bearing capacity of the voussoir beams. Then, it was used in combination with the weight of the mobilised rock mass by the anchor to assess the maximum anchoring resistance of the rock mass (anchoring capacity). The suggested method was calibrated by numerical modelling and relevant published pull-out test results. The technique developed in this paper shows the significance of rock block size, shear behaviour of rock discontinuities, Young's modulus of the rock mass, and uniaxial compressive and tensile strength of the intact rock in anchoring capacity of rock masses.publishedVersio

The Mechanism, Functioning and installation of Reinforced Ribs of Shotcrete in Norwegian Method of Tunnelling (NMT)

publishedVersio

Hydrogeological report

NGI is engaged by Norwegian water resources and energy directorate (NVE) to carry out stability analysis of the Åknes rock slide. This report is presenting and analysing the most important investigations and monitoring, such as packer tests and piezometer measurements in both open boreholes and boreholes with packers dividing the borehole in sections. A hydrogeological model is presented. The main conclusions are: - Water head on the sliding plane are small, only small peaks up to 3-4 meter over short periods of time are monitored. For most of the boreholes the water table is several meters below the main sliding plane. - The hydraulic communication in the slope is high. It seems like the hydraulic system is fed from a bigger area, not only the backscarp.Norges vassdrags- og energidirektorat (NVE

3D Numerical Modeling of Rock Mass Failure in an Uplift Test of a Rock Anchor with Focus on the Role of Rock Joints

This paper aims to numerically study the failure mechanisms of the rock mass during the uplift test of a rock anchor installed in limestone. The rock anchor consisted of a smooth tendon, with short threaded sections at both ends. At the distal end of this tendon, a washer and a nut were screwed together to form a plate-anchor. The annulus between the rock anchor and the borehole wall was filled with cementitious grout. The rock mass had three distinctive joint sets. One of them was vertical, i.e., parallel to the anchor axis, while the other two were inclined with a dip angle ≤ 40°. The uplift test was designed so that failure would occur in the rock mass. Geotechnical monitoring data collected from the field was used to calibrate a numerical model with the 3DEC (Three-Dimensional Distinct Element Code) software. This software was chosen because it can simulate the behavior of rock joints. The calibrated numerical model was able to replicate quite accurately the pull force over anchor displacement during testing until the ultimate load and the shallow rock cone failure observed in the field. The model indicates that rock joints start to open at 60% of the ultimate pull load and the key input parameters controlling the ultimate load capacity of the anchor are the tensile strengths of the inclined joint sets.publishedVersio

Numerical modelling of the longwall mining and the stress state in Svea Nord Coal Mine

This thesis presents numerical and analytical investigation of the geomechanics underlying longwall mining. It was tried out to study the disturbances induced by longwall mining in nearby rocks and their influence on the stability of the gates, pillars and main tunnels of longwall mines. The thesis consists of two major parts: numerical and analytical investigations. The study site is the Svea Nord coalmine, Svalbard, Norway. A novel algorithm was proposed for numerical simulation of the longwall mining process. In the proposed algorithm progressive cave-in and fracturing of the roof strata, consolidation of the cave-in materials and stress changes are simulated in detail. In order to outline the caved-in roof rocks a criterion based on maximum principal strain (in tension) was used. The critical tensile strain of roof cave-in was determined through back-calculation of the surface subsidence above a longwall panel at the mine. The results of the simulations were then used to analyse stress changes induced by longwall mining and the stability of gates. The simulations revealed that the stability of the gates and the loading to the rock bolts are closely related to the width of the chain pillars. With slender pillars, shear displacements along weak interlayers and bedding planes result in heavy loading to the rock bolts. Therefore, the locations of weakness zones should be taken into account in rock bolt design. The developed algorithm was implemented to study the loading and stability of the barrier pillar of the mine. The barrier pillars protect the main tunnels and border area of the mine from disturbances induced by longwall mining in the panels. The simulations show that the stresses in the barrier pillars fluctuate up and down during mining because of periodic cave-in events behind the longwall face. A failure zone of about 12 m exists in the wall of the barrier pillars. A large portion of the barrier pillar is still intact and is, thus, capable of protecting the border area. The results of the detailed simulations of longwall mining via the developed algorithm were, also, implemented in a large-scale numerical model. The model consists of all of the longwall panels and the border area of the mine. It is intended that the coal in the border area on the other side of the longwall panels will be mined after completion of the longwall mining. There is concern about how the longwall mining affects the stress state in the border area and how stress changes would affect future mining in the border area. A failure zone of about 20 m developed in the wall of the main tunnels on the side of the border area after all the longwall panels were mined out. The stress state in the remaining portion of the border area remains unchanged. Therefore, it will be possible to mine the border area in the future. In order to investigate the roof strata cave-in mechanism in detail a discontinuous numerical simulation of roof cave-in process was conducted by UDEC code. The block size in the roof strata and the mechanical parameters of the discontinuities were obtained through back-calculations. The back-calculations were conducted with a statistical method, Design of Experiment (DOE). Numerical simulations revealed that jointed voussoir beams formed in the roof strata before the first cave-in. Beam bending results in stress fluctuations in the roof strata. The maximum deflection of a roof stratum at the study site before the first cave-in is about 70% of the stratum thickness. The simulations and field measurements show no periodic weighting on the longwall shields in this mine. Numerical sensitivity analyses show, however, that periodic weighting may occur in strong roof strata. Roof strata with a high Young’s modulus and large joint spacing are not suitable for longwall mining. The maximum sustainable deflection of the roof strata before cave-in depends upon the horizontal in-situ stress state. It slightly increases with the in-situ horizontal stress in the stratum beams, but the horizontal stress would increase the possibility of rock-crushing in deflected roof beams. The implemented numerical method would be useful in assessment of the cavability of the roof strata and in selection of longwall shields with adequate load capacity. As shown through discontinuous numerical simulations, the roof strata above the underground opening constructed in the stratified rocks form voussoir beams. The stability of those beams is the major concern in the study of the gate stability and roof cave-in assessment in the longwall panels. Two different analytical methods were developed for cases with and without the in-situ horizontal stress acting along the beams. In the analytical model for the beams without horizontal stress a bilinear shape was assumed for the compression arch generated within the voussoir beams. The stability of the compression arch is governed by the energy method. The model requires an iterative procedure for convergence, and an algorithm was proposed for it. The analytical method was verified with numerical simulations by means of a discrete element code, UDEC. For the beams subjected to in-situ horizontal stress, the classic beam theory was employed to drive the analytical solution for it. The superposition method was used to obtain bending/deflection equations of the beam. The validity of both the assumptions and the developed method were, also, investigated by numerical simulations. The developed analytical method revealed that high Young’s modulus of a beam rock increases the stability of the beams against buckling but it causes higher stress within the compression arch which increases the probability of crushing failures in the beam abutments and midspan. In-situ horizontal stress along beams increases their stability against buckling and abutment sliding failure, but it raises the possibility of crushing failure at the abutments and the midspan