Abstract

The gob-side entry retaining plays an important role in improving working face ventilation, alleviating working face connection, and increasing mining revenue. According to the characteristics of the crossheading roof at the 2103 working face of a mine in Shanxi, a structural mechanics model of the roof was established to derive the theoretical formulae for the ultimate hanging arch length of the layered roof with anchors and the initial support resistance of the entry-side support. The influence factors of the ultimate hanging arch length were evaluated using local sensitivity analysis. Based on the theoretical study, the work proposed the collaborative support technology of the crossheading, collaborative support at the 2103 working face. The results showed that the ultimate hanging arch length was most influenced by the width of the plastic zone, followed by the width of the roadway, supporting strength, anchoring strength, layered thickness, and mining depth, while the ultimate tensile strength had little influence. The initial support resistance of the entry-side supports was closely related to the ultimate hanging arch length and the process of gob-side entry retaining. The improved entry-retaining supporting process could control the sharp surface convergence of the surrounding rocks of the entry retaining, the sinkage of the roof of the entry-retaining section was controlled below 100 mm, and that of the advanced section was controlled below 50 mm. The stability of the supports next to the entry is improved, and the needs of the site project are met.

1. Introduction

Gob-side entry retaining retains the haulage roadway of the upper-section working face for ventilation of the lower-section working face by strengthening the supports or using other effective methods after the upper section of the working face has been mined [14]. It has been highly valued by research departments and coal mining enterprises because of its advantages of easing mining succession tension, realizing coal pillar-free mining, improving coal resource recovery rate, isolating empty areas, avoiding air leakage, realizing Y-shaped ventilation, and solving corner gas accumulation. However, the gob-side entry retaining is subjected to the influence of two mining with strong mine pressure, and the rock destruction zone and plastic zone around the entry increases significantly, so maintenance is very difficult [58]. A large number of scholars have studied the interaction between the surrounding rock and supports, the law of surrounding rock activity along the gob-side entry retaining, and the supports in the entry, with useful results.

Li and Hua [9] studied the interaction between the key block and the surrounding rock along the gob-side entry retaining to determine the stability criterion of the key block, thus deriving the calculation formula of the support resistance next to the entry. Li [10] analyzed the process of rock movement and its deformation characteristics of the roof along the gob-side entry retaining. The design principles of controlling the support resistance of the filling body along the gob-side entry retaining at each stage are determined to establish the corresponding mathematical model of support resistance and reasonable compression. Xie [11] used numerical analysis software UDEC to establish the corresponding numerical analysis model and analyzed in detail the movement law of the surrounding rock along the integrated gob-side entry retaining. Han et al. [12] considered that the slow sinking of the hard roof plate with long-term overhang and the sudden fracture with violent incoming pressure are the root causes of disturbance along the gob-side entry retaining. Changing the fracture location, timing, and structure through active fracturing can realize the active control of the “big structure” of the roof. The dual active control mechanism of “decompression-anchoring” is described using theoretical analysis, numerical simulation, and industrial tests to address the problem of maintenance and control in the hard roof. Zhang et al. [13] focused on the problem of direct top stability control of the loose and broken top plate along the gob-side entry retaining, such as the released top coal and composite top plate in comprehensive mining. According to the spatial and temporal characteristics of the top plate activity along the gob-side entry retaining, the elastic damage mechanics model of the direct top in the filling area is established, with the expression of the direct top stress distribution in the filling area derived. The evolution law of direct top stress distribution in the filling area is studied to obtain the equations of direct top tensile stress range and horizontal misalignment range in the filling area and the corresponding evolution law.

Zheng et al. [14] proposed the in situ coal pillar containment technology for the technical difficulties faced in the traditional gob-side entry retaining. The “limited deformation” in situ coal column mechanics model is established, and the principle of structural synergistic support is proposed to determine the synergistic rock control technology of anchor beam net rope as basic support and -type steel beam + monolithic pillar as reinforced support. Kan et al. [15] clarified the movement law of the overlying rock layer of the secondary gob-side entry retaining for the problem of keeping it as a return wind passage for a long time after reusing it. A structural model of secondary gob-side entry retaining on the roof is established, with the concept of secondary along-air stay lane surrounding rock support structure proposed. Wu et al. [16], aiming at the control problem that the roof of the deep shaft along the gob-side entry retaining was difficult to form the bearing structure or even the topping disaster, integrated field research, theoretical analysis, and actual measurement methods and classified the deformation and damage forms in depth. A “multiprotection structure” control system, with the composition, control principle, support materials, and support timing studied.

However, the ultimate hanging arch length of the cantilever beam in the anchored state within the gob-side entry retaining after retrieval is usually neglected. Therefore, based on the research results of previous scholars, the work further studied the problem of gob-side entry retaining. By establishing the mechanical calculation model of the layered top slab in the anchored state, the mechanical expression of the ultimate hanging arch length of the top slab was derived and its influencing factors were evaluated primarily and secondarily. Based on the solved value of overhanging roof distance, the interaction model between entry-retaining support and the direct roof was analyzed to derive the mathematical expression for calculating the initial support resistance of entry-retaining fillers. The research results are of great practical guidance for the control of the top slab of the gob-side entry retaining and the understanding of the interaction mechanism of the surrounding rock and support. It has been applied to the field-engineering practice.

2. Mechanics Modeling and Calculation under Anchoring Condition

2.1. Mechanical Model of the Layered Roof with Anchors in the Entry Retaining

Both theory and practice show that due to the low tensile strength of the rock seam, most of the roof damage is in the form of compound flexural breakage under the action of self-weight and axial thrust. The breakage rule of the roof is generally to break after leaving the seam first (see Figure 1), and a fracture is produced in the overhanging coal wall. Then, the overhanging roof will be kept in equilibrium with the clamping structure. When the axial thrust reaches a certain limit, the bending deformation generated by it will be a vicious circle. The beam cannot maintain the new equilibrium state, leading to damage [1719].

Basic assumptions are as follows. (1) The surrounding rock satisfies all the assumptions of classical linear elasticity theory, that is, continuous, completely linear elastic, uniform, isotropic, small deformation, and mainly studies the plane deformation state. (2) The roof is in small deformation before breaking, and its overburden load does not change significantly with the advance of working face, and it is a uniform load. (3) The coal seam is strong and large, and the plastic deformation of the coal wall caused by retraction does not change its supporting effect on the roof plate. (4) The anchoring force of the anchor section is equal, and the support force provided by the anchor solid is evenly distributed on the rock layer of the roof. (5) Rock beam (plate) meets the basic conditions of elastic foundation beam (plate).

2.2. Calculation of the Ultimate Hanging Arch Length of the Layered Roof with Anchors

Based on the analysis of the lithology of the roof, the breaking process, and the characteristics of incoming pressure, we divided the roof motion into two processes: significant motion and relatively stable motion. The corresponding mechanical models were the rock beam model with solid support at both ends and the cantilever beam model with solid support at one end and cantilever at the other, and the work focused on its relative stability. Considering the unfilled wall section behind the lagging working face, it was approximated as a short-time cantilever beam model [2023]. The width of the stress limit equilibrium zone was taken as the elastic zone boundary, and the mechanical calculation equations were established to give the mathematical expressions of its limit equilibrium.

In the following, the formula for the ultimate hanging arch length is analyzed.

In Figure 2(b), the horizontal direction load is neglected under the action of overlying uniform load q, pillar initial support force P, coal pillar support force Fm in the plastic zone, and anchorage uniform load . From mechanical equilibrium equation , we obtain and . MA is the bending moment at the fixed end, FA is the y-directional binding force at the fixed end, x0 is the width of the stress limit balance zone, Fm acts at . The width of stress’s ultimate equilibrium zone x0 can be obtained from [22]. γ is the rock capacity; H is the mining depth; m is the thickness of coal seam mining; k is the stress concentration coefficient; A is the lateral pressure coefficient; Px is the support resistance of the bracket to the coal gang; C0 and φ are the cohesive force and internal friction angle at the intersection of the coal seam and the roof, respectively:

As measured, the length of x0 is generally 2–10 m.

There are three unknown quantities, and only two can be listed. This solution belongs to the primary super-stationary problem, for which another equilibrium equation needs to be sought.

Assuming that the rock beam is a rigid-like body, deflections and at the point of action of Fm and P approximately satisfy the compatibility equation of deformation (see Figure 3):

By the superposition method, and ; are the deflection values caused at action points Fm and P by q, q0, Fm, and P, respectively.

After substituting the result of the calculation into (2), . The solution is , where ; . Applying the knowledge of material mechanics, the bending moment equation of the beam is

In , are obtained by solving jointly.

In , let , the bending moment at takes maximum value ; in , the bending moment at takes maximum value ; in , the bending moment at takes maximum value ; in , the bending moment at takes maximum value . The bending moment values at each turning point can be found from Table 1.

Accordingly, we make the cantilever beam moment diagram roughly as follows.

By combining Figure 4 and comparing the extreme bending moments and values at each turning point, we can obtain the maximum value of bending moment (taking the absolute value) at the fixed end: . Considering the low tensile strength of the rock cantilever beam, the tensile stress generally exceeds the limit and causes the rock beam to be destabilized and broken. The strength condition was calculated by (4), and (5) can be obtained:

In (5), , , , , and ,where , , , , and . Considering the layered roof as a combined beam with unit width, n as the number of anchorage layers, and the thickness of each stratum as the combined beam of h, antibending section factor W is . After substituting into (5), (6) can be obtained:and is defined as the anchorage strength factor, and is the support strength factor. The resulting solution is as follows:

3. Parameter-Sensitivity Analysis of the Ultimate Hanging Arch Length of the Layered Top Plate with Anchors

3.1. Parameter Sensitivity Analysis

Parameter changes lead to numerical changes in model behaviors, and sensitivity analysis is a measure of the extent to which a particular parameter change affects the model output results. It can qualitatively or quantitatively analyze the characteristics of the parameters affecting the ultimate hanging arch length, the effect of a single parameter, or different combinations of parameters on the model performance [24]. The work used the local sensitivity analysis to analyze the response variation of the cantilever-beam mechanics’ model.

Local sensitivity analysis considers the effect of a change in a single parameter on the model response. The basic principle is to fix the other parameters and analyze only the model response results caused by the variation of a parameter in a small range [25, 26]. Assuming our objective function is with n parameters, xi is the ith parameter among them, and is the value of the sensitivity analysis of parameter as a function value of a fixed design point. In addition to parameters at this analysis point are fixed values. A small change in xi at this analysis point changes function :

Equation (8) is the sensitivity value of function to the variable parameter xi. The larger the value, the greater the influence of parameter xi on the function [27]. Ultimate hanging arch length has n influencing factors. When δ tends to be infinitesimal, the calculated expression for the variation of sensitivity is obtained as follows:where is the value of the hanging arch length at the fixed point, the variation factor of parameter xi, and is the variation value of the hanging arch length.

3.2. Evaluation of Local Sensitivity Analysis Results

Equation (7) was used to obtain the influencing parameter of ultimate hanging arch length l. The sensitivity of a given factor was analyzed using the factor change method; that is, a single factor was taken as the variable and the rest of the factors were taken as the central values. Assuming ±5%, ±10%, ±15%, ±20%, and ±25% fluctuations of each regulating factor, the value-added of the hanging arch length was obtained, and the sensitivity of each parameter change was obtained by substituting into (9).

Figure 5 shows that hanging arch length increases by ±1.41 m when the width of the plastic zone fluctuates by ±25%, and it is the most influential factor on among the parameters. The increased plastic zone makes the model boundary penetrate deeper into the coal mass. The coal wall provides more support to the overlying cantilever rock beam, resulting in the rock beam can have a larger overhead distance. When the width of the roadway fluctuates by ±25%, hanging arch length increases by ±0.87 m, which is the second most influential factor after the width of the plastic zone. As the width of the roadway increases, the range of the corresponding anchorage zone increases. For the layered roof slab, the effect of adding anchors [28] gives full play to the self-stabilizing ability of the rock mass and reduces the span of the roof. The ability of the roof to withstand deflection is further enhanced, and the ultimate suspension capacity is further increased. Analyzing the influence of anchorage strength and entry-side support strength, when the parameters fluctuate ±25%, increases by ±0.27 and ±0.60 m, respectively. The support strength is more capable of increasing the hanging arch length than the anchorage strength. When the parameters fluctuate within the same range, increased caused by mining depth, stratification thickness, and ultimate tensile strength are all below ±0.02 m, which means that the hanging arch length is less affected by these three factors.

Figure 6 reflects the degree of influence of each factor, and the sensitivity of each parameter is ranked as . From this, the plastic zone width has the greatest influence on the hanging arch length, followed by the width of the roadway, support strength, anchorage strength, stratification thickness, mining depth, and ultimate tensile strength successively.

4. Determination of the Support Resistance of the Entry-Retaining Support

4.1. Immediate-Roof Mechanical Model

During the gob-side entry retaining, the upper roof is in stable condition. After the working face advances a certain distance, timely strengthening of the support plays an important role in controlling the deformation of the entry retaining. The support near the entry acts on the immediate roof to stop the delamination of the immediate roof and the upper roof above the bearing structure. The support near the entry retaining should provide a certain initial support force or inhibit the resistance of rock beam sinking to maintain the structural stability of the entry roof and fillers. Meanwhile, the initial resistance of the entry should cut off the immediate roof on the side of the mining area and guide the broken collapse of the upper roof to maintain the stability of the coal gang on the nonrecovery side and the entry-side support as the bearing support of the immediate roof. Figure 7 shows the established mechanical model of the support near the entry retaining and the immediate-roof mechanical model.

4.2. Calculation of the Support Resistance of the Entry-Side Support

After the immediate roof is stabilized, rock mass A, B, C, and D in Figure 7 are temporarily stabilized under the support of the coal wall and the entry-side support. is used to obtain the calculation formula of the initial support resistance:where Pu is the support resistance of fillers, l is the ultimate hanging arch length, calculated by (7), q is the average load acting on the immediate roof, γ is the unit weight of the immediate roof, d is the width of the support zone near the entry, j is the distance of the support near the entry from the heading stope, u is the length of the entry-retaining every day, and q1 is the average support strength of the entry support to the immediate roof.

4.3. Analysis of Examples
4.3.1. Project Overview

Coal 2# was mainly mined in a mine in Shanxi with an average coal seam thickness of 1.62 m. The retrieval height along the roof of coal 2# was 1.8 m, with an average mining depth of 230 m. The track crossheading was dug along the roof of the coal seam, and the immediate roof of the coal seam was mudstone with a thickness of 1.6 m. The upper roof was fine sandstone with a thickness of 2.25 m. The width and height of the entry section were 4.1 × 2.2 m. The lithology and physical and mechanical parameters of the roof of the entry can be found from Table 2, and the geological comprehensive histogram of 2# coal seam is shown in Figure 8.

4.3.2. Calculation of the Hanging Arch Length

The formula for calculating the load on the upper floor is used to obtain q1 = 57.6 KPa, (q2)1 = 72.2 KPa, (q3)1 = 83.3 KPa, and (q4)1 = 31.0 KPa. The influences of the first, second, and third layers on the first layer should be considered, so the size of the load on the first layer of the upper floor is 83.3 KPa. Then, the load on the cantilever beam is q = 83.3 +  = 124.9 (KPa), anchor preload is 80 kN, and is calculated as 2.98 m by (1); values of each parameter can be found from Table 3.

The data in Table 4 is substituted into (7) to find that the ultimate hanging arch length l is (10.2–10.5) m, which is consistent with the actual observation in the field.

4.3.3. Calculation of Initial Support Resistance of the Gob-Side Entry Retaining

The lower crossheading of the 2103 working face is the gob-side entry retaining, and the reinforced support section of the entry lags behind the working face by 4 m. The parameters are as follows: l = (10.2–10.5) m; q = 124.9 KPa; γ = 25 kN/m3; d = 2.2 m; j = 4 m; u = 6 m; q1 = 125 KPa.

Each parameter is substituted into (11) to get initial support resistance Pu = (166.7–229.2) kN/m. Considering the safety factor taken as 1.2, the initial uniaxial compressive strength of the entry-side support should be not less than 0.46 MPa.

5. Cosupporting Technology

5.1. Roof Support

The support strength of the roof needs to be increased to stabilize the roof anchorage structure before mining and in the entry-retaining period. Reference [29] was studied for comprehensive support design, and the basic parameters of reinforcement are as follows:(1)A set of monolithic anchor ropes was supplemented in the middle of the roof, with anchor rope specification Φ17.8 × 5,200 mm and a rectangular tray of 300 × 300 × 15 mm. Each anchor cable adopted 1 CK2335 and 2 Z2350 anchoring agents with preload not less than 80–100 kN. The row distance is 1,000 mm, the same as anchor rods.(2)Both sides of the roof were 400–600 mm from the gang, and the outward angle of the anchor cable was 20–30°. Equip with channel steel 14 and arrange in the form of the anchor cable beam with a beam length of 2.8 m and a hole distance of 2.4 m. Each set of anchor cables was equipped with a 150 × 100 × 5 mm flat steel plate, and the anchor cable specification and anchoring performance and construction requirements were the same as those of the anchor cable in the middle of the floor.(3)The order of reinforcement was carried out from inside out to the stopping line at the location of the working face.

5.2. Auxiliary Strengthening Support in the Entry Retaining

The internal auxiliary reinforcement range of the entry was 60 m ahead of the working face and the space left in the entry retaining. Four monohydraulic pillars with cross-articulated roof beams were used to strengthen the support in the advanced area. Five monoblock pillars were arranged in the gob-side entry retaining, and the bottom beam was installed to control the floor heave of the entry actively in the recovery period (see Figure 9).

The roof should be reinforced by slurry spraying when the following three conditions occur. ① For the collapse column and continuous fault-zone-affected area, the roof sinks 500 mm or more or continuous netting. ② For the small fault affected area, the interior of the entry retaining shows steps, with anchor cables pulled off. ③ Drenching sections is the third condition.

5.3. In Situ Observation
5.3.1. Observation Program

The deformation of the roadway surface is an important index reflecting the magnitude of the roadway surface displacement and the degree of convergence of the roadway section. The main indices include roof subsidence, floor swelling, and roadway displacement. The observed values can be used to judge whether the deformation of the surrounding rock can ensure the normal use of the roadway and whether the support method can meet the requirements of safe production of the roadway. Therefore, the observation method of multisection rapid roadway surface convergence is essential to measure the safety of the roadway. The principle is to set up dense measuring stations at a certain distance in the excavated roadway with a considerable length. Each observation station conducts high-frequency observations from its initial setup. The space is exchanged for time, and the law of surrounding rock deformation is mastered in a relatively short period.

(1) Survey Station Arrangement. The station layout is based on the cross-point method (see Figure 10).

The rules of station layout are as follows. Taking the direction of advancement of 2103 working face as the positive direction, point A is located on the nonrecovery side of the entry; point B is located on the recovery side; point C is located on the roof of the entry; point D is located on the floor of the entry. When the measurement points are laid out, short anchor rods and reflective stickers are posted at points A, B, C, and D as permanent observation points.

(2) Spacing of Measuring Stations. The layout of the measuring station should take into account the speed of the working face and the law of mine pressure appearing in the entry. When the retrieval speed of the working face is faster, the measurement station should be arranged denser to increase the sampling frequency of deformation of the entry-surrounding rocks in time and space. According to the advance rate of the working face, the spacing of the measuring station arrangement is 10 m. The track crossheading of the 2103 working face is arranged in the whole range (see Figure 11).

According to the original data of each measurement point and the monitoring results, Figure 11 shows the analysis. The horizontal coordinate of the graph indicates the distance of the monitoring point from the current working face, with the direction of the working face advancement as the positive direction and the direction of the mining void area as the negative direction. The vertical coordinates indicate the displacement of each monitoring point.

It should be noted that in Figure 12, the range from 0 to −300 m has been mined, and the range from 0 to 300 m is the range that has not been mined. So, it is the deformation of goaf from 0 to −300 m about “distance of goaf” in Figure 12. On the other hand, it is the deformation of coal wall from 0 to 300 m about “distance of goaf” in Figure 12. The range from 0 to −300 on one side of the goaf in Figure 12 shows the deformation on one side of the goaf within the mined range.

5.3.2. Analysis of Observation Results

(1)Entry-retaining section is as follows: the average width of the entry in the entry-retaining section of 2103 working face was 3,620 mm, and the deformation of the two gangs was mainly for the solid coal gangs. The average height of the entry was 1,991 mm, and the roof and bottom plates were between 300 and 800 mm. Besides, the compression deformation of the coal mass and rock mass on the solid coal side was also large, generally around 300 mm.(2)Advanced-supporting section is as follows: the advanced-supporting distance was 60 m, and the average width and height of the entry in this section were 3,928 and 2,104 mm, respectively. The two gangs were shifted 100–150 mm, with the heaving floor above 200 mm and the roof sinking less than 50 mm.(3)The average width and height of the entry from the advanced-supported section to the 250 m entry in front of the working face were 4,043 and 2,110 mm, respectively. The convergence and deformation of the entry were not obvious, with the displacement of the two gangs less than 100 mm and the sinking of the roof less than 50 mm. There was no obvious separation between the anchor rod anchorage area and the anchor cable anchorage area.

In summary, the overall entry in the track crossheading of the 2103 working face is good (see Figure 13), but the monitoring of mine pressure should be further strengthened, especially the change of top plate sinking and monolithic support force. In addition to the collapse column and other structural areas that need to be observed, it does not need to reinforce the entry and the pre-mining roof. After retaining the entry according to the law of mine pressure monitoring, the retraction was performed gradually and orderly.

6. Conclusions

In this paper, we studied the ultimate hanging arch length of the layered roof with anchors and the gob-side entry retaining by theoretical analysis, numerical simulation, and field test. The following results were obtained:(1)A simplified mechanical model of the cantilever beam in the stable period of the layered roof with anchors was established, and the theoretical formula of its ultimate hanging arch length was derived accordingly. The local sensitivity analysis method was used to study the priority of the influencing factors of the ultimate hanging arch length. The width of the plastic zone had the greatest influence on the hanging arch length, followed by the width of the roadway, support strength, anchorage strength, layered thickness, mining depth, and the ultimate tensile strength orderly. Based on this formula, the ultimate hanging arch length of the roof of the intake airway at the 2103 working face of coal 2# in a mine in Shanxi was 10.2–10.5 m, which was consistent with the actual observation on site.(2)The mechanical model of the interaction between the entry-side support and the immediate roof was constructed when the immediate roof collapsed, and the calculated expression of the initial support resistance of the immediate roof was analyzed. Initial support resistance Pu = (166.7–229.2) kN/m was obtained from this formula for the intake airway of the 2103 working face of coal 2# in a mine in Shanxi, according to which the initial uniaxial compressive strength of the roadside support body should be not less than 0.46 MPa.(3)According to the theoretical calculation of the ultimate hanging arch length and the initial support resistance value of the entry-side support, guidance was given to the process of retaining the entry along the intake airway of the 2103 working face. The observation showed that the sinkage of the roof of the entry-retaining section was controlled below 100 mm, and that of the advanced section was controlled below 50 mm. The overall effect of the entry was good and met the requirements, which provides a strong guarantee for the stability control of roadway surrounding rock and mine safety production.

Data Availability

The data used in the field measurement can be obtained from the first author upon request.

Conflicts of Interest

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

Acknowledgments

The authors acknowledge the financial support for this work provided by the 2021 Year Postgraduate Research and Practice Innovation Project of Jiangsu Province (Grant no. KYCX21_2382) and the National Natural Science Foundation of China (Grant no. 51204159).