Next Article in Journal
Overview of the Role of Energy Resources in Algeria’s Energy Transition
Next Article in Special Issue
Good Practices to Counteract Epidemic Emergency in Mining Companies in Poland
Previous Article in Journal
A Multi-Agent Adaptive Architecture for Smart-Grid-Intrusion Detection and Prevention
Previous Article in Special Issue
Research on Calculation Method for Discharge Capacity of Draining Well in Tailing Ponds Based on “Simplification-Fitting” Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods

by
Karolina Adach-Pawelus
Faculty of Geoengineering, Mining and Geology, Wroclaw University of Science and Technology, 50-370 Wroclaw, Poland
Energies 2022, 15(13), 4729; https://doi.org/10.3390/en15134729
Submission received: 26 April 2022 / Revised: 7 June 2022 / Accepted: 9 June 2022 / Published: 28 June 2022
(This article belongs to the Special Issue Volume II: Mining Innovation)

Abstract

:
An important element of numerical modeling for specific mining issues is the selection of model parameters. The incorrect determination of geomechanical parameters can result in significant calculation errors carried throughout the entire problem. This paper presents a method for determining effective geomechanical parameters for technological and residual pillars through the use of numerical modeling, specifically, back-calculation. This is based on the results of numerical simulations, measurement data (e.g., excavation convergence measurements), and statistical methods (a non-linear regression model with “dummy” variables). The result is that appropriate parameters of pillars are set out iteratively so that the displacements of selected points in the numerical model correspond (with some approximation) to the results of mine measurements. The procedure of determining pillar parameters is presented using a case study of one mining field in an underground copper mine, where the deposit is mined using the room and pillar system. Numerical calculations were performed using a Phase2 v. 8.0 program (Rocscience, Toronto, Canada), while statistical calculations used a Statistica computer program. The results of excavation convergence measurements performed in the analyzed mine have been applied. This paper shows that for the presented method, the resulting matching of theoretical values of convergence determined numerically for specified pillar parameters to in-situ results of convergence measurements, is very good (R2 = 0.9896). This work exemplifies a set of the parameters of pillars for an elastic model of rock mass, but this method can also be applied to other models.

1. Introduction

In the past, rock mechanics was primarily described using empirical methods, based on experience and analytical methods giving a closed form of solution. As a result of the development of computer technology and computational methods, the most common methods used to solve geomechanics problems are numerical methods, which allow for the analysis of complex geometry (in close to real conditions) and different material behaviors. Both two- and three-dimensional tasks are solved. The most commonly used numerical calculation methods in geomechanics are the finite difference method, finite element method, and the boundary element method; though the discrete element method or hybrid methods (hybrid continuum/discrete methods) [1,2] are quite often used. Numerical modeling plays a very important role in the design of underground mines. It allows assessment of the current mining situation and predicts the behavior of rock mass along with the progress of work. Numerical methods in underground mining are used for, among others, stability analysis of mining excavations and the design of their support [3,4,5,6,7,8,9,10], to compare different types of support [11], for simulation of deposit excavation using various mining systems [12,13,14,15,16,17,18,19,20], for studies of the behavior and stability of pillars and remnants [21,22,23,24,25,26,27,28,29], for analysis of mining salt deposits with backfill and the determination of the backfilling influence on the dynamics of deformation of the undermined rock mass [30], as well as for the assessment of seismic and rockburst hazard [20,22,27,28,31,32,33,34,35], etc.
A key element of any rock mass behavior modeling is the correct determination of its parameters and verification of obtained results. The introduction of incorrect parameters into the model will result in incorrect results. These parameters should be chosen so that the values set out as a result of numerical calculations correspond to actual values obtained on the basis of observations and mine measurements. The main source of input data for a numerical model are laboratory tests of rock samples and field tests. The impact of the sample size on the mechanical characteristics of rock samples is referred to as the “size effect” and is a significant issue [36,37]. The values of strength and deformation parameters set out by laboratory tests are often reduced to take into account the “size effect” between the sample (micro scale) and the rock mass (macro scale). The parameters of rock mass for numerical calculations can be set out using rock mass weakening factors chosen depending on the structure of the rock mass and the compressive strength of the rock samples or by using Hoek–Brown classification [38,39,40].
In practice, however, it happens that the obtained calculated values (e.g., of rock mass displacements) differ significantly from the values measured in mining conditions. The choice of suitable geomechanical parameters for numerical calculations is a particularly important issue in the case of technological and residual pillars in room and pillar mining systems and when considering their progressive destruction, when the mining front is moved. The yielding of pillars in numerical calculations is introduced by reducing their strength and deformation parameters. In order to make the best projection of the actual working conditions of pillars, the strength and deformation parameters are chosen most often using a method called “back-calculation”. The required parameters are set out based on in-situ measurements in such a way, that a numerical model should fit best to the actual behavior of the object. Back-calculation based on measurements of displacements was initiated by Sakurai in 1981 [41] and is widely applied in geomechanics.
This paper presents a method for determining effective geomechanical parameters for technological and residual pillars used in numerical calculations for room and pillar mining systems using back-calculation, which is based on the results of numerical simulations, measurement data (e.g., excavation convergence measurements), and statistical methods. Numerical calculations have been performed as a plane strain problem using Phase2 v. 8.0 program (Rocscience, Toronto, Canada), which uses the finite element method, while the statistical analysis was done using Statistica. The above method is applied to one of the mining fields of the underground copper mines belonging to KGHM Polska Miedz SA. The analysis covered the D-IE mining field located in the Polkowice-Sieroszowice mine in south-western Poland. Copper ore deposit in the analyzed field was mined using the room and pillar system with a roof deflection. This work exemplifies how to set out the parameters of pillars for an elastic model of rock mass, but this method can also be used for other models.

2. The Method of Determination of Effective Geomechanical Parameters for Technological and Residual Pillars for Numerical Modeling Using Back-Calculation Based on the Results of Numerical Simulations, Measurement Data, and Statistical Methods

For the best projection of reality to determine the effective geomechanical parameters of technological and residual pillars applied to numerical calculations for room and pillar mining systems, back-calculation has been suggested based on the results of numerical simulations and measurement data, in which a non-linear regression model was applied with so-called dummy variables. This allows determining the values of the pillars parameters, for which the values of excavation convergence calculated in the numerical model will be most closely related to the results of convergence measured is-situ in this excavation. The model of non-linear regression with “dummy” variables is described by the formula:
Y = P1Z1 + P2Z2 + P3Z3 + … + PkZk,
where:
  • Y—explained variable (results of in-situ measurements of convergence for a chosen excavation),
  • Z1, …, Zk—so-called “dummy” variables (with a value of 1 or 0),
  • P1, …, Pk—quadratic functions acting as parameters of a regression model with “dummy” variables.
The established effective geomechanical parameters of technological and residual pillars are independent variables of quadratic functions P1, …, Pk. The number of variables in a quadratic function depends on the number of determined pillar parameters. Therefore, the actual determined parameters of a regression model with dummy variables are independent variables of the quadratic function being the parameters of pillars. In a case when all technological and residual pillars are characterized by one parameter such as E then functions P1, …, Pk take the form of a quadratic function of one variable:
P = β0 + β1E + β11E2
If we use more parameters to characterize the work of technological and residual pillars, the Formula (1) takes the form:
for two parameters (e.g., E, c):
P = β0 + β1E + β2c+ β11E2 + β12Ec + β22c2,
for three parameters (e.g., E1, E2, c):
P = β0 + β1 E1+ β2 E2 + β3 c + β11 E12 + β12 E1 E2 + β22 E22 + β13 E1 c + β33 c2 + β23 E2 c
The parameters of quadratic functions: β0, β1, β2, β3, β11, β12, β13, β11, β22, and β33 are set out by the least square method on the basis of convergence values defined in numerical simulations for different combinations of established values of the parameters characterizing the pillars (the matrix of the numerical experiment). The matrix of the numerical experiment is selected iteratively in such a way that the calculations are based on interpolation, i.e., the calculated parameters of pillars should be within a space defined by the matrix of the numerical experiment.
In order to match the results of the numerical experiment to measurement data, the regression model parameters are set out using backward stepwise regression. This involves eliminating from the model the effective geomechanical parameters of pillars, which had no significant effect on the optimal adjustment of convergence values (identified on the basis of numerical modeling) to the measurement data. The statistical significance of respective parameters is evaluated on the basis of confidence intervals or p-values for the assumed level of significance α. The regression model with dummy variables, which incorporates all the selected parameters of pillars, is individually deprived of those where p > α (and therefore those for which there was no reason to reject the hypothesis, that the established parameter is equal to 0). In the next step, a new regression model is set out with a number of parameters reduced by 1. Finally, when all parameters of the pillars are considered statistically significant (significantly different from 0), the process of elimination is interrupted and matching of the results of the predicted excavation convergence to measurement data is estimated based on the coefficient of determination R2. On this basis, a decision is made about the values of the geomechanical parameters for technological and residual pillars adopted for further numerical calculations.
The procedure of determining effective geomechanical parameters for technological and residual pillars based on the results of numerical simulations and data from mine measurements using statistical methods is presented in Figure 1. According to Figure 1, the first step (1) includes the choice of strength and deformation parameters of pillars, which have a decisive impact on the value of calculated vertical displacements for a given model of rock mass. Then, after construction of a numerical model of the selected mining field, (2) one should specify the number of quested effective geomechanical parameters characterizing the work of technological and residual pillars (3) and select the values of these parameters, for which numerical simulations will be carried out (construction of the numerical experiment matrix) (4). In the case of a numerical experiment matrix consisting of three parameters describing the work of the pillars (e.g., E1, E2, and c) and three different values adopted for each parameter, the number of numerical simulations will be 27. Simulations performed for a numerical experiment matrix are a basis for determining a convergence graphs model in a chosen excavation in the following steps of mining excavation (5). On the established convergence curves, one should choose representative points for the regression analysis (6) and determine the quadratic functions of an appropriate number of variables depending on the number of sought geomechanical parameters of the pillars (7). Knowing the quadratic functions serving as model parameters, one can specify the non-linear regression model with “dummy” variables (8) and determine the values of effective geomechanical parameters for the technological and residual pillars (9). The next step is to examine whether the received parameters have physical meaning (e.g., c ≥ 0, E ≥ 0, φ ≥ 0) (10). If not, calculations should be terminated (11), as this may mean that the numerical model has not been properly constructed: too small density of the grid, insufficient accuracy of calculations, etc. When the calculated parameters have physical meaning, it should be checked whether they belong to one type of parameter (12). If not, then one should check whether the set parameters are statistically significant (14) and if so, it should be checked whether the parameters belonging to one type are significantly different from each other (13). If among the established pillar parameters are those for which there are no grounds to reject the statistical hypothesis that they are equal, one should eliminate one of them, reduce the number of sought pillar parameters, (15) and repeat the procedure. If all set-out parameters are significantly different from each other and in the case of parameters not belonging to one type, one should check their statistical significance (14). In a case when among the established parameters are those not statistically significant, one of them must be eliminated (15), and following the steps of the procedure should be repeated. However, when all pillar parameters have been considered statistically significant, one should check whether they are in a space defined by the experiment matrix (16). If the set parameters do not belong to this space, one should return to the construction of the experiment matrix (4) and change the range of selected values. Otherwise, the set of effective geomechanical parameters for technological and residual pillars can be applied for further numerical modeling, evaluating the matching of numerical experiment data to measurement data on the basis of the coefficient of determination R2 (17). The procedure of selecting pillar parameters based on numerical simulations and measurement data can be performed repeatedly for different models of rock mass, searching for optimal parameters that will produce the best match of numerical experiment data to measurement data.

3. Case Study in Polkowice-Sieroszowice Mine

3.1. Characteristics of Research Area

The method of determining effective geomechanical method parameters for technological and residual pillars for numerical modeling using back-calculation based on the results of numerical simulations, measurement data, and statistical methods are presented in the example of one mining field (D-IE) in Polkowice-Sieroszowice underground copper mine. The Polkowice-Sieroszowice mine belongs to KGHM Polska Miedz SA. It is located in the south-western part of Poland and mines the copper ore deposit which covers the central part of a geological unit known as Sudetic Monocline. Sudetic Monocline falls gently towards the northeast. It is constructed of Permian and Triassic sediments, which have a base made from Proterozoic crystalline rocks and Carboniferous sedimentary rocks. The deposit occurs in Permian formations, contacted by a dolomite limestone series, red sandstone, and Lower Permian limestone. Shaped in the form of a pseudobed of variable thickness (from 0.4 to approx. 20 m) and low gradient (approx. 4°), it lies at great depth (from 600 to 1400 m). The copper ore deposit is formed by an accumulation of sulfides, mainly chalcocite, bornite, and chalcopyrite. Sulfide mineralization occurs at the contact of red sandstone and Permian limestone layers. It includes carbonate rocks (dolomites and limestones), copper-bearing shales in the bottom part of the Permian limestone, and white sandstones. The deposits of the Polkowice-Sieroszowice mine subject to ore mineralization are mainly carbonate rocks and shales. Mining of the deposit is performed using a variety of room and pillar systems, depending on geological and mining conditions in a given mining field [29].
D-IE mining field was located in the mining area Sieroszowice I, in the Polkowice-Sieroszowice mine. It was a closing field explored since March 2005, where mining works were carried out in the vicinity of gobs. In February 2008, due to problems with maintaining roof stability, the remnant was left behind on the right side of the D-IE mining field (Figure 2).
In the D-IE mining field, the deposit balance occurs in the lower part of a carbonate series of Permian limestones and the roof part of new red sandstone; it is comprised of grey quartz, fine-grained sandstone, loamy copper-bearing shale, and dolomite loamy shale, as well as streaked, dark grey, crypto-crystalline dolomite. The roof is made of rock layers, being part of a Permian limestone carbonate series, namely of calcareous dolomites with clear divisibility of bed (occurring in intervals of 0–2 m above the excavation roof), of concise calcareous dolomites with quite clear divisibility (occurring in intervals of 2–5 m above the excavation roof), and calcareous dolomites and dolomitic limestones with a bed structure (that occur in intervals of more than 5 m above the excavation roof). The carbonate series is directly covered by anhydrites. The direct floor is built of grey sandstones of red Permian sediment rock. These are fine-grained quartz sandstones with a loamy bond, carbonate-loam bond, and locally anhydrite bond (in the eastern part of the area). The roof part of the sandstones, due to the larger amount of carbonate bond, is harder and more concise. The deposit is oriented towards NW-SE and its decline (2–3°) towards NE. The rock formation has marginal tectonic sensitivity. The height of the mined deposit is 2.0–2.8 m.
Until 2008 mining of the deposit in the discussed area was conducted using a room and pillar system with roof deflection and closing pillar (J-UGR-PS), while in 2008 the closing pillar was liquidated, and further works were performed using a room and pillar system with roof deflection (J-UG-PS). Exploration using the room and pillar system consists in cutting the deposit with rooms and strips with separation of technological pillars of a certain geometry, which protect the roof over the working area. The size of the pillars is chosen to provide its work in the post-critical state. In the D-IE mining field, the cutting work was carried out using technological pillars situated perpendicular to the mining front, with basic dimensions of 6 × 8 m (J-UGR-PS and J-UG-PS). In the discussed mining systems, the height of excavation in the cutting phase depends on the thickness of the deposit and the requirements of working machines and is not more than 4.5 m. The width of excavations does not exceed 7 m. The minimum size of the opening face of the mining front is equal to the sum of two strips and the length of two rows of pillars into undisturbed rock. In the D-IE mining field, the width of the opening was 4–5 strips. Along with the progress of the mining front, the technological pillars from the last row before gobs, depending on the degree of their disintegration, are adjusted or cleft into smaller ones. The resultant support pillars are adjusted to residual dimensions in elementary plots and then left in gobs. They work as supports to mitigate the deflection of roof layers. For D-IE field size of residual pillars left in the gobs amounted to approx. 5.2 m2. In a room and pillar system with roof deflection and a closing pillar, the technological pillars are left in a separated part of the field, creating a gradually lengthening closing pillar. In one of the excavations of this pillar, a conveyor belt is assembled, which successively extends along with the progress of cutting. The width of a closing pillar depends on local geological and mining conditions and is generally 40–120 m [39]. The work of mining systems J-UGR-PS and J-UG-PS are shown in Figure 3 and Figure 4.

3.2. Characteristics of Numerical Modeling

Numerical calculations were performed in a plane strain state by a computer program Phase2 v. 8.0 (Rocscience, Toronto, ON, Canada) [43]. The computational model was a plate, which comprises the rock layers creating the rock mass (Figure 5). Construction of the mass resulted from geological recognition conducted in the analyzed field. The upper edge of the model was loaded with vertical pressure, replacing the influence of overlaying rocks. It was assumed that at the upper edge of the plate the stress should be equal to 17.657 MPa, corresponding to the value of the vertical stress set out for the D-IE mining field on the basis of data from the borehole S-294. The calculations considered the deadweight of rock layers. Horizontal stress values were determined on the basis of Poisson’s ratio υ of a given rock layer. For the edges of the plate, the displacement edge conditions were assumed. At the lower edge of the model—lack of vertical displacements while at the side edges—no displacements in directions perpendicular to the surface of the edge. An applied grid of finite elements was composed of three nodal elements of triangular shape. In the central part of the plate, adjacent to the excavations, the grid was densified to improve the accuracy of numerical calculations.
The calculations were performed stepwise, simulating the mining carried out using a room and pillar system with parameters characteristic of an analyzed mining field (64 calculation steps). The first step covered the situation in the rock mass before the creation of mine excavations (Figure 6a). The second step consisted of the cutting of undisturbed rock into technological pillars (8 m in width) (Figure 6b). In the next steps, the size of the technological pillars was reduced to residual size (3 m width) and further technological pillars were cut out (Figure 6c). Cutting of the deposit was carried out with strips having dimensions 6 m width under the roof. Numerical simulations considered the width of the working area opening consisting of 5 strips.

3.3. Determination of Parameters for Rocks and Rock Mass

The parameters of rock mass, which were assumed for numerical modeling, were determined based on Hoek–Brown classification. The results of laboratory tests of rock samples taken from geotechnical boreholes located in the analyzed mining field were applied. The averaged parameters of the rocks, which were designated in the laboratory for the D-IE field, are shown in Table 1, while the parameters of the rock mass are shown in Table 2.

3.4. Determination of Effective Geomechanical Parameters for Technological and Residual Pillars by Numerical Modeling Using Back-Calculation Based on the Results of Numerical Simulations, Measurement Data and Statistical Methods

The parameters of the technological and residual pillars described by an elastic model were determined using the procedure presented in the article. On the basis of in-situ tests carried out in the ZG Polkowice-Sieroszowice mine in 2002–2007, the course of the excavation convergence in time was determined for room and pillar mining systems with roof deflection and closing pillar. The pattern convergence curve, which was adopted as a reference for numerical calculations performed for the D-IE field, is shown in Figure 7. Each measurement of convergence made on a chosen test post was referred to as the mining step and the position of the working front line on a given day of measurement. This enabled the construction of a model that would fit the actual situation in the field and allow the determination of the effective geomechanical parameters of the technological and residual pillars.
The values of longitudinal elasticity modulus E were determined for an elastic model, (the parameter having a decisive impact on the numerically calculated values of displacements). Depending on the distance from the face of the mining front, the degree of pillar disintegration is varied. In the first step, the pillars were divided into three groups. Three technological pillars (located close to undisturbed rock) were characterized by longitudinal elasticity modulus E1, the other two technological pillars with modulus E2, and the residual pillars with modulus E3 (Figure 8). It was assumed that E1 = E2 + a. The following values of E1, E2, and E3 were adopted for the matrix of the numerical experiment:
E1E2E3
2000 MPa500 MPa100 MPa
6000 MPa1500 MPa150 MPa
10,000 MPa2500 MPa200 MPa
27 numerical simulations were performed for all combinations of assumed values of Young’s modulus E. The convergence of the selected excavation in subsequent steps of the executed mining for a few chosen cases is presented in Figure 9.
Using the method of surface regression with Statistica v. 10 program, the quadratic function parameters of three variables were determined for selected points that are included in a non-linear regression model with dummy variables:
Y = Z1(−145.75967 + 0.01(E2 + a) − 0.00000046(E2 + a)2 + 0.05718E2 − 0.00001043E22 − 0.1659E3 + 0.00054E32 0.00000100(E2 + a)E2 + 0.00000087(E2 + a)E3 + 0.00000142E2E3) + … + Z9(−4620.2027 + 0.0405(E2 + a) − 0.00000188(E2 + a)2 + 0.36882E2 − 0.00006627E22 + 22.9864E3 − 0.04617E32 − 0.00000146(E2 + a)E2 + 0.00002395(E2 + a)E3 − 0.00019174E2E3)
The values of a, E2, and E3 calculated using Statistica v. 10 are shown in Table 3.
E1 = E2 + a = 725.066 + 5541.145 = 6266.211 MPa
Based on confidence intervals, statistical inference was performed to check whether E1 and E2 are significantly different from each other, namely if the added value of a is different from zero. The assumed level of significance α = 5%.
H0: a = 0
H1: a 0
The a = 0 parameter is within 95% of the confidence interval, which indicates that there is no basis to reject the zero hypothesis H0 (Figure 10). Therefore, it cannot be stated that at the level of significance α = 5% E1E2. In addition, statistical inference based on test probability of p-value (p = 0.893210 > α = 0.05) also indicates that there is no basis to reject H0: a = 0 hypothesis. In such a situation, it is assumed that E1 = E2 and the number of sought parameters of pillars was reduced to two. Technological pillars were characterized by longitudinal elasticity modulus E1, while residual pillars with longitudinal elasticity modulus E2 (Figure 11). It was assumed, that E1 = E2 + a and the following values of E1 and E2 were applied to the matrix of the numerical experiment:
E 1 E 2 200   MPa 80   MPa 400   MPa 120   MPa 600   MPa 160   MPa
There were nine numerical simulations. Convergence in the selected excavation is shown in Figure 12 in subsequent steps.
Using the method of surface regression with Statistica v. 10 program, the quadratic function parameters of two variables were determined for selected points that are included in the non-linear regression model with dummy variables:
Y = Z1(−220.077 + 0.424(E2 + a) − 0.000322(E2 + a)2 − 0.0178E2 + 0.00005993E22 + 0.00000610(E2 + a)E2) + … + Z16(−4746.12 + 1.84(E2 + a) − 0.000983(E2 + a)2 + 28.539E2 − 0.06987090E22 − 0.00252617(E2 + a)E2)
The values of a and E2 were determined using Statistica v.10 presented in Table 4.
E1 = E2 + a = 136.931 + 174.424 = 311.355 MPa
Statistical inference carried out based on confidence intervals was performed to check whether E1 and E2 are significantly different from each other, namely if the added value of a is different from zero. The assumed level of confidence α = 5%.
H0: a = 0
H1: a 0
The parameter a = 0 is located outside the 95% of confidence interval, which leads to rejection of zero hypothesis H0 (Figure 13). It can therefore be stated that at significance level α = 5% E1E2. Statistical inference based on test probability of p-value (p = 0.010470 < α = 0.05) also suggests rejection of hypothesis H0: a = 0.
The final values of E1 and E2 parameters determined using Statistica v.10 are shown in Table 5.
On the basis of p-values it can be concluded that E1 and E2 parameters are statistically significant (E1: p = 0.0000364562557 < α = 0.05; E2: p = 0.0000000000119 < α = 0.05). The effective parameters of pillars E1 = 311.355 MPa and E2 = 136.931 MPa are also located in the space defined by the matrix of the numerical experiment. The resulting matching of theoretical values of convergence determined for specified pillar parameters, to in-situ results of convergence measurements is shown in Figure 14.
Reduced values of elasticity modulus for technological and residual pillars were adopted for further numerical modeling of the geomechanical situation in the D-IE mining field. The final validation of the numerical models was based on the results of convergence measurements of excavations carried out in this field.

4. Conclusions

In this paper, the innovative method for determining the geomechanical parameters of technological and residual pillars by numerical calculations using back-calculation based on the results of numerical simulations, in-situ measurement data, and a non-linear regression model with dummy variables, has been presented. The method has the potential to be useful in numerical modeling for specific mining issues due to the importance of model parameter selection. Incorrect determination of the geomechanical parameters results in significant calculation errors throughout the entire analyzed problem. These parameters should be chosen so that the values set out as a result of numerical calculations correspond to actual values obtained on the basis of observations and mine measurements.
In the case of room and pillar mining system modeling, a particularly difficult issue is the choice of parameters for technological and residual pillars (considering the progress of their destruction when a mining front is moved). Incorrect determination of the geomechanical parameters of pillars results in significant calculation errors throughout the whole analyzed mining field. The presented method enables the determination of pillar parameter values, for which the calculated numerical values of convergence of a chosen excavation in the model would be most closely related to the results of excavation convergence measurements carried out in mining conditions. This method facilitates estimation of the matching of predicted excavation convergence results to measurement data (based on the coefficient of determination, R2. The case study presented in this paper proved that by using the described method, the matching between results of in-situ convergence measurements in mines to convergence values obtained for effective parameters of pillars determined using statistical methods is very good, with R2 is equal to 0.9896.
The method has a lot of advantages and can be used to calibrate the numerical models of other engineering problems, e.g., for the determination of gob parameters for different mining systems. The method allows the determination of pillar parameters for different models of rock mass. In the case of an elastic-plastic model, it enables the description of the post-destructive part of stress-strain characteristics using the appropriate values of strength parameters such as the Coulomb–Mohr hypothesis. This is of vital importance for the accuracy of obtained results of numerical modeling and thus for the degree of real situation representation (e.g., a real situation in a modeled mining area). An appropriately verified numerical model allows for accurate analysis of the current situation in an analyzed region, for forecasting rock mass behavior during the progress of work, and for identification of dangerous phenomena which might create a threat to people working in a particular area. Future research will include the application of the method described in this manuscript to the numerical analysis of different mining situations.

Funding

The research work was co-founded with the research subsidy of the Polish Ministry of Science and Higher Education granted for 2022.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Jing, L.; Hudson, J.A. Numerical methods in rock mechanics. Int. J. Rock Mech. Min. Sci. 2002, 39, 409–427. [Google Scholar] [CrossRef]
  2. Jing, L. A review of technics, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering. Int. J. Rock Mech. Min. Sci. 2003, 40, 283–353. [Google Scholar] [CrossRef]
  3. Beer, G.; Pousen, A. Efficient numerical modelling of faulted rock using the boundary element method. Int. J. Rock Mech. Min. Sci. 1995, 32, 117A. [Google Scholar] [CrossRef]
  4. Kayupov, M.A.; Kuriyagawa, M. DDM modelling of narrow excavations and/or cracks in anisotropic rock mass. In Proceedings of the ISRM International Symposium-EUROCK 96, Torino, Italy, 2–5 September 1996; Balkema: Rotterdam, The Netherlands, 1996; pp. 351–358. [Google Scholar]
  5. Cerrolaza, M.; Garcia, R. Boundary elements and damage mechanics to analyze excavations in rock mass. Eng. Anal. Boundary Elements 1997, 20, 1–16. [Google Scholar] [CrossRef]
  6. Pan, E.; Amadei, B.; Kim, Y.I. 2D BEM analysis of anisotropic half- plane problems—Application to rock mechanics. Int. J. Rock Mech. Min. Sci. 1998, 35, 69–74. [Google Scholar] [CrossRef]
  7. Shou, K.J. A three-dimensional hybrid boundary element method for non-linear analysis of a weak plane near an underground excavation. Tunn. Undergr. Sp. Technol. 2000, 16, 215–226. [Google Scholar] [CrossRef]
  8. Mark, C.; Molinda, G.; Dolinar, D. Analysis of roof bolt stability. In Proceedings of the 20th International Conference on Ground Control in Mining, Morgantown, WV, USA, 7–9 August 2001; Peng, S.S., Mark, C., Khair, A.W., Eds.; West Virginia University: Morgantown, WV, USA, 2001; pp. 218–225. [Google Scholar]
  9. Pawelus, D. Evaluation of corridor excavation stability in the area of R-XI shaft using elastic-plastic model of rock formation and Coulomb- Mohr criterion. “Cuprum” Sci. Tech. Pap. Ore Min. 2013, 4, 21–40. [Google Scholar]
  10. Adach-Pawelus, K.; Pawelus, D. Influence of driving direction on the stability of a group of headings located in a field of high horizontal stresses in the Polish underground copper mines. Energies 2021, 14, 5955. [Google Scholar] [CrossRef]
  11. Yassien, A.M.; Zhang, Y.; Han, J.; Peng, S.S. Comparison of Some Aspects of Bolting Mechanisms between Fully-Grouted Resin and Tensioned Bolts in Underground Mine Entries. In Proceedings of the 21st International Conference on Ground Control in Mining, Morgantown, WV, USA, 6–8 August 2002; Peng, S.S., Mark, C., Khair, A.W., Heasley, K.A., Eds.; West Virginia University: Morgantown, WV, USA, 2002; pp. 114–125. [Google Scholar]
  12. Ghanmi, A.; Planeta, S.; Djelloud, J. Numerical modelling of mining sequences for steep orebodies. In Mine Planning and Equipment Selection; Balkema: Rotterdam, The Netherlands, 1998; pp. 53–57. [Google Scholar]
  13. Board, M.; Brummer, R.; Seldon, S. Use of Numerical Modelling for Mine Design and Evaluation. In Underground Mining Methods: Engineering Fundamentals and International Case Studies; Society for Mining, Metallurgy and Exploration: Englewood, CO, USA, 2001; pp. 483–492. [Google Scholar]
  14. Yasitli, N.E.; Unver, B. 3D numerical modelling of longwall mining with top-coal caving. Int. J. Rock. Mech. Min. Sci. 2005, 42, 219–235. [Google Scholar] [CrossRef]
  15. Peng, S. Coal Mine Ground Control, 3rd ed.; Society for Mining Metallurgy: Englewood, CO, USA, 2008. [Google Scholar]
  16. Villegas, T.; Nordlund, E. Numerical Analysis of the Hangingwall Failure at the Kiirunavaara Mine. In Proceedings of the InInternational Conference & Exhibition on Mass Mining, Luleå, Sweden, 9–11 June 2008; Schunnesson, H., Nordlund, E., Eds.; Luleå Tekniska Universitet: Luleå, Sweden, 2008. [Google Scholar]
  17. Blachowski, J.; Ellefmo, S. Numerical Modelling of Rock Mass Deformation in Sublevel Caving Mining System. Acta Geodyn. Geomater. 2012, 9, 379–388. [Google Scholar]
  18. Shabanimashcool, M.; Li, C.C. Numerical modelling of longwall mining and stability analysis of the gates in a coal seam. Int. J. Rock. Mech. Min. Sci. 2012, 51, 24–34. [Google Scholar] [CrossRef]
  19. Suchowerska, A.M.; Merifield, R.S.; Carter, J.P. Vertical stress changes in multiseam mining under supercritical longwall panels. Int. J. Rock. Mech. Min. Sci. 2013, 61, 306–320. [Google Scholar] [CrossRef]
  20. Lafont, G.; Gunzburger, Y.; Mitri, H.S.; Al Heib, M.; Didier, C.; Piguet, J.P. Influence of mining excavation on energy redistribution and rockburst potential. In Proceedings of the 23 World Mining Congress, Montreal, QC, Canada, 11 August 2013. [Google Scholar]
  21. Ozbay, M.U.; Ryder, J.A. The effect of foundation damage on the performance of stabilizing pillars. J. S. Afr. Inst. Min. Metall. 1990, 90, 29–35. [Google Scholar]
  22. Kaiser, P.K.; Tang, C.A. Numerical simulation of damage accumulation and seismic energy release during brittle rock failure part II: Rib pillar collapse. Int. J. Rock. Mech. Min. Sci. 1998, 35, 123–134. [Google Scholar] [CrossRef]
  23. York, G. Numerical modelling of the yielding of a stabilizing pillar/foundation system and a new design consideration for stabilizing pillar foundations. J. S. Afr. Inst. Min. Metall. 1998, 98, 281–297. [Google Scholar]
  24. Vieira, F.M.C.C.; Dede, T.; Maccelari, M.J.; Spottiswoode, S.M. Stabilizing and bracket pillar design to reduce seismicity. J. S. Afr. Inst. Min. Metall. 1998, 98, 11–17. [Google Scholar]
  25. Yavuz, H. Yielding Pillar Concept and its Design. In Proceedings of the 17th International Mining Congress and Exhibition of Turkey–IMCET, Antalya, Turkey, 19–22 June 2001. [Google Scholar]
  26. Watson, B.P.; Roberts, M.K.C.; Nkwana, M.M.; Kuijpers, J.; van Aswegen, L. The stress-strain behaviour of in-slope pillars in the Bushveld Platinum deposits in South Africa. J. S. Afr. Inst. Min. Metall. 2007, 107, 187–194. [Google Scholar]
  27. Jia, P.; Zhu, W.-c. Dynamic−static coupling analysis on rockburst mechanism in jointed rock mass. J. Cent. South Univ. 2012, 19, 3285–3290. [Google Scholar] [CrossRef]
  28. Adach, K. Impact of the Size of Remnants on Seismic and Rockburst Hazard in LGOM Mines. Ph.D. Dissertation, Wroclaw University of Science and Technology, Wroclaw, Poland, 2015. (In Polish). [Google Scholar]
  29. Adach-Pawelus, K.; Pawelus, D. Application of hydraulic backfill for rockburst prevention in the mining field with remnant in the Polish underground copper mines. Energies 2021, 14, 3869. [Google Scholar] [CrossRef]
  30. Rybak, J.; Khayrutdinov, M.; Kuziev, D.; Kongar-Syuryun, C.B.; Babyr, V.N. Prediction of the geomechanical state of the rock mass when mining salt deposits with stowing. J. Min. Inst. 2022, 253, 61–70. [Google Scholar] [CrossRef]
  31. Tang, C. Numerical Simulation of Progressive Failure and Associated Seismicity. Int. J. Rock. Mech. Min. Sci. 1997, 34, 249–261. [Google Scholar] [CrossRef]
  32. Wiles, T.D.; Lachenicht, R.; van Aswegen, G. Integration of deterministic modelling with seismic monitoring for the assessment of rockmass response to mining: Part I Theory. In Proceedings of the Fifth Intermational Symposium on Rockburst and Seismicity in Mines, Johannesburg, South Africa, 17–19 September 2001; pp. 379–387. [Google Scholar]
  33. Zhu, W.C.; Li, Z.H.; Zhu, L.; Tang, C.A. Numerical simulation on rockburst of underground opening triggered by dynamic disturbance. Tunn. Undergr. Sp. Technol. 2010, 25, 587–599. [Google Scholar] [CrossRef]
  34. Adach-Pawelus, K.; Butra, J.; Pawelus, D. An attempt at evaluation of the remnant influence on the occurrence of seismic phenomena in a room-and-pillar mining system with roof deflection. Studia Geotech. Mech. 2017, 39, 3–16. [Google Scholar] [CrossRef] [Green Version]
  35. Adach-Pawelus, K. Application of seismic monitoring and numerical modelling in the assessment of the possibility of seismic event occurrence in the vicinity of ore remnant. In IOP Conference Series—Earth & Environmental Science, Proceedings of the World Multidisciplinary Earth Sciences Symposium, Prague, Czech Republic, 3–7 September 2018; IOP Publishing Ltd.: Bristol, UK, 2019; Volume 221, pp. 1–10. [Google Scholar] [CrossRef]
  36. Wang, S.; Masoumi, H.; Oh, J.; Zhang, S. Scale-Size and Structural Effects of Rock Materials; Woodhead Publishing: Sawston, UK, 2020. [Google Scholar]
  37. Kashfia, M.; Shada, S.; Zivarba, D. Evaluation of sample scale effect on geomechanical tests. Pet. Res. 2022. [Google Scholar] [CrossRef]
  38. Hoek, E.; Brown, E.T. Practical estimates of rock mass strength. Int. J. Rock. Mech. Min. Sci. 1997, 34, 1165–1186. [Google Scholar] [CrossRef]
  39. Hoek, E.; Marinos, P. GSI: A geologically friendly tool for rock mass strength estimation. Presented at the ISRM International Symposium, Melbourne, Australia, 19–24 November 2000. [Google Scholar]
  40. Hoek, E.; Carranza-Torres, C.T.; Corkum, B. Hoek-Brown failure criterion—2002 edition. In Proceedings of the North American Rock Mechanics Society Meeting in Toronto, Toronto, ON, Canada, 7–10 July 2002. [Google Scholar]
  41. Sakurai, S. Direct strain evaluation technique in construction of underground openings. In Proceedings of the 22nd US Symposium on Rock Mechanics, Cambridge, MA, USA, 29 June–2 July 1981; MIT Press: Cambridge, MA, USA, 1981; pp. 278–282. [Google Scholar]
  42. Catalogue of Copper Deposit Mining Systems in KGHM Polska Miedz SA Mines; KGHM Polska Miedz SA Mines: Lubin, Poland, 2007. (In Polish)
  43. RocScience Phase2 Finite Element Analysis for Excavations and Slopes v. 8.0 Documentation; RocScience Inc.: Toronto, ON, Canada, 2015.
Figure 1. Method of determination of effective geomechanical pillar parameters on the basis of numerical simulations and measurement data using statistical methods.
Figure 1. Method of determination of effective geomechanical pillar parameters on the basis of numerical simulations and measurement data using statistical methods.
Energies 15 04729 g001
Figure 2. Mining field in the Polkowice-Sieroszowice mine where remnant has been left.
Figure 2. Mining field in the Polkowice-Sieroszowice mine where remnant has been left.
Energies 15 04729 g002
Figure 3. Room and pillar mining system with closing pillar J-UGR-PS [42].
Figure 3. Room and pillar mining system with closing pillar J-UGR-PS [42].
Energies 15 04729 g003
Figure 4. Room and pillar mining system with roof deflection J-UG-PS [42].
Figure 4. Room and pillar mining system with roof deflection J-UG-PS [42].
Energies 15 04729 g004
Figure 5. Calculation scheme for D-IE mining field.
Figure 5. Calculation scheme for D-IE mining field.
Energies 15 04729 g005
Figure 6. Simulation of room and pillar mining in the analyzed field: (a) step 1, (b) step 2, (c) step 4.
Figure 6. Simulation of room and pillar mining in the analyzed field: (a) step 1, (b) step 2, (c) step 4.
Energies 15 04729 g006
Figure 7. The results of excavation convergence measurements for J-UGR-PS mining system in Polkowice-Sieroszowice mine.
Figure 7. The results of excavation convergence measurements for J-UGR-PS mining system in Polkowice-Sieroszowice mine.
Energies 15 04729 g007
Figure 8. Fragment of the model of analyzed exploration field with marked division of pillars.
Figure 8. Fragment of the model of analyzed exploration field with marked division of pillars.
Energies 15 04729 g008
Figure 9. Chart of convergence determined numerically for chosen excavation in D-IE mining field for different combinations of parameter values assumed for pillars.
Figure 9. Chart of convergence determined numerically for chosen excavation in D-IE mining field for different combinations of parameter values assumed for pillars.
Energies 15 04729 g009
Figure 10. Intervals of confidence for a.
Figure 10. Intervals of confidence for a.
Energies 15 04729 g010
Figure 11. Fragment of the model of analyzed mining field with marked division of pillars.
Figure 11. Fragment of the model of analyzed mining field with marked division of pillars.
Energies 15 04729 g011
Figure 12. Results of convergence calculated numerically for the matrix of the experiment.
Figure 12. Results of convergence calculated numerically for the matrix of the experiment.
Energies 15 04729 g012
Figure 13. Confidence intervals for a.
Figure 13. Confidence intervals for a.
Energies 15 04729 g013
Figure 14. Matching of results of in-situ convergence measurements in Polkowice-Sieroszowice mine to convergence values obtained for effective parameters of pillars determined using statistical methods (elastic model).
Figure 14. Matching of results of in-situ convergence measurements in Polkowice-Sieroszowice mine to convergence values obtained for effective parameters of pillars determined using statistical methods (elastic model).
Energies 15 04729 g014
Table 1. Averaged geomechanical parameters of rock.
Table 1. Averaged geomechanical parameters of rock.
LocationName of Rockh [m]ρ [kg/dm3]Rc [MPa]Rr [MPa]Es [MPa]v [-]
ROOFMain anhydrite100.02.9093.16.456,1000.24
Loamy anhydrite breccia10.02.2536.01.713,6500.18
Basic anhydrite73.02.9095.55.554,6000.25
Calcareous dolomite I15.02.53132.58.351,0900.24
Calcareous dolomite II2.02.74213.016.099,3200.27
MINED DEPOSITMined deposit 2.72.63110.97.434,4500.21
BOTTOMQuartz sandstone I8.22.1222.11.481900.15
Quartz sandstone II194.51.9516.70.761900.13
Table 2. Parameters of rock mass assumed for numerical modeling.
Table 2. Parameters of rock mass assumed for numerical modeling.
LocationName of Rock h [m]Es [MPa]v [-]σt [MPa]c [MPa]φ [°]
ROOFMain anhydrite100.041,1100.240.7466.96738.66
Loamy anhydrite breccia10.071000.180.0932.50739.06
Basic anhydrite73.040,0100.250.7657.14638.66
Calcareous dolomite I15.044,9800.242.93312.08539.00
Calcareous dolomite II2.087,4400.274.71519.89539.00
MINED DEPOSITMined deposit 2.725,2400.210.8258.42439.31
BOTTOMQuartz sandstone I8.242600.150.0571.53839.06
Quartz sandstone II194.532200.130.0431.16039.06
Marking in the above tables: h—thickness of rock layers, ρ—volume density, Rc—rock sample uniaxial compression strength, Rr—rock sample tensile strength, Es—longitudinal elasticity modulus, v—Poisson coefficient σt—tensile strength of rock mass, c—coefficient of cohesion, ϕ—angle of internal friction.
Table 3. Parameters of pillars determined for the elastic model using Statistica v. 10.
Table 3. Parameters of pillars determined for the elastic model using Statistica v. 10.
ParameterEstimationpLower Confidence Limit Upper Confidence Limit
E2725.0660.739776−4582.56032.6
a5541.1450.893210−95,323.3106,405.6
E3135.8080.000182122.6209.0
Table 4. Pillar parameters set for the elastic model using the Statistica v. 10 program.
Table 4. Pillar parameters set for the elastic model using the Statistica v. 10 program.
ParameterEstimationpLower Confidence LimitUpper Confidence Limit
E2136.9310.0000000000119122.134151.729
a174.4240.0000162422638238.441391.832
Table 5. Pillar parameters set for the elastic model using the Statistica v. 10 program.
Table 5. Pillar parameters set for the elastic model using the Statistica v. 10 program.
ParameterEstimationpLower Confidence LimitUpper Confidence Limit
E1311.3550.0000364562557198.823423.887
E2136.9310.0000000000119122.134151.729
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adach-Pawelus, K. Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods. Energies 2022, 15, 4729. https://doi.org/10.3390/en15134729

AMA Style

Adach-Pawelus K. Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods. Energies. 2022; 15(13):4729. https://doi.org/10.3390/en15134729

Chicago/Turabian Style

Adach-Pawelus, Karolina. 2022. "Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods" Energies 15, no. 13: 4729. https://doi.org/10.3390/en15134729

APA Style

Adach-Pawelus, K. (2022). Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods. Energies, 15(13), 4729. https://doi.org/10.3390/en15134729

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop