New Method to Design Coal Pillar for Lateral Roof Roadway Based on Mining-Induced Stress : A Case Study

State Key Laboratory of Mining Disaster Prevention and Control Co-Founded by Shandong Provence and the Ministry of Science and Technology, Shandong University of Science and Technology, Qingdao 266590, China College of Mining and Safety Engineering, Shandong University of Science and Technology, Qingdao 266590, China State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China


Introduction
Multiseam longwall coal mining is widely applied both here in China and abroad [1,2], for example, Germany has a history of multiseam mining over the past 100 years.It is very common to mine multiple coal seams simultaneously due to the close distance coal seams.e roadway in the upper panel is located in the overlying strata of the lower coal seam, that is, lateral roof roadway.
e lateral roof roadway influenced by the lower coal seam mining is prone to play an important role in transportation, ventilation, etc., during close distance coal seams repeat mining [3][4][5].However, the multiseam longwall coal mining will lead to strata movement, stress redistribution, roof caving, and overburden fracturing, resulting in lateral roof roadway convergence.Nevertheless, the protection of lateral roof roadway requires the determination of the delimitation of the risk of collapse.
e traditional strata movement contour is foundation and precondition of the prediction of strata movement and deformation as well as the protection of roadway outside the mining area [6][7][8][9].An appropriate coal pillar is designed based on the evolution of strata movement contour [10,11].Traditionally, the strata movement contour is considered to be a straight line based on the mining-induced movement angle of the overlying strata [6,12].According to the definition of the mining-induced movement angle, only one angle of the rock strata is used to represent the characteristics of overlying strata movement, while the internal deformation mechanism as well as the stress distribution of the overlying strata is ignored.erefore, the lateral roof roadway would mostly be damaged as the coal pillar determined by the mining-induced movement angle of the rock strata was not appropriate.For example, based on a 61 °mining-induced movement angle [11], coal pillar in 90 m wide was reserved to protect the lateral roof roadway in Pingdingshan No. 1 coal mine, which makes it outside of the overlying strata movement area of longwall panel 31010.However, during the mining process, serious heave was still occurred in the roadway. is paper analyzed the deformation mechanism of the roadway based on the mining-induced stress and proposed a method to design coal pillar for lateral roof roadway.
e contents are as follows: (1) Deformation mechanism of the lateral roof roadway was studied by field measurements and numerical simulations.
(2) Distribution of mining-induced stress in the overlying strata was studied by theoretical analysis and numerical and physical simulation.(3) e method was proposed to design coal pillar for the lateral roof roadway based on the mininginduced stress contour.(4) According to the method, when the width of coal pillar was 150 m, the maximum deformation of the roadway was 270 mm, which ensured the normal usage of this roadway.In order to avoid the influence of mining, coal pillar for lateral roof roadway was designed via the displacement angle of rock strata.According to the displacement angle of rock strata (β � 61 °) [11] in the coal mine, the location of the roof roadway was originally determined as shown in Figure 1.

Mechanism of Lateral Roof Roadway Failure
e average vertical distance between the roadway and the #10 coal seam is 85 m, and the width of the coal pillar on the side of panel 31010 is 90 m.Before Oct 2011, the width of the coal pillar between panel 21210 and the roadway was gradually decreased.From Oct 2011, the width of the coal pillar between the panel and the roadway was maintained at 165 m.Although coal pillar was reserved for the lateral roof roadway, severe deformations still occurred in partial positions of roadway during the mining period of panel 31010.

Field Observation.
Measuring points were installed to monitor deformation in the roadway at intervals of 50 m on Jan 21, 2012.After the points were arranged, horizontal convergence (rib-to-rib convergence) and vertical convergence (roof-to-floor convergence) of the lateral roof roadway were also monitored on Feb As shown in Figure 2, after the mining of panel 21210, there is no horizontal and vertical convergence before the panel 31010 excavation, indicating that the mining of panel 21210 has almost no effect on the lateral roof roadway.During the mining of panel 31010 until July 2012, severe deformation occurred, so it is caused by panel 31010 excavation.e maximum horizontal convergence and vertical convergence were 488 mm and 587 mm.

Numerical Simulation.
UDEC 2D was used to study the deformation mechanism of the roadway.Section A-A in Figure 1 was selected to build the model with the dimensions of 1500 m wide and 300 m height shown in Figure 3. e hierarchies of the overlying strata were correspondingly simplified according to the field drilling columnar.Field measurement showed that the vertical stress was 14.67 MPa.
e maximum principal stress was 17.44 MPa, and the minimum principal stress was 15.64 MPa. is in situ stress state was incorporated into the model.During the simulating process, the elastic-plastic model [13,14] was applied to the model, which can represent shear failure in rock and is widely used in underground engineering.e mechanical properties of different rock layers are shown in Table 1. Figure 3 shows that good agreement was achieved between the numerical results and the laboratory test for stress-strain curve of the silty mudstone sample.ree plans were designed in this experiment: mining panel 31010; mining panel 21210; and mining both panels.Ten measuring lines separately located in the positions of Y � 26 m, 66 m, 108 m, 166 m, 216 m, 266 m, 316 m, 366 m, 416 m, and 466 m were arranged in the coal seam and overlying strata to monitor the mining-induced stress.
In order to determine the influence distance of mininginduced stress, the curve formed by connecting the points, where the mining-induced stress is equal to 1.05 times the in situ stress [15,16], is defined as mining-induced stress contour.
erefore, the stress inside the mining-induced stress contour is redistributed, while the outside strata are still located in the in situ stress zone.Mining-induced stress on the overlying strata is distributed as shown in Figure 4.
(1) After mining of panel 21210, there is no deformation because the roadway is located outside the mininginduced stress contour.(2) After excavation of panel 31010, the width of the coal pillar in the roadway of panel 31010 is significantly smaller than that of panel 21210.Considering that the roadway is located inside the mining-induced stress contour, there are slight deformations of roadway.

Advances in Civil Engineering
(3) After panel 21210 and panel 31010 are both mined, the roadway is located in the cumulated in uence area of the two panels.e roadway shows severe deformations due to the signi cant increase in vertical stress compared to the case of mining panel 31010.e results of numerical simulation reveal the deformation mechanism of the roadway.

Evolution of Mining-Induced Stress Contour
e aforementioned results indicate that the mininginduced stress contour presents an external expansion type from the coal seam to the overlying strata in the numerical simulation.In the following section, the distribution of the mining-induced stress contour within the overlying   Advances in Civil Engineering strata is developed by means of theoretical analysis and physical simulation experiments.

eoretical Analysis.
e excavation of the panel causes movement of the overlying strata and results in a caving zone, fractured zone, and continuous deformation zone in the vertical direction [15,17], which leads to stress redistribution.In fact, there is an inseparable relationship between overlying strata movement and mining-induced stress.On the basis of the strata subsidence law at di erent distances away from the coal seam, the distribution of mining-induced stress in overlying strata is analyzed below.Advances in Civil Engineering 3.1.1.Mining-Induced Subsidence. Figure 5 shows the mining-induced subsidence of the overlying strata in different heights of the overlying strata.Similar to the law of surface subsidence [10], Equation (1) can be used to express the distribution of subsidence at di erent heights of overlying strata: where x represents the coordinate in horizontal direction (m); y represents the coordinate of vertical direction (m); L represents the half width of the model (m); ω(x, y) represents the subsidence at di erent heights of overlying strata (m); ω max (y) represents the maximum subsidence at the same height (m); and C(y) is an undetermined coe cient.
In Equation (1), there are three undetermined coe cients, L, ω max (y), and C(y), among which L can be determined according to the width of the panel, while the other two undetermined coe cients cannot be obtained through theoretical calculation.erefore, the subsidence in di erent heights of the overlying strata is tted using Equation (1), and the results are as shown in Figures 5 and 6 and can be expressed as follows:  Advances in Civil Engineering Equation ( 2) is substituted into Equation ( 1) to obtain the subsidence of the overlying strata at different heights, expressed as follows: (3)

Relation of Mining-Induced Subsidence to Stress
Distribution.According to Equation (3), the subsidence at any height in the overlying strata can be obtained.e load distribution can be obtained according to the relation between load distribution and deflection on an arbitrary beam [18]: where E represents the elasticity modulus (GPa) and I represents the inertia moment (m 4 ).e schematic diagram of the overlying strata subsidence is shown in Figure 7.By taking the first and second stratum under the surface as an example, the stress on the rock stratum is shown in Figure 7.It can be obtained from the equilibrium of first and second stratum that where c 1 and c 2 represent the volume weight of the first and second strata (kN/m 3 ) and q 0 represents the top boundary load on the stratum (Pa).Because the first stratum is the ground rock stratum, q 0 equals to 0, q 1 (x, y) represents the stress distribution on the floor boundary of the first rock stratum (Pa), q 2 (x, y) represents the stress distribution on the floor boundary of second rock stratum (Pa), h 1 and h 2 represent the thickness of the first and second strata (m), E 1 and E 2 represent the elasticity modulus of the first and second strata (GPa), and I 1 and I 2 represent the inertia moment of the first and second strata (m 4 ).e equilibriums of i . . .n overlying strata are shown in Equation (6): where c i andc n represent the volume weight of the i-th and n-th strata (kN/m 3 ); q i−1 (x, y) represents the stress distribution on the floor boundary of the (i − 1)th rock stratum (Pa); q n−1 (x, y) represents the stress distribution on the floor boundary of (n − 1)th rock stratum (Pa); h i and h n represent the thickness of the i-th and n-th strata (m); E i and E n represent the elasticity modulus of the i-th and n-th strata (GPa); and I i and I n represent the inertia moment of the i-th and n-th strata (m 4 ).By adding each left item and each right item of Equations ( 5) and ( 6), the stress distribution of an arbitrarily overlying stratum can be obtained by the following: According to the definition of mining-induced stress contour, that is, q n (x, y) � 1.05 •  n i�1 c i h i , the previous equation can be simplified and taken as q n (x, y) �  n i�1 c i h i .erefore, the distribution characteristics of the mininginduced stress contour at any height of overlying strata can be obtained according to Equation ( 8): According to Equation ( 8) and numerical simulation results, the distribution characteristic of the mining-induced stress contour is depicted in Figure 8. e results of the numerical simulation correlate well with theoretical results.

Physical Simulation.
e physical simulation model shown in Figure 9 was constructed in order to study the mining-induced stress distribution in the overlying strata.
e model composed of seven rock layers has dimensions of 2500 mm × 1300 mm × 200 mm (length × height × width).A scaled physical model has to be designed in accordance with the laws of the similarity theory [19,20].As per these laws, the parameters should satisfy the fundamental condition of the similarity theory.
where C σ is the similarity constant for strength, C ρ is the similarity constant for density, and C L is the similarity constant for geometry; and the subscript P denotes the prototype based on field conditions, while m represents the physical model.e geometric similarity ratio was 1 : 100, weight similarity ratio was 1 : 1.25, rigidity similarity ratio was 1 : 1.67, 8 Advances in Civil Engineering and stress similarity ratio was 1 : 1.67.River sand was used as aggregate and gypsum, and calcium carbonate was used as to allow for the cohesion of each layer.e loose bed was composed of the river sand and dry sawdust.e mixture material ratio and the physical and mechanical parameters of the strata are presented in Table 2. Four vertical survey lines including four points in each line were installed in the model to monitor the mininginduced stress.BW-type pressure cells embedded in each rock layer in advance were adopted for stress monitoring during the entire mining process (Figure 9).In the experiment, the panel was advanced for an accumulative total length of 1,700 mm.At both boundaries of the model, 400 mm coal pillar remained to eliminate the boundary in uence.e panel was advanced by 50 mm at an interval of 15 min.e distance between the mining-induced stress contour and the mining boundary is shown in Table 3.
ω n (x,y) First stratum Second stratum  Advances in Civil Engineering stress enlarges gradually.For example, the stress in measuring point 2-1 starts to increase as the distance between the measuring point and the panel is 30 m.While the distance is 60 m for measuring point 2-4.When the panel is continuously advanced forward, there is also the same stress development law on the other measuring points.erefore, when the panel is mined, as the vertical distance between the measuring point and the coal seam increases, the in uence ranges of mining-induced stress are 28 m, 43 m, 54 m, and 60 m on average.e in uence distance of mining-induced stress gradually enlarges.e mining-induced stress contour presents a concave-upward type from the coal seam to the overlying strata.ese conclusions are consistent with the numerical simulation results.

Method to Design Coal Pillar for Lateral Roof Roadway
Mining-induced stress contour is the curve to judge whether the overlying strata is a ected by coal mining.Inside of it, stress inside the contour is redistributed, whereas outside of it, the overlying strata are still within the in situ stress area.As shown in Figure 10, the mining-induced stress contour extends in a concave way from the coal seam to surface.e higher the coal seam, the larger the horizontal distance between stress contour and longwall mining edge.erefore, the width of coal pillar should be designed according to stress contour and the position of protectors to ensure that the lateral roof roadway is not affected by excavation.e detailed steps of this method are as follows: determining mining-induced stress contour based on actual geological conditions; calculating the intersection of the level of lateral roof roadway and mining-induced stress contour and making its projection point on the coal seam; the distance between projection point and longwall mining edge is the width of the coal pillar.
A schematic diagram of mining-induced stress distribution in the overlying strata is shown in Figure 10.When the vertical distance between the roadway and the coal seam is L, the width of coal pillar designed according to the traditional displacement angle of rock strata is a.It is obvious that the roadway located at position 1 would be damaged because it is still located within the mininginduced stress concentration zone.us, the coal pillar according to the traditional method cannot ensure that no deformation would occur to the roadway.According to the method outlined in this paper, the width of the coal pillar for roadway would be b, and the roadway would be located in position 2. e roadway is outside the mining-induced stress contour and is not subject to the influence by excavation, which meets the design requirements.Compared with the traditional method, the modified method takes into account the distribution law of mining-induced stress within the overlying strata, so it can better meet engineering demands.

Case 1: Longwall Panel 31010 of Pingdingshan Coal Mine.
e lateral roof roadway was still located in the overlapping influence area of panels 21210 and 31010, and the roadway was damaged.e traditional design method ignores both the strata deformation and the mining-induced stress distribution of the overlying strata.Now, the coal pillar was redesigned based on the new methodology.According to the aforementioned results, the average distance between the mining-induced stress contour and the edge of panel 21210 and 31010 was shown in Table 4.
As shown in Figure 1, the average vertical distance between the lateral roof roadway and the #10 coal seam is 85 m. e coal pillar is 143 m and 158 m in width for panel 21210 and 31010 based on the new method according to Table 4. e lateral roof roadway is not influenced by mining of the panel 21210 because the actual size of coal pillar is wider than the design value.Because the width of the coal pillar in the roadway of panel 31010 is significantly smaller than that of panel 21210, the lateral roof roadway had undergone serious deformation.According to the method, the coal pillar's width is 158 m for the lateral roof roadway of longwall panel 31010.

Mining Conditions. By taking Pingdingshan No. 1 Coal
Mine as an example for designing coal pillar for lateral roof roadway, the deformation law of the roadway was also studied when mining panel 31100.e average burial depth of panel 31100 is 750 m. e thickness of the coal seam is 2.4 m, and the average dip angle is 6.5 °.With a section of 4.2 m wide and 3 m high, the roadway is supported by the combination of bolt and anchor.e bolt is 20 mm in diameter and 2.2 m in length with raw and line space of 700 mm × 700 mm. e anchor is 20 mm in diameter and 7 m in length.
e lateral roof roadway is adjacent to panel 31100 passenger roadway.e vertical distance between the lateral roof roadway and the coal seam is 100 m.According to the method, the coal pillar's width should not be less than 160 m, and the actual width of coal pillar for roadway is 150∼175 m.

Field Measurement and eoretical Analysis.
In order to observe the deformation characteristics of the lateral roof roadway under the mining impact of panel 31100, a total of 12 measurement points at intervals of 100 m were arranged in the roadway.During the later observation, the data of displacement of roof-to-oor and two sides are measured by the cross-distributing point method and the initial data are obtained on Oct. 10, 2013.Field measurement has been carried out for 14 times.e results of displacement of roof-to-oor and two sides are shown in Figure 11.After the mining of panel 31100, slight deformations were observed.e maximum deformation of roadway was 270 mm, which could satisfy normal use without any repair.erefore, the coal pillar designed according to the method can de nitely meet engineering practice.

Conclusion
(1) Evolution of mining-induced stress was simulated by numerical and physical simulations.e stress of the overlying strata was redistributed during excavation.A contour line of 1.05 times the in situ stress is used to de ne the mining-induced stress contour.Mining-induced stress contour was used to determine whether the stress of the overlying strata was redistributed or not.Stress inside the contour is redistributed while that outside the contour the    Advances in Civil Engineering overlying strata is still within the in situ stress area.Mining-induced stress contour presents a concaveupward type from coal seam to the overlying strata that cannot be merged into a straight one.(2) A method was proposed to design coal pillar for the lateral roof roadway based on the evolution of mining-induced contour.Compared with the conventional method, the method outlined in this paper made up for the deficiency which ignored the variations of the internal strata deformation law as well as stress distribution characteristics.(3) When the vertical distance between the lateral roof roadway and the coal seam was 100 m in longwall panel 31100 of Pingdingshan No. 1 coal mine, the width of the coal pillar was 160 m based on the aforementioned method, and the maximum deformation of roadway was 270 mm.e new method can definitely meet engineering demands.

2. 1 .
Mining Conditions.As shown in Figure 1, Pingdingshan No. 1 Coal Mine is located in Pingdingshan, Henan Province, China.e average depth of the roadway studied in this paper is 470 m, with a section size (width × height) of 4.2 m × 3.0 m. e lateral roof roadway is located in the overlying strata of the longwall panels 21210 and 31010 in #10 coal seam.e width of panel 21210 and 31010 is 125 m and 180 m on average.e thickness of #10 coal seam ranges from 2.1 m to 3.8 m, and the inclination of the coal seam ranges from 4 °to 12 °.e drilling columnar and lithological parameters of the roof are shown in Figure 1.

Figure 1 :
Figure 1: Layout sketch of longwall panel and lateral roof roadway: (a) location of longwall panel and lateral roof roadway; (b) pro le of longwall panel and lateral roof roadway; (c) borehole columnar section of mining area.※Location of the lateral roof roadway in our study.

Figure 3 :
Figure 3: Numerical simulation: (a) numerical simulation model; (b) stress-strain curves in laboratory test and numerical simulation of the silty mudstone UCS test.

Figure 8 :
Figure 8: Comparison between numerical simulation and theoretical results.

Figure 9 :
Figure 9: Physical model and stress monitoring system: (a) plan of physical simulation model (unit: m); (b) front view of physical model; (c) stress monitoring system.

Figure 10 :
Figure 10: Method to design coal pillar for lateral roof roadway.β represents the displacement angle, 1 and 2 represent the roadways in different positions, a and b represent the width of the coal pillar, and L represents the vertical distance between the roadway and coal seam.

Table 2 :
Mixture material ratio and the physical and mechanical parameters of the strata.

Table 3 :
Distance between the mining-induced stress contour and the edge of panel (unit: m).

Table 4 :
Distance between mining-induced stress contour and edge of panel.