sustainability
Article Mine Size Effects on Coal Pillar Stress and Their Application for Partial Extraction
Yang Yu 1 ID , Ka-Zhong Deng 1,* ID and Shen-En Chen 2 ID
1 School of Environment Science and Spatial Information, China University of Mining and Technology, Xuzhou 221116, China; [email protected] 2 Department of Civil and Environmental Engineering, University of North Carolina at Charlotte, 9201 University City Blvd., Charlotte, NC 28223, USA; [email protected] * Correspondence: [email protected]; Tel.: +86-0516-8388-5986
Received: 30 January 2018; Accepted: 9 March 2018; Published: 13 March 2018
Abstract: Coal is a nonrenewable resource. Hence, it is important to improve the coal recovery ratio and ensure the stability of coal mines for sustainable development of mining cities. Partial extraction techniques, such as strip pillar mining or room-and-pillar mining, are efficient methods to extract coal. Pillar stress is a critical property for pillar design and for the assessment of mine stability after partial extraction. Current pillar stress calculation methods can sometimes overestimate the pillar stress and unnecessarily large coal pillars may be left underground, which leads to a waste of coal resources. In this paper, the size effects of mining activity on the maximum vertical pillar stress were investigated using numerical simulations. Both strip pillar mining and room-and-pillar mining were considered as possible mining scenarios at different mining depths. The results show that the maximum pillar stress of a mine is primarily controlled by four factors: the mine size to mining depth ratio, the mining width to pillar width ratio, the overburden elastic modulus, and the mining depth. The maximum pillar stress of a mine gradually increases to an ultimate value as the mine size increases. Simplified formulas and methodology have been derived for stress calculations under consideration of mine size effects and, therefore, can reduce the waste of coal resources from the overestimation of pillar stress.
Keywords: strip pillar mining; room-and-pillar mining; pillar stress; pillar design
1. Introduction Underground mining activities can result in severe ground subsidence and damage man-made structures (buildings, railways, or other infrastructure) on the surface [1–5]. Mining subsidence is an important factor that restricts the sustainable development of mining cities [5]. To reduce mining subsidence and protect surface and ground structures, partial extraction methods such as strip pillar mining or room-and-pillar mining are widely adopted [5–7]. Figure1 shows typical layouts of room-and-pillar mining and strip pillar mining techniques. In a partial extraction operation, massive coal pillars will be left underground to support the overburden, where the effectiveness of subsidence control for the partial extraction will depend on the stability of these pillars. Additionally, the failure of a coal pillar may result in violent ground movements [8,9]. Therefore, the stability of the coal pillars in a partial extraction is very important during and after mining activities. The pillar stress is a critical factor affecting pillar stability, and the pillar stress calculation is an important aspect of mine design. The stress calculation for room-and-pillar mining is a three-dimensional problem because both the pillar length and the pillar width are relatively small (Figure1a), whereas the stress calculation for strip pillar mining can be simplified to a two-dimensional problem because the pillar length is usually far longer than the pillar width (Figure1b). In coal mine
Sustainability 2018, 10, 792; doi:10.3390/su10030792 www.mdpi.com/journal/sustainability SustainabilitySustainability2018 2018,,10 10,, 792x FOR PEER REVIEW 22 of 12 11
design and pillar stability evaluation, the pillar stress is usually determined by using TAT (tributary designarea theory), and pillar PAT stability (pressure evaluation, arch theory), the pillaror numerical stress is simulations usually determined [5,10–16]. by TAT using assumes TAT (tributary that the areaoverburden theory), load PAT will (pressure be evenly arch distributed theory), or numericalamong coal simulations pillars, and [ 5the,10 –overburden16]. TAT assumes load above that thethe overburdenmined room load is equally will be shared evenly by distributed neighboring among pillars coal [5]. pillars, On the and other the overburdenhand, PAT loadassumes above that the a minedpressure room arch is will equally form shared between by large neighboring barrier pillarspillars [and5]. Oncarry the most other of hand, overburden PAT assumes weight, that the aabutment pressure angle arch willis utilized form betweento describe large the barriershape of pillars the pressure and carry arch most [5,15]. of overburdenFor narrow mine weight, panels, the abutmentTAT overestimates angle is utilized the pillar to describestresses since the shape pressure of the arches pressure may archcarry [5 the,15]. weight For narrow of the mineoverburden. panels, TATIn such overestimates cases, PAT is the usually pillar adop stressested sinceto calculate pressure pillar arches stress may [5,14–16]. carry the However, weight the of the abutment overburden. angle Inused such by cases, PAT is PAT mostly is usually summarized adopted from to calculate longwall pillar mini stressng [5,10–20], [5,14–16 which]. However, is not thealways abutment accurate angle for usedpartial by extraction. PAT is mostly On the summarized other hand, from the longwall shape and mining height [5 ,of10 –the20 ],pressure which is arch not alwaysmay change accurate during for partialexcavation extraction. [21], indicating On the other that hand,the portions the shape of andthe overburden height of the load pressure that transfer arch may to changethe pillars during are excavationaffected by [the21], mine indicating size. Therefore, that the portions the existing of the methods overburden fail to loadprovide that a transfer reliable tostress the pillarsestimation are affectedfor partial by theextraction. mine size. Inaccurate Therefore, pillar the existingstress calculations methods fail usually to provide result a reliablein unnecessarily stress estimation large coal for partialpillars extraction.that are permanently Inaccurate left pillar underg stressround; calculations the coal usually resource result is wasted. in unnecessarily large coal pillars that areTo permanentlyaddress this leftissue, underground; empirically theestablishe coal resourced techniques is wasted. and methodologies which estimate pillarTo stress address have this been issue, suggested empirically [22]. established He et al. techniques[23] observed and that methodologies the mine subsidence which estimate and pillarstrata stressmovement have beenare affected suggested by [the22]. mine He et size al. [23. ]Additionally, observed that they the minefound subsidence that the mine and size strata movement to mining aredepth affected ratio by the(Figure mine 1c, size D⁄ .) Additionally, is an important they parameter found that to the det mineermine size whetherD to mining the subsidence depth H ratio can (Figurereach its1c, peakD/H value.) is an This important observation parameter suggests to determine that the deformation whether the subsidenceand stress of can strata reach may its peakbe a value.function This of observationthe mine size suggests . Roberts that et the al. deformation [11] reported and similar stress observations of strata may that be the a function pillar stress of the is mineaffected size byD . Roberts⁄ and et proposed al. [11] reported a stress similarcalculation observations formula for that room-and-pillar the pillar stress mines is affected that included by D/H andmine proposed size effects. a stress calculation formula for room-and-pillar mines that included mine size effects. InIn this paper, paper, the the effects effects of of mine mine size size variation variation on on coal coal pillar pillar stress stress are investigated are investigated for both for bothstrip strippillar pillar mining mining and androom-and-pillar room-and-pillar mining mining using using numerical numerical software software UDEC UDEC (Universal DistinctDistinct ElementElement Code)Code) and and 3DEC 3DEC (3 Dimension(3 Dimension Distinct Distinct Element Element Code). Code). Simulation Simulation of mines of includes mines differentincludes pillardifferent sizes pillar and differentsizes and minedifferent depths. mine The depths. mine sizeThe refersmine size to the refers mined to areathe mined of the area coal seamof the with coal coalseam pillars with coal and pillars mined and rooms mined (Figure rooms1c). (Figure Simple 1c). formulas Simple to formulas calculate to the calculate maximum the maximum pillar stress pillar for differentstress for mine different sizes mine are derived sizes are and derived proposed and for pr mineoposed design. for mine Based design. on the proposedBased on formulas,the proposed the recoveryformulas, ratio the recovery may be improved ratio may ifbe the improved mine size if isthe not mi toone size large, is not and too the large, waste and of coal the resourcewaste of coal can beresource reduced. can be reduced.
FigureFigure 1. TypicalTypical layout layout of of partial partial extraction extraction sites sites [5]: [ 5(a]:) Planar (a) Planar view view of the of working the working face of face a room- of a room-and-pillarand-pillar mine; mine; (b) Planar (b) Planar view view of the of theworking working face face of ofa astrip strip pillar pillar mine; mine; ( c) Two-dimensionalTwo-dimensional schemescheme ofof aa partialpartial extractionextraction sitesite (front(front viewview ofof aa partialpartial extractionextraction mine). mine).
Sustainability 2018, 10, 792 3 of 12
2. Mine Modeling and Simulation The discrete element modeling program, UDEC, is used in the current study. The scheme of the model is shown in Figure1c. The ground surface and the coal seam are modeled as horizontal strata, where the thicknesses of the coal seam and the mine floor are modeled as 7 m and 20 m, respectively. The mining depth H is varied from 100 m to 800 m; thus, the models included both shallow and deep mines. The pillar width wp is assumed to range from 7 m to 134 m. The mining width wc of strip pillar mining is usually 0.1 H–0.25 H to control the ground subsidence, thus, the pillars with very large sizes are proposed mainly for the theoretical analysis of strip pillar mining at deep depths. The ratio of mining width to pillar width rm is varied from 0.75 to 2. The geometric parameters of the models are listed in Table1. Two boundary pillars with width of 200 m are reserved at both ends of the model. Other element sizes include 1 m for pillars, 2 m for the floor, 2 m for the roof (thickness is 20 m), and 4 m for other parts.
Table 1. The geometric parameters of models: H is mining depth, wp is pillar width, wc is mining width, D is mine size, and rm is the ratio of mining width to pillar width. All lengths are in meters.
H wp wc D rm H wp wc D rm H wp wc D rm 14 10 1066 0.71 17 20 1093 1.18 300 75 1575 1 12 10 516 0.83 14 20 1006 1.43 54 40 2108 0.74 10 1510 1 10 20 1190 2 40 1560 1 9 10 656 1.11 7 20 1073 2.86 20 40 2080 2 5 10 670 2 23 1587 45 1575 12 1620 200 25 1525 50 1550 100 13 1573 29 1537 1 400 57 1539 1 15 1575 1 34 1530 67 1541 17 1547 41 1517 80 1520 20 1500 67 50 1454 0.75 50 100 2050 2 34 25 1146 0.74 50 1550 1 100 1500 1 25 1525 1 25 50 1475 2 134 100 2206 0.75 13 25 1127 1.92 15 30 1650 2 50 1550 7 1127 20 30 1680 1.5 56 1512 8 1288 25 30 1735 1.2 63 1575 9 1449 30 1530 1500 72 1512 10 1510 1 35 30 1785 0.86 84 1596 1 11 1507 300 40 30 1780 0.75 100 1500 200 13 1521 34 1530 125 1625 15 1515 38 1558 600 60 3060 1 27 20 1195 0.74 43 1505 700 70 3570 24 20 1076 0.83 50 1550 800 80 4080 20 1500 1 60 1500
In the model, all non-coal strata are assumed to consist of one type of rock. The properties of non-coal strata are chosen to mimic sandstone except in Section 3.3. In Section 3.3, the elastic modulus and Poisson’s ratio are changed while the other properties of non-coal strata refer to sandstone to study the effect of different rock properties on pillar stress. Rock and coal properties are taken from the Yulin mine area in China and are listed in Table2[ 24]. The Mohr–Coulomb failure criterion is assumed for the models. Since the joint stiffness only affects the joint deformation, a large stiffness (5 × 109 Pa) is assigned to the interface between non-coal strata and the coal seam, so that the maximum stresses in pillars can be calculated. The strength parameters of the interfaces are matched to the coal strength. The model base is fixed in the vertical direction and the lateral boundary of the model is fixed in the horizontal direction. Static equilibrium analysis is conducted with consideration of gravitational loads on the system. Every pillar contained one vertical observation line located at the pillar center. Sustainability 2018, 10, 792 4 of 12
Sustainability 2018, 10, x FOR PEER REVIEW 4 of 11 The model is mined one tunnel at a time and from the left-hand side until the mine reached the design size (Table1). For each time the tunnel number (or mine size D) increases, the maximum vertical pillar pillar center. The model is mined one tunnel at a time and from the left-hand side until the mine stress of the whole mine is recorded. reached the design size (Table 1). For each time the tunnel number (or mine size ) increases, the maximumTable 2.verticalMechanical pillar properties stress of ofthe strata: wholeσT ismine the cohesion,is recorded.σC is the tensile strength, ϕ is the internal friction angle, E is the elastic modulus, υ is the Poisson’s ratio, and ρ is the density. Table 2. Mechanical properties of strata: is the cohesion, is the tensile strength, is the ◦ 3 internalStrata friction angle,σT (MPa) is theσ elasticC (MPa) modulus,ϕ ( is) the Poisson’sE (GPa) ratio, and υ is the density.ρ (kg/m )
Sandstone 6.36 11.6 34.6 15.1 0.153 2370 Strata (MPa) (MPa) (°) (GPa) (kg/m ) Coal 0.72 2.61 42 2.2 0.27 1290 Sandstone 6.36 11.6 34.6 15.1 0.15 2370 Interface 0.72 2.61 42 Coal 0.72 2.61 42 2.2 0.27 1290 Interface 0.72 2.61 42 3. Analysis of Simulation Results
3. AnalysisMine failures of Simulation did not occurResults in all of the numerical simulations. The numerical simulation results indicate that the maximum vertical stress of a single pillar appeared at the pillar bottom, and the Mine failures did not occur in all of the numerical simulations. The numerical simulation results maximum vertical pillar stress of the mine is dominated by the mine size D, the mining width to pillar indicate that the maximum vertical stress of a single pillar appeared at the pillar bottom, and the width ratio r , and the overburden elastic modulus E . In this study, the maximum vertical pillar maximum verticalm pillar stress of the mine is dominatedo by the mine size , the mining width to stress of the mine is expressed by the stress concentration factor K: pillar width ratio , and the overburden elastic modulus . In this study, the maximum vertical pillar stress of the mine is expressed by the stress concentrationS factor : K = , (1) ρ gH + ρ gh =1 2 , (1) where g is the acceleration due to gravity, ρ and ρ are the density of non-coal strata and the density where is the acceleration due to gravity,1 and2 are the density of non-coal strata and the of the coal seam, respectively; H and h are mining depth and pillar height, respectively, and S is the density of the coal seam, respectively; and are mining depth and pillar height, respectively, and maximum vertical pillar stress of the mine from the simulation. S is the maximum vertical pillar stress of the mine from the simulation. 3.1. Effects of the Mine Size on the Maximum Vertical Pillar Stress of the Mine 3.1. Effects of the Mine Size on the Maximum Vertical Pillar Stress of the Mine By fixing rm to unity, the mine size effect on K was investigated. The results (Figure2) show that thereBy is anfixing exponential to unity, relationship the mine between size effectK and on D/ wasH: investigated. The results (Figure 2) show that there is an exponential relationship between and ⁄ : D − K = KS + ae H , (2) = + e , (2) wherewhere K S is is the the ultimate ultimate stress stress concentration concentration factor factor when wh theen the pillar pillar stress stress is stable is stable and reachesand reaches its peak its valuepeak value (Figure (Figure2); e is 2); the e natural is the natural base, and base,a is and a constant is a constant dependent dependent on H, rm on, and , E o., and .
Figure 2.2. Example of relationships between the the stress concentration factor factor K and the mine size toto mining depth ratio D ⁄/H for the mining width width to to pillar pillar width width ratio ratio r m = 1: ( (aa)) K ofof modelsmodels withwith a mining heightheight ofof 100100 m m and and a a pillar pillar width width from from 10 10 m tom 25to m;25 (m;b) the(b) regressionthe regression curve curve for the for case the withcase miningwith mining depth depth 200 m 200 and m pillar and pillar width width 10 m. 10 m.
Figure2 2 shows shows examplesexamples ofof K distribution distribution (Figure(Figure2 2a)a) and and the the regression regression results results (Figure (Figure2 b)2b) as as a function of D ⁄/H.. Initially, Initially, K increasesincreases when D /⁄H increases, and and gradually, gradually, K become become stablestable andand reaches itsits peakpeak valuevalue whenwhen D /⁄H ≥≥3 to 4. Therefore, the parameter D ⁄/H can be an indicationindication to determine if the pillar stress in the mine reaches its maximum value; if ⁄ <3 to 4, the pillar stress of the mine is not fully developed and will increase with further increase of ⁄ . Otherwise, reaches its maximum value and will no longer increase with the increase of ⁄ . Sustainability 2018, 10, 792 5 of 12 determine if the pillar stress in the mine reaches its maximum value; if D/H <3 to 4, the pillar stress ofSustainability the mine 2018 is, not 10, x fully FOR PEER developed REVIEW and K will increase with further increase of D/H. Otherwise, 5 of 11K reaches its maximum value and will no longer increase with the increase of D/H. The parametersparameters inin Equation Equation (2) (2) are are listed listed in Tablein Table3, which 3, which indicates indicates that a thatis dependent is dependent on mining on depthmining but depth is independent but is independent of the mining of the and mining pillar and widths. pillar For widths. each miningFor each depth, mining the depth, mean valuethe mean of a canvalue be of used can toestimate be used to the estimate maximum the pillarmaximum stress pillarK (Table stress4). It is (Table shown 4). in It Figure is shown3 that in theFigure average 3 that coefficientthe averagea coefficientincreases as the increases mining as depth the mining increases. depth increases.
Table 3. 3.Original Original parameters parameters of Equationof Equation (2) when (2) when the mining the mining width towidth pillar to width pillar ratio, widthrm =ratio, 1. R indicates = 1. R adjusted indicates R-square adjusted values. R-square values.