Field Measurement and Study on Overburden Fracture and Surface Subsidence Law of Solid Filling Mining under Buildings

College of Mining, Liaoning Technical University, Fuxin 123000, Liaoning Province, China State Key Laboratory of Mining Disaster Prevention and Control Co-founded By Shandong Province and the Ministry of Science and Technology, Shandong University of Science and Technology, Qingdao 266590, China College of Safety and Environmental Engineering, Shandong University of Science and Technology, Qingdao 266590, China College of Energy and Mining Engineering, Shandong University of Science and Technology, Qingdao 266590, China Tiandi (Changzhou) Automation Co., Ltd., Changzhou 213000, China


Introduction
In recent years, with the development of coal mining science and technology, filling mining has become one of the important technical ways to solve the problem of coal pressure under buildings, water bodies, railways, and confined water bodies [1][2][3][4][5]. Compared with the traditional caving method, the overburden of filling mining can maintain a certain degree of integrity and continuity and will not trigger greater deformation and destruction, thus reducing and slowing down the surface subsidence over stope [6][7][8][9]. e roof can be controlled in filling mining by filling the mining area with solid waste such as gangue, which can not only effectively control surface subsidence and improve the coal recovery rate but also rationally utilize solid waste [10][11][12]. e waste is not transported and stored, which is beneficial to reduce the cost and obviously reduce the surface deformation, thus avoiding the impact on the surface environment [13][14][15].
In order to study the law of rock pressure appearance and the law of overlying strata change under the condition of solid waste filling mining in mines, the surface movement and deformation of backfilling mining are predicted and analyzed through the field measurement of surface subsidence in backfilling mining, and then, the control effect of backfilling mining on overlying strata and surface is evaluated [16][17][18][19]. By making full use of the existing engineering conditions of a mine in the eastern mining area of China, the monitoring system of the mining hydraulic support system was established to analyze the variation law of working resistance of filling hydraulic support with the advancing of working face. e monitoring system of surrounding rock deformation and stress was established to study the law of rock pressure appearing in two roadways. e real-time monitoring system of borehole peeping for the stress in the filling body of gob and the crack development in the roof of gob is established to get the movement law of overlying strata after filling with solid waste in gob. At the same time, the observation station of surface subsidence is arranged to analyze the measured data of surface subsidence [20,21], so as to provide data support for the prediction of surface movement and deformation model.

Coal Seam and Roof and Floor Conditions.
e coal seam of 8103 filling face is 8# coal seam, which belongs to complex structure coal seam. e inclination angle of coal seam is 7-10°, the average inclination angle is 8.5°, the bulk density of coal is 1.49 t/m3, the average thickness of coal seam is 3.7 m, the mining face elevation is from -730m to -680m, and the corresponding surface elevation is about +19m. Coal seam contains 1-2 layers of intercalated rock with thickness of 0.1-0.2 m, mainly carbonaceous mudstone and mudstone. e roof and floor of 8103 filling face are shown in Figure 2.

Mining Face Parameters.
e reserves of 8103 filling face are about 220,000 t, the mining strike is 128m long, the inclined mining width is 120m, and the recovery rate is 95%. e actual mining-filling ratio of 8103 working face is 1 : 1, which means 1 ton of coal is mined with 1 ton of gangue filled. is working face adopts a new type of comprehensive mechanized gangue filling mining system, which can ensure full contact between the filling body and roof, that is, filling underjointed roof quantity h q � 0. According to the compression test of waste rock in laboratory, the compression ratio of waste filling body is η � 15%.

Stope Pressure Monitoring
In order to reveal the characteristics and regularities of rock pressure in filling face and mining roadway and the movement law of overburden strata in stope, on-site monitoring of rock pressure in 8103 filling face during mining process is carried out. e contents include stress monitoring of mining hydraulic support system, stress monitoring of surrounding rock affected by mining, stress monitoring of the filling body in gob, and borehole peek monitoring of roof crack development in overlying strata of gob.

Monitoring and Analysis of Support of Hydraulic Support in Coal Mining
3.1.1. Layout Scheme of Monitoring Equipment. In order to study the variation of working resistance of hydraulic support in 8103 filling face with the advance of the working face, a monitoring instrument of working resistance is arranged along the hydraulic support in the working face, as shown in Figure 3.

Monitoring Results and Analysis.
rough the analysis of monitoring results of working resistance of backfilling support in 8103 working face, the measured curve of working resistance distribution is shown in Figures 4-7.
In the solid filling mining process of 8103 working face, the average initial support force is between 6600 KN and 7000 KN, which is little different from the average working resistance of the support 7200 KN to 7800 KN. When the working face is pushed to 15m, 35m, and 60m positions, respectively, the average working resistance of the support increases slightly to 7400 KN, 7700 KN, and 7500 KN, which are not up to the rated working resistance of the hydraulic support 8800 KN, and there is no sharp increase of working resistance, which indicates that the working resistance of the support increases sharply. ere are no initial and periodic weighting phenomena in the mining process of filling coal face. e hydraulic support of filling coal mining effectively controls the roof subsidence before filling.

Monitoring and Analysis of Surrounding Rock
Deformation in the Working Face

Layout Scheme of Monitoring Equipment in Roadways of the Working Face.
e layout equipment of two roadways in the working face is mainly to observe the rock pressure appearance, advance abutment pressure, and surrounding rock deformation of the two roadways.
ere are three stations in the two lanes, each of which is 10m in length, as shown in Figure 8.
Layout of mine pressure monitoring instruments in two roadways: in order to observe the movement of roof and floor and the approximation of two sides of roadway in 8103 filling face, the longitudinal force of bolt (cable), and the displacement of deep surrounding rock, six monitoring points of roadway surface displacement, twelve multipoint displacement meters, and twelve bolt (cable) stress monitors are arranged in the two roadways, which are distributed evenly in each station to monitor various data.
Layout of monitoring equipment for advance abutment pressure: in order to monitor the advance abutment pressure of two roadways in the mining process of the working face, six supporting pillars are installed in each roadway of the working face, and pressure monitors are arranged. e advanced pressure monitor is arranged at the positions of 20m, 80m, and 140m from each roadway, respectively. Before the start of mining, six advanced pressure monitors are installed, which will not move along with the advance of the working face and will be removed when the working face reaches the position where they were installed.

Monitoring Results and Analysis.
Based on the analysis of the actual measurement of the advance stress of the coal body in two roadways of 8103 working face, the measured curve of the advance monomer working resistance in the gangue haulage roadway is obtained, as shown in Figure 9.
e measured curve of the advance monomer working resistance in the coal haulage roadway, the measured curve of the anchor force in the coal haulage roadway, and the measured curve of the anchor force in the gangue haulage roadway are similar to those shown in Figure 9.
Based on the analysis of the measured deformation of surrounding rock in two roadways of 8103 working face, the measured curve of surrounding rock deformation in coal haulage roadway of 8103 filling face is obtained, as shown in Figure 10. At about 18m away from the working face, the closeness of the roof and the floor and the closeness speed of two sides of rock mass are the fastest, the closeness speed of the roof and the floor is about 5.5 m/d, and the closeness speed of rock mass is about 6.4 mm/d; the closer to the working face it is, the faster the displacement of the roof and floor is. e larger the proximity is, the smaller the change is. It can be divided into slow growth stage, rapid growth stage, and stable stage. e change trend of working resistance of advanced monomer in gangue roadway is basically the same. In the same position, the displacement of the roof and the floor is greater than that of two sides of rock mass.

Layout Scheme of Monitoring Equipment in the Filling
Body. In order to monitor the pressure change of the filling body in 8103 filling face under the action of roof, 25 filling body stress monitors were installed in the filling body. e instrument is arranged in five rows, five in each row and symmetrically arranged. e distance from the roadway side of 8103 coal haulage roadway is 5 m, 30 m, and 60 m, respectively. When the self-cutting eye is pushed to 10m, the first row of equipment is installed. When the working face is pushed to 30m, 50m, 80m, and 100m, the other rows of equipment are installed, respectively, as shown in Figure 11.

Analysis of Stress Change in the Filling
Body. e stress monitoring and analysis of the filling body in the working face is carried out by using the stress monitor. e change of the stress with the working advancing distance is shown in Figure 12.  From the analysis, we can see that, with the advancing of the working face, the stress variation law of the filling body can be divided into three stages: initial stress zone (0-15m from working face), stress increasing zone (15-33m), and stress stable zone (after 33m), which are consistent with the variation law of roof subsidence with advancing distance of the working face.
Within the initial stress range of 0-15m, the internal stress of the filling body varies between 1.6 and 2 MPa, and there is almost no stress increase, which indicates that the roof bending subsidence in the initial stress area after filling is small.
In the stress increasing area of 15-33m, the stress inside the filling body increases gradually, but the growth rate is        Shock and Vibration 5 slow, which indicates that the roof of the goaf is gradually bent and subsided under the action of overlying strata, and the filling body is gradually compacted. After 33m away from the working face, with the advancing of the working face, the internal stress of the filling body gradually stabilizes, with the maximum value reaching 12 MPa, which is close to the original rock stress, indicating that the overall bending subsidence of the overburden strata has become stable under the support of the filling body at this time.

Shock and Vibration
According to the results of borehole detection, with the advancing of working face, the overburden strata rupture in gob develops spatially from shallow to deep; the shallow is mainly in the form of rupture, and the deep is mainly in the form of fissure. In terms of time, due to the lag of filling operation and overburden movement, deep overburden fracture lags behind mining process in time. During the process of filling mining, the overburden strata are mainly bent and subsided, and there is no obvious separation between the direct roof and the basic roof. In the vertical direction, the direct roof within 2m distance from the coal seam appears to be fractured and cracks are developed, but the integrity is good, and the basic roof does not appear to be fractured. is indicates that, in the process of solid filling  mining, the direct roof is locally fractured. ere was no collapse zone, but the main roof was bending and sinking, and no breakage occurred. e supporting resistance of hydraulic support in backfilling mining of 8103 working face, the surrounding rock of working face, the change of internal stress in backfilling body, and the crack development of roof are monitored and analyzed. Comprehensive analysis shows that the measured data analysis of 8103 filling face is consistent with the change rule of roof subsidence with the advance of the working face. ere is no phenomenon of initial and periodic weighting in the filling mining process. e roof bends and subsides, but no fracture occurs. Monitoring and analysis of roof fracture development characteristics proves that, in the filling process, the direct roof breaks locally without the caving zone, while the basic roof mainly bends and subsides without breaking. e filling body in gob plays an effective role in controlling the movement of overlying strata. With the advancing of the working face, the internal stress of the filling body is gradually stable, which is close to the original rock stress. At this time, the overall bending and subsidence of the overlying strata has become stable under the support of the filling body.

Analysis of Surface Subsidence in Filling
Mining

Analysis of Surface Subsidence Data.
e subsidence curves along lines 1 to 3 are obtained by processing the observed data of lines 1, 2, and 3, as shown in Figures 18-20. From the graph, it can be seen that the subsidence values of observation points N10-N15 and W4-W7 are on the high side of 8103 working face, and the cumulative maximum subsidence value of observation point N15 is -38 mm. e actual measurement shows that the mining of 8103 working face has little influence on the surface.

Prediction Methods and Determination of Parameters.
e probability integral method is used to predict surface subsidence in 8103 working face. When the probability integral method based on equivalent mining thickness is used to predict the surface subsidence of solid filling mining such as gangue, the corresponding parameters of caving mining in this mining area are selected. According to the regularity and parameters of strata and surface movement, the prediction parameters of mining subsidence are listed in Table 1. e filling mining of 8103 working face is a new type of comprehensive mechanized gangue filling mining system, which can ensure the full contact between the gangue filling body and the roof, that is, the underfilling roof quantity hq � 0. In filling mining, the actual mining height minus the height of filling material is equivalent mining height [22], and then, the movement and deformation characteristics of overburden rock in stope can be analyzed by the traditional theory of strata control in nonfilling mining. According to the equivalent mining height theory, the equivalent mining thickness HD of solid filling mining such as gangue can be expressed as where L denotes roof control distance of the filling face, m, L 0 denotes limit span of the roof, m, h m denotes actual mining thickness of coal seam, m, and η is compressibility of backfill.
According to the compression test of waste rock in laboratory, η, compressibility of backfill, is chosen to be equal to 15%. According to the geological and mining data of the strata in the mining area, it is known that the maximum roof control distance L � 9.3 m, limit span of roof L 0 ＝30m, and actual mining thickness of coal seam h m ＝3.5 m. From

Prediction and Analysis of Surface Subsidence.
is paper calculates the surface subsidence caused by filling mining in 8103 working face with the method of outline integral and obtains the maximum predicted value of surface movement and deformation, as listed in Table 2. e prediction contours of surface subsidence, horizontal movement, and horizontal deformation along strike and dip directions in filling mining of 8103 working face are obtained, as shown in  Since the compaction process of filling material in 8103 working face is gradual, the stability time of surface subsidence in filling mining is generally 3-6 months, and the prediction result of surface subsidence is the final subsidence value, so the prediction value of surface subsidence is larger than the observation value. For example, the predicted maximum subsidence (65 mm) is larger than the observed maximum subsidence (38 mm).
After filling mining in 8103 working face, the predicted value of maximum subsidence is 65 mm, the predicted value of maximum inclination is 0.2 mm/m, and the predicted value of maximum horizontal deformation is 0.32 mm/m, which has little impact on the surface and buildings.

Conclusions
is study reveals the characteristics of rock pressure in solid filling mining face. e average working resistance of the support does not reach the rated working resistance of the hydraulic support of 8800 KN, and there is no sharp increase in working resistance. e stress variation law of filling body can be divided into three stages: initial stress zone, stress increasing zone, and stress stable zone without stress surge.
is shows that there are no initial and periodic weighting phenomena in the mining process of filling coal face.
It reveals the evolution characteristics of overburden fracture in solid filling mining. e overburden morphology in solid filling mining is mainly characterized by the bending zone and the fracture zone with certain height, without the caving zone. In the process of solid filling mining, the direct roof is broken locally without the caving zone, while the basic roof is mainly bending subsidence without breaking. e research shows that the probabilistic integral method model based on equivalent mining thickness is unified, and the surface subsidence of 8103 working face is effectively predicted by the probabilistic integral method. In order to make full use of the observation results of surface movement in caving mining, it is suggested that the probabilistic integral model should be used to predict surface subsidence in solid filling mining. Furthermore, a method of selecting parameters system of the probability integral model based on equivalent mining height is proposed.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.