Abstract

Rock-like materials can be used as a filling material to improve the stability of underground engineering substantially. And the optimization of the rock-like material is an effective way to improve the performance of the filling material. Firstly, the AHP-FUZZY comprehensive optimization model was constructed for the complexity and fuzziness of the rock-like filling material ratio optimization. Secondly, according to the mechanical properties of rock-like filling materials under special environment (such as high temperature, high humidity, high stress, and high airtight) in old goaf, mechanical properties of rock-like filling materials were studied with the method of field core and laboratory test, which revealed the variation law of mechanical properties with time, and the regression equation between mechanical parameters and time was established with the method of least squares. Finally, the strength and deformation characteristics of rock-like filling materials were monitored by the monitoring and early warning technology. The results show that the optimal ratio of the rock-like filling materials is “E3”; that is, the cement content is 9%, the ash ratio is 2 : 5, and the mass concentration is 74%. The mass concentration is the main factor that affects the slump of the slurry, and the proportion of fly ash and coal gangue content directly affects the stratification and bleeding rate of the slurry. Reasonably increasing the ash and gangue ratio can significantly improve the workability and water retention of the rock-like filling materials. Also, the amount of composite cementitious material is the main factor that affects the setting time and the strength of the rock-like material. What is more, the special environment in gob is good for each chemical reaction fully in rock-like filling materials and strengthens the gelling property of deformation resisting capability, which can be a benefit for disaster prevention of underground engineering.

1. Introduction

After a decade of golden development, the coal industry has gradually entered a difficult period from 2012. On the one hand, extensive coal mining has caused destructive damage to the mine environment, such as large-scale surface subsidence, increasingly serious environmental pollution, and frequent occurrence of haze [1]. On the other hand, the development and utilization technology of various new clean energy, such as wind energy, solar energy, and nuclear energy, has led to overcapacity and difficult operation of coal mines. How to find a way of coal resources development and utilization and sustainable development of ecological environment is the focus of domestic and foreign scholars [2].

As early as the 1950s, Canada began to develop cemented filling. In China, cemented filling mining began in the mid-1960s, mainly for metal ore filling. In the 1960s, filling mining began to develop in China [3]. As a key technology of green coal mining, it can effectively control the surface subsidence and improve the ecological environment of the mining area, which has been the research hotspot of domestic mining scholars [4]. Filling mining is mainly used in noncoal mines at first and has been widely used in coal mines in recent years. In the 1980s and 1990s, new technologies such as paste-like filling, rock-like filling, waste rock cemented filling, full-tailings cemented filling, and high water filling have been developed rapidly in coal mines [5]. The mechanical properties of rock-like filling body are the core of rock-like filling technology, and it is the key to effectively control surface deformation and prevent surface subsidence [6]. In 2003, Qian Ming, a senior academician [7], put forward a green mining technology system for underground coal mining, which laid a foundation for the development of green mining technology in China's coal mines. The optimization of the filling material ratio is an effective way to improve the performance of filling material and an important guarantee for controlling surface subsidence.

Mu et al. conducted orthogonal experiments on the preparation of mine filling materials with calcium carbide slag, blast furnace slag, fly ash, and flue gas desulfurization gypsum [8]. Gardner et al. studied the properties of magnesium potassium phosphate cement mixed with fly ash and ground blast furnace slag [9]. Du et al. studied the influence of fly ash characteristics on geopolymer properties [10]. Fall et al. contributed to understanding the effect of curing temperature on the mechanical properties of the colloidal end sand filling of mines [11]. Hao et al. studied the surface modification of fly ash by calcium carbide slag and its effect on the compressive strength and autoshrinkage of the blended cement paste [12]. Zhou et al. studied the diffusion behavior of respirable dust in fully mechanized coal mining face: CFD numerical simulation and engineering application [13]. Liu et al. did a study on the effects of the installation position of a multiradial swirling air-curtain generator on dust diffusion and pollution rules in a fully mechanized excavation face [14]. Zhang et al. studied and practiced the local filling technology of fly ash high water material short wall [15]. Zhang et al. studied the bearing capacity and roof stability of the filling body during the strip coal pillar paste recovery process [6]. Sobhi et al. conducted a numerical study on the earth pressure coefficient along the centerline of the backfill stope [16]. Zhang et al. studied the hydraulic characteristics of coal gangue-fly ash paste-filling and its application in underground coal mines [17].

The stability of the filling system is an important guarantee for the successful implementation of the filling mining technology and is the key to effectively controlling the surface deformation and preventing surface collapse. The filling system is a complex system, affected by multilevel and multidimensional uncertain factors. At present, scholars at home and abroad have conducted a lot of exploration and research on the stability evaluation of the filling system.

Macías et al. proposed anomalous metal-rich phosphogypsum: characterization and classification according to international regulations [18]. de Rezende et al. presented a laboratory study on phosphogypsum, stabilizers, and tropical soil mixtures [19]. Dong et al. used multiple indicators to comprehensively consider the uncertain factors of the mine to quantitatively evaluate the risk degree of underground mined-out areas and conduct case studies [20]. Wang et al. put forward a multiplane detection optimization algorithm for the design of the interval charge structure of large-diameter deep-hole blasting based on rock crushing [21]. Dong et al. indicated a quantitative evaluation and case study of multi-index clean mining considering uncertain factors for phosphate rock [22]. Sun et al. designed X-ray CT three-dimensional reconstruction and discrete element analysis of the pore structure of cement paste under uniaxial compression [23].

It can be seen that most of the domestic research on rock-like filling materials and mechanical properties is based on laboratory proportioning to prepare filling material test blocks and carry out routine mechanical property tests on them under laboratory curing conditions. Most of them have no field sampling to study the mechanical properties of the rock-like filling body in a specific underground environment (high temperature, high humidity, high stress, and high sealing), which leads to the measured data. It is not very consistent with the field measurement, and cannot reflect the actual stress and strain of the rock-like filling body. On the basis of previous research, the author explores the combination of the traditional orthogonal test method and hierarchical fuzzy mathematics method, and the combination of qualitative and quantitative methods, which not only avoids the error of human factors in the traditional filling material test but also overcomes the lack of a large amount of data in traditional mathematical statistics method. At the same time, combining with the experience and opinions of experts, the proportion of filling materials is optimized. Through field coring, the mechanical properties of the rock-like filling body and its relationship with time and the long-term stability of goaf filling body and overlying strata structure under mining disturbance are studied. It provides a reliable basis for guiding the field practice.

2. Experimental Study on Material Properties of Rock-Like Filling Materials in Coal Mine

2.1. Test Materials and Preparation

Coal gangue is selected as coarse aggregate, fly ash as fine aggregate, and 325 ordinary Portland cement as cementing material. The size of the test specimen is 150 mm × 150 mm × 150 mm. The setting time, stratification, slump, and other parameters of the filling slurry were measured during the sample making process. After the specimen was made, the curing box was placed. The temperature was set at about 22°C and the humidity was adjusted to 90%, as shown in Figures 1 and 2. Then, core and grind the large module with test sample processing equipment (grinding machine (AHM-200); cutting machine (DQ-4)). As the standard rock sample is cylindrical with a diameter of 50 mm 100 mm, a 50 mm drill bit (inner diameter of 50 mm) is used. Finally, the uniaxial compressive strength of the samples at different ages was measured by the pressure testing machine, and the experimental results were listed in the L9 (33) orthogonal analysis table, as shown in Table 1.

2.2. Experimental Scheme Design

Three factors and three levels of orthogonal table were used to design the orthogonal experiment scheme. Three factors are selected in this filling material ratio test: cement quality/mixture quality is recorded as cement content a; fly ash mass/coal gangue mass is recorded as ash gangue ratio B; and slurry mass concentration C. Three levels of factor a cement content are selected, which are A(a): 3%, A(a): 9%, A(c): 16%; factor B: ash gangue ratio: B (a): 5 : 2, B (b): 1 : 1, B (c): 2 : 5; and factor C: slurry mass concentration, C (a): 74%, C (b): 63%, C (c): 68%; the test scheme is designed according to the design principle of rock-like filling materials of Jinchuan company. The relevant indexes meeting the requirements of pumpability, fluidity, and mechanical properties of filling paste were selected, including slump D1, setting time D2, stratification D3, bleeding rate D4, and compressive strength (28 d) D5. Through range analysis and visual analysis, the average of each test factor was compared, and the corresponding indexes D1, D2, D3, D4, and D5 were obtained. The optimal proportion combination of rock-like filling materials was A (b), B (a), and C (b); A (c), B (b), and C (c); A (c), B (c), and C (c); A (c), B (c), and C (c); A (c), B (c), and C (1). Therefore, according to each index D1 ∼ D5, the L9 (33) orthogonal test scheme of nine groups of rock-like filling materials in Table 1 was optimized, and finally five groups of optimal proportion combinations were determined, which were, respectively, recorded as schemes E1, E2, E3, E4, and E5. It can be seen that, corresponding to indexes D3 and D4, the optimal mix proportion combination of rock-like filling materials, namely, E3 and E4 test schemes, are all A (c), B (c), and C (a). Since the range values of factors A (b) and A (c) corresponding to index D3 are basically the same, so the optimal ratio combinations of E3 are A (b), B (c), and C (a), as shown in Table 2.

2.3. Construction of Optimization Index System of Filling Material Ratio

According to the practice of mine filling and the previous research in literature [2224], the main factors affecting the optimization of filling material ratio can be divided into internal and external factors. The internal factors include the pumpability and fluidity of the filling material and the strength and mechanical stability of the filling body. The external factors include the stability of the filling system and the economy of filling materials. Therefore, based on APH theory, a three-level optimization index system of filling material ratio is constructed, which takes the optimization of filling material ratio as the target layer, four factors as the criterion layer, four factors as the criterion layer, and seventeen factors as the index layer, such as the strength of filling body and the compactness of filling body, to form a three-level optimization index system of filling material proportion, as shown in Figure 3.

Stability of filling body: the stability of the filling body is directly related to the filling quality, that is, whether the filling body can effectively support the overlying strata and control the surface subsidence. The strength of the filling body includes initial strength, namely, self-supporting strength and long-term strength [25]. It is found that the content of cementitious material is the main factor affecting the strength of filling slurry, and the strength of filling slurry can be effectively improved by increasing the content of cementitious material. The compactness and cementation of the backfill are also the main factors affecting the stability [26]. The better the compactness and cementation are, the less likely the backfill is to be compressed and more stable. In addition, the top connection rate of the filling body will also affect the stability of the filling body [27]. The higher the top connection rate is, the earlier the deformation and movement of surrounding rock can be limited, and the support system of roof filling body floor can be formed, so as to improve its own stability, effectively control the movement and deformation of surrounding rock, and reduce the surface subsidence.

Pumpability of filling materials: the pumpability of filling materials is the main index to measure the performance of filling materials. It refers to that the filling materials are pumped into the goaf for filling through filling pumps, filling pipelines, and sometimes expressed by mobility. The pumpability or fluidity of filling materials will directly affect the performance of filling materials and even cause the plugging of filling pipelines and paralyze the whole filling system. The flow properties of filling materials are directly affected by filling mass concentration, setting time, stratification, workability, and bleeding. It is found that [2831], the mass concentration is the main factor affecting the workability of filling materials, and the specific gravity of fly ash and coal gangue directly affects the layering degree and bleeding rate of filling materials. The amount of cement is the main factor affecting the setting time and strength of the filling material. Therefore, the pumpability of filling materials can be improved by adjusting the mass concentration of filling materials and the content of fly ash.

Stability of filling system: the complete filling system consists of five parts: the supply system of filling raw materials, the preparation system of filling slurry, the filling pumping system, the filling pipeline system, and the filling stope system. The proportion and performance of filling materials are related to the composition and stability of the filling system and then affect the filling efficiency [32]. The stability of the filling system mainly includes the working environment of stope, stability of filling pipeline, stability of filling equipment and facilities, occupational health, and hygiene of personnel.

Economy of filling materials: the economy of filling materials is an important condition for optimizing the proportion of filling materials and an important part of filling mining costs, including the transportation economy of filling materials, the economy of cementing materials and aggregates, the economy of supporting filling systems, and the recycling of filling materials [33].

2.4. Determination of Index Weight for Optimization of Filling Material Ratio

Based on the analysis of the factors influencing the optimization of the filling material ratio, the relative importance of the stability of the filling body C1, the pumpability of the filling material C2, the stability of the filling system C3, and the economic efficiency of the filling material C4 are analyzed by the general 1–9 scaling method [34], and the judgment matrix is obtained through the investigation and analysis of the actual situation of the site. Finally, the expert's opinions were obtained. The specific steps are as follows .(1)The general 1–9 scale method is used to determine the judgment matrix A, as shown in Table 3.(2)The weight vector of influencing factors of the stability of the filling body, the pumpability of the filling material, the economy of the filling material, and the stability of the filling system on the optimization of the filling material proportion of the target layer are obtained by the summation method (approximate method). The column vectors of the judgment matrix are normalized, then the row vector is summed to get the weight matrix, and after that the column vector is normalized. The result is the criterion layer weight of each factor.(3)Calculate the maximum eigenvalue of judgment matrix: it is solved by.(4)The maximum eigenvalue is obtained as follows:(5)The consistency index is obtained as follows:(6)Since the judgment matrix is a fourth-order matrix, the random consistency index table [35] shows that(7)Consistency ratio:It is obvious that the consistency is satisfied.The weight vector of other index layers can be obtained by the same method, as shown in Table 1.

2.5. Fuzzy Comprehensive Optimization of Filling Material Ratio
2.5.1. The Optimization Factor Set and Weight Set of Filling Material Ratio Were Established

According to the index system Figure 3 of rock-like filling material ratio optimization, the influence factors of filling material proportion optimization are divided into various factors, and the factor sets are established, respectively. The factor set of criterion layer includes the stability of the filling body, the pumpability of filling material, the stability of the filling system, and the economy of filling material.

2.5.2. Construction of Fuzzy Evaluation Matrix

By using the Delphi method, the relevant persons in charge of filling mining in coal mining enterprises, construction directors of labor engineering companies, and relevant experts from universities and industries are invited to evaluate the index system for the optimization of filling material proportioning, combining with the practice of mine filling mining, and the influence of each criterion layer and index layer factor on the optimization of filling material ratio under different test schemes (E1 ∼ E5) is obtained at the sound level, as shown in Table 4.

For example, R1: for the index filling body strength, 4 out of 40 people think that the strength of filling body in scheme E1 has a great influence on the optimization of filling material proportion, and 2 people think that the strength of filling body in scheme E2 has a great influence on the optimization of filling material proportion (the same as others), as shown in Figures 4(a)–4(d).

2.5.3. First-Order Fuzzy Comprehensive Optimization

Operator is used to make a fuzzy comprehensive evaluation on the optimization scheme of coal mine filling material proportion, and the evaluation result of filling body stability factor index is A1. In the same way, the evaluation result vector of other factors affecting the optimization of coal mine filling material ratio can be obtained as B:

Therefore, according to the principle of maximum membership degree, the fuzzy comprehensive optimization scheme of gangue gypsum filling material proportion is ranked as follows: E3 > E2 > E4 > E5 > E1. The optimal proportioning scheme is “E3.” They are A (b), B (c), and C (a). The cement content is 9%, the ratio of ash to gangue is 2 : 5, and the mass concentration is 74%. It can be seen that the sensitivity of “ash to gangue ratio” is the highest in the process of filling material proportion optimization, followed by “cement content” and “mass concentration,” and the two are basically the same sensitivity to the optimization of filling material proportion. Therefore, the proportion of fly ash and gangue, that is, the proportion of fly ash and coal gangue, should be paid more attention to in the optimization of filling material proportion. The test results show that a proper amount of fly ash has a positive effect on improving the pumpability of filling material, but excessive fly ash will affect the performance of filling material.

3. Experimental Study on Mechanical Properties of Rock-Like Filling Materials in Coal Mine

3.1. Site Coring Scheme

In order to study the mechanical properties of backfill behind the wall, the paste samples were drilled at 2301 working face in Daizhuang coal mine of Zibo Mining Group. The filling body with filling time t ≤ 30d is selected as the research object. According to the filling technology of 2301 working face in Daizhuang coal mine, a cycle of three days is conducted, and the drilling depth should be 40 m. The drilling depth of No. 3 hole in this sampling is 50 m, which is located in the middle of the working face; the drilling depth of No. 4 hole is 42 m, which is 20 m away from the track trough of 2301 working face, and the distance between No. 3 hole and No. 3 hole is 27 m. The positions of hole 3 and hole 4 on the 2301 working face are shown in Figure 5.

3.2. Preparation of Standard Specimen

This time, a 50 mm drill bit was used for coring, which was sealed after coring on-site, so as to prevent weathering and keep the original as much as possible. The drilling depth is 4 m, 8 m, 12 m, ..., 44m (i.e., the age is 1 d, 4 d, 7 d, …, 31 d); the cores were processed again in the laboratory. The processed specimen is a neat cylinder, the ratio of height to diameter is 2.5–3.0, the diameter is about 50 mm, the end face of the specimen is flat to 0.02 mm, and the perpendicularity of the specimen axis is not more than 0.05 mm, which can meet the standard of rock test specification. According to the different ages (1 d, 4 d, 7 d, …, 31 d) of the specimens made from hole 3 and hole 4, respectively, they were divided into 11 groups and numbered as Ai (hole 3) and Bi (hole 4), i = 1,2 11 [36]. Six standard specimens were selected for each group, numbered Ai − j (hole 3) and Bi − j (hole 4), i = 1, 2 11, j = 1, 2, …, 6. Finally, three standard specimens were selected from each group for uniaxial tensile test, and the other three specimens for uniaxial compression test. The average value of each group was taken as the test results. In order to reduce the influence of drilling position on the mechanical property test of rock-like filling body and reduce the experimental error as much as possible, the average value of the results obtained from the samples of No. 3 hole and No. 4 hole in the same group is taken. Some standard samples obtained in this test are shown in Figure 6.

3.3. Mechanical Properties Analysis of Rock-Like Filling Body
3.3.1. Relationship between Tensile Strength and Time of Backfill

By using the least square method, the regression curves 1 and 2 of the tensile strength of the filling body with time are obtained, as shown in Figure 7.

It can be seen from Figure 7 that the tensile strength of the filling body increases with time. The change of tensile strength of the filling body with time can be divided into three stages: the first stage (1∼18 d) is a rapid growth period, the tensile strength of filling body increases rapidly with time, and the increase range is 86%, the compressive strength of filling body reaches 89.8% of the final strength after 18 d; the second stage (18∼30 d), which is a slow growth period, increases slowly with time; the third stage (≥30 d) is a stable period. With the increase of time, the compressive strength of backfill almost does not increase, reaching the ultimate tensile strength of 1.198 MPa.

3.3.2. Relationship between Uniaxial Compressive Strength and Time of Backfill

By using the least square method, the regression curves of the uniaxial compressive strength of the filling body with time are obtained, as shown in Figure 8. It can be seen from the figure that, overall, with the increase of time, the uniaxial compressive strength of the filling body gradually increases, the increase range decreases from large to small, and finally tends to be stable, and the compressive strength of the filling body is almost unchanged. Locally, with the increase of time, the uniaxial compressive strength of the filling body is almost unchanged. The uniaxial compressive strength of the backfill specimens showed a leap forward increase. There are two peaks in the process of “leap forward.” The first peak appears at the age of 6 days, which is 3.06 times the compressive strength of 1D filling body. The second peak appears at the age of 18 days, which is 4.25 times of the compressive strength of 1D filling body. It shows that when the age is 18 days, the chemical reaction of various compounds in the filling body is the strongest, forming C-S-H gel and Ca(OH)2, C3AH6, AFt, and other products. These products have strong gelling properties and in a certain way rely on a variety of gravitational links to form a cement stone structure to produce greater strength. The uniaxial compressive strength of the filling body tends to be stable after 21 days, and the compressive strength of the filling body is basically stable at 6.0 MPa. This is because the internal chemical action and mechanical properties of the filling body tend to be stable, and the effect of various factors on the strength of the filling body will also tend to be stable.

It can be seen from Figure 9 that the uniaxial compressive strength of the 28d test piece is 4.7 MPa under the optimal proportioning conditions of coal gangue gypsum filling material and laboratory curing conditions (curing temperature is 21 ± 1°C, humidity is 83%), and the compressive strength of filling body after 21 days measured by field coring is 6.0 MPa, which is 1.3 times of the former. This shows that the specific environmental conditions in coal mine, namely, high temperature and humidity and confined space (the actual measured internal temperature of filling paste is about 50°C), are conducive to the full reaction of various chemical substances in the filling body, so that the cementitious property of the filling body is enhanced, so that it is stable and has high strength before 28 days.

4. Field Monitoring of Long-Term Stability of Rock-Like Filling Materials

The online monitoring system for the long-term stability of the rock-like filling body in the coal mine is mainly composed of monitoring server, communication interface, underground monitoring master station, multifunctional monitoring substation, mine data optical terminal, and supporting explosion-proof system. The well monitoring server uses an industrial PC to extend RDS or optical fiber interface. Office LAN users access the well monitoring server through telephone line, optical fiber, or/and Ethernet ring network, and the monitoring information of the well monitoring server is shared by client mode or IE mode. The underground monitoring master station is responsible for patrolling each underground monitoring substation, and the patrolling data is transmitted to the underground monitoring server. The system adopts a multilevel bus distributed structure and intrinsically safe design. The maximum monitoring capacity of the system can reach 1000 measuring points. The monitoring server can access the mine LAN through OPC or FTP. The monitoring server software adopts C/S + B/S structure. The users of office LAN can share monitoring information through client mode or IE mode. The underground monitoring master station has built-in RDS and NPORT interfaces. Users can choose a telephone line, optical fiber, and ring network to connect with the underground monitoring server. The underground monitoring substation has the functions of data display, bus communication, and alarm. The substation is connected with the substation through the RS485 bus. The monitoring system is shown in Figure 10.

4.1. Subsidence Monitoring System of Filling Body Top and Bottom Plate

The special KBU101-200 roof subsidence position state instrument is used to measure the roof and floor subsidence of goaf. The measuring range is extended to 500 mm. The gear rack structure is used to convert the linear displacement into angular displacement. The displacement signal is converted into a voltage signal through the internal angular displacement sensor. The digital signal is collected by a single chip microcomputer. The monitor has a built-in 485 communication interface to transmit the data to RS48 5 bus.

KBU101 roof subsidence position instrument adopts elastic energy storage piston rod displacement structure, which can make the sensor waterproof and sealed, and the position instrument has zero calibration function. The length of the instrument itself is 1.7 m, and the measuring range of the extension rod can be up to 3.0 m. If the height to be measured is larger, the supporting rod needs to be added. The simple way is to fix the metal tube on the floor to a certain height and install the position meter on the upper part of the metal tube to contact with the roof. The structure diagram of the roof subsidence monitoring instrument is shown in Figure 11.

4.2. Stress Monitoring of Filling Body

The load on the filling body can directly act on the pressure sensor. The sensor is horizontally placed on the stope floor. When the floor is irregular, the sensor can be fixed to the steel plate and then placed on the bottom plate, as shown in Figure 12.

The structure of the pressure sensor is shown in Figure 13. The pressure of the filling medium directly acts on the strain of the sensor, which makes the strain produce elastic deformation. The strain gauge outputs a voltage signal proportional to the force applied, and the transmitter amplifies the voltage signal and outputs it.

4.3. Monitoring Scheme

Four survey lines are set up in the stope. The first measuring line was arranged 138.7 m away from the filling cut of the working face. Due to the unclear understanding of the field temperature and equipment, the equipment failed about 8 hours after filling. The second measuring line is installed at 158.1 m of the total advance of the working face, in which 011Y # pressure sensor fails at the beginning of installation, and other equipment operates normally. The third measuring line is installed at 212.3 m of the working face, and the equipment operates normally. The fourth measuring line is installed at 456 meters ahead of the working face. The location and layout of each measuring line are shown in Figure 14, and the actual installation process is shown in Figure 15.

4.4. Monitoring Results and Analysis

Due to the high temperature of the water injection box, the connection box may be damaged due to the high temperature.

The following can be seen from the monitoring data Table 5:(1)The maximum compression amount of the filling body occurs in the middle of the filling body, which is 104.4 mm. The maximum compression amount at both ends of the filling body is only 52.7 mm, and the maximum compression rate is 1.5∼3.2% (the thickness of coal seam is 3.1 m). The maximum pressure on the filling body also occurs in the middle of the filling body, which is 5.2 MPa, and the maximum pressure on both ends is only 3.8 MPa. It can be seen that the stress in the middle of the filling body is obviously greater than that at both ends, and the compression in the middle of the filling body is the largest.The author thinks that this is mainly due to the fact that the roof in the middle of the filling body completely collapses on its upper part under pressure, and the weight of the roof and overlying strata is completely applied to the filling body, while the roof at both ends of the filling body plays a supporting role for the two ends of the roof due to the role of surrounding rock at both ends, and part of the top plate is in the “hanging” state due to the support; that is, a part of the weight is transferred to the surrounding rock All applied to the filling body. Therefore, in the later filling process, the strength of the middle part of the filling body should be increased as much as possible, and the strength at both ends of the filling body can be appropriately reduced. If the partition type filling technology is adopted and a different filling material ratio is adopted for different areas, the better filling effect can be achieved, the movement of the roof and overlying strata can be well controlled, the surface deformation can be limited, and the filling can be reduced. The cost of materials has achieved good economic benefits.(2)The maximum pressure on the top of the filling body is 5.2 MPa, the maximum pressure at the bottom is 3.7 MPa, and the difference is 1.5 MPa, which is about 30% of the maximum pressure on the top of the filling body. It can be seen that the stress on the top of the filling body is obviously greater than that on the bottom. It shows that the rock-like filling has a good function of absorbing and transferring stress and can absorb and transfer about 30% of the pressure when the roof is pressed. It also verifies that the filling proposed in [7] has good plastic strengthening characteristics, and it needs to absorb and consume more energy when it is damaged, so as to overcome the plastic deformation. This kind of performance of the filling body plays a good supporting role on the roof and reduces the risk of damage to the floor by the filling body.

5. Discussions

(1)Generally speaking, with the advance of working face, the stress of filling body increases gradually. When the working face is advanced about 40 m, the stress of filling body increases linearly with the advance of working face. With the continuous advance of working face, the increasing range of stress gradually decreases, and the roof is basically stable. What is more noteworthy is that during the advancing process of the whole working face, there is no sudden change of stress on the filling body. This is mainly because the filling body gradually forms a common support system with the surrounding rock after compression deformation, which prevents the further subsidence of the roof rock and makes the roof rock gradually tend to be stable. Therefore, the stress of the filling body is also a gradually increasing and gradually stable process.(2)The reading of the pressure box at the roof of the goaf is greater than that of the bottom, that is, the pressure on the top of the filling paste is obviously greater than that on the bottom, and the difference is about 2∼2.5 MPa, which is about 1/4∼1/2 of the maximum pressure on the roof of the filling body. This shows that the rock-like filling body has a good function of absorbing and transferring stress and can absorb and transfer the stress. At the same time, it further verifies that the rock-like filling has good plastic strengthening characteristics, and it needs to absorb and consume more energy to overcome the plastic deformation. This kind of performance of the filling body plays a good supporting role on the roof and reduces the risk of damage to the floor by filling body.

6. Conclusions

(1)Through the hierarchical fuzzy decision-making method, it can be seen that the optimal proportion scheme of gangue gypsum filling material is “E3,” that is, the cement content is 9%, the ratio of cement to gangue is 2 : 5, and the mass concentration is 74%. The fuzzy decision-making results are consistent with the results of the orthogonal test. At the same time, the data show that the sensitivity degree of “ash to gangue ratio” in the optimization of gangue gypsum filling material proportion is the highest, followed by “cement content” and “mass concentration.” Therefore, in the design of gangue gypsum filling material proportion, we should pay special attention to the “ash to gangue ratio” parameter and strictly control the amount of fly ash in the actual production.(2)AHP and the fuzzy decision-making method are used to optimize the proportion of filling materials, and the fuzzy and complex problems are mathematically and systematically combined with qualitative and quantitative methods. The mathematical model of filling material ratio optimization is constructed, which not only avoids the error of the traditional filling material test but also overcomes the disadvantage of a large amount of data of the traditional mathematical statistics method. At the same time, combined with expert experience and opinions, the optimization result is scientific, reasonable, and consistent with the actual situation, which provides a new decision-making method for the selection of filling material proportion optimization scheme.(3)The mass concentration is the main factor affecting the slump of slurry, and the specific gravity of fly ash and gangue directly affects the stratification and bleeding rate of slurry. The fluidity of the slurry can be improved by properly reducing the mass concentration of the slurry; the workability and water retention of the slurry can be significantly improved by reasonably increasing the ash gangue ratio, that is, the amount of fine aggregate in the slurry, and the stratification and bleeding rate of the slurry can be reduced. The amount of composite cementitious materials is the main factor affecting the setting time of grouting materials and the strength of filling paste. Adding a small amount of desulfurized gypsum into cement clinker can promote the hydration reaction of cement clinker, shorten the setting time of slurry, and improve the strength of the filling body.(4)It is found that the compressive strength, tensile strength, and maximum elastic modulus of backfill increase gradually with the increase of age, while Poisson's ratio first increases and then decreases with the increase of age, and the mechanical parameters tend to be stable after 28 days of age. The 28 d compressive strength is 5.89 MPa, the tensile strength is 1.17 MPa, the elastic modulus is 842.45 MPa, and the maximum Poisson’s ratio is 0.0698.(5)It is found that the stress in the middle of the filling body is greater than that at both ends, and the stress on the top is greater than that on the bottom. The maximum compression of the filling body is 104.4 mm, and the maximum compression ratio is 1.5∼3.2%. The maximum pressure on the filling body is 5.2 MPa, which is less than the filling body strength of 5.3 MPa under the optimal ratio in 3.3 above and the measured compressive strength of the filling body measured on site is 5.89 MPa, which verifies that the optimization scheme of rock-like filling material performance is feasible and the long-term stability of the filling body in goaf is good.

Data Availability

The data used in this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was partially funded by the National Natural Science Foundation of China (Grant no. 51804178), Shandong Provincial Natural Science Foundation (ZR2020QE125), and 2020 Shandong Province Key R&D Program [Soft Science Project] (ZR 2020RKB01083 and 2020RKB01167).