Three-Dimensional Physical Simulation and Control Technology of Roof Movement Characteristics in Non-Pillar Gob-Side Entry Retaining by Roof Cutting

An overlying rock structure plays a key role in controlling the roof deformation of nonpillar gob-side entry retaining by roof cutting. On the bases of the actual geological conditions of II 632 Haulage Roadway at the Hengyuan coal mine, a similar threedimensional simulation experiment of roof precutting is conducted.,ereafter, the caving characteristics andmigration law of the roof strata in the strike and dip directions are obtained. Moreover, the roof of the retained roadway and key strata of the goaf can form a hinge structure of the key blocks. By monitoring the deformation of the surrounding rock and stress distribution of the roof, the skew deformation characteristics of roadway roof are obtained. By observing the borehole peeping technology, the roof subsidence near the goaf is determined to be greater than that of the solid coal side, and the roof subsidence of the gob-side entry retained by roof cutting is greater than that of the floor heave and two sides approaching. Results of the three-dimensional similar simulation experiment indicate that the mechanical structure model of the key block of the retained roadway roof is constructed, and the mechanical analytical solution of the required support resistance of the retained roadway roof is obtained. ,is study proposes the constant resistance and large deformation anchor cable reinforcement support method to control the roof deformation of the retaining roadway. ,rough engineering application, the maximum value of the roof and floor movement of the retained roadway is stable at approximately 650mm. ,e retained roadway can meet the demand of the next mining face.


Introduction
Nonpillar gob-side entry retaining by roof cutting mining is one of the key technologies for safe and efficient coal production in China. e application of this technology can effectively improve the recovery rate of coal resources and alleviate the shortage of coal resources [1][2][3][4]. In actual mining activities, overburden movement has a substantial effect on the roof of roadway; the disturbance range of the overburden breaking and falling is large, thereby resulting in considerable roof deformation in the process of retaining roadway [5][6][7].
Numerous Chinese scholars have conducted research to solve the problem of roof stability control of the gob-side entry retaining by roof cutting without coal pillar. Zhang et al. [8] analyzed the roof fracture structure and its movement and deformation law of gob-side entry retaining by roof cutting and established the fracture mechanics model of a bilateral cantilever beam of the lower cutting slit and upper bending crack. Su et al. [9] studied the mechanism of roof cutting and pressure relief in the gob-side roadway (GR) in deep coal mine, and the overburden structure and migration law of GR in deep coal mine have been studied. Wang et al. [10] studied the entire process of overburden movement of self-formed roadway by roof cutting and pressure releasing through physical simulation test and determined that, through roof presplitting and cutting joints, the surrounding rock pressure of the roadway can be effectively relieved, thereby controlling the surrounding rock deformation. Zhang et al. [11] analyzed the disturbance factors during the period of "mining-retaining-using," through numerical simulation and on-site measurement; the stress and deformation characteristics of the surrounding rock during the reuse of the self-formed roadway were analyzed comprehensively. He et al. [12][13][14] used the mechanism of roof presplitting and pressure releasing at the goaf side as basis to establish the mechanical model of "surrounding rock structure-roadway side support body" under different roof positions. ey also proposed surrounding rock control technologies, such as shaped-charge presplitting blasting, constant resistance large deformation anchor cable, and dense single pillar beside roadway.
Hua et al. [15] analyzed the fracture mechanism of basic roof presplitting blasting, established the fracture criterion based on tensile strength, and determined the minimum quantitative relationship between charge length and hole spacing when basic roof presplitting blasting was formed. Yang et al. [16] took coal mining in Shen Fu and Dongsheng coal mining areas as the engineering background; the roof deformation laws and mechanism in the process of gob-side entry by roof cutting were comprehensively studied by means of theoretical analysis, numerical simulation, and field test. Wang et al. [17] explored the mine pressure behavior characteristics of roadway under the condition of the gob-side entry retaining by roof cutting, a structural mechanics model, where "load-bearing structure" and "loadreleasing structure" are established and analyzed. Zhang and Han [18] proposed the concept of composite load-bearing structure of gob-side entry retaining, in which roof pressure is shared by the roof, wall, and floor. Gao and Guo [3] and Ma and He [19] analyzed the evolution law and mechanical mechanism of the surrounding rock structure of gob-side entry retaining based on the mechanical mechanism of the nonpillar gob-side entry retaining by roof cutting. To study the deformation mechanism of the filling body in the goaf side, Sun et al. [20] used key block theory as basis to establish the mechanical model of the basic roof fracture structure of the retaining roadway. e preceding research also analyzed the influence of the cutting angle and cutting height on the roof support resistance of the retaining roadway. Sun and Zhang [21] analyzed the deformation and failure characteristics of the surrounding rock and established the mechanical damage model of goaf support based on Coulomb Earth pressure theory. Wang et al. [22] introduced the roof presplitting blasting technology of gob-side entry retaining by roof cutting, the technology of forming support column by caving gangue along the goaf, and the concept of automatically forming the gob-side entry by roof collapse.
Guo et al. [23] deduced the relationship between the cutting angle and cutting height by establishing the mechanical model of the overburden structure. ey also determined the reasonable parameters of the precrack roof cutting using a numerical calculation method. Guo et al. [24] studied the interaction relationship among roof bolt, tray, and bolt through numerical simulation and field combination. Sun et al. [25] determined that the key parameters affecting gob-side entry retaining in thin coal seam are presplitting cutting height, presplitting cutting angle, and the distance between presplitting blasting holes by analyzing the roof stress state in the process of thin coal seam mining. Based on the principle of gob-side entry retaining technology, Chen et al. [26] deeply studied the pressure relief effect, surrounding rock stress evolution law, and control technology of deep gob-side entry retaining technology. Yuan et al. [27] studied the reasonable parameters of gobside entry retaining by roof cutting in the face of fully mechanized mining with thin direct roof and large mining height by combining theoretical analysis, numerical simulation, and field monitoring. e preceding studies have indicated the extensive achievements in the research of roof structure, roof cutting parameters, and supporting force of gob-side entry retaining by roof cutting. However, only a few studies have been conducted on the interaction between the roof structure of roof gob-side entry retaining by roof cutting and the overall movement of the overlying strata. In addition, a two-dimensional plane model is often used in physical experiments, and the actual underground activities are in three-dimensional space. A three-dimensional model is used to study a physical simulation experiment platform, which can objectively reflect the movement characteristics of the overlying rock of the retaining roadway. erefore, this research takes the II 632 haulage drift of the Hengyuan coal mine as an example. By conducting a three-dimensional similar simulation physical experiment, the roof structure characteristics and overlying rock movement law under the two dimensions of strike and dip are studied. On the bases of the deformation characteristics of the retaining roadway roof, the mechanical structure model of the key block of the retaining roadway roof is constructed, and the mechanical analytical solution of the retaining roadway roof supporting resistance is obtained. is study proposes the method of constant resistance and large deformation anchor cable reinforcement support to control the roof deformation of retaining roadway. Field application shows that the effect of roadway retaining is good. Moreover, the research results can provide a theoretical basis and reference value for the roof stability control of gob-side entry retaining by roof cutting.  Table 1, and the comprehensive histogram is shown in Figure 1.

Model Experiment Scheme Design.
e three-dimensional stress model is adopted in this similarity simulation test. And, the experimental process is as follows.
(1) e model size is 1100 mm × 1000 mm × 1000 mm, simulation cutting angle is 15°, cutting height is 8 m, and the iron wire with a similarity ratio of 1 : 100 is used to simulate the bolt and anchor cable (see Figure 2) (2) e simulated bolt diameter is 20 mm, length is 2400 mm, anchor cable diameter is 17.8 mm, and length is 6200 mm (3) e long wood bar is used to simulate the coal seam mining, and one wood bar is drawn out each time to simulate coal mining (4) Rock layers with different lithology in the model are made and laid by water, gypsum powder, fine river sand, and barite powder with different proportions; moreover, the rock specimens with different proportions meet the specified strength and weight requirements through indoor rock mechanics tests (5) During the experiment, according to the rock stratum size of the experimental model, the floor rock stratum was laid first (6) e laying of coal seam is replaced by wooden bars, and the coal seam excavation is carried out by extracting the wooden bars; the specification of each wooden bar is L × B × H � 1000 mm × 50 mm × 30 mm (7) When laying the test model, the thin plastic sheet is embedded in advance to cut the seam in the roadway roof, the angle between the plastic sheet and the vertical direction is 15°, and the size of plastic sheet is L × B � 1000 mm × 80 mm (8) e strata above the coal seam are laid in order according to the size of the model, during the laying process, the angle of the plastic shall not change, and this method is used until the end of the laying of the overlying strata (9) When the model is stationary for about 10 days, the outer baffle of the model will be removed when the excavation requirements are met (10) e model surface is coated with white coating, the displacement observation line and displacement observation point are arranged, the strain gauge and strain gauge are connected, and the data terminal is debugged, and then the next excavation experiment is carried out (11) Each time along the strike of the working face, a wooden bar is drawn out to simulate the mining of a knife of coal, a total of 20 knives are mined, that is, the working face is pushed forward for 100 m, and each knife is allowed to stand for 2 h (12) A group of strain gauges are embedded at the direct roof 40 mm above the roadway, and the spacing of each strain gauge is 200 mm, as shown in Figure 3; the stress change of the roof of the retained roadway during the mining process is observed

Movement Law of Overburden in the Working Face along
Incline. When the working face advances 10 m, as shown in Figure 4(a), a small amount of roof subsidence began to appear above the roadway, bed separation fissures appeared at 4 m and 8 m above the direct roof of the coal seam, and the basic roof was in a stable state. When the working face is advanced to 20 m, as shown in Figure 4(b), the direct roof strata over the goaf collapses, and the roof of the key block at the cutting slit side of the roadway collapses along the reserved roof cutting line, thereby cutting off the physical connection between the roof above the roadway and the roof of the goaf. Furthermore, the roof of the roadway and roof of the goaf lose the mechanical constraints. When the working face is advanced to 45 m, as shown in Figure 4(c), the reserved cutting joint developed and expanded upward. Vertical microcracks appear near the cutting slit, and the gangue in the goaf is gradually compacted in the movement of the overlying strata. e goaf tends to produce a large area of arch fracture, and evident separation layer appears above the basic roof. A relatively stable rock block hinge structure is formed above the retained roadway, and the stable hinge structure plays an important role in the process of controlling the roof stability of the retained roadway.
When the working face is advanced to 80 m, as shown in Figure 4(d), the goaf is further compacted under the action of the dead weight of the roof rock and the bearing pressure. Moreover, the separation value of the upper strata decreases, and the roof shows evident bending and subsidence deformation, thereby forming a bending zone. Lastly, the vertical cracks continue to develop at the upper strata. e process of coal seam mining is simulated. In the process of coal seam mining, close range photography is used to monitor the activity of the overlying strata. Four groups of measuring points are set at heights of 4 m, 8 m, 13 m, and 18 m of the coal seam roof. e roof subsidence curves of the different layers are shown in Figure 5.
When the working face is advanced by 20 m, the direct roof collapses for the first time. At the initial stage of the working face, the roof subsidence is small, mining activities have minimal influence on the movement of the overlying strata, displacement changes gradually, and the displacement curve trend is trough-shaped. At the level of 4 m height of the coal seam roof, the roof collapses between 10 m and 50 m away from the open-off cut, and the maximum subsidence is 2.5 m.
When the working face is advanced 45 m, as shown in Figure 6, the direct roof of the working face completely collapses. From the roof subsidence change of the working face, the roof of the No. 6 coal seam collapses between 10 m Shock and Vibration 3 and 50 m away from the cut hole at 4 m height. Furthermore, the maximum subsidence is approximately 2.8 m.

Movement Law of the Overlying Rock along the Strike of the
Working Face. Figure 7 shows that when the working face is advanced by 20 m, the roof of the goaf is suspended, and the area of the direct roof above the working face is gradually increased. When the working face is advanced to 30 m, cracks appear in the roof 5 m and 8 m above the coal seam. Moreover, the roof of the goaf lagging behind the temporary support area begins to sink gradually. When the working face is advanced 45 m, the immediate roof above the goaf collapses for the first time, the roof caving height is 7 m, and the caving step distance is 40 m. Cracks in the upper strata of the basic roof further develop and expand, and these cracks evidently show a trend of expansion. When the working face advances 80 m, the roof shows a second large-area collapse phenomenon, and the collapse occurs 20 m above the basic roof, which is an arch fracture. As the working face continues to advance, square collapse occurs, the roof subsidence deformation is severe, and the horizontal roof falling strata are hinged with each other, thereby forming a masonry beam structure. When the working face advances 100 m, the broken gangue in the goaf is further compacted, and the caving layer of the overlying strata in the goaf roof gradually increases. Figure 8 shows that when the coal seam is advanced to 45 m along its strike direction, the roof is broken at

Deformation Characteristics of the Roadway Surrounding
Rock.
e borehole peeper is placed 20 m away from the open-off cut in the roadway to observe the internal deformation of the surrounding rock of the retained roadway (see Figure 9). When the working face is advanced to 20 m, roof subsidence occurs near the goaf, and the floor heave is negligible. When the working face is advanced to 45 m, the roof subsidence near the goaf is evident. When the working face is advanced to 80 m, the maximum roof subsidence is approximately 10 mm, showing skew deformation, and the roof subsidence near the goaf is greater than that of the coal pillar side. During the mining period, the roof subsidence of the gob-side entry retaining by roof cutting is greater than that of the floor heave and roadway wall approaching.

Distribution Law of the Roadway Roof Stress.
rough the embedded strain gauge, the roof stress change above the roadway is obtained, as shown in Figure 10. When it is 20 m to 40 m ahead of the working face, stress increases

Bolt
Anchor cable Roadway

Mechanical Analysis of the Hinge Key Block for Roof Cutting and Pressure Releasing
According to the law of overburden breaking and migration analyzed by the three-dimensional similar simulation experiment, the roof of the retained roadway is fractured on the side of the solid coal, and the roof of the gob-side entry retaining overall produced rotary subsidence. Hence, the hinge structure of the key block is formed, which jointly carries the load of the overlying strata. After the roof of the retaining roadway rotates and sinks, the hinge structure of the key block formed by the roof of the retaining roadway and support resistance of the surrounding rock of the retaining roadway form a stable structure together, thereby forming a stable structure, as shown in Figure 11. In this figure, F s is the supporting force of the solid coal side, kN; Fm is the supporting resistance of the roadway, kN; F g is the supporting force of the gangue, kN; x 0 is the plastic zone range of solid coal side, m; B is the width of the roadway, m; a is the action range of gangue supporting force, m; h is the height of the coal seam, m; h f is the height of immediate roof, m; h d is the cutting height, m; β is the cutting angle,°; L i is the breaking length of the key block B, m; F h is the weight of the direct roof and key block B; and θ is the turning angle of the key block B.

Ultimate Breaking Length of Key Block B L i .
eory of broken line states that when the goaf beside the retained roadway is periodically collapsed, the broken length of the key block above the retained roadway roof L i is as follows [28]: where S is the length of the working face, m, and L is the basic jacking step, m.

Supporting Force of Goaf of the Gangue on Key Block F g .
When the long cantilever is cut off to form a short cantilever, in the process of the key block B of the retained roadway roof turning and sinking and with the continuous compaction of the gangue contacting with the retaining roadway roof, the compression of gangue in the goaf under the rotation of key block B is △u. e entire key block B forms the hinged pressure bearing structure of the key block from the start of sinking to the release of sinking. e supporting force F g of the gangue per unit area in the goaf is as follows [28]: Δu where K g is the support coefficient of the caving gangue, MPa/m; K c is the residual crushing expansion coefficient of the collapsed gangue after compaction, which can be determined according to the simulation test or field measurement; K s is the crushing expansion coefficient of the roof rock; h f is the direct roof height, m; and θ is the rotation angle of the retained roadway roof. ereafter, according to formula (1)-(3), the supporting force F g of the goaf gangue to the key block B is as follows: e expression of the limit equilibrium zone width x 0 of the coal body in the retaining roadway is as follows [12]:     Shock and Vibration where A is the lateral pressure coefficient; p x is the support strength of the coal wall, MPa; h is the buried depth of the roadway, m; c 0 is the bonding strength of coal and rock mass, MPa; δ is the friction angle in the coal body,°; c is the average density of the overlying strata, kN/m 3 ; k is the stress concentration factor; and h is the height of roadway, m.    Shock and Vibration

Supporting Force of the Solid Coal
where φ 0 is the internal friction angle of the interface between the top and bottom slate,°.

Supporting
Force of the Roadway Side Support. e horizontal thrust of key block B is as follows: where q 1 is the weight per unit length of the basic roof and its load, kN; h f is the thickness of the basic roof, m; and △u is the subsidence of rock block B at B′, m. e static balance method is used to analyze key block C. From the formula F x � 0, From the formula F y � 0, From the formula M B′ � 0, e results are as follows: From the formula M A′ � 0, e results show that the support resistance of the roadway side is as follows: where F m is the support resistance of the roadway. Formula (13) shows that when the retaining roadway roof breaks at the solid coal side, the roof of retaining roadway is in a stable state under the action of its own structure and supporting force of solid coal wall side and active supporting resistance. Under the condition that the supporting force of retaining roadway roof structure and solid coal side is fixed, the roof deformation control is mainly realized by increasing the supporting resistance provided by anchor cable. erefore, constant resistance and large deformation anchor cable is used to reinforce the roof of the retained roadway.

Deformation Control of the Surrounding
Rock in the Gob-Side Entry Retaining by Roof Cutting

Support Parameters and Control Technology of the Gob-Side Entry Retaining by Roof Cutting.
According to the results of three-dimensional similar simulation experiment and theoretical analysis, combined with the engineering background of II 632 working face in Hengyuan coal mine, the key technologies of active and passive strengthening support during roadway retaining are proposed, as shown in Figure 12. e roof is composed of a bolt with diameter of Φ20 mm, length of 2400 mm, ribbed steel bar, spacing of 750 mm, row spacing of 800 mm, anchor cable diameter of Φ17.8 mm, length of 6200 mm, spacing of 1300 mm, row spacing of 1600 mm, M3 steel strip, and 5.0 m × 1.0 m steel mesh. In the mining advance affected area, constant resistance large deformation anchor cable support and single pillar support are added. e diameter of the constant resistance large deformation anchor cable is Φ21.8 mm and length is 9500 mm. ree rows are arranged. e distance between the constant resistance anchor cable on the cutting slit side and roof cutting line is 500 mm and the row distance is 1000 mm. e distance between the second row of constant resistance anchor cable and roof cutting line is 1500 mm and row distance is 1600 mm. e distance between the constant resistance anchor cable on the solid coal side and the coal wall is 1500 mm and row distance is 3200 mm. ree rows of single pillars are arranged along the roadway section, and the advance support range is 40 m. Two rows of single pillars are added 150 m behind the working face, and U-shaped steel and metal mesh are used as gangue retaining support at the side of the goaf. After the retaining roadway is stable, the single passive support is removed.

Analysis of the Roof Control Effect of the Gob-Side Entry
Retaining by Roof Cutting. To verify the roof control effect of nonpillar gob-side entry retaining by roof cutting, the deformation of the surrounding rock of the roadway retaining along the goaf 500 m in II 632 transportation chute was observed. Moreover, a survey station is arranged every 50 m in advance working face 100 m; the roof and floor displacements were measured using the cross-shaped measuring point method. And five observation stations were arranged in the transportation chute to continuously observe the roadway deformation, as shown in Figure 13. e influence range of one mining advance is 42 m, in which the deformation of roof and floor is considerably severe, the roof subsidence increases from 32 mm to 114 mm, the floor heave increases from 30 mm to 102 mm, and the roof and floor convergence increases from 52 mm to 216 mm. e influence range of the roadway retaining by one time mining is 96 m behind the working face. e deformation of the roof and floor is the most severe in this range, and mine pressure behavior is evident. e roof and floor approach increases from 216 mm to 650 mm. e deformation rate of roof and floor is small in the influence stage of roadway retaining by one mining, and the deformation rate is from 16 mm to 28 mm during severe deformation.
Under the condition of grading support in different stages, the maximum roof and floor movement of the II 632 transport gateway is stable at approximately 650 mm. rough field photography, as shown in Figure 14, the stability of surrounding rock of roadway is good. According to statistics, the grading support method in different stages can meet the needs of safety production of working face and also achieve significant economic benefits.

Conclusions
(1) rough the three-dimensional physical simulation experiment, the following results are obtained. In the dip direction, vertical microcracks appear at the overburden cutting joints; when the basic roof collapse and fall, the overlying strata of the roadway and the strata on the side of the goaf form a relatively stable hinge key block structure. In the strike direction, the caving layers gradually increase; the horizontal roof caving strata are hinged with each other to form a stable masonry structure. (2) e deformation monitoring of the roadway surrounding rock indicates that the roof of the roadway deformation is skewed, the roof subsidence near the goaf is greater than that of the coal pillar side, and the roof subsidence is greater than the floor heave and the two sides approach. e roof of the roadway stress reaches the peak when it is approximately 20 m away from the working face, and the distance of the lagging working face affected by one-time mining is from 20 m to 60 m.
(3) rough a similar simulation experiment, the stable hinged key block structure is formed after the basic roof of the overburden is broken. By constructing the mechanical model of the hinged key block structure, the mechanical analytical solution of the retaining roadway roof support resistance is obtained.  Combined with the actual geological conditions of the Hengyuan coal mine II 632 and the constant resistance and large deformation anchor cable proposed, the method of roof reinforcement and support effectively controls the deformation of the retained roadway roof. (4) According to the field measurement, the maximum subsidence of the roof of the retained roadway roof is approximately 360 mm, the maximum floor heave of the floor of retained roadway is approximately 290 mm, and the retained roadway without coal pillar is between the canceled working faces; moreover, it can meet the use requirements of the next mining face.
Data Availability e data used in the field measurement can be obtained from the corresponding author upon request.

Conflicts of Interest
e authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.