Control Technology of Surface Movement Scope with Directional Hydraulic Fracturing Technology in Longwall Mining: A Case Study

: Mining-induced surface subsidence causes a series of environmental hazards and social problems, including farmland destruction, waterlogging and building damage in the subsidence area. To reduce mining damages, an innovative method of controlling the surface movement scope via artiﬁcial weak planes generated by hydraulic fracturing technology was proposed in this paper. Numerical models were built to analyze the inﬂuence of weak planes with di ﬀ erent heights and dips on the overlying strata movement. The numerical simulation results showed that the weak planes structure cut o ﬀ the development of the overlying strata displacement to the surface and a ﬀ ected the surface movement scope. When the weak planes’ dips were bigger than the angle of critical deformation, with the increase of the weak planes’ heights (0–120 m) the advance angle of inﬂuence changed from 53.61 ◦ to 59.15 ◦ , and the advance distance of inﬂuence changed from 173.31 m to 140.27 m which decreased by 30.04 m. In applications at Sihe coal mine in China, directional hydraulic fracturing technology was used in panel 5304 to form artiﬁcial weak planes in overlying strata. The measured surface subsidence and deformation value met the numerical simulation results and the mining-induced surface movement scope reduced. Moreover, no damage occurred to the surface buildings which were predicted to be in the a ﬀ ected area after extraction. This technology provided a new method to protect the surface structures from damages and had great beneﬁts for the sustainable development of coal mines.


Introduction
Surface subsidence is a global problem [1][2][3] and mostly results from human activity such as groundwater and mineral extraction and the construction of underground facilities. This problem is particularly significant within the underground longwall mining [4,5] industry due to the severe disturbance to the overlying strata during the mining process ( Figure 1). Developing effective measures to control mining-induced surface subsidence has been the focus of researchers over the past few decades.
At present, there are two main methods of mining under buildings in China.
(1) Management methods, that is, the relocation of surface buildings.
However, these methods have their shortcomings as follows: (1) The relocation mode has problems such as land acquisition, contradiction between workers and residents, and one-time investment in relocation funds. (2) Strip mining [11,12] causes serious loss of resources, and its recovery efficiency is generally low, which is unfavorable to the development of the enterprise, moreover, it is difficult to ensure the long term stability of coal pillars, and the coal pillars in the gob are prone to failure and cause serious damage to the surface. (3) For filling mining, the demand for filling materials is large, and the initial input and production costs are relatively high; the effect of reducing subsidence of grouting bed separation is poor [13], because the location of the separation layer is difficult to determine. (4) The effect of coordinated mining [14] dealing with reducing subsidence is not good, and the resulting surface deformation position will cause seasonal water accumulation.
With the development of technology, some scholars propose the coupling of mining technology of the above. "Strip mining-gob grouting-filling-strip pillar mining" method [15,16] and "miningbackfilling-keeping" method [17] combined the advantages of strip mining and backfilling mining organically.
Previously, all of the technologies above were aimed at reducing the degree of overburden deformation and subsidence to protect surface buildings and constructions. In the course of its formation and long-term geological activities, many randomly oriented strength weakening planes had formed in the rock mass, which were also called rock mass fracture surfaces or rock mass discontinuities. With the development of mining subsidence monitoring and the in-depth analysis of overlying strata movement data, Xie et al. [18] found that the development degree of mining subsidence would be completely different under the condition of the same mining factors and different geological factors. The surface environmental disasters in the coal mining area related to Developing effective measures to control mining-induced surface subsidence has been the focus of researchers over the past few decades.
At present, there are two main methods of mining under buildings in China.
(1) Management methods, that is, the relocation of surface buildings.
However, these methods have their shortcomings as follows: (1) The relocation mode has problems such as land acquisition, contradiction between workers and residents, and one-time investment in relocation funds. (2) Strip mining [11,12] causes serious loss of resources, and its recovery efficiency is generally low, which is unfavorable to the development of the enterprise, moreover, it is difficult to ensure the long term stability of coal pillars, and the coal pillars in the gob are prone to failure and cause serious damage to the surface. (3) For filling mining, the demand for filling materials is large, and the initial input and production costs are relatively high; the effect of reducing subsidence of grouting bed separation is poor [13], because the location of the separation layer is difficult to determine. (4) The effect of coordinated mining [14] dealing with reducing subsidence is not good, and the resulting surface deformation position will cause seasonal water accumulation.
With the development of technology, some scholars propose the coupling of mining technology of the above.
Previously, all of the technologies above were aimed at reducing the degree of overburden deformation and subsidence to protect surface buildings and constructions. In the course of its formation and long-term geological activities, many randomly oriented strength weakening planes had formed in the rock mass, which were also called rock mass fracture surfaces or rock mass discontinuities. With the development of mining subsidence monitoring and the in-depth analysis of overlying strata movement data, Xie et al. [18] found that the development degree of mining subsidence would be completely different under the condition of the same mining factors and different geological factors. The surface environmental disasters in the coal mining area related to underground mining, although derived from the mining of coal seams, were essentially controlled by these rock mass discontinuities. Under different rock mass discontinuities, the extent of ground subsidence and cracking in coal mining areas caused by the same intensity of underground mining activities were significantly different. Furthermore, the surface subsidence law was significantly different from that of surface subsidence under ordinary geological conditions [19].
This paper studied the relationship between the movement characteristic of overlying strata and weak planes by means of theoretical analysis, numerical simulation and field measurement. Moreover, the technology of "mining-induced surface movement scope control by hydraulic fracturing" was put forward for the first time. The field test verified the effect of surface movement scope controlled by overlying strata hydraulic fracturing during longwall mining, providing technical support and acting as reference for similar trials investigating the mitigation of surface subsidence and building damages in coal mining districts, which is of great significance for the sustainable development of aging coal mine enterprises.

Mechanism of Mining-Induced Surface Movement Scope Control by Weak Planes
Xia et al. [20,21] systematically studied the control mechanism of tectonic environment (tectonic medium, tectonic interface, tectonic stress, etc.) on coal mining subsidence from the geological point of view and put forward the theory of "tectonic control of hazards" (TCH). Tectonic weak planes are one of the important factors determining the geological environment bearing capacity of coal mining areas. Weak planes destroy the continuity of the rock strata and weaken the mechanical strength of the rock strata, controlling the deformation, failure mode and stability of the rock strata.
Taking faults as an example [22], during the advancing of the working face, the overlying rock mass above the gob is destroyed by tensile stress, which leads to the movement of rock strata since faults belong to the weak structure planes, and have no resistance to the tensile stress. When the movement of the rock strata reaches faults, the overlying strata on both sides of faults are easy to shear and slip along the fault planes. Within the influence area by underground excavation, most of the extrusion displacement and rebound deformation pointing to the gob are limited to the fault zones, which leads to different characteristics of rock strata movement and deformation on both sides of the fault planes. Guo [5] made an analysis with a profile perpendicular to the strike of the fault and studied the influence of faults on the law of strata and surface movement. When the faults dip is greater than 2 • and the drop of the fault is greater than 10 m, the law of surface subsidence will get significantly different from that of surface subsidence under ordinary geological conditions. Especially, when the faults dip α is larger than the angle of critical deformation δ, the scope of surface movement above gob will reduce (Figure 2a), while the angle of critical deformation δ is larger than the faults dip α, the scope of surface movement above gob will enlarge ( Figure 2b). Therefore, the weak fault planes may aggravate the mining damage, and may slow down the mining damage.
Energies 2019, 12, x FOR PEER REVIEW  3 of 19 underground mining, although derived from the mining of coal seams, were essentially controlled by these rock mass discontinuities. Under different rock mass discontinuities, the extent of ground subsidence and cracking in coal mining areas caused by the same intensity of underground mining activities were significantly different. Furthermore, the surface subsidence law was significantly different from that of surface subsidence under ordinary geological conditions [19]. This paper studied the relationship between the movement characteristic of overlying strata and weak planes by means of theoretical analysis, numerical simulation and field measurement. Moreover, the technology of "mining-induced surface movement scope control by hydraulic fracturing" was put forward for the first time. The field test verified the effect of surface movement scope controlled by overlying strata hydraulic fracturing during longwall mining, providing technical support and acting as reference for similar trials investigating the mitigation of surface subsidence and building damages in coal mining districts, which is of great significance for the sustainable development of aging coal mine enterprises.

Mechanism of Mining-Induced Surface Movement Scope Control by Weak Planes
Xia et al. [20,21] systematically studied the control mechanism of tectonic environment (tectonic medium, tectonic interface, tectonic stress, etc.) on coal mining subsidence from the geological point of view and put forward the theory of "tectonic control of hazards" (TCH). Tectonic weak planes are one of the important factors determining the geological environment bearing capacity of coal mining areas. Weak planes destroy the continuity of the rock strata and weaken the mechanical strength of the rock strata, controlling the deformation, failure mode and stability of the rock strata.
Taking faults as an example [22], during the advancing of the working face, the overlying rock mass above the gob is destroyed by tensile stress, which leads to the movement of rock strata since faults belong to the weak structure planes, and have no resistance to the tensile stress. When the movement of the rock strata reaches faults, the overlying strata on both sides of faults are easy to shear and slip along the fault planes. Within the influence area by underground excavation, most of the extrusion displacement and rebound deformation pointing to the gob are limited to the fault zones, which leads to different characteristics of rock strata movement and deformation on both sides of the fault planes. Guo [5] made an analysis with a profile perpendicular to the strike of the fault and studied the influence of faults on the law of strata and surface movement. When the faults dip is greater than 2° and the drop of the fault is greater than 10 m, the law of surface subsidence will get significantly different from that of surface subsidence under ordinary geological conditions. Especially, when the faults dip α is larger than the angle of critical deformation δ, the scope of surface movement above gob will reduce (Figure 2a), while the angle of critical deformation δ is larger than the faults dip α, the scope of surface movement above gob will enlarge ( Figure 2b). Therefore, the weak fault planes may aggravate the mining damage, and may slow down the mining damage. The overlying strata movement and surface subsidence caused by underground mining is a complicated process of mechanical change. Underground coal mining, roof collapse, overburden subsidence and surface subsidence are an overall process of excavation, stress redistribution, deformation and destruction. When solving the coal mining problem under the buildings, the angle The overlying strata movement and surface subsidence caused by underground mining is a complicated process of mechanical change. Underground coal mining, roof collapse, overburden subsidence and surface subsidence are an overall process of excavation, stress redistribution, deformation and destruction. When solving the coal mining problem under the buildings, the angle of critical deformation δ and the loose layer deformation angle ϕ are usually used to define N (the boundary of surface subsidence basin) from M (the boundary of gob) to define the mining influence scope. When the coal seam is mined out, the overlying strata can be divided into three different mining influence zones in the horizontal direction according to their vertical displacement [23]: (A) the stacking zone of broken rock block, (B) the triangular slip zone and (C) the undisturbed zone, as shown in Figure 3. of critical deformation δ and the loose layer deformation angle φ are usually used to define N (the boundary of surface subsidence basin) from M (the boundary of gob) to define the mining influence scope. When the coal seam is mined out, the overlying strata can be divided into three different mining influence zones in the horizontal direction according to their vertical displacement [23]: (A) the stacking zone of broken rock block, (B) the triangular slip zone and (C) the undisturbed zone, as shown in Figure 3. (1) The stacking zone of broken rock block. When the underground coal seam mines out, the immediate roof of gob moves downward and bends under the action of gravity stress and overlying strata. When the internal tensile stress of the immediate roof exceeds the tensile strength of it, it breaks and collapses one after another. Then, the main roof moves and bends along the normal direction of the bedding plane in the form of "cantilever beam" structure, resulting in fracture and separation. Because of the stress redistribution in the rock mass around the mining area caused by strata movement, the layered bent strata subsides, and the collapsed and fractured rock blocks are compacted gradually.
In this area, the rock strata are severely broken, and there is a general lack of force between the broken rock blocks in the horizontal direction. When the upper rock blocks break, they simply pile up on the lower broken rock blocks and fill the gob. Therefore, this area is named as "the stacking zone of broken rock block".
(2) The triangular slip zone (also known as "laminated beam" structure zone with dislocation end).
With the advancing of the working face and the continuing development of the overburden movement, the gob gradually compacts. The bedrock strata besides the gob gradually dislocate upward from the coal seam, forming a "laminated beam" structure. Due to the uncompacted space in the gob, the "laminated beam" structure rotates and moves toward the gob, forming "the triangular slip zone" [24]. According to the morphological characteristics of laminated beams in this area, Yan [23] also vividly named the composite beam as the "laminated beam structure zone with dislocation end". The boundary between "the triangular slip zone" and "the undisturbed zone" is the mechanical boundary. The rock strata in the undisturbed zone bear primary rock stress only, while the rock strata in the "triangular slip zone" bear both primary rock stress and side abutment stress. When the mining area is large enough, the bottom-up rock movement develops to the surface gradually, and the overburden subsides in an inverted trapezoid, forming a subsidence deformation basin on the surface with a much larger size than the gob area [1].
(3) The undisturbed zone. The rock strata in this area is not disturbed. There is no change in the stress and displacement, and surface is not disturbed in this area.
In the field of coal mining, hydraulic fracturing technology was proposed as an effective method to fracture hard rock mass, achieving strength weakening, stress release, crack expansion and slowing of mine pressure. Li et al. [25] proposed hydraulic fracturing for hard top-coal and optimization of the caving process to lower the integrity and strength of the top-coal, and upgrade its crushing effect to weaken the support load effectively during the weight period, which reduces the likelihood of chock support accidents occurring. Yan et al. [26] explored the mechanism of hydraulic fracturing (1) The stacking zone of broken rock block. When the underground coal seam mines out, the immediate roof of gob moves downward and bends under the action of gravity stress and overlying strata. When the internal tensile stress of the immediate roof exceeds the tensile strength of it, it breaks and collapses one after another. Then, the main roof moves and bends along the normal direction of the bedding plane in the form of "cantilever beam" structure, resulting in fracture and separation. Because of the stress redistribution in the rock mass around the mining area caused by strata movement, the layered bent strata subsides, and the collapsed and fractured rock blocks are compacted gradually.
In this area, the rock strata are severely broken, and there is a general lack of force between the broken rock blocks in the horizontal direction. When the upper rock blocks break, they simply pile up on the lower broken rock blocks and fill the gob. Therefore, this area is named as "the stacking zone of broken rock block".
(2) The triangular slip zone (also known as "laminated beam" structure zone with dislocation end).
With the advancing of the working face and the continuing development of the overburden movement, the gob gradually compacts. The bedrock strata besides the gob gradually dislocate upward from the coal seam, forming a "laminated beam" structure. Due to the uncompacted space in the gob, the "laminated beam" structure rotates and moves toward the gob, forming "the triangular slip zone" [24]. According to the morphological characteristics of laminated beams in this area, Yan [23] also vividly named the composite beam as the "laminated beam structure zone with dislocation end". The boundary between "the triangular slip zone" and "the undisturbed zone" is the mechanical boundary. The rock strata in the undisturbed zone bear primary rock stress only, while the rock strata in the "triangular slip zone" bear both primary rock stress and side abutment stress. When the mining area is large enough, the bottom-up rock movement develops to the surface gradually, and the overburden subsides in an inverted trapezoid, forming a subsidence deformation basin on the surface with a much larger size than the gob area [1].
(3) The undisturbed zone. The rock strata in this area is not disturbed. There is no change in the stress and displacement, and surface is not disturbed in this area.
In the field of coal mining, hydraulic fracturing technology was proposed as an effective method to fracture hard rock mass, achieving strength weakening, stress release, crack expansion and slowing of mine pressure. Li et al. [25] proposed hydraulic fracturing for hard top-coal and optimization of the caving process to lower the integrity and strength of the top-coal, and upgrade its crushing effect to weaken the support load effectively during the weight period, which reduces the likelihood of chock support accidents occurring. Yan et al. [26] explored the mechanism of hydraulic fracturing dealing with a hard roof. Huang et al. [27] constructed the hydraulic fracturing theory and technological framework of coal and rock mass. Successful field tests in the Tongxin coal mine (Datong mine area, China) showed that hydraulic fracturing in both a hanging roof over an adjacent gob area and in the gob area behind the advancing working face controlled the behavior of strong strata material on the gob-side in longwall mining and also guaranteed safe extraction at the working face. Yu et al. [28] took Tashan coal mine (Datong mine area, China) as an example and performed ground hydraulic action experiments on the high-level hard rock strata. The experiments changed the fracturing structures of the high-level rock strata and controlled the strong strata behavior successfully, thereby relieving the stress concentration around the working face. Feng et al. [29] used the Polish-like drilling and grooving method to control the expansion of the fracture surface. In his study, the hydraulic fracturing technology was successfully applied to the top of the hard roof of the working face for the first time in Wangtaipu coal mine (Jincheng Anthracite Mining Group, China), avoiding the safety hazard caused by a large-scale roof collapse. The test results showed that the transverse notch was introduced successfully in the roof by means of KZ54 drill and was capable of reducing the pressure required for crack breakdown during the directional hydraulic fracturing. The inflatable straddle packers adopted for borehole sealing ensures the stepwise fracturing in the hard strata, the stepwise fracturing could create multiple cracks in the roof. This test proved the feasibility of hydraulic fracturing in a hard rock formation. However, no field trials of this technique had been implemented in controlling mining-induced surface subsidence scope in China, and thus, ability to control subsidence in mining influenced area has yet to be demonstrated. In this study, field trials using this technique were carried out during longwall mining under buildings in the Sihe coal mine, China.
Within the influence area by underground excavation, the existence of weak planes would change the normal distribution law of displacement field and stress field of surrounding rock strata. The broken zone formed by crisscrossed weak planes played the role of a barrier affecting the deformation and stress propagation of overlying strata, as shown in Figure 4. dealing with a hard roof. Huang et al. [27] constructed the hydraulic fracturing theory and technological framework of coal and rock mass. Successful field tests in the Tongxin coal mine (Datong mine area, China) showed that hydraulic fracturing in both a hanging roof over an adjacent gob area and in the gob area behind the advancing working face controlled the behavior of strong strata material on the gob-side in longwall mining and also guaranteed safe extraction at the working face. Yu et al. [28] took Tashan coal mine (Datong mine area, China) as an example and performed ground hydraulic action experiments on the high-level hard rock strata. The experiments changed the fracturing structures of the high-level rock strata and controlled the strong strata behavior successfully, thereby relieving the stress concentration around the working face. Feng et al. [29] used the Polish-like drilling and grooving method to control the expansion of the fracture surface. In his study, the hydraulic fracturing technology was successfully applied to the top of the hard roof of the working face for the first time in Wangtaipu coal mine (Jincheng Anthracite Mining Group, China), avoiding the safety hazard caused by a large-scale roof collapse. The test results showed that the transverse notch was introduced successfully in the roof by means of KZ54 drill and was capable of reducing the pressure required for crack breakdown during the directional hydraulic fracturing. The inflatable straddle packers adopted for borehole sealing ensures the stepwise fracturing in the hard strata, the stepwise fracturing could create multiple cracks in the roof. This test proved the feasibility of hydraulic fracturing in a hard rock formation. However, no field trials of this technique had been implemented in controlling mining-induced surface subsidence scope in China, and thus, ability to control subsidence in mining influenced area has yet to be demonstrated. In this study, field trials using this technique were carried out during longwall mining under buildings in the Sihe coal mine, China. Within the influence area by underground excavation, the existence of weak planes would change the normal distribution law of displacement field and stress field of surrounding rock strata. The broken zone formed by crisscrossed weak planes played the role of a barrier affecting the deformation and stress propagation of overlying strata, as shown in Figure 4. The weak planes destroyed the cohesion of rock mass, and the rock stratum was easy to slide along the fault under the action of high-strength shear stress. Due to the barrier effect of the weak planes on the movement of the overlying strata, the boundary of the surface movement basin redrew by the broken zone dip α and the loose layer deformation angle φ and moves from N to N′ ( Figure  4). Obviously, the scope of the triangular slip zone decreased while the scope of the undisturbed zone expanded, achieving the purpose of controlling the surface movement scope and protecting the surface buildings.

Numerical Simulation Analysis
To verify the effectiveness of the proposed method, numerical simulations were performed to study the influence of weak planes on overlying strata movement. The discontinuity characteristics in 3DEC software could well express the faults, joints and contact surfaces encountered in mining, and the sliding, opening and closing of the above structural planes could be well simulated. The

Tensile crack
Stepped crack The weak planes destroyed the cohesion of rock mass, and the rock stratum was easy to slide along the fault under the action of high-strength shear stress. Due to the barrier effect of the weak planes on the movement of the overlying strata, the boundary of the surface movement basin redrew by the broken zone dip α and the loose layer deformation angle ϕ and moves from N to N ( Figure 4). Obviously, the scope of the triangular slip zone decreased while the scope of the undisturbed zone expanded, achieving the purpose of controlling the surface movement scope and protecting the surface buildings.

Numerical Simulation Analysis
To verify the effectiveness of the proposed method, numerical simulations were performed to study the influence of weak planes on overlying strata movement. The discontinuity characteristics in 3DEC software could well express the faults, joints and contact surfaces encountered in mining, and Energies 2019, 12, 3480 6 of 18 the sliding, opening and closing of the above structural planes could be well simulated. The program could also better calculate the movement and deformation of overlying strata caused by coal seam excavation and revealed the variation law of displacement field of overlying strata.

Description of Mining and Geological Conditions
Sihe coal mine is located in Jincheng City, Shanxi Province, China ( Figure 5). The coalfield area is 114.5 km 2 with 15 coal seams, of which 3# coal seam and 5# coal seam are the main workable seams. The panel 5304 is 245.7 m wide and 3212.7 m long on average. The coal seam designed for panel 5304 belongs to the 3# coal seam with an average thickness, depth and dip angle of 6.1 m, 240 m and 3 • (nearly horizontal), respectively. The simplified stratigraphic column of the panel 5304 is shown in Figure 6. program could also better calculate the movement and deformation of overlying strata caused by coal seam excavation and revealed the variation law of displacement field of overlying strata.

Description of Mining and Geological Conditions
Sihe coal mine is located in Jincheng City, Shanxi Province, China ( Figure 5). The coalfield area is 114.5 km 2 with 15 coal seams, of which 3# coal seam and 5# coal seam are the main workable seams. The panel 5304 is 245.7 m wide and 3212.7 m long on average. The coal seam designed for panel 5304 belongs to the 3# coal seam with an average thickness, depth and dip angle of 6.1 m, 240 m and 3° (nearly horizontal), respectively. The simplified stratigraphic column of the panel 5304 is shown in Figure 6.    program could also better calculate the movement and deformation of overlying strata caused by coal seam excavation and revealed the variation law of displacement field of overlying strata.

Description of Mining and Geological Conditions
Sihe coal mine is located in Jincheng City, Shanxi Province, China ( Figure 5). The coalfield area is 114.5 km 2 with 15 coal seams, of which 3# coal seam and 5# coal seam are the main workable seams. The panel 5304 is 245.7 m wide and 3212.7 m long on average. The coal seam designed for panel 5304 belongs to the 3# coal seam with an average thickness, depth and dip angle of 6.1 m, 240 m and 3° (nearly horizontal), respectively. The simplified stratigraphic column of the panel 5304 is shown in Figure 6.

Introduction of 3DEC Software
The 3DEC 5.20 software is based on the theory of UDEC software. According to Newton's second law and stress-displacement law, it deals with the mechanical behavior of rock blocks and joints. 3DEC regards the research object as a set of two basic elements, continuity (such as rock mass) and discontinuity (such as structural plane) and defines the mechanical behavior of these basic elements by mechanical laws. Firstly, the motion of blocks is calculated by Newton's second law, and the velocity and displacement of blocks are calculated by known forces. Then, combined with the stress-displacement law, according to the displacement of blocks, the stress between discontinuous planes in rock mass is calculated as the initial boundary condition for the next time order calculation.

Establishment of Model and Selection of Parameters
(1) Determination of Model Size In the construction of the model, there was an influence area from the excavation area to the boundary of the model, and the dip angle of the model coal seam was taken as 0 • according to the horizontal coal seam. The simulated size designed to be 800 m × 250 m × 10 m (length × height × width) along the incline of the working face and the roof to surface, as shown in Figure 7. In the model, the rock strata with a width of 10 m and a low mechanical property was used to simulate the broken zone formed by hydraulic fracturing [30]. The 3DEC analysis software has the function of dividing joints, and the reasonable division of joints can improve the quality of simulation analysis. The principles of joint division in this paper are as follows: 1.
The division of joints in the surface soil layer is smaller than that in general bedrock, and the size of its fragmentation is relatively small.

2.
The fragmentation of the joints of the key layers in the bedrock is longer than that of other strata, and the division is mainly based on the fracture distance.

Introduction of 3DEC Software
The 3DEC 5.20 software is based on the theory of UDEC software. According to Newton's second law and stress-displacement law, it deals with the mechanical behavior of rock blocks and joints. 3DEC regards the research object as a set of two basic elements, continuity (such as rock mass) and discontinuity (such as structural plane) and defines the mechanical behavior of these basic elements by mechanical laws. Firstly, the motion of blocks is calculated by Newton's second law, and the velocity and displacement of blocks are calculated by known forces. Then, combined with the stressdisplacement law, according to the displacement of blocks, the stress between discontinuous planes in rock mass is calculated as the initial boundary condition for the next time order calculation.

Establishment of Model and Selection of Parameters
(1) Determination of Model Size In the construction of the model, there was an influence area from the excavation area to the boundary of the model, and the dip angle of the model coal seam was taken as 0° according to the horizontal coal seam. The simulated size designed to be 800 m × 250 m × 10 m (length × height × width) along the incline of the working face and the roof to surface, as shown in Figure 7. In the model, the rock strata with a width of 10 m and a low mechanical property was used to simulate the broken zone formed by hydraulic fracturing [30]. The 3DEC analysis software has the function of dividing joints, and the reasonable division of joints can improve the quality of simulation analysis. The principles of joint division in this paper are as follows: 1. The division of joints in the surface soil layer is smaller than that in general bedrock, and the size of its fragmentation is relatively small. 2. The fragmentation of the joints of the key layers in the bedrock is longer than that of other strata, and the division is mainly based on the fracture distance. (2) Determination of Boundary Conditions  The horizontal displacement was fixed at the lateral boundary.  The vertical displacement was fixed at the bottom boundary.  The model is directly simulated to the surface, so the top is a free boundary, and the model is only affected by its own gravity.
(3) Selection of parameters The selection of mechanical parameters comprehensively considered the geological conditions and overburden conditions of the coal seam ( Table 1). The Mohr-Coulomb failure criterion was adopted. The properties required for the model are density, bulk modulus, shear modulus, internal (2) Determination of Boundary Conditions

•
The horizontal displacement was fixed at the lateral boundary.

•
The vertical displacement was fixed at the bottom boundary.

•
The model is directly simulated to the surface, so the top is a free boundary, and the model is only affected by its own gravity. (3) Selection of parameters The selection of mechanical parameters comprehensively considered the geological conditions and overburden conditions of the coal seam ( Table 1). The Mohr-Coulomb failure criterion was adopted. The properties required for the model are density, bulk modulus, shear modulus, internal friction angle, cohesion and tensile strength. The mathematical relationships of these properties are as follows: where K is bulk modulus, E is Young's modulus, G is shear modulus and ν is Poisson's ratio.

Numerical Simulation Results and Analysis
Four numerical models (Plan I) with the same weak planes dip of 90 • but different weak planes heights (0, 60 m, 90 m and 120 m) as changing conditions (a~d); four numerical models (Plan II) with the same weak planes height of 60 m but different weak planes dips (0, 30 • , 60 • and 90 • ) as changing conditions (e~h) were established as shown in Figure 8. friction angle, cohesion and tensile strength. The mathematical relationships of these properties are as follows: where K is bulk modulus, E is Young's modulus, G is shear modulus and ν is Poisson's ratio.

Numerical Simulation Results and Analysis
Four numerical models (Plan I) with the same weak planes dip of 90° but different weak planes heights (0, 60 m, 90 m and 120 m) as changing conditions (a~d); four numerical models (Plan II) with the same weak planes height of 60 m but different weak planes dips (0, 30°, 60° and 90°) as changing conditions (e~h) were established as shown in Figure 8.  Plan II Numerical simulation models of different weak planes dips The coal seam was excavated and the gob was stable, the displacement fields of overburden strata in the direction of inclination under different conditions were as shown in Figure 9.

Plan I Distribution of displacement field under different weak planes heights
Plan II Distribution of displacement field under different weak planes dips The broken block in the lower key stratum at the end could not be able to articulate with the incompletely fractured rock stratum in front of it, and formed a "cantilever beam" structure together, while the high key stratum showed the articulation characteristics of the "voussoir beam" structure. According to the distribution of vertical displacement field of overburden (Figure 9), the overburden were divided into three different mining influence zones in the horizontal direction: (A) the stacking The coal seam was excavated and the gob was stable, the displacement fields of overburden strata in the direction of inclination under different conditions were as shown in Figure 9. Plan II Numerical simulation models of different weak planes dips The coal seam was excavated and the gob was stable, the displacement fields of overburden strata in the direction of inclination under different conditions were as shown in Figure 9.  The broken block in the lower key stratum at the end could not be able to articulate with the incompletely fractured rock stratum in front of it, and formed a "cantilever beam" structure together, while the high key stratum showed the articulation characteristics of the "voussoir beam" structure. According to the distribution of vertical displacement field of overburden (Figure 9), the overburden were divided into three different mining influence zones in the horizontal direction: (A) the stacking The broken block in the lower key stratum at the end could not be able to articulate with the incompletely fractured rock stratum in front of it, and formed a "cantilever beam" structure together, while the high key stratum showed the articulation characteristics of the "voussoir beam" structure. According to the distribution of vertical displacement field of overburden (Figure 9), the overburden were divided into three different mining influence zones in the horizontal direction: (A) the stacking zone of broken rock block, (B) the triangular slip zone and (C) the undisturbed zone, as mentioned in Section 2.
In the case of a normal caving working face (Figure 9a,e), due to the geometric symmetry of the excavation area itself, the distribution regularity and symmetry of the vertical displacement field of surrounding rock in the two sides of the gob were obvious after excavation. When there had been a broken zone in the mining affected area, the vertical displacement distribution of excavated surrounding rock were as shown in Figure 9b-d,f-h with weak planes heights of 60 m, 90 m, 120 m and weak planes dips of 30 • , 60 • and 90 • , respectively. Owing to the strength of the weak planes, they were generally lower than that of the rock mass on both sides of it, therefore it was easier for the rock mass to deform and destroy at the broken zone [31]. With the change of rock mass potential energy in the mining affected area, the work done by the gravity stress was mainly consumed by the frictional action and the slip deformation of the weak planes, which made it difficult to move across the broken zone and spread farther outward. The movement and deformation of rock mass directing to gob were limited to the broken zone. Therefore, the displacement of rock mass near the broken zone was obviously large, while the displacement of surrounding rock on the other side of the broken zone was relatively small, showing the barrier effect of the broken zone on the displacement field, and the scope of the triangular slip zone (B) reduced. Under the same burial depth and mining height, the barrier effect became more obvious as the weak planes heights increased.
According to the new "Rules for the mine extraction and coal pillars establishment under buildings, water bodies, railways and main laneways" [32] promulgated and implemented on July 1, 2017, the location where the surface deformation value reaches 10 mm, is usually considered the influencing boundary by mining. The influencing boundaries under different conditions were drawn, as shown in Figure 10. 10.

Plan I Surface movement boundary under different weak planes heights.
Plan II Surface movement boundary under different weak planes dips. In Figure 10, the angle between the connecting line (which connects the line that affects the position of the boundary point and the working face on the side of the coal pillar) and the horizontal line is called the advance angle of influence, which is expressed in ω. The horizontal distance from the influencing boundary point to the working face is called advance distance of influence, which is expressed in l. Finally the numerical simulation results under different fault heights (Table 2) and different fault dips (Table 3) were obtained. The tables above show that the movement boundary of the rock layer had the following rules due to the influence of weak planes.
(1) After the coal seam excavation was completed and the model was stable, the advance angle of influence changed from 53.61 • to 59.15 • with weak planes heights ranging from 0 m to 120 m, and the advance distance of influence changing from 173.31 m to 140.27 m which decreased by 30.04 m. (2) Compared with the condition without weak planes, when the weak planes dips were 30 • and 60 • , the advance distance of influence increased to more than 250 m and 209.83 m, respectively. Meanwhile, the advance angle of influence decreased to less than 45 • and 48.26 • , respectively. However, when the weak planes dip was 90 • , the advance distance of influence decreased to 167.13 m, and the advance angle of influence increased to 54.57 • .
To show the influence of different types of weak planes on surface subsidence characteristics vividly, the comparison diagram of surface subsidence curves was drawn based on the simulation results, as shown in Figure 11.
The surface subsidence curves of Plan I and Plan II were compared in Figure 11. When there were no weak planes, the surface subsidence curve (the black line) symmetrically distributed after mining, and the maximum subsidence point (4537.88 mm) of the surface appeared in the middle of the working face. When there were weak planes with height of 120 m and dip of 90 • at the end of the working face, the maximum subsidence curve (the magenta line) became asymmetric obviously, and the maximum value reached 5204.33 mm. The boundary point of the surface subsidence basin both shifted to the broken zone, and the surface mining influence scope decreased. With the increase of weak planes height, the mining influence scope showed a decreased trend. In Plan II, compared with the angle of critical deformation (57 • ) in Sihe coal mine, when there were weak planes of which dips were smaller than the angle of critical deformation in overburden, the surface mining influence scope increased obviously. When there were weak planes which dips were far greater than the angle of critical deformation in overburden, the surface mining influence scope would decrease.  The surface subsidence curves of Plan I and Plan II were compared in Figure 11. When there were no weak planes, the surface subsidence curve (the black line) symmetrically distributed after mining, and the maximum subsidence point (4537.88 mm) of the surface appeared in the middle of the working face. When there were weak planes with height of 120 m and dip of 90°at the end of the working face, the maximum subsidence curve (the magenta line) became asymmetric obviously, and the maximum value reached 5204.33 mm. The boundary point of the surface subsidence basin both shifted to the broken zone, and the surface mining influence scope decreased. With the increase of weak planes height, the mining influence scope showed a decreased trend. In Plan II, compared with the angle of critical deformation (57°) in Sihe coal mine, when there were weak planes of which dips were smaller than the angle of critical deformation in overburden, the surface mining influence scope increased obviously. When there were weak planes which dips were far greater than the angle of critical deformation in overburden, the surface mining influence scope would decrease.
In the numerical models established in this paper, the broken zone created by directional hydraulic fracturing was regarded as homogeneous weakened rock strata, and the actual fracturing zone was composed of a complex fracture network. Future research work will be conducted to establish a numerical model considering the complex fracture network formed in overburden, and further study the effects of weak planes on mining-induced overlying strata deformation and surface movement.  In the numerical models established in this paper, the broken zone created by directional hydraulic fracturing was regarded as homogeneous weakened rock strata, and the actual fracturing zone was composed of a complex fracture network. Future research work will be conducted to establish a numerical model considering the complex fracture network formed in overburden, and further study the effects of weak planes on mining-induced overlying strata deformation and surface movement.

Study Area Overview
There are extensive buildings nearby besides the main remove gateway of panel 5304, as shown in Figure 12.
Mining subsidence prediction was used to assess possible subsidence influences on surface structures based on the probability-integral method. According to the specific geological conditions of panel 5304, the calculated parameters were as follows: subsidence factor q = 0.83; displacement factor b = 0.21; tangent of major influence angle tan β = 2.4; influence transference angle θ = 86.6 • . Taking these parameters into account, the predicted results were visualized by Surfer software, and the surface subsidence contour map of the 5304 working face after mining was fitted, as shown in Figure 13.

Study Area Overview
There are extensive buildings nearby besides the main remove gateway of panel 5304, as shown in Figure 12. Mining subsidence prediction was used to assess possible subsidence influences on surface structures based on the probability-integral method. According to the specific geological conditions of panel 5304, the calculated parameters were as follows: subsidence factor q = 0.83; displacement factor b = 0.21; tangent of major influence angle tan β = 2.4; influence transference angle θ = 86.6°. Taking these parameters into account, the predicted results were visualized by Surfer software, and the surface subsidence contour map of the 5304 working face after mining was fitted, as shown in Figure 13. Taking the subsidence of 10 mm as the surface moving basin boundary, the calculated contour lines of surface subsidence and the most external influence boundary (the surface subsidence 10 mm contour line) are shown by the red line in Figure 13. It was obvious that some surface buildings (structures) were in the mining influence scope of panel 5304.
To control the mining damage effectively, overburden hydraulic fracturing technology was applied to create artificial weak planes in the tailgate of panel 5304 according to the spatial position relationship between surface buildings and working face, as shown in Figure 14.

Panel 5304
Subsidence basin boundary  Mining subsidence prediction was used to assess possible subsidence influences on surface structures based on the probability-integral method. According to the specific geological conditions of panel 5304, the calculated parameters were as follows: subsidence factor q = 0.83; displacement factor b = 0.21; tangent of major influence angle tan β = 2.4; influence transference angle θ = 86.6°. Taking these parameters into account, the predicted results were visualized by Surfer software, and the surface subsidence contour map of the 5304 working face after mining was fitted, as shown in Figure 13. Taking the subsidence of 10 mm as the surface moving basin boundary, the calculated contour lines of surface subsidence and the most external influence boundary (the surface subsidence 10 mm contour line) are shown by the red line in Figure 13. It was obvious that some surface buildings (structures) were in the mining influence scope of panel 5304.
To control the mining damage effectively, overburden hydraulic fracturing technology was applied to create artificial weak planes in the tailgate of panel 5304 according to the spatial position relationship between surface buildings and working face, as shown in Figure 14.

Panel 5304
Subsidence basin boundary Taking the subsidence of 10 mm as the surface moving basin boundary, the calculated contour lines of surface subsidence and the most external influence boundary (the surface subsidence 10 mm contour line) are shown by the red line in Figure 13. It was obvious that some surface buildings (structures) were in the mining influence scope of panel 5304.
To control the mining damage effectively, overburden hydraulic fracturing technology was applied to create artificial weak planes in the tailgate of panel 5304 according to the spatial position relationship between surface buildings and working face, as shown in Figure 14.

Design of Overburden Hydraulic Fracturing
Hydraulic fracturing weakening should be coordinated with mining work in time and space. The hydraulic fracturing at the face end generally was carried out before mining, and artificial hydraulic cracks generated in the target rock layer in advance to weaken the overburden, based on the bearing pressure distribution and the hydraulic fracturing weakening characteristics, the hydraulic fracturing construction site is determined. At the same time, considering the safety of underground hydraulic fracturing, measures should be taken to strengthen support in fracturing area to avoid water fracturing destroying the stability of roof support body, thus causing roof fall and other accidents. During field implementation, the effect of hydraulic fracturing on subsidence scope controlling of surface mining could be judged by the distribution of surface cracks and the observed surface subsidence value.
According to the wetting radius R of the water injection hole [33], the spacing is based on the design criteria of forming through cracks. To achieve a good crushing effect, the adjacent fractured Open-off cut Main remove gateway

Design of Overburden Hydraulic Fracturing
Hydraulic fracturing weakening should be coordinated with mining work in time and space. The hydraulic fracturing at the face end generally was carried out before mining, and artificial hydraulic cracks generated in the target rock layer in advance to weaken the overburden, based on the bearing pressure distribution and the hydraulic fracturing weakening characteristics, the hydraulic fracturing construction site is determined. At the same time, considering the safety of underground hydraulic fracturing, measures should be taken to strengthen support in fracturing area to avoid water fracturing destroying the stability of roof support body, thus causing roof fall and other accidents. During field implementation, the effect of hydraulic fracturing on subsidence scope controlling of surface mining could be judged by the distribution of surface cracks and the observed surface subsidence value.
According to the wetting radius R of the water injection hole [33], the spacing is based on the design criteria of forming through cracks. To achieve a good crushing effect, the adjacent fractured area should coincide. The maximum distance d between hydraulic fracturing boreholes is twice that of the wetting radius of the water hole: where d is the hydraulic-fracture drilling spacing, R is the water injection hole wetting radius, t is water injection time, V is the water injection speed, γ is the rock bulk density, n is the fractional rock water absorption, k is the nonuniform coefficient (0.08-0.2), l is the length of the water injection hole. According to the numerical simulation results and site's operational conditions, the borehole positions were chosen to be 1 m from the coal wall. The borehole length was 60 m, the angle between the borehole and the axis of the tunnel was 5 • , the angle of elevation was 70 • , and the borehole spacing was 10 m. During the hydraulic fracturing, we adopted the backward fracturing method and used a KZ54 type cutting drill bit to prefabricate multiple radial cutting grooves in the rock layer to increase the number of fractures at the surface. A 3ZSB80/62-90 type high pressure water pump was used in the same drill to form a hydraulic-fracture zone from the bottom of the hole to the bottom every 1.5 m-2 m. A technical diagram of the hydraulic fracturing progress is shown in Figure 15.
(1) Before mining, drill holes were arranged according to the set drilling parameters on the roof of the tailgate. When the bit drilled to the slotting position, the slotting bit replaced and grooved at the bottom of the borehole. After grooving, the slotting bit withdrew, and the ordinary bit replaced it to continue drilling. Then the grooving process was carried out at the next slotting position, so that the drilling cycle would be completed at the end of the drilling, as shown in Figure 15c. (2) Then, the inflatable straddle packer used to seal the fracturing section before and after grooving, and the backward fracturing method used for sealing and water injection fracturing, that is, from the groove at the bottom of the borehole to the outward sequential fracturing. The high-pressure water pump was connected with water first and then electrified, and then slowly pressurized. At the same time, the pressure gauge data of the pump and manual pump recorded. The pressure would continue to pressurize until the pre-crack cracks occurred. At this time, the pressure would suddenly drop, and the pressure-retaining water injection would make the cracks continue to expand. When water seeped out from the roof, coal wall or borehole of the roadway, the fracturing stopped immediately, and the main hydraulic crack and branch airfoil crack was created in the overburden rock at the fracturing site, as shown in Figure 15b. (1) Before mining, drill holes were arranged according to the set drilling parameters on the roof of the tailgate. When the bit drilled to the slotting position, the slotting bit replaced and grooved at the bottom of the borehole. After grooving, the slotting bit withdrew, and the ordinary bit replaced it to continue drilling. Then the grooving process was carried out at the next slotting position, so that the drilling cycle would be completed at the end of the drilling, as shown in Figure 15c.
(2) Then, the inflatable straddle packer used to seal the fracturing section before and after grooving, and the backward fracturing method used for sealing and water injection fracturing, that is, from the groove at the bottom of the borehole to the outward sequential fracturing. The highpressure water pump was connected with water first and then electrified, and then slowly pressurized. At the same time, the pressure gauge data of the pump and manual pump recorded. The pressure would continue to pressurize until the pre-crack cracks occurred. At this time, the pressure would suddenly drop, and the pressure-retaining water injection would make the cracks continue to expand. When water seeped out from the roof, coal wall or borehole of the roadway, the fracturing stopped immediately, and the main hydraulic crack and branch airfoil crack was created in the overburden rock at the fracturing site, as shown in Figure 15b.
(3) The above operations were repeated in turn until a hydraulic fracturing broken zone intersecting airfoil branching cracks formed between boreholes as shown in Figure 15a. When the face mined near the hydraulic fracturing section, the overburden at face end slid along the broken zone under the action of overburden stratum pressure and self-gravity.

Application Effects and Field Measurements
A subsidence observation station was established over the working surface to analyze the effect of mining on the nearby buildings after hydraulic fracturing. The setting of the displacement observation station and the layout of the working face are shown in Figure 16.

Application Effects and Field Measurements
A subsidence observation station was established over the working surface to analyze the effect of mining on the nearby buildings after hydraulic fracturing. The setting of the displacement observation station and the layout of the working face are shown in Figure 16. Line X was laid at the hydraulic fracturing area, as shown in Figure 17. In order to analyze and compare any surface movement before and after hydraulic fracturing, Origin Pro.8.6 software was used to draw the measured and the predicted maximum subsidence curve of line X, as shown in Figure 17.  Line X was laid at the hydraulic fracturing area, as shown in Figure 17. In order to analyze and compare any surface movement before and after hydraulic fracturing, Origin Pro.8.6 software was used to draw the measured and the predicted maximum subsidence curve of line X, as shown in Figure 17.
According to Figure 17, the results showed that the line X subsidence curve got steeper after hydraulic fracturing, and the measured maximum subsidence value was greater than the predicted value. The mean values of the predicted and the measured subsidence were −857.99 mm and −962.75 mm, and the standard deviation of them were 1098.29 and 1262.03, respectively. The boundary point of subsidence basin moved 31 m to the middle of the gob after hydraulic fracturing and the mining subsidence scope reduced. Line X was laid at the hydraulic fracturing area, as shown in Figure 17. In order to analyze and compare any surface movement before and after hydraulic fracturing, Origin Pro.8.6 software was used to draw the measured and the predicted maximum subsidence curve of line X, as shown in Figure 17. According to Figure 17, the results showed that the line X subsidence curve got steeper after hydraulic fracturing, and the measured maximum subsidence value was greater than the predicted value. The mean values of the predicted and the measured subsidence were −857.99 mm and −962.75 mm, and the standard deviation of them were 1098.29 and 1262.03, respectively. The boundary point of subsidence basin moved 31 m to the middle of the gob after hydraulic fracturing and the mining subsidence scope reduced.
During the field surface investigation after the end of mining of panel 5304, it was found that the collapse cracks (step drop 0.3-0.7 m) shown in Figure 18 appeared near the corresponding surface position of hydraulic fracturing section along the tailgate. However, combined with household investigation, no cracks or damage occurred after mining with hydraulic fracturing on overburden, which verified that hydraulic fracturing in overburden reduced the mining subsidence scope and protected the constructions above the surface. Moreover, control measures were proposed such as strengthening the observation of cracks evolution and cracks landfill compaction.

Conclusions
(1) After the coal seam mined out, the overburden bedrock was divided into three different mining influence zones in the horizontal direction according to vertical displacement of overburden: the undisturbed zone, the triangular slip zone and the stacking zone of broken rock block. Among them, the characteristics of the triangular slip zone was the key to controlling the surface movement scope.
(2) Numerical simulation models were established to study the influence of weak planes on overlying strata movement. When the weak planes dip was bigger than the angle of critical deformation, the advance angle of influence changed from 53.61° to 59.15° with weak plane heights ranging from 0 m to 120 m, and the advance distance of influence changed from 173.31 m to 140.27 m which decreased by 30.04 m. With the increase of the weak planes' height, the mining influence scope showed a decreasing trend.
(3) The field test results of panel 5304 in Sihe coal mine showed that buildings (structures) formerly in mining-affected areas did not appear to have cracks or mining damages after mining by using hydraulic fracturing technology in overburden, which proved the feasibility of this technology. This work is of important theoretical value and practical significance to surface subsidence control However, combined with household investigation, no cracks or damage occurred after mining with hydraulic fracturing on overburden, which verified that hydraulic fracturing in overburden reduced the mining subsidence scope and protected the constructions above the surface. Moreover, control measures were proposed such as strengthening the observation of cracks evolution and cracks landfill compaction.

Conclusions
(1) After the coal seam mined out, the overburden bedrock was divided into three different mining influence zones in the horizontal direction according to vertical displacement of overburden: the undisturbed zone, the triangular slip zone and the stacking zone of broken rock block. Among them, the characteristics of the triangular slip zone was the key to controlling the surface movement scope. (2) Numerical simulation models were established to study the influence of weak planes on overlying strata movement. When the weak planes dip was bigger than the angle of critical deformation, the advance angle of influence changed from 53.61 • to 59.15 • with weak plane heights ranging from 0 m to 120 m, and the advance distance of influence changed from 173.31 m to 140.27 m which decreased by 30.04 m. With the increase of the weak planes' height, the mining influence scope showed a decreasing trend.
(3) The field test results of panel 5304 in Sihe coal mine showed that buildings (structures) formerly in mining-affected areas did not appear to have cracks or mining damages after mining by using hydraulic fracturing technology in overburden, which proved the feasibility of this technology. This work is of important theoretical value and practical significance to surface subsidence control and sustainable coal mining with similar geological mining conditions, and provides a new method for mining damage and protection under buildings (constructions), railways, and water bodies.
Author Contributions: W.G. conceived this paper and made modifications several times; Z.F., D.Y. and F.X. analyzed the data and calculated the results; W.Y. conducted the numerical simulation; Z.F. wrote the paper.

Funding:
The authors gratefully acknowledge the financial support for this work provided by the National Natural Science Foundation of China (grant no. 51774111).