CC BY
62
Vladimir N. Tyupin DOI: 10.31897/PMI.2019.2.167
Estimation of Critical Depth of Deposits ...
UDC 622.831.325:550:34
ESTIMATION OF CRITICAL DEPTH OF DEPOSITS BY ROCK BUMP
HAZARD CONDITION
Vladimir N. TYUPIN
Belgorod State National Research University, Belgorod, Russia
During the development of minerals by the underground method, dynamic manifestations of rock pressure occur at a certain depth, which significantly reduces the safety of mining operations. Regulatory documents prescribe at the exploration and design stages to establish the critical depth for classifying a deposit as liable to rock bumps. Currently, there are a number, mainly instrumental, methods for determining the liability of rock mass to rock bumps and methods based on the determination of physical and technical properties and the stress-strain state of rock massifs.
The paper proposes a theoretical method for determining the critical depth for classifying a deposit as liable to rock bumps. A formula for determining the critical depth of the rock bump hazard condition is obtained. A mathematical analysis of the influence of the physical and technical parameters of the formula on the critical depth is carried out. Its physical and mathematical validity is substantiated. The numerical calculations of the critical depth for 17 developed fields were carried out using a simplified formula. It also provides a comparison of calculated and actual critical depth values.
It is established that the variation of the actual and calculated critical depth is due to the lack of actual data on the value of the friction coefficient and parameters of fracturing of the rock mass in the simplified formula. A simplified calculation formula can be used to estimate the critical depth of a field at the survey and design stages.
More accurate results can be obtained if there are actual data on fracture parameters, friction coefficients and stress concentration near the working areas.
Key words: rock bump hazard; survey and design; critical depth; physical and technical properties of rocks; depth determination formula; numerical calculations
How to cite this article: Tyupin V.N. Estimation of Critical Depth of Deposits by Rock Bump Hazard Condition. Journal of Mining Institute. 2019. Vol. 236, p. 167-171. DOI: 10.31897/PMI.2019.2.167
Relevance and coverage of problem studies. Issues of rock-bump hazard arise at deep horizons during underground mining. Regulatory documents* state that the liability to rock bumps at certain mining depth must be determined established at the survey and design stages [1].
The analysis of Russian and foreign literature [1-8, 12-17] shows that liability of rock mass to rock bumps is determined by elastic and strength parameters of rocks, peculiarities of micro and natural fracturing (size of structure, crack width, number of crack systems, the angle of inclination to the outcrop), the parameters of tectonic faults and their location. The liability of rock mass to rock bumps is also determined by the magnitude of the gravitational pressure of rocks and the stress caused by the movement of the lithospheric plates of the Earth [8]. In addition, during the formation of workings, a redistribution of the natural stress state occurs with a significant increase in its proximity to the open surfaces of the rock mass. In addition, in the formation of cavities in the rock mass, some blasting operations lead to the formation of an impact and explosion residual stresses zones [10].
An important parameter of a deposit liable to rock bumps is the critical depth, i.e. the depth from the Earth's surface, at which rock bumps occur during mining, or it is labeled as «dangerous», which corresponds to the hazardous state of the rock mass in the marginal reservoir area, at which rock bump may occur.
The aim of the work is the theoretical determination of the critical depth of the deposit at which the rock bumps can occur during heading in a single goaf.
* Provision for safe mining operations in deposits liable and dangerous for rock bumps. Federal rules and regulations in the field of industrial safety. Order of Rostekhnadzor N 576. 02.12.2013 // URL: http: // www.consultant.ru/document/cons_doc_LAW_1616871
Journal of Mining Institute. 2019. Vol. 236. P. 167-171 • Mining
Vladimir N. Tyupin
Estimation of Critical Depth of Deposits
The initial mining of single (or parallel) workings, as a rule, is carried out during mining and exploration operations. It should be noted that in case of great angularity of rock mass the critical depth must be reduced.
Theoretical definition and numerical calculation of the critical depth of the rock mass depending on rock bump hazardousness. In [10], based on the energy conservation law, a theoretical formula for calculating the radius of the loose zone was obtained during drivage with blasting technique in the stressed fissured rock mass. The formula is obtained based on an experimentally proven mechanism of impact from the explosion of a group of charges of explosives in the stressed fissured rock mass, which consists of two stages.
Stage I. The explosion generates stress waves in the jointings around the charges, which ensures the destruction of cleavages with radial cracks. Further, under the action of the quasistatic pressure of the detonation products (DP), the broken particles begin to shift in the radial direction from the explosive charges, which is accompanied by deformations and friction between the faces of the remote cleavages, as well as their elastic deformations. In this case, the magnitude of deformations and stresses in separate cases decreases with distance from explosive charges, i.e. the strained array at this stage can be represented as an unevenly compressed spring with maximum compression near the charges and decreased compression at a distance.
Stage II. After the pressure drops in the charging planes, the reaction of elastically deformed (maximally stressed by the explosion) cleavages in combination with rock pressure leads to the movement of the massif towards the developed space with the formation of artificial or natural cracks approximately parallel to the goaf outline, i.e. the loose zone is formed.
The distance from the baring outlining to the far edge of the loose zone is determined from the following condition: the energy of the elastic response of the array, accumulated under the action of successive blasting from the center of the working to the periphery and the rock pressure, is equal to the energy required to overcome the friction forces between the cleavages in the marginal working zone:
where D, pb, d3 - detonation rate, powder-loading density and explosive charge diameter; c, v, ^ -longitudinal wave velocity, Poisson's ratio, separate massifs, friction coefficient between cleavages, respectively; Kn±, K±(N), Kop - indicators that take into account the interaction of explosive charges and open surface; E - cleavage elasticity modulus; A, de, k, P - accordingly, the amount of displacement of cleavages, the size of cleavages, the number of fracture systems, the angle of inclination of the fracture system to the outcropping; O - rock fracture index; P - rock pressure,
K - stress concentration factor created due to the formation of workings and tectonic plate movement; p - a specific weight of rock mass; g - acceleration of gravity; H - depth of the Earth surface.
The reliability of formula (1) was proved by industrial experiments at the mines of PJSC «Priargunskoye Proizvodstvennoe Ob'edinenie» (PPGHO) at a depth of 300-600 m from the Earth's surface. The results of calculations and practical data are shown in Fig. 1.
The critical depth is determined by the rock-bump hazard condition. It is known that with an increase in the distance from the Earth's surface, the stressed state of the rock mass increases and rock
(1)
P = KpgH,
(2)
168 -
Journal of Mining Institute. 2019. Vol. 236. P. 167-171 • MINING
ê Vladimir N. Tyupin
Estimation of Critical Depth of Deposits
bumps occur during blasting (about 60 %) and drilling (about 25 %) operations, i.e. with dynamic impact on rock mass. During a rock bump, as a rule, there is a pulsating displacement of the mass at depth [6, 7] and pressure bump directed towards the developed space. At the same time, «the energy of rock impact is formed from the energy of elastic deformations of a collapsing part of the rock mass and energy released from adjacent rocks» [6, 7]. During rock bump the impact zone is at a radial distance from the working outline, it significantly exceeds the determined size of the loose zone (Fig.1). Then it can be mathematically assumed that in this case, the denominator in (1) tends to zero and R3 to infinity.
Substituting (2) into (1) and equating the denominator to zero and solving the equation, we obtain the formula for determining the critical depth of the rock bump impact hazard for drivage of a single goaf:
"cr k V-V
Kpg (1 - v)(1 -|)de Z sin2 P,
¿=i
Analysis of dependence (3) shows that with an increase in the transverse strain coefficient, the compressive load is compensated by lateral expansion of the rock, i.e. with an increase in v, Hcr increases, which is mathematically reflected in (3). The friction between the faces of the cleavages holds the mass from bumping and, with increasing | the depth should be greater. The greater the volume of the overlying mass of rocks, the less is Hrf. Field observations on horizon 10 of mine No. 8 of PJSC PPGHO [11] found that the most frequent rock-bump hazard situations (dynamic roof rock slip) are observed when there are systems of cracks perpendicular to the outcrops of mine workings, which is mathematically obtained in formula (3). According to the data of [6, 7], in a rock mass with large cleavages, de the magnitude of rock pressure is greater, which
means that the probability of rock bump is higher. This mathematical analysis indicates the validity of the obtained formula for determining the depth of the dynamic manifestations of rock pressure when driving single workings. For further proof of the validity of the obtained formula (3), its numerical calculations and comparison with the data* were carried out for several fields in the Russian Federation, developed by underground mining. The physical-mechanical properties of rocks for calculations were selected in [9] and presented in the table as average values.
Numerical average values of parameters for calculations in the formula (3) are equal to the following: p = 2.6-103 kg/m3; g = 9.8 m/s2; de = 0.4-1.0 m (0.7 m in average); | = 0.3-0.6 (0.45 in aver-i k
age); A = 10-3 m; Z sin2 P, =1.75 at p, = 0; 60; 90°; K = 2.
i=1
* Provision for safe mining operations in deposits liable and dangerous for rock bumps. ...
- 169
Journal of Mining Institute. 2019. Vol. 236. P. 167-171 • Mining
Fig.1. Dependence of the distance between the working outline and the boundary of the loose zone (maximum of stresses) on the size of the cleavage, determined by the goaf walls in various rock massifs instrumentally (core disking, ultrasonic measurements, unloading method, parallel well method) (1) and theoretically (2) ( Formula 1)
1 - granites; 2 - felsites; 3 - conglomerates; 4 - trachydacites;
5 - tufogravellites
2EpvA
Vladimir N. Tyupin
Estimation of Critical Depth of Deposits
Assumed and calculated critical depth for rock bump hazard condition
Deposit Rocks and ores prone to brittle fracture ^ m Е, 1010 Pa v Нcr, m
Goroblagodatskoe Microsyenites, skarns, garnet-magnetite rocks 300 2.2 0.21 152
Korobkovskoe Iron Quartzites 600 7.0 0.24 576
Tashtagolskoe Syenites, skarns, tuffs, iron ore 400 8.0 0.2 588
Shereshegskoe Syenites, granites, hornfels 600 5.8 0.24 644
Krivorozhskoe Granites 500 6.0 0.23 468
Dzhezkazganskoe Sandstones 400 6.9 0.20 450
Zyryanovskoe Microquartzite 600 7.5 0.23 585
Mirgalimsayskoe Limestone, Dolomites 400 4.0 0.3 448
Tekeliyskoe Quartzits 600 6.0 0.24 495
Oktyabrskoe and Talnakhskoe Hornbacks, gabbro-dolerites 700 9.4 0.22 692
Khaydorkasnkoe Yuzhnoe Limestone 215 3.97 0.25 345
Tyrnauzskoe Skarns 800 6.0 0.25 523
Apatitoviy tsirk Iyolite-urtites, poor and rich ores 200 5.4 0.25 470
Partomchorrskoe Iyolite-urtites, poor and rich ores 400 5.0 0.25 435
Rasvumchorrskoe Iyolite-urtites, poor and rich ores 400 7.5 0.26 688
Kukisvumchorrskoe Iyolite-urtites, poor and rich ores 300 7.0 0.22 516
Streltsovskoe ore deposit Granites 500 5.3 0.23 413
(PJSC PPGHO)
Note. Нa - assumes depth*; ^ - calculated critical depth.
The simplified calculation formula is
#cr = 2,62 -10-8-^-. (4)
cr (1 - v)
Numerical calculations for (4) are presented in the table and in Fig.2.
The analysis of Fig. 2 shows that with an increase in the value of Ev(1 - v)-1, the critical depth increases, while the calculated values are practically on the straight line. The spread of actual values relative to the calculated curve is on average ± 150 m, which is quite significant. This is explained by the
fact that assumed unknown mean values of the parameters are taken: the coefficient of friction between the edges of cleavages, the magnitude of the displacement and the size of the cleavages, the number of fracture systems and their angle of inclination to the outcrop plane, the coefficient of stress concentration. Therefore, formula (4) can approximately estimate the critical depth of rock bump hazard. In this case, the elastic modulus and the Pois-son's ratio is determined in laboratory conditions by the core of exploration wells. The coefficient of friction correlates with the coefficient of the strength according to M.M.Protodyakonov [10].
* Provision for safe mining operations in deposits liable and dangerous for rock bumps. ...
170 -
Journal of Mining Institute. 2019. Vol. 236. P. 167-171 • MINING
Н„r, m
800 -
600
400
200 -
O- 1 • - 2
0,5
—1—1—1—1—1—1—1—1—1—1—1—1—1—1—1—r-
1 15 t C..//1 ,.\-1 m10i
2 Ev/(1-v)-1, 1010 Pa
Fig.2. Dependence of critical depth and rock bump hazard from the elastic properties of rocks
1 - actual values ; 2 - calculated values
0
é Vladimir N. Tyupin
Estimation of Critical Depth of Deposits
Conclusions
Based on the analysis of technical literature, theoretical studies, numerical calculations and comparing them with the development of hazardous deposits, we can draw the following conclusions:
1. A theoretical formula has been developed for calculating the critical depth of a deposit according to the rock bump hazard condition when driving a single working. The formula considers the physical-technical properties of the rock mass, the parameters of its fracturing, and stress concentration near the working.
2. Mathematical analysis and numerical calculations using the formula prove its validity. The scatter of the actual and calculated critical depth in the simplified calculation formula (4) is explained by the lack of actual data on the values of the friction coefficient and parameters of the rock mass fracture, and stress concentration the working.
3. A simplified calculation formula can approximate the critical depth of a deposit at the survey and design stages.
4. More accurate results can be obtained if there are actual data on fracture parameters, friction coefficients and stress concentration near the working.
REFERENCES
1. Aksenov A.A. Improving the practice of classifying deposits as liable to rock bursts. Bezopasnost' truda v promyshlennosti. 2018. N 1, p. 58-60. DOI: 10.24000/0409-2961-218-1-58-60 (in Russian).
2. Rasskazov I.Yu., Saksin B.G., Usikov V.I., Potapchuk M.I. Geodynamic state of the rock massif of the Nikolaev polymetal-lic deposit and the characteristics of the manifestation of rock bump hazard during its development. Gornyi zhurnal. 2016. N 12, p. 23-25. DOI: 10.17580/gzh.2016.12.03(in Russian).
3. Freidin A.M., Neverov S.A., Neverov A.A., Konurin A.I. Geomechanical evaluation of geotechnologies for underground mining of ores at the design stage. Gornyi zhurnal. 2016. N 2, p. 39-44. DOI: 10.17580/gzh.2016.02.08 (in Russian).
4. Kuranov A.D. The use of numerical modeling to select the safe parameters of ore deposits development systems in the highly stressed rock mass. Zapiski Gornogo instituta. 2013. Vol. 206, p. 60-64 (in Russian).
5. Lomakin V.S., Grigorovich S.V., Potekhin R.P., Khalevin N.I. On the relationship of the volume of the focal destruction zone with the seismic energy of rock bursts. Geologiya i geofizika. 1989. N 5, p. 129-132 (in Russian).
6. Petukhov I.M. Classification of rock bumps. Bezopasnost' truda v promyshlennosti. 1987. N 12, p. 41-43 (in Russian).
7. Petukhov I.M. On the nature of pulsing impact rock mass deformation. Gornyi zhurnal. 1989. N 7, p.45-48 (in Russian).
8. Rasskazov I.Yu. Control and management of rock pressure in the mines of the Far Eastern region. Moscow: Gornaya kniga, 2008, p. 329 (in Russian).
9. Reference (inventory) of the physical properties of rocks. Pod red. N.V.Mel'nikova, V.V.Rzhevskogo, M.M.Protod'yakonova. Moscow: Nedra, 1975, p. 279 (in Russian).
10. Tyupin V.N. Explosive and geomechanical processes in the fractured strained rock mass. Belgorodskii natsional'nyi issle-dovatel'skii universitet. Belgorod, 2017. p. 192 (in Russian).
11. Wang N., Wan B.H., Zhang P., Du X.L. Analysis on deformation development of open-pit slope under the influence of underground mining. Proceedings of International Symposium on Land Reclamation and Ecological Restoration. Beijing. China. 2015, p. 53-58.
12. Braun L.G. Seismic hazard evaluation using apparent stress ratio for mining-induced seismic events. Ph. D. Thesis, Laurentian University, 2015, p.257.
13. Ingraham M.D., Issen K.A., Holcomb D.J. Use of acoustic emission to investigate localization in high-porosity sandstone subjected to true triaxial stresses. Acta Geotechnica. 2013. Vol. 8. N 6, p. 646-663.
14. Marcak M., Mutke G. Seismic activation of tectonic stresses by mining. Journal of Seismology. 2013.Vol. 17. N 4, p. 1139-1148.
15. Snelling P.E., Godin L., McKinnon S.D. The role of geologic structure and stress in triggering remote seismicity in Creigh-ton Mine, Sudbury, Canada. International Journal of Rock Mechanics and Mining Sciences. 2013. Vol. 58, p. 166-179.
16. Young D.P. Energy variations in mining-induced seismic events using apparent stress. MASc Thesis. Laurentian University, 2012, p. 85.
17. Wesseloo J., Woodward K., Pereira J. Grid-based analysis of seismic data. The Journal of Southern African Institute of Mining and Metallurgy. 2014. Vol. 114, p. 815-822.
Author: Vladimir N. Tyupin, Doctor of Engineering, Professor, [email protected] (Belgorod State National Research University, Belgorod, Russia).
The paper was received on 18 March, 2018.
The paper was accepted for publication on 25 January, 2019.
Journal of Mining Institute. 2019. Vol. 236. P. 167-171 • Mining
