Experience gained until now underground mining worldwide and in Poland indicates that remnants may have an impact on the occurrence of seismic phenomena. Remnants are stress concentration sites encompassing both the deposit and the layers of rock mass located above and below the undisturbed rock. In the case where stresses in the remnant exceed its strength, it may collapse, and under unfavourable geomechanical conditions, stress-induced rockburst may occur. Remnants may also cause breaking of strong roof layers above their edges, which results in the occurrence of high-energy shocks (Salustowicz [30], Adach [3], Adach and Butra [4]). This article presents the possibility of utilizing numerical modeling to evaluate the influence of remnant upon the occurrence of seismic phenomena. The results of numerical calculations performed for a model room-and-pillar mining system with roof deflection under the conditions of copper ore mines in the Legnica-Głogów Copper District (LGOM) are presented. Numerical calculations in a plane strain state were performed by means of Phase2 v. 8.0 software for the analyzed mining system in which remnant was left behind. The results of numerical modeling showed that sudden fracturing of roof layers above the mined out space may occur on the edge of the remnant. This may cause a shock with very high energy, and under the appropriate conditions, this may lead to the rockburst phenomenon.

The values of primary stresses are not allowed for as a criterion in the selection of roof bolting systems in mining excavations located at various depths in Polish copper ore mines. Therefore, in order to ensure enduring and safe operation of excavations, in particular, those driven in unfavourable geological and mining conditions, this problem has required solutions based on numerical methods. This article presents an example of applying numerical simulations to the evaluation of the stability of headings in Polish copper ore mines. The analyses included mining excavations located at various depths in the rock mass. This issue is of great importance, as safety regulations are prioritised in mining excavations which remain in operation even for several decades. The stability of the headings was evaluated with the use of the RS2 specialist numerical simulation software. This computer program uses the finite element method (FEM) for calculations. The rock parameters used in the numerical models have been determined on the basis of the Hoek–Brown classification. For that purpose, the RocLab 1.0 software was used. The parameters of the stress field were identified from the profile of the GG-1 shaft with the assumed hydrostatic state of stress. The numerical modelling was performed in a triaxial stress state and in a plane strain state. The numerical analyses were based on the Mohr–Coulomb failure criterion. The rock medium was described with the elastic-plastic model with softening (roof and walls) and with the elastic-plastic model (floor). The results of the numerical analyses served to provide an example of the application of a roof bolting system to protect headings located at the depths of 1000 m b.g.l. and 1300 m b.g.l.

This paper concerns the issue of selecting appropriate stress field parameters for predicting the stability of headings driven under the geological and mining conditions of Polish underground copper mines. The problem is of key importance due to strict safety requirements in mine workings that serve ventilation and transport functions. Numerical analyses were carried out for four stress field variants: the stress state determined based on Bulin’s formulas (variant 1), the hydrostatic stress state (variant 2), and stress states determined from in situ measurements conducted in the Rudna mine (variant 3 and variant 4). Numerical simulations were performed for a group of four headings, supported with fully grouted rock bolts, in the geological and mining conditions of the Rudna mine. Stability assessment was performed using the finite element method (FEM). Rock mass input parameters for the modeling were obtained with RocLab 1.0, applying the Hoek–Brown classification, while numerical analyses employed the Mohr–Coulomb failure criterion. The elastic–plastic model with softening was used to describe the rock mass behaviour. Numerical calculations were conducted in the RS2 computer program in a triaxial stress state and in a plane strain state. The range of the yielded rock mass zone in the roof of the headings was assumed as the optimal measure of the headings stability. The obtained simulation results provided a basis for recommending suitable rock bolting systems to protect the stability of headings developed under various initial stress field conditions.