Characteristics and Control of Mining Induced Fractures above Longwall Mines Using Backﬁlling

: Water conservation in mining is the key to solving the conﬂict between coal resource exploitation and ecological environment protection, especially in arid and semi-arid mining areas. Continuous excavation and continuous backﬁlling (CECB) in longwall mining is an important method to realize water conservation mining. Considering the di ﬀ erent boundary conditions of the main roof stress in di ﬀ erent mining phases, the mechanical models of clamped–clamped beam, continuous beam, and elastic foundation beam among ﬁlling body, main roof, and strata are established. Furthermore, the spatio-temporal evolution mechanisms of mining-induced fractures (MIF) are studied. It is found that there is a hyperbolic function relationship between MIF and the mining roadway (MR) ﬁlling percentages. Based on mining the XV coal seam under CECB in the Wangtaipu Coal Mine, the distribution patterns of MIF are studied. It is concluded that the distribution pattern is an isosceles trapezoid with the moving angle of overlying strata as the bottom angle, and the upper and lower boundary of MIF as the two parallel sides. Based on the inﬂuence coe ﬃ cient of MR ﬁlling percentages on MIF, the curve of the MIF height is divided into three ranges, which include the stability control range, the critical range, and the lost control range. The controlling e ﬀ ects of MR ﬁlling percentages are studied, and the calculation expression of the MIF height in the stability control range is given. In engineering practice, 90% MR ﬁlling percentage is used for CECB. The MIF height is about 3.0 times of mining height, and the main roof beam is not broken. The water-resisting property of aquiclude III is not destroyed, thus, the mining does not adversely impact the water. The results provide theories and practices for controlling MIF under CECB in the conditions of extremely close distance aquifers.


Introduction
China is the largest coal producing nation in the world. For the past decade, the average coal production in China is approximately 3.48 billion tons a year, which represents 43.81% of the total coal production in the world [1]. Coal has been a tremendous booster for China's economic development, but has also brought considerable negative impacts on ecological environment, especially on water resource. On the average, 2.04 m 3 water should be discharged for producing a ton of coal [2][3][4]. Coal mining can cause the deterioration of regional ecological environment especially in arid and semi-arid mining areas which have vulnerable ecosystems [5,6]. Achieving the conservation of water resource in coal mining process is of great significance to maintaining ecological balance in the mines [7][8][9][10].
As early as the 1960s, with the aim of attempting the conservation of water resource in the coal mining process, Austrian scholars began to use tracer agents to examine the relation between subsidence and underground water [11]. American researchers investigated the strata disturbance region while performing coal mining under a surface water, and established related mathematical models to predict water diversion performances of the mining-induced fractures [12,13]. Chinese scholars proposed the concept of water conservation mining from the end of the 20th century. After the development over the past 30 years, systematical water conservation mining technical and theoretical system has been initially established in China [14][15][16][17]. At present, the water conservation mining methods mainly include drift mining, slice mining (height-limit mining), room and pillar mining, and backfilling. Specifically, drift mining exhibits certain limitations, including low coal recovery rate and difficulty in recovery of stranded coal pillars; using slice mining (height-limit mining), the mining time and the effects on aquifer can hardly be determined; room and pillar mining also encounters some difficulties, such as low coal recovery rate and complex production processes. By contrast, the backfilling mining (also referred to mining with backfilling) becomes the most effective water conservation mining method; nevertheless, traditional longwall backfilling technology also shows some problems, such as insufficient constraints of backfilling space and great difficulties in mining activity coordination, thereby seriously restricting the development of water conservation mining [18][19][20].
In recent years, by combining the advantages and concepts of the Wongawilli rapid mining method [21,22] and pillar backfilling mining [23], the writers proposed continuous excavation and continuous backfilling (CECB) in longwall mining and explored the preparation of filling materials [24,25]. This method adopts the interval excavated MRs for coal mining. Once a MR is excavated, it will be filled immediately. At the same time, the next MR will be excavated until all MRs of the whole working face are excavated and filled, as shown in Figure 1. This method can successfully overcome the limitations in traditional longwall backfilling mining and provide the backfilling face with sufficient time and effective space. Accordingly, the filling bodies can be solidified, and the requirements on bearing strength can be satisfied to effectively control mining-induced fractures (MIF). subsidence and underground water [11]. American researchers investigated the strata disturbance region while performing coal mining under a surface water, and established related mathematical models to predict water diversion performances of the mining-induced fractures [12,13]. Chinese scholars proposed the concept of water conservation mining from the end of the 20th century. After the development over the past 30 years, systematical water conservation mining technical and theoretical system has been initially established in China [14][15][16][17]. At present, the water conservation mining methods mainly include drift mining, slice mining (height-limit mining), room and pillar mining, and backfilling. Specifically, drift mining exhibits certain limitations, including low coal recovery rate and difficulty in recovery of stranded coal pillars; using slice mining (height-limit mining), the mining time and the effects on aquifer can hardly be determined; room and pillar mining also encounters some difficulties, such as low coal recovery rate and complex production processes. By contrast, the backfilling mining (also referred to mining with backfilling) becomes the most effective water conservation mining method; nevertheless, traditional longwall backfilling technology also shows some problems, such as insufficient constraints of backfilling space and great difficulties in mining activity coordination, thereby seriously restricting the development of water conservation mining [18][19][20]. In recent years, by combining the advantages and concepts of the Wongawilli rapid mining method [21,22] and pillar backfilling mining [23], the writers proposed continuous excavation and continuous backfilling (CECB) in longwall mining and explored the preparation of filling materials [24,25]. This method adopts the interval excavated MRs for coal mining. Once a MR is excavated, it will be filled immediately. At the same time, the next MR will be excavated until all MRs of the whole working face are excavated and filled, as shown in Figure 1. This method can successfully overcome the limitations in traditional longwall backfilling mining and provide the backfilling face with sufficient time and effective space. Accordingly, the filling bodies can be solidified, and the requirements on bearing strength can be satisfied to effectively control mining-induced fractures (MIF). It is one of the important scientific problems in water conservation mining that avoiding MIF conducting through aquifers [26]. Scholars have undertaken considerable researches on the law of MIF and its influence on water resources when longwall mining. American scholars have studied the importance of aquicludes between aquifers and coal seam by studying MIF in longwall working face [27,28]. Ukrainian scholars studied the distribution pattern and evolution law of MIF in a longwall working face of the Donetsk coalfield, and divided the overlying strata into caving zone, fracture zone and continuous deformation zone [29]. At the same time, the caving zone and fissure zone were defined as MIF connected area of roof water. Chinese scholars have detected MIF distribution patterns of hundreds of longwall working faces under different geological conditions, and concluded It is one of the important scientific problems in water conservation mining that avoiding MIF conducting through aquifers [26]. Scholars have undertaken considerable researches on the law of MIF and its influence on water resources when longwall mining. American scholars have studied the importance of aquicludes between aquifers and coal seam by studying MIF in longwall working face [27,28]. Ukrainian scholars studied the distribution pattern and evolution law of MIF in a longwall working face of the Donetsk coalfield, and divided the overlying strata into caving zone, fracture zone and continuous deformation zone [29]. At the same time, the caving zone and fissure zone were defined as MIF connected area of roof water. Chinese scholars have detected MIF distribution patterns of hundreds of longwall working faces under different geological conditions, and concluded that MIF in overlying strata are saddle shaped [30,31]. In addition, scholars have laid a theoretical foundation for backfilling water conservation mining by studying the law of MIF under longwall backfilling and drift backfilling conditions [32,33].
However, there are relatively few studies on the control of MIF in roadway backfilling, especially the characteristics and control of MIF under CECB are still unclear. The interval mining technology under CECB leads to alternate distribution between coal pillar (or filling body) and mined roadway [34]. The supporting structures of strata are discontinuous, and the distribution patterns of MIF are also different. In addition, CECB is essentially a backfilling mining method. The main control principle of MIF is to fill before the goaf forms and occupy the subsidence space of strata, reduce the degree of strata migration, and control the development of MIF. The filling percentage is an important parameter to measure filling qualities and filling effects. Therefore, the control of MIF through the mining roadway (MR) filling percentages is one of the important scientific problems of CECB. At present, there is relatively little research in this area.
Based on the engineering background of CECB in the Wangtaipu Coal Mine, and on the basis of different boundary conditions of the main roof stress, this paper establishes the mechanical models of clamped-clamped beam, continuous beam and elastic foundation beam between the filling body, main roof, and strata. Moreover, the mechanism and characteristics of MIF are analyzed, and the controlling effects of filling rates on MIF are studied. The present results should provide theoretical and practical guidance for water conservation mining under extremely close distance aquifers and other practices of water conservation mining with extremely thin barriers.

Profile of the Study Area
The Wangtaipu Coal Mine is located at the south of Shanxi, China, and its mining area is about 33.79 km 2 . The residual recoverable reserves in this mine are approximately 57.46 Mt which are mostly deposited below buildings or water resources. To be specific, the XV23 mining segment has a mean depth of 220 m, a strike length of 330 m, a dip length of 156 m and the reserves of approximately 0.19 Mt. The under mining coal seam, namely XV coal seam, has a mean thickness of 2.38 m and a dip angle of 1~2 • . Figure 2 displays the study area. that MIF in overlying strata are saddle shaped [30,31]. In addition, scholars have laid a theoretical foundation for backfilling water conservation mining by studying the law of MIF under longwall backfilling and drift backfilling conditions [32,33]. However, there are relatively few studies on the control of MIF in roadway backfilling, especially the characteristics and control of MIF under CECB are still unclear. The interval mining technology under CECB leads to alternate distribution between coal pillar (or filling body) and mined roadway [34]. The supporting structures of strata are discontinuous, and the distribution patterns of MIF are also different. In addition, CECB is essentially a backfilling mining method. The main control principle of MIF is to fill before the goaf forms and occupy the subsidence space of strata, reduce the degree of strata migration, and control the development of MIF. The filling percentage is an important parameter to measure filling qualities and filling effects. Therefore, the control of MIF through the mining roadway (MR) filling percentages is one of the important scientific problems of CECB. At present, there is relatively little research in this area.
Based on the engineering background of CECB in the Wangtaipu Coal Mine, and on the basis of different boundary conditions of the main roof stress, this paper establishes the mechanical models of clamped-clamped beam, continuous beam and elastic foundation beam between the filling body, main roof, and strata. Moreover, the mechanism and characteristics of MIF are analyzed, and the controlling effects of filling rates on MIF are studied. The present results should provide theoretical and practical guidance for water conservation mining under extremely close distance aquifers and other practices of water conservation mining with extremely thin barriers.

Profile of the Study Area
The Wangtaipu Coal Mine is located at the south of Shanxi, China, and its mining area is about 33.79 km 2 . The residual recoverable reserves in this mine are approximately 57.46 Mt which are mostly deposited below buildings or water resources. To be specific, the XV23 mining segment has a mean depth of 220 m, a strike length of 330 m, a dip length of 156 m and the reserves of approximately 0.19 Mt. The under mining coal seam, namely XV coal seam, has a mean thickness of 2.38 m and a dip angle of 1~2°. Figure 2 displays the study area.

Stratigraphic Structure
The XV coal seam is located in upper the Taiyuan formation of carboniferous system (C3t), as the stratigraphic structure shown in Figure 3. The main roof of the coal seam is composed of limestone

Stratigraphic Structure
The XV coal seam is located in upper the Taiyuan formation of carboniferous system (C 3 t), as the stratigraphic structure shown in Figure 3. The main roof of the coal seam is composed of limestone  In the stratigraphic structure, there exist three aquifers. Aquifer I is situated below the Quaternary yellow clay overburden with a thickness of 2~5 m. It is an aquifer with water supply significance and ecological value. Aquifer II is located near the boundary between Permian and Carboniferous, mainly consisting of 4~5 aquifer sandstones with a total thickness of about 10.35 m. Aquifer III locates approximately 13 m above the XV coal seam, with a mean thickness of 3.89 m, a normal water inflow of 1~3 m 3 /h, and a maximum water inflow of 10 m 3 /h.

Mining Method of CECB
Aquifer III is within the influence range of XV coal mining. If the caving mining is adopted, the aquifer will be within the range of overburden caving zone (8 times of mining height), which will inevitably cause the damage of aquifer, leading to the loss of water resources and water inrush into the working face. CECB ( Figure 4) is an effective method to solve the problem of water conservation mining under extremely close distance aquifers.  In the stratigraphic structure, there exist three aquifers. Aquifer I is situated below the Quaternary yellow clay overburden with a thickness of 2~5 m. It is an aquifer with water supply significance and ecological value. Aquifer II is located near the boundary between Permian and Carboniferous, mainly consisting of 4~5 aquifer sandstones with a total thickness of about 10.35 m. Aquifer III locates approximately 13 m above the XV coal seam, with a mean thickness of 3.89 m, a normal water inflow of 1~3 m 3 /h, and a maximum water inflow of 10 m 3 /h.

Mining Method of CECB
Aquifer III is within the influence range of XV coal mining. If the caving mining is adopted, the aquifer will be within the range of overburden caving zone (8 times of mining height), which will inevitably cause the damage of aquifer, leading to the loss of water resources and water inrush into the working face. CECB ( Figure 4) is an effective method to solve the problem of water conservation mining under extremely close distance aquifers.
CECB is a mining method of combining wall system mining technology with continuous and rapid room and pillar mining. According to longwall mining method, head entry, tail entry, and open-off cut are laid out [25]. Before excavating, the CECB working face is divided into n mining roadways (MRs), and MRs are divided into m (usually m = 2~5) mining phases. During excavating, MRs are excavated in turn, and the excavated MRs are filled immediately. In the CECB working face, the "parallel excavating and filling" operation mode is always maintained, until all the MRs are excavated and filled.  CECB is a mining method of combining wall system mining technology with continuous and rapid room and pillar mining. According to longwall mining method, head entry, tail entry, and open-off cut are laid out [25]. Before excavating, the CECB working face is divided into n mining roadways (MRs), and MRs are divided into m (usually m = 2~5) mining phases. During excavating, MRs are excavated in turn, and the excavated MRs are filled immediately. In the CECB working face, the "parallel excavating and filling" operation mode is always maintained, until all the MRs are excavated and filled.

Development Mechanisms of MIF
Based on the stress characteristics of the main roof in different mining phases, the mechanical models of the clamped-clamped beam, continuous beam, and elastic foundation beam among the filling body, main roof, and strata are established. The deformation characteristics of the main roof in each mining phase are analyzed, and the expressions of the MIF height are given.

Mechanical Models
After the excavation of MRs, the main roof is subjected to the load of strata and the supporting effect of coal pillars on both sides of MR, forming a rock beam structure. After backfilling, the filling bodies act as cushion and supports the main roof synergistically, providing support for the rock beam structure. During excavation, the exposed area of the main roof and the length of beam in mining space are gradually enlarged. Although the rock beam structure of the main roof is continuous in space, the stress and subsidence of the main roof above the MRs of adjacent mining area are not continuous in time. The deformation of the main roof appears a cusp (non-differentiable point) at the junction of the adjacent MRs.
During the whole CECB, the main roof is a kind of rock beam structure which is always supported by the backfilling and the support resistance increases as mining progresses. However,

Development Mechanisms of MIF
Based on the stress characteristics of the main roof in different mining phases, the mechanical models of the clamped-clamped beam, continuous beam, and elastic foundation beam among the filling body, main roof, and strata are established. The deformation characteristics of the main roof in each mining phase are analyzed, and the expressions of the MIF height are given.

Mechanical Models
After the excavation of MRs, the main roof is subjected to the load of strata and the supporting effect of coal pillars on both sides of MR, forming a rock beam structure. After backfilling, the filling bodies act as cushion and supports the main roof synergistically, providing support for the rock beam structure. During excavation, the exposed area of the main roof and the length of beam in mining space are gradually enlarged. Although the rock beam structure of the main roof is continuous in space, the stress and subsidence of the main roof above the MRs of adjacent mining area are not continuous in time. The deformation of the main roof appears a cusp (non-differentiable point) at the junction of the adjacent MRs.
During the whole CECB, the main roof is a kind of rock beam structure which is always supported by the backfilling and the support resistance increases as mining progresses. However, with the increase of mining phase, the rock beam structure evolves gradually from clamped-clamped beam to continuous beam and then to elastic foundation beam. Figure 5 shows the mechanics models of clamped-clamped beam, continuous beam, and elastic foundation beam among the filling body, main roof, and strata in different mining phases.

Clamped-Clamped Beam Model
After the first mining phase, the coal pillars with a width three times that of the MR are isolated by MRs. According to the engineering practice experience, the interaction between the two MRs should be ignored. It is considered that the influence of the excavating and filling MRs on the strata is independent [35]. In this phase, the coal pillar can provide fixed support constraint for the main roof, which can not only transfer stress, but also transfer bending moment at the constraint end. The stress characteristics of the main roof after the first mining phase are simplified to a clamped-clamped beam model, as shown in Figure 5a. In the figure, L 0 is the length of the beam which is equal to the width of the MR, q is the load on the strata which is simplified as a uniform distribution, and f is the supporting force of the filling body. The connecting joints marked 1 and 2 represent the fixed support constraint of the coal pillars on the main roof.

Elastic Foundation Beam Model
After all MRs are excavated and filled, the main roof is in contact with the filling body continuously. All the coal pillars are replaced by filling, resulting in the removal of stress isolation between continuous beams, and all spacing continuous beams are connected as a whole [37]. In order to establish the mechanical model of elastic foundation beam, the following assumptions are made: (1) the boundary of mining block imposes a constraint on the main roof, and there exists a supporting stress and bending moment on the filling body near the boundary to provide the main roof with a static balance. (2) When the main roof reaches a rebalancing state, the anchorage restraint of the connecting joints are relieved, and all points on the main roof are continuous and differentiable. (3) The strength of filling bodies meets the requirements of supporting and bearing, and there is no failure, so the pillar foundation bearing the main roof can be regarded as the elastic foundation. (4) The internal part of filling bodies can be transferred by interaction, and the elastic foundation can provide a uniformly-distributed support.
The mechanical model of beam on elastic foundation is shown in Figure 5d. Strata act on elastic foundation beam with uniformly distributed load q, and the filling bodies provide uniformly distributed support stress f. In this figure, the elastic foundation beam of 1-4 segment represents that the main roof has reached the rebalance state, the 5-6 segment represents that the basic roof is in the L0 L0 f f f

Continuous Beam Model
With the increase of mining, the exposed area of basic roof increases, and the structure of rock beam will change. In order to establish the mechanical model of the continuous beam, the following basic assumptions are made: (1) The mining technology of CECB can protect the main roof from breaking, which can always transmit stress and bending moment during the mining and filling. (2) Before the excavation of MR, the adjacent filling body has been fully compacted and in a stable state, which can replace the coal pillar to effectively support the main roof. (3) The connecting joint between the two MRs act as a fixing and restraining function to the main roof, and the angles of the main roofs on both sides of the point match each other (rotation angles are equal in size and reverse in direction) [36]. (4) Coal pillars can play a role of stress isolation, and the mechanical characteristics of the main roof can be simplified as a continuous beam structure with multiple intervals and independence. The mechanical model of continuous beam is shown in Figure 5b,c. Figure 4b corresponds to the stress characteristics of the main roof after the second mining phase. The continuous beam model is connected by two equal-length clamped-clamped beams through connecting joint 2. The connecting joint 1 and connecting joint 3 are two fixed-supported ends and are on the same level. Due to the sinking of the main roof, the connecting joint 2 is lower than this level. In Figure 5b, the continuous beam model is asymmetrically distributed at connecting joint 2, and the 1-2 and 2-3 segments rotate θ 2 clockwise and counterclockwise, respectively. Figure 5c corresponds to the stress characteristics of the main roof after the third mining phase, which is an extension of Figure 5b. The continuous beam model is composed of three equal-length clamped-clamped beams through connecting joint 2 and connecting joint 3. Connecting joint 1 and connecting joint 4 are two fixed-supported ends which are on the same level. In the process of forming the 3-4 segment, connecting joint 3 sinks with the main roof. According to the assumption of continuous beam model, it is advisable to simplify the 2-3 segment into horizontal ones. In Figure 5c, the continuous beam model is asymmetrically at the midpoint of the 2-3 segment, and the 1-2 and 3-4 segments rotate θ 3 clockwise and counterclockwise, respectively.

Elastic Foundation Beam Model
After all MRs are excavated and filled, the main roof is in contact with the filling body continuously. All the coal pillars are replaced by filling, resulting in the removal of stress isolation between continuous beams, and all spacing continuous beams are connected as a whole [37]. In order to establish the mechanical model of elastic foundation beam, the following assumptions are made: (1) the boundary of mining block imposes a constraint on the main roof, and there exists a supporting stress and bending moment on the filling body near the boundary to provide the main roof with a static balance. (2) When the main roof reaches a rebalancing state, the anchorage restraint of the connecting joints are relieved, and all points on the main roof are continuous and differentiable. (3) The strength of filling bodies meets the requirements of supporting and bearing, and there is no failure, so the pillar foundation bearing the main roof can be regarded as the elastic foundation. (4) The internal part of filling bodies can be transferred by interaction, and the elastic foundation can provide a uniformly-distributed support.
The mechanical model of beam on elastic foundation is shown in Figure 5d. Strata act on elastic foundation beam with uniformly distributed load q, and the filling bodies provide uniformly distributed support stress f. In this figure, the elastic foundation beam of 1-4 segment represents that the main roof has reached the rebalance state, the 5-6 segment represents that the basic roof is in the original rock stress state, and the 4-5 segment represents the transition stage from the original rock stress state to the rebalance state. Connecting joints 1, 2, and 3 represent the joints after removing the fixed support and constraints, which only provide support to the elastic foundation beam and no longer transfer the bending moment. Connecting joint 4 represents the critical point of the elastic foundation beam reaching the rebalancing state. Connecting joint 5 represents the boundary of the CECB working face, which is the fixed support end. Connecting joint 6 represents the original rock stress state.

Deformation Characteristics of the Beam Models
Based on the above analysis, taking the main roof as an example, the deformation characteristics of strata after each mining phase are studied.

The Structure of the Continuous Beam
After the third mining phase, the static equilibrium equation of the beam structure in the main roof can be written as: The maximum bending moment of the main roof at the connecting joints can be acquired and, therefore, the bending moments at different connecting joints can be calculated as: where F i and N i , the supporting force and the supporting force components at the i-th connecting joint (i = 1, 2, . . . , 5), respectively, kN; M ij, the bending moments at the i-th connecting joint at the j-th phase (i = 1, 2, . . . , 5; j = 1, 2, 3, 4), kN·m; I 0, the inertia moment of the main roof, m 4 ; L 0, the width of the MR, m; f, the supporting force of the filling body, MPa; q, the load on the strata, MPa; θ j , the inclination angle of the beam at the j-th phase (j = 1, 2, 3, 4), • ; and E, the elastic modulus of the main roof, GPa. The conditions at the third mining phase are evolved from those at the first two mining phases. The bending moment equations at the first two mining phases can thus be derived according to Equation (4) [17]. The bending moments at each connecting joint after the completion of the first two mining phases can be written as:

The Structure of the Elastic Foundation Beam
Based on Winkler elastic foundation beam theories, the supporting strength of the bearing filling bodies in the main roof can be calculated as: where w(x) is the subsidence of the main roof, and m and k are the foundation coefficients. The subsidence of the main roof, the load in the strata, and the supporting strength of the filling bodies satisfy the following relation: As shown in Figure 4d, the initial subsidence of the main roof, denoted as h 0 (with a unit of m) can be calculated as: where x is the distance between the sinking position and coal wall, m.

Fractures in the Overlying Strata
The critical condition for the generation of fractures in the overlying strata is that the tensile strain produced in bending deflection of the strata reaches the ultimate strain. Assuming that fractures in the strata develop towards the n-th stratum above the coal seam, a unit with a length of dx and a length of 2y was selected on the edge of fracture development. The distances between the upper/lower boundary of the unit and the neutral axis O-O are denoted as y. As shown in Figure 6, during the deformation of left and right cross-sections, the unit is rotated by dϕ around the P-P axis, with a curvature of ρ.  Under critical conditions, the strain of the side AB of the unit in the stratum can be calculated as: The ultimate strain (εm) of the stratum can be calculated as: where hR is the thickness of the n-th stratum, m; γ is the unit weight of the n-th stratum, kN/m 3   Under critical conditions, the strain of the side AB of the unit in the stratum can be calculated as: The ultimate strain (ε m ) of the stratum can be calculated as: where h R is the thickness of the n-th stratum, m; γ is the unit weight of the n-th stratum, kN/m 3 ; and h up is the development height of the fractures in the n-th stratum, m. The subsidence of the stratum, the length of beam and the radius of curvature satisfy the following geometrical relation: where w, the subsidence of the stratum, m, L 0, the beam length, m, ρ, of the curvature radius, rad. It can be calculated from the above that: Accordingly, the MIF height is denoted as H up (with a unit of m) in the stratum can be calculated as: where h i denote the thickness of the i-th stratum above the coal seam (i = 1, 2, 3 . . . ), m. The MR filling percentage is the ratio of the volume of fully compacted filling bodies to the volume of coal mined out, which expressed in η. It is one of the key factors controlling the development of MIF. From Equation (21), the relationship between the MIF height and the MR filling percentage can be expressed in the form of hyperbolic function. Definition coefficient λ is used to characterize the influence of MR filling percentages on MIF. The relationship between H up and η can be expressed as: where M is the thickness of the coal seam, m.

The Numerical Model and the Scheme
The universal distinct element code (UDCE) produced by Itasca Consulting Group, Inc. (Minneapolis, MN, USA) was used to analyze the characteristics of MIF under CECB. After rock failure, the parameters such as cohesion, internal friction angle, and tensile strength decrease rapidly with the increasing strain, which reflects strain-softening characteristics. However, the main characteristics of CECB include: (a) mining in the form of excavating the MRs, the exposed area of roof is smaller, and the deformation space of strata is minimized; (b) filling the MRs immediately after excavation to minimize the deformation time of strata; (c) interval excavated MRs can provide mutual lateral restraint between coal and filling to control the deformation of strata further. Although the parameters of fractured rock mass, such as cohesion are decreased, the reduction is limited in the conditions of confined deformation space and limited deformation time. In addition, CECB is put forward to solve the water conservation mining problems of low damage strata control. If there is a large deformation in the strata, this method will not be applicable.  Tables 1 and 2 [24,25,38]. the parameters of fractured rock mass, such as cohesion are decreased, the reduction is limited in the conditions of confined deformation space and limited deformation time. In addition, CECB is put forward to solve the water conservation mining problems of low damage strata control. If there is a large deformation in the strata, this method will not be applicable. Based on the characteristics and suitable conditions of CECB, a Mohr-Coulomb model with a size of 300 m × 230 m (length × height) is established to study the characteristics of MIF. The width of mining area is 120 m and 90 m boundary pillars are left on both sides of the mining area. Displacement constraints are imposed on the left, right, and bottom boundaries of the model. Surveying lines are laid at the top and bottom of each aquifer and each surveying line has 150 points, as shown in Figure 7. The parameter values of each stratum are shown in Table 1 and Table 2 [24,25,38].    The mining area is divided into four mining phases, and each mining phase contains 5 MRs, the size of which is 6 m (width) × 2.5m (height). Equivalent mining height theories hold that filling materials occupy part of overburden subsidence space, which is equivalent to reducing mining height [39,40]. The MR filling percentage can be expressed as: where M e is equivalent to the mining height, m.
In the numerical calculation, the method of a deleting grid is used to realize the different MR filling percentages. We delete the height of M e from top to bottom in the model of a coal seam, and set the rest as the elastic model. Then, we set enough elastic modulus in the elastic model to ensure that the deformation of elastic model is far less than M e in the whole mining process, so as to realize CECB in different MR filling percentages. For example, a 90% MR filling percentage is achieved by deleting 10% from the height of the coal seam model in each mining phase from top to bottom. According to the above method, the characteristics of MIF in the excavating of MRs are calculated under the conditions of 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90 MR filling percentages respectively.

The Elastic Foundation Beam Structure
The CECB can protect the main roof from breaking or caving by filling the excavated MR in time to reduce the development of MIF, and ultimately realize water conservation mining. Therefore, the higher the MR filling percentage is, the better the controlling effect of the MIF is, and the more advantageous to water conservation mining. Taking a 90% MR filling percentage as an example, the formation and evolution process of MIF under CECB is further studied. The "plot pen block open" command is used to display the distribution of MIF in each mining phase under CECB, as shown in Figure 8.
The first two mining phases have interval distribution of MIF. At the end of the first mining phase, MIF only exist above the mining area. The highest point of MIF is located in the middle of the mining area, and gradually decreases to both sides of the mining area, the shape of which is approximately triangular. MIF have discontinuous distribution because of the isolation of the coal pillars. At the end of the mining phase, MIF form a saddle shape with high sides and low middle sides over the mining area. Micro-cracks appear in aquiclude III and aquifer III. In this phase, MIF still have a discontinuous distribution. phase, MIF only exist above the mining area. The highest point of MIF is located in the middle of the mining area, and gradually decreases to both sides of the mining area, the shape of which is approximately triangular. MIF have discontinuous distribution because of the isolation of the coal pillars. At the end of the mining phase, MIF form a saddle shape with high sides and low middle sides over the mining area. Micro-cracks appear in aquiclude III and aquifer III. In this phase, MIF still have a discontinuous distribution. However, MIF in the last two mining phases are continuously distributed. At the end of the third phase, MIF are connected with each other, which include interlayer fractures formed by failure along the upper and lower surface of the grids and longitudinal fractures formed by failure along the left and right surfaces of the grids. Among them, there are only longitudinal fractures along the left and right surfaces of the grids above the coal pillar. This shows that the excavating MRs lead to the subsidence of overlying strata, and the asynchronism of the subsidence of each strata produces interlayer fractures, while the strata are also subjected to horizontal tensile stress to produce fractures. The strata above the coal pillar are mainly caused by horizontal tension stress and longitudinal fractures. At the end of the fourth phase, MIF are connected with each other, but the density of the fractures decreases compared with the third phase. After recovering the last part of the coal pillar, all the cushions bearing the main roof become the same kind of medium, and the main roof acts on the elastic foundation under a uniform force, and the horizontal tensile stress decreases. In the last phase, the longitudinal fractures are closed due to the horizontal tensile stress. During the whole CECB, the MIF height in each phase is almost unchanged, which is 7.0 m (2.8 times the mining height). Figure 8 shows the distribution patterns of MIF after CECB. By comparing and analyzing the characteristics of MIF in different MR filling percentages, it can be concluded that the displacement angle of overlying strata is about 66 degrees, and this hardly changes by the changed MR filling percentages. Combining with Figures 8d and 9, when the filling percentage (of the mined voids) is not less than 30%, there is no obvious difference between the middle and the edge of the working face, which is approximately the same height. Therefore, unlike the saddle-shaped distribution pattern of traditional longwall mining, the MIF under CECB are approximately isosceles trapezoidal distribution.

The Elastic Foundation Beam Structure
In the conditions of different MR filling percentages, the distribution patterns of MIF are a group of isosceles trapezoids with the same bottom edge, the same inclination angle of two waists and different heights. The above analysis shows that the parameters such as fracture distribution pattern and displacement angle of overlying strata are not sensitive to the change of MR filling percentages. The MR filling percentage, however, is an important index affecting the MIF height. However, MIF in the last two mining phases are continuously distributed. At the end of the third phase, MIF are connected with each other, which include interlayer fractures formed by failure along the upper and lower surface of the grids and longitudinal fractures formed by failure along the left and right surfaces of the grids. Among them, there are only longitudinal fractures along the left and right surfaces of the grids above the coal pillar. This shows that the excavating MRs lead to the subsidence of overlying strata, and the asynchronism of the subsidence of each strata produces interlayer fractures, while the strata are also subjected to horizontal tensile stress to produce fractures. The strata above the coal pillar are mainly caused by horizontal tension stress and longitudinal fractures. At the end of the fourth phase, MIF are connected with each other, but the density of the fractures decreases compared with the third phase. After recovering the last part of the coal pillar, all the cushions bearing the main roof become the same kind of medium, and the main roof acts on the elastic foundation under a uniform force, and the horizontal tensile stress decreases. In the last phase, the longitudinal fractures are closed due to the horizontal tensile stress. During the whole CECB, the MIF height in each phase is almost unchanged, which is 7.0 m (2.8 times the mining height). Figure 8 shows the distribution patterns of MIF after CECB. By comparing and analyzing the characteristics of MIF in different MR filling percentages, it can be concluded that the displacement angle of overlying strata is about 66 degrees, and this hardly changes by the changed MR filling percentages. Combining with Figure 8d and Figure 9, when the filling percentage (of the mined voids) is not less than 30%, there is no obvious difference between the middle and the edge of the working face, which is approximately the same height. Therefore, unlike the saddle-shaped distribution pattern of traditional longwall mining, the MIF under CECB are approximately isosceles trapezoidal distribution.  In the conditions of different MR filling percentages, the distribution patterns of MIF are a group of isosceles trapezoids with the same bottom edge, the same inclination angle of two waists and different heights. The above analysis shows that the parameters such as fracture distribution pattern and displacement angle of overlying strata are not sensitive to the change of MR filling percentages. The MR filling percentage, however, is an important index affecting the MIF height.

Influences of MR Filling Percentages on MIF
According to Equations (22) and (23), the effects of MR filling percentages on the MIF height can be expressed by two parameters: the influence of MR filling percentages on MIF (λ) and the ratio of the MIF height to equivalent mining height (μ). In order to display the movement of grids more intuitively, the numerical results are redisplayed by using "boundary plot" command. The curves of Hup-η, λ-η and μ-η are plotted under different MR filling percentages, respectively, as shown in Figure 10.
When the MR filling percentage decreases from 90% to 30%, the growth rate of Hup is relatively slow, from 2.8 times to 8.4 times, and the value of λ keeps around 133.1. When the MR filling percentage decreases from 30% to 10%, Hup increases from 8.4 times to 23.8 times, and the value of λ increases from 126.3 to 376.1. Therefore, the MIF height (Hup) is correlated with λ. In addition, when the MR filling percentage decreases from 90% to 30%, the trend of the increasing height of MIF is slowed down and the value of μ gradually decreases. When the MR filling percentage is less than 30%, the MIF height increases rapidly and the value of μ increases gradually. Therefore, the growth rate of the MIF height (ΔHup) is correlated with μ. The comparison between Hup-λ and ΔHup-μ in different MR filling percentages is shown in Table 3.
The Pearson correlation coefficients of Hup-λ and ΔHup-μ are calculated using the following formula: where Cov (x, y), the covariance of two sets of data, x and y; Var(x), variance of data x; and Var(y), variance of data y. The calculated results show that the Pearson correlation coefficients of Hup-λ and ΔHup-μ are 0.98 and 0.61, respectively, and the correlation grades are very strong and moderate, respectively.

Influences of MR Filling Percentages on MIF
According to Equations (22) and (23), the effects of MR filling percentages on the MIF height can be expressed by two parameters: the influence of MR filling percentages on MIF (λ) and the ratio of the MIF height to equivalent mining height (µ). In order to display the movement of grids more intuitively, the numerical results are redisplayed by using "boundary plot" command. The curves of H up -η, λ-η and µ-η are plotted under different MR filling percentages, respectively, as shown in Figure 10.  By using parameter λ as the main index, the controlling effects of MR filling percentages on MIF can be described more accurately. The Hup-η curve is divided into three ranges according to the change of Hup and λ in different MR filling percentages. They are the stability control range, the When the MR filling percentage decreases from 90% to 30%, the growth rate of H up is relatively slow, from 2.8 times to 8.4 times, and the value of λ keeps around 133.1. When the MR filling percentage decreases from 30% to 10%, H up increases from 8.4 times to 23.8 times, and the value of λ increases from 126.3 to 376.1. Therefore, the MIF height (H up ) is correlated with λ. In addition, when the MR filling percentage decreases from 90% to 30%, the trend of the increasing height of MIF is slowed down and the value of µ gradually decreases. When the MR filling percentage is less than 30%, the MIF height increases rapidly and the value of µ increases gradually. Therefore, the growth rate of the MIF height (∆H up ) is correlated with µ. The comparison between H up -λ and ∆H up -µ in different MR filling percentages is shown in Table 3. The Pearson correlation coefficients of H up -λ and ∆H up -µ are calculated using the following formula: where Cov(x, y), the covariance of two sets of data, x and y; Var(x), variance of data x; and Var(y), variance of data y. The calculated results show that the Pearson correlation coefficients of H up -λ and ∆H up -µ are 0.98 and 0.61, respectively, and the correlation grades are very strong and moderate, respectively. By using parameter λ as the main index, the controlling effects of MR filling percentages on MIF can be described more accurately. The H up -η curve is divided into three ranges according to the change of H up and λ in different MR filling percentages. They are the stability control range, the critical range and the lost control range. When the MR filling percentage is more than 50%, the MIF height increases steadily, and the overall stability of λ is maintained. Filling bodies can effectively support the strata, and the MIF height increases steadily with the decrease of the MR filling percentage, which reflects the stable control of the filling bodies on the development of MIF. When the MR filling percentage is between 30% and 50%, although the MIF height increases steadily, the amplitude of λ fluctuates and decreases. The gradual decrease of the MR filling percentage leads to the intensification of overlying strata migration, the separation between soft rock and hard rock, and the temporary slowing down the MIF development. Abscission layer is the sign of unstable state of overlying strata migration. Therefore, this range is the critical state of filling bodies' transition from stable control to lost control of MIF development. When the MR filling percentage is less than 30%, the MIF height increases sharply, which indicates the control of filling bodied is invalid.
According to the influence of MR filling percentages on MIF, the necessary condition for filling bodies to play a stable control role is that the MR filling percentage is not less than 50%. When the MR filling percentage is more than 50%, the MIF height can be expressed as follows:

Arrangement of MRs
In the study area using two wings to arrange MRs, which opens at 50 degrees on both sides of the head entry, and the width of the MR is 6 m. The MRs are divided into four mining phases. In the same mining phase, the width of the isolation structure is 18 m (that is, three MRs per interval). All the MRs are excavated in turn, and filling of the excavated MRs is done immediately. The operation mode of "parallel excavating and filling" is always maintained until the end of mining, as shown in Figure 11.
In the study area using two wings to arrange MRs, which opens at 50 degrees on both sides of the head entry, and the width of the MR is 6 m. The MRs are divided into four mining phases. In the same mining phase, the width of the isolation structure is 18 m (that is, three MRs per interval). All the MRs are excavated in turn, and filling of the excavated MRs is done immediately. The operation mode of "parallel excavating and filling" is always maintained until the end of mining, as shown in Figure 11.

The Critical MR Filling Percentage
The filling materials are high-water swelling materials. The aggregate is fly ash, and the auxiliary materials are lime, gypsum, cement, additives, etc. The filling material is mixed with water to make high-water swelling materials with a mass ratio of 1.40:1 and a density of 1.55 g/cm 3 , which is filled into the MRs. After 28-day solidification, the strength of filling body was approximately 5.0 MPa [25]. By using the condition of whether MIF cut through the aquiclude as the judgement criterion of the aquiclude failure, the critical MR filling percentage for achieving water conservation mining in the study area is determined as 83.2%.

Effects of Water Conservation
By taking into account safety coefficient and construction factors, the MR filling percentage was set as 90%. The theoretical calculation result of the MIF height is 7.55 m (3.0 times mining height) and the numerical results is 7.0m (2.8 times mining height). The above two methods calculate results that the MIF height are less than the main roof thickness, so MIF do not penetrate the main roof. However, affected by mining, water resources in aquifer III might still be lost.
In order to study the influence of water level in aquifer above MIF under the condition of CECB, a fluid-mechanical coupling model using FLAC3D was established. The size of CECB working face is 330 m × 210 m (strike × dip), and the distance between the working face and the model boundary is considered to be no less than 1.2 times of coal seam buried depth. The physical and mechanical parameters of each strata are consistent with those in Table 1. The size of the fluid-mechanical coupling model is determined as 900 m × 800 m × 230 m (length × width × height), as shown in Figure 12.

The Critical MR Filling Percentage
The filling materials are high-water swelling materials. The aggregate is fly ash, and the auxiliary materials are lime, gypsum, cement, additives, etc. The filling material is mixed with water to make high-water swelling materials with a mass ratio of 1.40:1 and a density of 1.55 g/cm 3 , which is filled into the MRs. After 28-day solidification, the strength of filling body was approximately 5.0 MPa [25]. By using the condition of whether MIF cut through the aquiclude as the judgement criterion of the aquiclude failure, the critical MR filling percentage for achieving water conservation mining in the study area is determined as 83.2%.

Effects of Water Conservation
By taking into account safety coefficient and construction factors, the MR filling percentage was set as 90%. The theoretical calculation result of the MIF height is 7.55 m (3.0 times mining height) and the numerical results is 7.0m (2.8 times mining height). The above two methods calculate results that the MIF height are less than the main roof thickness, so MIF do not penetrate the main roof. However, affected by mining, water resources in aquifer III might still be lost.
In order to study the influence of water level in aquifer above MIF under the condition of CECB, a fluid-mechanical coupling model using FLAC3D was established. The size of CECB working face is 330 m × 210 m (strike × dip), and the distance between the working face and the model boundary is considered to be no less than 1.2 times of coal seam buried depth. The physical and mechanical parameters of each strata are consistent with those in Table 1. The size of the fluid-mechanical coupling model is determined as 900 m × 800 m × 230 m (length × width × height), as shown in Figure 12. At the beginning of numerical calculation, use the code "CONFIG Fluid" to set the seepage calculation mode. The Mohr-Coulomb model is the same as the previous one. We set the water level in aquifer III, and the vertical distance to the bottom boundary of aquifer III is 3 m. The initial value and pore pressure gradient in the code "INITIAL pp" can be calculated. The default porosity of At the beginning of numerical calculation, use the code "CONFIG Fluid" to set the seepage calculation mode. The Mohr-Coulomb model is the same as the previous one. We set the water level in aquifer III, and the vertical distance to the bottom boundary of aquifer III is 3 m. The initial value and pore pressure gradient in the code "INITIAL pp" can be calculated. The default porosity of FLAC3D is 0.5. In the fluid-mechanical coupling model, the porosities of limestone and other hard rocks are slightly higher than that of mudstone and other soft rocks. In order to achieve the calculated equilibrium state, the porosities are not less than 0.2. It is divided into four phases to excavate and fill the MRs, and calculate the pore pressure in the strata. After data processing, the contour map of water level drawdown in aquifer III is obtained (Figure 13).
Energies 2019, xx, x FOR PEER REVIEW 20 of 24 Figure 13. The contour map of water level drawdown in the study area after CECB.

Discussions
(1) Advantages of CECB The background of CECB is to solve the limitation of traditional longwall backfilling mining, which concludes: (a) the exposed area of the roof in longwall mining is larger, and the overburden migration is relatively significant, which is not conducive to the support of the roof by the filling body; (b) the filling body is lack of lateral constraints after filling into the goaf, and its deformation resistance is weak; and (c) the mining activities and filling operations are mutually constrained, which not only affects the mining efficiency, but is difficult to ensure the full compaction of backfill.
The advantages of CECB are as follows: (a) mining in the form of excavating MRs, the exposed roof area is smaller than longwall mining, and the roof is supported by coal pillars or filling bodies on both sides of MRs, so as to minimize the migration degree of overburden; (b) the filling body is filled into the confined space, which is always in the state of three-dimensional stress, and its deformation is relatively small; (c) the mining activities and filling operations are independent in space and do not interfere with each other, so the efficiency of mining and filling is high; and (d) both sides of the MRs just filled are not excavated immediately, which It is determined that the maximum drawdown of aquifer III is 0.61 m, which occurs in the middle of CECB working face. The section lines are made along the middle strike and dip of the working face. The result is that on both sides of the strike direction of CECB working face, the range of water level drawdown is 130 m and 180 m, respectively, and on both sides of the dip direction, the range of water level drawdown is 100 m and 140 m, respectively. Considering the maximum influence range, the average recovery of water level is 3.39 × 10 −3 m along the strike direction for every 1m away from CECB working face, and 4.36 × 10 −3 m along the dip direction. Although the water level in the vicinity of the mining area is decreased due to mining activities, the water resources in aquifer III did not completely leak. As the distance to the CECB working face increases, the water level gradually returns to the initial state. It is considered that aquiclude III is not destroyed, and water conservation mining is realized in the study area [41,42].

Discussions
(1) Advantages of CECB The background of CECB is to solve the limitation of traditional longwall backfilling mining, which concludes: (a) the exposed area of the roof in longwall mining is larger, and the overburden migration is relatively significant, which is not conducive to the support of the roof by the filling body; (b) the filling body is lack of lateral constraints after filling into the goaf, and its deformation resistance is weak; and (c) the mining activities and filling operations are mutually constrained, which not only affects the mining efficiency, but is difficult to ensure the full compaction of backfill.
The advantages of CECB are as follows: (a) mining in the form of excavating MRs, the exposed roof area is smaller than longwall mining, and the roof is supported by coal pillars or filling bodies on both sides of MRs, so as to minimize the migration degree of overburden; (b) the filling body is filled into the confined space, which is always in the state of three-dimensional stress, and its deformation is relatively small; (c) the mining activities and filling operations are independent in space and do not interfere with each other, so the efficiency of mining and filling is high; and (d) both sides of the MRs just filled are not excavated immediately, which provide sufficient time and effective space for the solidification and full bearing of the filling body.
(2) Comparison the results by using strain-softening model and Mohr-Coulomb model To check the feasibility by using Mohr-Coulomb model and rationality of the numerical results, a strain-softening model with 90% MR filling percentage is established and studied. The general parameters of each stratum are the same as those in Mohr-Coulomb model, such as bulk modulus, shear modulus, density, internal friction angle, cohesion, and tensile strength. Then the additional parameters, such as ctable, ftable, and ttable, are determined by the plastic zone characteristics of the total stress-strain curves through uniaxial compression tests. The numerical results are shown in Figure 14. parameters, such as ctable, ftable, and ttable, are determined by the plastic zone characteristics of the total stress-strain curves through uniaxial compression tests. The numerical results are shown in Figure 14. It is determined that the density and the range of MIF obtained by the Mohr-Coulomb model and strain-softening model are similar. In strain-softening model, the height of MIF is 7.2 m (2.88 times the mining height), which is slightly higher than the 7.0m (2.8 times of mining height) calculated by using Mohr-Coulomb model. The differences between the two results are about 3%, but both are smaller than the theoretical calculation value. Therefore, it is also feasible to use the Mohr-Coulomb model to study the characteristics of MIF under CECB. In addition, the numerical results of the two models are smaller than the theoretical results. From the perspective of water conservation mining, the height of MIF is predicted by the theoretical calculation value, which can also avoid the influence caused by the selection of different numerical models. It is determined that the density and the range of MIF obtained by the Mohr-Coulomb model and strain-softening model are similar. In strain-softening model, the height of MIF is 7.2 m (2.88 times the mining height), which is slightly higher than the 7.0m (2.8 times of mining height) calculated by using Mohr-Coulomb model. The differences between the two results are about 3%, but both are smaller than the theoretical calculation value. Therefore, it is also feasible to use the Mohr-Coulomb model to study the characteristics of MIF under CECB. In addition, the numerical results of the two models are smaller than the theoretical results. From the perspective of water conservation mining, the height of MIF is predicted by the theoretical calculation value, which can also avoid the influence caused by the selection of different numerical models.

(3) Selection of mining and filling parameters
The influence of MR filling percentages on MIF is affected by factors, such as mining technology, backfilling technology, and geometrical conditions. The influencing factors mainly include the strength of the backfilling materials, the MR filling percentage, the MR height-to-width ratio, and the interval MR width. This study performed single-variable analysis for examining the effect of the filing percentage on the MIF height, during which the effects of other factors were not taken into account.
A smaller MR height-to-width ratio, more mining phases, and slighter roof subsidence can contribute to controlling the MIF height, but simultaneously reduce the mining efficiency. Higher strength of the filling materials and higher MR filling percentages can more remarkably control the MIF height; however, the utilization rate of the backfilling materials dropped. A reasonable arrangement of the MR and setting of backfilling parameters can not only achieve the function of water conservation mining in extremely close distance aquifers, but can also ensure the effectiveness of mining and increase the utilization rate of materials.
(4) Methods of improving coal mining efficiency In actual production, the MR filling percentage is generally greater than that of excavation. In order to achieve coordinated operation of mining and filling, parallel excavation of multiple MRs can be carried out, which is similar to the simultaneous cutting of multiple shearers in the same working face. On the premise of meeting the requirement of water conservation mining, filling materials are not limited to one kind. High water expansion materials, aeolian sand paste-like materials, and industrial solid waste can be selected.

(5) Advantages of promotion
Filling cost is an important issue affecting the promotion and application of CECB. One advantage of this method is that the filling of MR can be adapted to local conditions and the comprehensive benefits of coal mining can be maximized. For example, in the northwest mining area of China, filling materials are difficult to obtain in situ materials. On the premise of water conservation mining, the filling quantity of MRs can be reduced appropriately, so as to reduce the filling cost. In the eastern mining area of China, the filling materials of industrial solid waste residue mainly consisting of gangue are sufficient and cheap, and MRs can be fully filled, which not only realizes the water conservation mining in the conditions of extremely close distance aquifers, but also solves the pollution problem of industrial solid waste residue to the surface environment.

(6) Scientific and social significance
The scientific significance of this study lies in providing theoretical guidance for water conservation mining under the extremely thin barrier layers, including mining under the extremely close-distance aquifers, mining close-distance coal seams, and other geological conditions. At present, CECB has been widely used in Ordos, Jincheng, Xinwen, and Zaozhuang mining areas in China. On the premise of water conservation mining, more than 2.6 million tons of recovered coal has been increased in total in the above mining areas, which increased profits by more than 67 million dollars, and achieved good economic and environmental benefits. With a view of realizing efficient, safe, and green mining of coal resources, CECB will be gradually extended to the core coal production bases in Western China.

Conclusions
According to the different boundary conditions of gradual replacement of the main roof stress from coal pillars to filling bodies, the mechanical models of clamped-clamped beam, continuous beam, and elastic foundation beam are established. The horizontal tensile strains of the overlying The distribution patterns of MIF under CECB are studied. It is concluded that the distribution pattern is an isosceles trapezoid with the moving angle of the overlying strata as the bottom angle, and the upper and lower boundary of MIF as two parallel sides. Furthermore, the moving angle of the overlying strata does not change with MR filling percentages.
The controlling effects of MR filling percentages on MIF are studied. Taking the value of λ as the main index, the curve of the MIF height is divided into three ranges, including the stability control range, the critical range, and the lost control range. In addition, the calculating expression of the MIF height in the stability control range is given.
In the study area, the 90% MR filling percentage is used for CECB. The MIF height is 7.55 m (3.0 times mining height), and the main roof is not penetrated. The aquiclude III is effectively protected by the main roof, and its water-resisting property is not destroyed, thus realizing water conservation mining.

Conflicts of Interest:
The authors declare that they have no competing financial interests in connection with the work submitted.