Displacement of mined ground as a consequence of the exploitation of Pb-Zn ore in the mine Crnac Pomik podkopanega terena kot posledica pridobivanja rude Pb-Zn v jami Crnac Blagoje Nedeljkovic1, *, IvicA RiSTovic2, Gabriel Fedorko3 University of Priština, Faculty of Technical Sciences, Kneza Miloša 7, 38220 Kosovska Mitrovica, Serbia ^University of Belgrade, Faculty of Mining and Geology, Djušina 7, 11000 Belgrade, Serbia 3Technical University of Košice, BERG Faculty, Letna 9, 04200 Košice, Slovakia Corresponding author. E-mail: blagojenedeljkovic@yahoo.com Received: August 05, 2009 Accepted: August 19, 2009 Abstract: Underground exploitation of mineral raw materials and the construction of sub-surface facilities cause changes on the ground surface. These changes on the ground surface are practically manifested as specific displacements and deformations on the surface of the mined ground. The objective of this Paper is to determine a safe working depth and to make an assessment of a prospective deformation hazard and subsidence of the ground over the mining zones of Pb-Zn ore veins in the Mine Crnac. A high-quality determination, forecasting of deformations and subsidence of a terrain is very important for the safety of underground mine workings and the protection of facilities on the surface of the ground. Izvleček: Podzemno pridobivanje mineralnih surovin, kot tudi gradnja podzemnih objektov povzročata spremembe na površini. Praktično se te spremembe na površini izražajo kot pomiki in deformacije na območju odkopa. Cilj članka je določanje varne globine odkopa-vanja ter ocena potencialne nevarnosti zaradi deformacij in ugrezanja nad območji pridobivanja Pb-Zn rudnih žil v jami Crnac. Kakovostno določanje, napovedovanje deformacij in ugreznin terena je zelo pomembno za varnost podzemnih jamskih prostorov ter zaščito objektov na površini. 450 Nedeljkovic, B., Rlstovic, I., Fedorko, G. Key words: working environment, underground works, displacement of the ground, deformation, angular parameters, analytic dependence Ključne besede: delovno okolje, podzemna dela, pomik terena, deformacija, vplivni koti, analitična odvisnost Introduction The extraction of lead-zinc ore (Pb-Zn) in the mine Crnac has been made intensively since 1967, year which is taken as the beginning of the modern exploitation. By investigation mining at the height of 862 m there were over 20 Pb-Zn ore veins identified, with thickness from 1 m to 3.5 m. Depending on the physical and mechanical characteristics of ore and accompanying rocks of ore veins, extraction methods for each ore vein are applied separately. For the extraction of ore veins the following methods are applied: extraction with back-filling of cavities, caving method and sublevel open stope mining method.[6] This paper shows the impact of Pb-Zn ore exploitation in the mine Crnac on the displacement of mined ground in the exploitation by open stope method on the example of the ore vein No. 3. Characteristics of a working environment In a geological environment of a deposit the following lithologic members are involved: amphibolites, Palaeozoic shales, serpentinites, diabases, diabasehornstone series and tertiary effusives with their pyroclastites. Tertiary mag-matic activity was manifested by the formation of significant masses of effusive rocks and a number of changes in the active rocks. The same magma-tism gave rise to dumping of the Pb-Zn mineralization in the form of ore veins with 1-3.5 m thickness and 60-900 angles of occurrence (dipping) and with coefficients of solid ore and accompanying rocks (/ 5-15) according to Pro-todjakonov. Impact of underground works on the surface of the ground As a consequence of the underground mineral deposit exploitation, it comes to roof caving above the working cavity, which can frequently be manifested on the surface of the ground.[9] The first signs of ground subsidence are manifested by deflections of the terrain, subsidence or caving of a stope roof. Moreover, the subsidence is reflected by an increased pressure in the roof and sidewall.[8] The underground mining leads to vertical and significantly mere horizontal displacements and deformations on the surface of the terrain. The shape, size and the process of terrain deformations depend on a number of factors where the main are the following:[5] • Structure of the characteristic of a rock massif (fissuring, bedding), • Physical and mechanical characteristics of rocks which constitute the ore massif, • Dipping conditions of ore bodies and accompanying rocks, • The shape and the size of ore bodies, their thickness, ratio between the size of the working cavity and the depth of works, • Mining system, • Damage degree of rocks, • Terrain relief. The above factors in each individual case define particularities of the displacement and deformation processes of the massif, by enabling the application of analytic, graphic and numerical and analogue methods for defining and studying the processes of deformation of massif and the terrain surface.[10] The Figure 1 shows the diagram of displacement and deformation zones on the surface of the mined ground and angular parameters and displacement assessment. Impact zone Displacement zone Deformation zone Caving zone Figure 1. Diagram of displacement and deformation zones on the surface of the mined terrain The first researches related to the im- Impact of the inclination angle of an pact of underground works on the sur- ore body a face of the terrain were from the late For the open stope method (caving 19th century. These issues were dealt by method and sublevel open stope mina number of researchers who gave their ing method) it is found that the inclina-indisputable contribution to the scien- tion angle a of an ore body is the factor tific thought in this field and served as which influences the most on the angle an instrument for developing new theo- fi of the deformation impact zone: ries in the contemporary mining, based on the mathematical processing of a fi = -37.27 + 1.37 a (n = 0.881; p = large number of collected data.[1] [2] [11] 9.68 ), fi = 0.002 a2 + 0.296 a + 35.572 (n = Criteria for determining angular 0.682; p^ =6.89), PARAMETERS B ACCORDING TO I. N. KISI- mov fi" = -0.052 a2 + 8.53 a - 264.898 (n = 0.596; pn = 6.47 ), Applying the theory of probability and mathematical statistics methods, and fi\ = 0.004 a2 - 0.070a + 57.572 (n = on the basis of measurements complet- 0.675; p^ = 10.55 ). ed, analytical dependencies and criteria for determining angular parameters Impact of the strength of accompa-were proposed (Figure 1.) when resolv- nying rocks / ing problems related to displacements Impact of the coefficient f on the an-of a mined terrain due to the impact of gle fi is significant and is considered as mining works. The basic criteria for a the approximate to the impact of an ore forecast calculation of displacements body dip a, is given in the form of the and deformations of a mined terrain expression for open stope methods: and angular parameters of the displacement process were given by the author fi = -17.07 + 0.97 a + 0.93 f (R = 0.896; I. N. Kisimov:[3] pR = 11.23), • Impact of the inclination angle of an ore body a, fi= -79.03 + 1.63 a + 1.62 f (R = • Impact of the strength of accompa- 0.847; pR = 7.30), nying rocks f, • Impact of the actual thickness of an fi" = 0.024 f2 - 0.086 f + 76.955 (n = ore body m, 0.484; p= 39.2 with f = 6-18), • Impact of the depth of mining works H. = -0.087 f2 - 2.456 f + 63.031(n = 0.589; un = 7.96 with f = 8-18). Impact of the actual thickness of an ore body m Impact of the thickness m on the angle ( is certain to the itself impact of the coefficient f: ( = 67 - 0,42 m + 0,81 f (with R = 0.375 and uR = 3,59 for m = (1:3m) Impact m to the factor a: ( = 28.01 + 0.55 a + 0.29 (R = 0.509; Ur = 4.28) Impact m to the factor H: ( = 66.52 + 0.04 H + 0.83 m (R = 0.529 and uR = 5.04 ( = 77.6 + 0.19 m + 0.38 f (R = 0.329 and uR = 3.09). Impact of the depth of mining works H The change in the angle ( depending on the depth H is typical for the decrease in value ( at depths from 50 m to 150 m and the increase in the angle ( in the interval from 150 m to 300 m; therefore the dependency of the impact H and f on the angle ( is determined: ( = 7.148 + 0.33 f - 0.004 H with f = 6-18 and H = 100-400 m ( = 0.0003 H2 - 0.094 H + 76.71 (n = 0.784 and un = 13 63) and value H = 15-400m ' ( = 74.25 + 0.74 f - 0.02 H with R = 0.368 and uR = 3.38 (T = - 0.0005 H2 + 0.152 H + 68.162 Determining the stability of the TERRAIN surface due To THE impact of mining works in Mine crnac Previous mining activities in the pit of Crnac Mine, more ore veins were excavated above the horizon H = 862 m (Figure 2). Extraction methods which were applied to ore mining in particular ore veins of different thickness (from 1-3.5 m) and dipping angle of ore veins from 60-900, were not adapted to mining conditions, and therefore it led to displacements of terrain surface above working cavities. Figure 3 shows a longitude cross-section of the stope in the ore vein No. 3 and its relation to the terrain surface where there was the displacement of the terrain surface (Figure 4.). In order to prevent displacements of the terrain surface above the working cavity it is necessary to determine the impact of mining on the surface. In case that mining works have impact on the surface, it is necessary to define Figure 3. Longitude cross-section of a stope in the ore vein No. 3 and its relation to the terrain surface Figure 4. Displacement of the terrain surface zones of displacements and caving and a possibility of creating a safe depth, and when there are not any, it is important to know if there are stable areas directly above mining works.[4] [7] The stability of surfaces directly above mining works depends on the rock excavation method and natural and technical conditions, such as: occurrence angle of the ore body a, and the strength coefficient of accompanying rocks f, mining lengths at the dip L, with mining thickness m, depth to the upper border of exploitation H, lengths of excavation in the direction of strike N, with the thickness of covering detritus hp the thickness of covering main rocks h2. According to KISIMOV there is a directly proportional interdependence between the surface stability and the parameters a, H, f, hv h2 and L, m, N respectively. The impact of these parameters may be determined when necessary data are analyzed and expressed through empirical mathematical dependencies. Certainly, the stability of the surface in mining is influenced also by fissuration of accompanying rocks, considering that all the rocks are fissured to a small or large extent, it can be accepted that the fissuration takes part through the value of empiric coefficients. Patterns for minimal extraction depths where there are no displacements on the surface of a working cavity, bearing in mind the thickness of the ore body m, strength coefficient of roof beds f and excavated length at strike N, for the conditions of the Mine Crnac are given as follows: • Impact of the excavated thickness: XT ^ 25,4-(L-cosa + /n-sina) Hm > -,-^--— • m [L • cosa + m • sina)+ 3,6 • m 25,4-(125-0,174 + 3-0,985) > (125 • 0,174 + 3 • 0,985)+ 3,6 • 3 •3 Hm > 53,02m m ' • Impact of the strength coefficient: jj > 5.6• (L • cosa + m-sin«) ^ / (L • cos a + m • sin a)+1.9 • / > 5,6-(125-0,174 +3-0,985) m ~ (125 • 0,174 + 3 • 0,985)+1,9 -8 Hf > 21,1m • Impact of the excavated length at strike: 9.3 • (L • cosa + m • sina) Hn> H„ > {L- cosa + m -sina )+5- 9,3-(125-0,174+ 3-0,985) (125-0,174 + 3-0,985)+5-600 Hm > 45,51m On the basis of these calculations, we obtained three different values for a minimal extraction depth, the maximum value is adapted, which is the dislocation due to the impact of the strength of the ore vein No. 3 of the Crnac Mine. Conclusion Applying the KISIMOV's pattern along with introducing local parameters of the deposit, gave real, approximate values for determining a safe exploitation depth or the possibility of the appearance of falling-in on the surface of the terrain. Under the conditions in the Mine Crnac we obtained the results which indisputably show that a safe mining depth had not been defined, so the works on the ore extraction on the example of the ore vein No. 3 had caused the displacement of the terrain surface. Bearing in mind that the ore exploitation above horizon N0 : 862 m in the mine Crnac and further open stope mining method, it is necessary to define the impact of extraction on the surface of the terrain for each ore vein separately. References [1] Djordjevic, D. (2006): Methods to forecast and calculating movements and deformations Displacement of mined ground. Faculty of Mining and Geology, Belgrade. [2] Elezovic, D., Nedeljkovic, B., Jaksic, M. (2001): Mechanics of rocks and soil. University of Pristina, Faculty of Technical Sciences, Kosovska Mitrovica. [3] Kisimov, I. (1971): Prognozirane ustolicivosta na zemnata povrh-nost namiranate se nad minite raboti na D. M. P. Rudodobiu i metalurgija. Sofija. [4] Milentijevic, G., Nedeljkovic, B., Jaksic, M. (2005): The protection of the surface flows from the impact of the mine waters of the ore deposit Koporic-Zuta Prlina-Jelakce Satorica. Underground Mining Engineering, XII, No. 14, pp. 49-57, Faculty of Mining and Geology, Belgrade. [5] Nedeljkovic, B., Milentijevic, G., Jaksic, M. (2008): Determination of static and dynamic module of the elasticity of rock samples. Makedonsko rudarstvo i geologija, No. 11, FYR of Macedonia. [6] Nedeljkovic, B., Milentijevic, G. (2006): Estimation of endanger-ment of surface and ground waters of the Ibar's middle river basin as a result of RMHK Trepca activity. Underground Mining Engineering, XIII, No. 15, pp. 61-68, Faculty of Mining and Geology, Belgrade. [7] NEDELJKOVIC, B., MlLENTIJEVIC, G., Lazic, M. (2007): Protection of the environment in the inactive industrial areas. Proceedings of 1st Conference Environmental protection in industrial areas, pp. 8-25, University of Pristina, Faculty of Technical Sciences, Kosovska Mi-trovica. [8] Ristovic, I., Slepcev G., Colic, M., BoGDANovic, D. (2001): Environmental Protection - Dependence of Level in Mining Production. Proceedings of 3rd international Symposium Mining and Environmental Protection, MEP 01, pp. 69-73, Faculty of Mining and Geology, Vrdnik. [9] VUKELIČ, Ž., SpORIN, J., VlŽINTIN, G. (2004): Pore pressure. RMZ-ma-ter. geoenviron., 2004, vol. 51, no. 4, pp. 2117-2125, Ljubljana. [10] Vižintin, G., stevanovič, L., vukelič, Ž. (2008): Development of environmental criteria for estimation of land development using GIS = Uporaba GIS-a za določitev naravnih kriterijev možnosti okoljske-ga razvoja. RMZ-mater. geoenviron., Vol. 55, No. 2, pp. 237-258, Ljubljana. [11] Vulic, M., Korelc, J. (2008): The use of the logistic function for forecasting vertical movements of surface. RMZ-mater. geoenviron., 2008, Vol. 55, No. 3, pp. 389-395, Ljubljana.