Geotech Geol Eng DOI 10.1007/s10706-015-9851-2

ORIGINAL PAPER

The Height of Water-Conducting Fractured Zones in Longwall Mining of Shallow Coal Seams Xuesheng Liu • Yunliang Tan • Jianguo Ning Chenlin Tian • Jun Wang

•

Received: 7 October 2014 / Accepted: 24 January 2015 Ó Springer International Publishing Switzerland 2015

Abstract The height of the water-conducting fractured zone is of significant importance for both the prevention of water inrush induced disasters and the design of longwall mines. To ensure the mining safety of shallow coal seams under water-rich aquifers, an approach for forecasting the height of the waterconducting fractured zone is proposed by analyzing the damage zone distribution in FLAC3D model. First, 30 numerical models with different face lengths and mining thicknesses of no. 1101 coal face at Luxin Mine, China, a typical shallow coal face, were built using FLAC3D. Then, the heights of water-conducting fractured zones and the thicknesses of protective layers with different face lengths and mining thicknesses were obtained from numerical simulations and empirical formulae. Results show that the height of the water-conducting fractured zone increases with both the face length and mining thickness. And there is a critical value of face length, which is about 300 m in no. 11 coal seam at Luxin Mine. To mine no. 11 coal seam, a slicing mining method is suggested and the face length and mining thickness of the first slice are recommended as 200 and 3.0 m, respectively.

X. Liu (&)  Y. Tan  J. Ning  C. Tian  J. Wang State Key Laboratory of Mining Disaster Prevention and Control Co-founded by Shandong Province and the Ministry of Science and Technology, Shandong University of Science and Technology, Qingdao 266590, China e-mail: [email protected]

Keywords Water-conducting fractured zone  Shallow coal seam  Protective layer  FLAC3D

1 Introduction The overlying strata of coal face are damaged to different degrees during longwall mining (Tan et al. 2012, 2013). The area where water can flow through is named water-conducting fractured zone (Du and Wang 2013; Tan et al. 2011). The development law of the water-conducting fractured zone, especially its height, is of significant importance for the safety of underground excavations. Quantitatively forecasting the height of the water-conducting fractured zone is not only helpful to prevent water inrush induced disasters but also instructive to the design of longwall mining (Liu et al. 2011; Zhang et al. 2011, 2014). Researchers have made great contributions on forecasting and monitoring the height of the waterconducting fractured zone in various conditions (Xu et al. 2010; Khanal et al. 2012; Booth and Bertsch 1999; Cherubini 2008). For example, Miao et al. (2011) monitored the height of the water-conducting fractured zone in undermined rock strata of Bulianta Coal Mine and proposed that the thick and strong rock layer termed key stratum plays an important role in controlling the height. Guo et al. (2012) studied the height of the water-conducting fractured zone in

123

Geotech Geol Eng

longwall mining of deep coal mines by using piezometers and compared it with that obtained from empirical formulae used in China. However, because of the differences in buried depth and depositional age, the damage law of overlying strata in shallow coal seams differs a lot from that of others (Huang 2002; Liu et al. 2014; Lokhande et al. 2012). Compared with other conditions, less research is reported in the literature on the development of water-conducting fractured zones above shallow coal seams. Furthermore, because the height of the waterconducting fractured zone is significantly influenced by mining method, geometric condition, roof type and many other factors, the analytic solution can hardly be obtained (Wu et al. 2000; Xu et al. 2012; Adhikary and Guo 2014). The numerical simulation software FLAC3D is increasingly used to estimate the development of water-conducting fractured zones (Zhang 2011). This article aims to study the variations of the height of water-conducting fractured zones. It begins with the establishment of numerical models of no. 1101 coal face at Luxin Mine, China, using FLAC3D. We then calculate the damage zone using numerical models and obtain the heights of water-conducting fractured zones for different face lengths and mining thicknesses. At last, the mining method, face length and mining thickness of no. 11 coal seam are determined.

2 Numerical Modelling 2.1 Geological and Mining Conditions Numerical models were built based on no. 11 coal seam at Luxin Mine, China. No. 1101 coal face is the first face of no. 11 coal seam in Luxin Mine. The depth of no. 1101 coal face ranges from 168 to 256 m, with an average of 218 m. The thickness of the coal seam ranges from 0.6 to 20.9 m, with an average of 10.76 m. The dip angle ranges from 3° to 5°, with an average of 4.1°. The coal seam is simple in geological formation and rich in joint development. The main roof is sandy mudstone with thickness ranging from 29.8 to 48.2 m, and the immediate roof is also sandy mudstone with thickness ranging from 18.9 to 31.7 m. Characteristics of roof and floor are shown in Table 1. There is a water-rich aquifer in the strata overlying the

123

coal seam, and its distance from the coal seam ranges from 70 to 93.6 m. 2.2 Model Building 2.2.1 Model Size and Grid Based on ‘‘Pillar design and mining regulations under buildings, water, rails and major roadways’’ in China, the height of the water-conducting fractured zone can be obtained from the following empirical formula (National Bureau of Coal Industry of China 2000). P 100 M P Hli ¼  5:6 ð1Þ 1:6 M þ 3:6 where Hli is the height of the water-conducting fractured zone, m; M is the mining thickness, m. The calculated heights of the water-conducting fractured zone of no. 1101 coal face range from 46.1 to 57.3 m. According to the mining practices, the height of the water-conducting fractured zone does not increase with excavating any more when the working face advances more than 300 m. To ensure that there are no significant boundary effects, two groups of numerical models are built, group A and group B. The size of the numerical model in group A (length 9 width 9 height) is set to 2l (l = face length) 9 300 9 200 m and the vertical distance below the seam is set to 30 m. The depths of no. 11 coal seam in both groups are set to 218 m. The size of the numerical model in group B is set to 3l 9 300 9 220 m and the vertical distance below the seam is set to 50 m. Both the haulage way and tail roadway are excavated in the coal seam along the floor of the coal face, with 4 m in width and 3 m in height. The brick shape is used to create the FLAC3D grid. To improve the accuracy, the elements of coal, immediate roof, main roof and floor should be small, the side length of which is less than 2.0 m. The side length of other elements ranges from 3.0 to 5.0 m. The model grid of group A is shown in Fig. 1. 2.2.2 Boundary Conditions Displacement boundary condition is applied on the boundaries of the model (Itasca Consulting Group, Inc., 1997). X-velocity is set to be zero at both left and right boundaries. Y-velocity is set to zero at both front

Geotech Geol Eng Table 1 Characteristics of roof and floor of no. 1101 coal face Appellation

Lithology

Thickness (m)

Description

Main roof

Sandy mudstone

29.8–48.2

Poorly cemented, with some hard interlayers

Immediate roof

Sandy mudstone

18.9–31.7

Poorly cemented, soft, with some coal granules

Floor

Sandy mudstone

29.37–38.92

Argillaceous structure and easy sliming

yield of zones complies with the Mohr–Coulomb criterion, as shown in Eq. (4). The tensile yield is estimated by tensile strength, as shown in Eq. (5). sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 þ sin u 1 þ sin u fs ¼ r1  r3  2c [0 ð4Þ 1  sin u 1  sin u r1 ¼ ½r

ð5Þ

where r1 is the maximum principal stress, MPa; r3 is the minimum principal stress, MPa; u is the internal friction angle, °; c is the cohesion, MPa; [r] is the tensile strength, MPa. 2.2.5 Material Properties

Fig. 1 Grid for no. 1101 coal face model in group A

and back boundaries. X-velocity, y-velocity and z-velocity are all set to be zero at bottom boundary. The top boundary is free.

There are ten rock strata in the model. The physical and mechanical parameters of coal and rocks are listed in Table 2.

2.2.3 Initial Conditions

2.2.6 Calculating Procedures

To reach the equilibrium state rapidly, the vertical and horizontal stress are initialized to rv and rh, respectively, where rv is obtained from Eq. (2) and rh is obtained from Eq. (3) (Itasca Consulting Group, Inc., 1997).

The calculating procedures are as follows:

rv ¼ ch rh ¼

l rv 1l

1.

ð2Þ ð3Þ

where rv is the vertical stress based on gravity loading of overlying strata, Pa; rh is the horizontal stress due to self-loading following Poisson’s effect, Pa; c is bulk density of rock masses, N/m3; h is the depth of coal seam, m; l is the Poisson’s ratio.

2.

2.2.4 Constitutive Model

3.

Coal and rock in model are represented using Mohr– Coulomb plasticity model. The roadway and goaf excavated are represented using null model. The shear

Calculating the model to its initial stress equilibrium state without any excavation. Results show that the vertical stress, the maximum horizontal stress and the minimum horizontal stress are 4.58, 1.93 and 1.12 MPa, respectively, as shown in Fig. 2. Comparing with the measured results (4.6, 1.9 and 1.1 MPa) by the hydrofracturing method, we can reach the conclusion that the numerical model can describe the real stress state of no. 11 coal seam in Luxin Mine. For each group, excavating no. 1101 coal face with different face lengths and mining thicknesses as shown in Table 3 and iterating to the stress equilibrium state. Comparing vertical stress, horizontal stress, vertical displacement and damage zone distribution of model in group A with those of corresponding model in group B to find out whether there are significant boundary effects, as shown in Fig. 3. It

123

Geotech Geol Eng Table 2 Physical and mechanical parameters of coal and rocks in model Lithology

Thickness (m)

Sandy mudstone

30

No. 11 coal seam

10.76

Sandy mudstone Sandy mudstone

43.67 7.92

No. 9 coal seam Siltstone No. 8 coal seam Sandstone No. 6 coal seam

Bulk modulus (GPa)

Shear modulus (GPa)

Cohesion (MPa)

Friction angle (°)

3.12

2.54

2.5

35.0

2.38

1.16

1.1

22.9

3.57 6.67

2.46 2.22

2.5 2.5

35.0 35.0

6.40

2.38

1.16

1.1

22.9

21.83

2.86

2.61

2.5

35.0

3.52

2.38

1.16

1.1

22.9

22.50

10.00

2.14

2.5

35.0 22.9

6.35

2.38

1.16

1.1

Fine sandstone

11.42

10.00

2.14

2.5

35.0

Conglomerate

35.63

2.14

1.67

2.2

34.0

3 Results and Discussions 3.1 Effects of Face Length on the Height of the Water-Conducting Fractured Zone

Fig. 2 The resulting horizontal stress distribution without excavation: a vertical stress, b the maximum horizontal stress

4.

can be easily concluded from the comparisons that any boundary effects are indeed negligible. By analyzing the damage zone distribution in models of group A, the area consisting of damage zones that directly connect the goaf can be seen as water-conducting fractured zone, as shown in Fig. 4. Then the heights of water-conducting fractured zones of different models can be obtained.

123

The results of Models 1–1 to 1–6 are shown in Fig. 5. Results of other models are similar to Fig. 5 and not listed. By analyzing the damage zone distribution in models, the heights of the water-conducting fractured zone were obtained. The variations of the height of the water-conducting fractured zone with different face lengths are shown in Fig. 6. As shown in Fig. 6, for the same mining thickness, the heights of water-conducting fractured zones increase with the face length. Many curves of the heights hit inflection points at the face length of 250 m. The heights increase slowly before the inflection points and rapidly thereafter. Then the heights reach peak values at the face length of 300 m and remain nearly constant. The largest height of the water-conducting fractured zone ranges from 140 to 180 m, about 3 times larger than that at the face length of 100 m. It can be concluded that in shallow coal seams the height of the water-conducting fractured zone generally increases with the face length. And there is a critical value of face length where the height of the water-conducting fractured zone remains nearly constant after the face length increases to the critical value. The critical value of face length of no. 11 coal seam is about 300 m.

Geotech Geol Eng Table 3 Face lengths and mining thicknesses of different models Appellation

Model 1–1

Model 1–2

Model 1–3

Model 1–4

Model 1–5

Model 1–6

Face length (m)

100

150

200

250

300

350

Mining thickness (m)

3.0

3.0

3.0

3.0

3.0

3.0

Appellation

Model 2–1

Model 2–2

Model 2–3

Model 2–4

Model 2–5

Model 2–6

Face length (m)

100

150

200

250

300

350

Mining thickness (m)

4.5

4.5

4.5

4.5

4.5

4.5

Appellation

Model 3–1

Model 3–2

Model 3–3

Model 3–4

Model 3–5

Model 3–6

Face length (m)

100

150

200

250

300

350

Mining thickness (m)

6.0

6.0

6.0

6.0

6.0

6.0

Appellation

Model 4–1

Model 4–2

Model 4–3

Model 4–4

Model 4–5

Model 4–6

Face length (m)

100

150

200

250

300

350

Mining thickness (m)

7.5

7.5

7.5

7.5

7.5

7.5

Appellation

Model 5–1

Model 5–2

Model 5–3

Model 5–4

Model 5–5

Model 5–6

Face length (m)

100

150

200

250

300

350

Mining thickness (m)

9.0

9.0

9.0

9.0

9.0

9.0

Fig. 4 The water-conducting fractured zone obtained from the damage zone distribution in model (paralleling the face, at y = 150 m). Blue is no-failure zone; Red is failure in tension; Light magenta is failure in shear; Green is failure in the past

Fig. 3 Some numerical results of models (face length 150 m, mining thickness 3.0 m) in group A and group B: a vertical stress, b damage zone distribution, c horizontal stress, d vertical displacement (paralleling the face, at y = 150 m)

Fig. 5 The damage zone distribution in models: a Model 1–1, b Model 1–2, c Model 1–3, d Model 1–4, e Model 1–5, f Model 1–6 (paralleling the face, at y = 150 m)

123

Geotech Geol Eng

layer between the water-rich aquifer and the waterconducting fractured zone. According to ‘‘Pillar design and mining regulations under buildings, water, rails and major roadways’’ in China, the minimum thickness of the protective layer needed can be obtained from Eq. (6). The thicknesses of the protective layer in Models 1–1 to 5–6 can be obtained from Eq. (7). The minimum thicknesses needed and the thicknesses of protective layers in Models 1–1 to 5–6 are shown in Table 4. Hmin ¼ 3M Fig. 6 Heights of water-conducting fractured zones for different face lengths. The different curves are for different mining thicknesses, 3.0, 4.5, 6.0, 7.5 and 9.0 m, respectively

3.2 Effects of Mining Thickness on the Height of the Water-Conducting Fractured Zone The variations of the height of the water-conducting fractured zone with mining thickness are shown in Fig. 7. The height of the water-conducting fractured zone generally increases steadily with the mining thickness. 3.3 Determination of the Face Length and Mining Thickness To ensure the safety of the shallow coal mine, the mining induced water-conducting fractured zone cannot reach the water-rich aquifer in the overburden. And there must be a certain thickness of the protective

ð6Þ

where Hmin is the minimum thickness of the protective layer needed, m; M is the mining thickness, m.  H  Hli H [ Hli HN ¼ ð7Þ 0 H  Hli where HN is the thickness of the protective layer in the numerical model, m; H is the minimum distance between coal seam and the water-rich aquifer, m; Hli is the height of the water-conducting fractured zone obtained from numerical calculations, m. As shown in Table 4, only Models 1–1, 1–2, 1–3, 2–1 and 2–2 can ensure enough thickness of the protective layer during mining. Therefore, no. 11 coal seam had to adopt slicing mining method. Comparing Models 1–1, 1–2, 1–3, 2–1 and 2–2, it was easy to exclude Models 1–1, 1–2 and 2–1 for their relatively inefficiency in mining. Considering that the abutment pressure of coal face with larger mining thickness is much higher, we selected Model 1–3 as the mining model in the first slice of no. 11 coal seam. Of course, to mine the remaining resources of no. 11 coal seam, some measures such as borehole dewatering and draining should be adopted to dewater the water-rich aquifers in the overburden. Or we can use some special mining methods such as backfill mining to greatly reduce the height of the water-conducting fractured zone. 4 Conclusions 1.

Fig. 7 Heights of water-conducting fractured zones with different mining thicknesses. The different curves are for different face lengths, 100, 150, 200, 250, 300 and 350 m, respectively

123

Numerical simulations show that the height of the water-conducting fractured zone increases with both the face length and mining thickness of coal face in a shallow coal seam. And there is a critical value of face length, which is about 300 m in no. 11 coal seam at Luxin Mine, China.

Geotech Geol Eng Table 4 The minimum thicknesses needed and the thicknesses of protective layers in numerical models

Models Hmin (m)

Model 1–3

Model 1–4

Model 1–5

Model 1–6

9

9

9

9

9

9

33.6

26.2

11.3

6.1

0

0

Models

Model 2–1

Model 2–2

Model 2–3

Model 2–4

Model 2–5

Model 2–6

Hmin (m)

13.56

13.56

13.56

13.56

13.56

13.56

HN (m)

28.4

23.1

Models

Model 3–1

Model 3–2

Model 3–3

Model 3–4

Model 3–5

Model 3–6

Hmin (m)

18

18

18

18

18

18

HN (m)

14.7

0

0

0

0

Models

Model 4–1

Model 4–2

Model 4–3

Model 4–4

Model 4–5

Model 4–6

Hmin (m)

22.5

22.5

22.5

22.5

22.5

22.5

2.5

0

0

0

0

0

8.6

1.3

2.8

0

0

Models

Model 5–1

Model 5–2

Model 5–3

Model 5–4

Model 5–5

Model 5–6

Hmin (m)

27

27

27

27

27

27

0

0

0

0

0

0

HN (m)

3.

Model 1–2

HN (m)

HN (m)

2.

Model 1–1

The slicing mining method is suggested to be adopted in no. 11 coal seam, and the face length and mining thickness of the first slice are recommended as 200 and 3.0 m, respectively. The thickness of the protective layer is expected to be 11.3 m. To mine of next slice of no. 11 coal seam, measures such as borehole pumping should be adopted to dewater the water-rich aquifer in overburden. Or some special mining methods such as backfill mining should be used.

Acknowledgments This study was supported by National Natural Science Foundation of China (Nos. 51274133, 51174129), Shandong Province Natural Science Fund (No. ZR2010EEZ002), Doctoral Scientific Fund Project of the Ministry of Education of China (No. 20123718110013) and Open Fund of State Key Laboratory of Mining Disaster Prevention and Control Co-founded by Shandong Province and the Ministry of Science and Technology (No. MDPC2013 KF12).

References Adhikary DP, Guo H (2014) Measurement of longwall mining induced strata permeability. Geotech Geol Eng 32(3): 617–626 Booth CJ, Bertsch L (1999) Groundwater geochemistry in shallow aquifers above longwall mines in Illinois, USA. Hydrogeol J 7:561–575

Cherubini C (2008) A modeling approach for the study of contamination in a fractured aquifer. Geotech Geol Eng 26(5):519–533 Du JP, Wang LQ (2013) Special mining method. China University of Mining and Technology Press, Xuzhou (in Chinese) Guo WB, Zou YF, Hou QL (2012) Fractured zone height of longwall mining and its effects on the overburden aquifers. Int J Min Sci Technol 22:603–606 Huang QX (2002) Ground pressure behavior and definition of shallow seams. Chin J Rock Mech Eng 21(8):1174–1177 (in Chinese) Itasca Consulting Group, Inc (1997) FLAC 3D (version 2.0) User’s manual. Minneapolis Khanal M, Adhikary D, Balusu R (2012) Numerical analysis and geotechnical assessment of mine scale model. Int J Min Sci Technol 22(5):693–698 Liu YD, Zhang DS, Fan GW et al (2011) Applicable conditions for a classification system of aquifer-protective mining in shallow coal seams. Min Sci Technol (China) 21:381–387 Liu XS, Ning JG, Tan YL (2014) Study on roof breaking model of near shallow seam. Chin J Rock Mech Eng 31(2): 214–219 (in Chinese) Lokhande RD, Murthy VMSR, Singh KB (2012) Pot-hole subsidence in underground coal mining: some Indian experiences. Geotech Geol Eng 31(2):793–799 Miao XX, Cui XM, Wang JA et al (2011) The height of fractured water-conducting zone in undermined rock strata. Eng Geol 120:32–39 National Bureau of Coal Industry of China (2000) Pillar design and mining regulations under buildings, water, rails and major roadways. China Coal Industry Publishing House, Beijing (in Chinese) Tan YL, Ning JG, Gu ST et al (2011) Ground pressure and strata control (Revised version). China Coal Industry Publishing House, Beijing (in Chinese)

123

Geotech Geol Eng Tan YL, Ning JG, li HT (2012) In situ explorations on zonal disintegration of roof strata in deep coalmines. Int J Rock Mech Min Sci 49:113–124 Tan YL, Yu FH, Chen L (2013) A new approach for predicting bedding separation of roof strata in underground coalmines. Int J Rock Mech Min Sci 61:183–188 Wu Q, Huang XL, Dong DL et al (2000) ‘‘Three maps-two predictions’’ method to evaluate water bursting conditions on roof coal. J China Coal Soc 25(1):60–65 Xu ZM, Sun YJ, Dong QH et al (2010) Predicting the height of water-flow fractured zone during coal mining under the Xiaolangdi Reservoir. Min Sci Technol 20:434–438 Xu JL, Zhu WB, Wang XZ (2012) New method to predict the height of fractured water-conducting zone by location of key strata. J China Coal Soc 37(5):762–769

123

Zhang SJ (2011) Reasonable workface Mining thickness under shallow cover and thin strata geological condition in Sanbula Coal Mine. Proc Eng 26:411–417 Zhang DS, Fan GW, Ma LQ et al (2011) Aquifer protection during longwall mining of shallow coal seams: a case study in the Shendong Coalfield of China. Int J Coal Geol 86:190–196 Zhang JG, Miao XX, Huang YL et al (2014) Fracture mechanics model of fully mechanized top coal caving of shallow coal seams and its application. Int J Min Sci Technol 24:349–352