A New Gob-Side Entry Layout for Longwall Top Coal Caving

: In China, gob-side entry is typically located one pillar-width (less than 5 m) away from the previously mined-out panel in a conventional longwall panel using top coal caving (LTCC). Design of gob-side entries is a challenge due to the complex dynamic loading process during their service life. A new gob-side entry (NGE) design practice is presented here for Zhenchengdi Colliery, which has many advantages. A theoretical analysis is presented followed by numerical modeling. The modeling included the double-yield constitutive model for gob behavior to analyze the stress environment for the gob-side entry with validation through ﬁeld observations. The results indicate that pre-mining stress within gob occurred 51 m away from the gob edge. The NGE is located within a destressed zone in the entire panel system, which results in lower side abutment pressure for the adjacent panel to be mined. The stress concentration around the gob edge near the gob-side entry is relatively low (less than 0.1) compared with the other side of the gob. Field observations indicate that: steel sets with wire mesh on the top, can maintain ground control and ventilation during active mining; (2) the periodic weighting interval is 9–12 m; (3) the length of the block “B” which “protects” the gob-side entry is 10–13 m; (4) stress distribution formula for the elevating section is derived; and (5) the roof pressure of the gob-side entry is much smaller than non-gob side entry and is smaller than the pre-mining stress. Deformation data of the gob-side entry shows that both roof-to-ﬂoor and rib-to-rib convergences are smaller than for the non-gob-side entry with an improved overall stress environment. Theoretical analysis, numerical modeling, and ﬁeld observation are consistent, which validates the scientiﬁc foundations of the new technology.


Introduction
Mechanized, single-cut retreat longwall mining with top coal caving (LTCC) is widely used for underground coal extraction around the world. Development entries, or gateroads, play an important safety and productivity role in mining operations in LTCC. A wide pillar is left unmined between adjacent panels to ensure stability of the entry, as shown in Figure 1. This, however, can lead to low coal recovery. To solve this problem, the gob-side entry approach has been proposed by researchers [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. As shown in Figure 2, it is a combination of a gateroad, pillar, and gob, which is a typical panel layout with slender pillar or artificial pillar or wall. Gob-side entry are of two types [2]: (i) gob-side entry driven in a new panel with only a slender (often less than 5 m) pillar away from the gob, as shown in Figure 2a; and (ii) retained gob-side entry (Figure 2b), which is an entry of a previous panel retained as the entry for a future panel by constructing an artificial wall along the gob edge of the previous mined-out panel [3,4]. retained as the entry for a future panel by constructing an artificial wall along the gob edge of the previous mined-out panel [3,4]. After excavation of a panel, the stress on the solid coal on the right (Figure 3) can be divided into four zones [2,5]: I-de-stressed yield zone, II-over-stressed plastic zone, III-over-stressed elastic zone, and IV-pre-mining vertical stress zone. Gob-side entry is located within zone I to avoid a high stress concentration in zones II and III, and to avoid a large size pillar in zone IV, which would result in low recovery. Gob-side entry makes use of the side abutment pressure redistribution to locate gobside entry in the de-stressed zone within the solid coal, i.e., zone I in Figure 3. Yan [3] stated that the de-stressed zone occurs 0-7 m and the main roof will form the stable masonry structure to provide a low-stress environment for gob-side entry. He suggested driving the gob-side entry along the gob with a 5 m slender pillar. Li [6] noted that a gob-side entry with a 6 m wide coal pillar based on numerical modeling was employed in LTCC in mining a 15 Figure 1. A conventional development entry design in LTCC.
Energies 2018, 11, x FOR PEER REVIEW 2 of 24 retained as the entry for a future panel by constructing an artificial wall along the gob edge of the previous mined-out panel [3,4]. After excavation of a panel, the stress on the solid coal on the right ( Figure 3) can be divided into four zones [2,5]: I-de-stressed yield zone, II-over-stressed plastic zone, III-over-stressed elastic zone, and IV-pre-mining vertical stress zone. Gob-side entry is located within zone I to avoid a high stress concentration in zones II and III, and to avoid a large size pillar in zone IV, which would result in low recovery. Gob-side entry makes use of the side abutment pressure redistribution to locate gobside entry in the de-stressed zone within the solid coal, i.e., zone I in Figure 3. Yan [3] stated that the de-stressed zone occurs 0-7 m and the main roof will form the stable masonry structure to provide a low-stress environment for gob-side entry. He suggested driving the gob-side entry along the gob with a 5 m slender pillar. Li [6] noted that a gob-side entry with a 6 m wide coal pillar based on numerical modeling was employed in LTCC in mining a 15 m thick coal After excavation of a panel, the stress on the solid coal on the right ( Figure 3) can be divided into four zones [2,5]: I-de-stressed yield zone, II-over-stressed plastic zone, III-over-stressed elastic zone, and IV-pre-mining vertical stress zone. Gob-side entry is located within zone I to avoid a high stress concentration in zones II and III, and to avoid a large size pillar in zone IV, which would result in low recovery. Gob-side entry makes use of the side abutment pressure redistribution to locate gob-side entry in the de-stressed zone within the solid coal, i.e., zone I in Figure 3.
Energies 2018, 11, x FOR PEER REVIEW 2 of 24 retained as the entry for a future panel by constructing an artificial wall along the gob edge of the previous mined-out panel [3,4]. After excavation of a panel, the stress on the solid coal on the right ( Figure 3) can be divided into four zones [2,5]: I-de-stressed yield zone, II-over-stressed plastic zone, III-over-stressed elastic zone, and IV-pre-mining vertical stress zone. Gob-side entry is located within zone I to avoid a high stress concentration in zones II and III, and to avoid a large size pillar in zone IV, which would result in low recovery. Gob-side entry makes use of the side abutment pressure redistribution to locate gobside entry in the de-stressed zone within the solid coal, i.e., zone I in Figure 3. Yan [3] stated that the de-stressed zone occurs 0-7 m and the main roof will form the stable masonry structure to provide a low-stress environment for gob-side entry. He suggested driving the gob-side entry along the gob with a 5 m slender pillar. Li [6] noted that a gob-side entry with a 6 m wide coal pillar based on numerical modeling was employed in LTCC in mining a 15 Figure 3. Stress distribution after extraction of an LTCC panel [2].
Yan [3] stated that the de-stressed zone occurs 0-7 m and the main roof will form the stable masonry structure to provide a low-stress environment for gob-side entry. He suggested driving the gob-side entry along the gob with a 5 m slender pillar. Li [6] noted that a gob-side entry with a 6 m wide coal pillar based on numerical modeling was employed in LTCC in mining a 15 m thick coal seam in the Tashan mine of Datong Coal Group in China. Li [7], on the other hand, argued that high pressure would also occur on the edge of the solid coal and a structural model was built, as shown in Figure 4. A formula for computing the width of "inner stress field" S 1 was derived: where L is the panel width, (m); h i is the thickness of the ith roof main stratum, (m); γ is the unit weight of the main roof, (N/m 3 ); σ t is the tensile strength, (Pa); G 0 is the stiffness of coal on the left of caving line, (Pa); y 0 is the average compression of the slender pillar; (q n ) 1 is the loading of the n th roof strata in ascending order above the key stratum (including the dead weight of the lowest key stratum), and (q n ) 1 is expressed as: where E i is modulus of elasticity of individual layers, [Pa].  [7], on the other hand, argued that high pressure would also occur on the edge of the solid coal and a structural model was built, as shown in Figure 4. A formula for computing the width of "inner stress field" S1 was derived: where L is the panel width, (m); hi is the thickness of the ith roof main stratum, (m); γ is the unit weight of the main roof, (N/m 3 ); σt is the tensile strength, (Pa); G0 is the stiffness of coal on the left of caving line, (Pa); y0 is the average compression of the slender pillar; (qn)1 is the loading of the n th roof strata in ascending order above the key stratum (including the dead weight of the lowest key stratum), and (qn)1 is expressed as: where Ei is modulus of elasticity of individual layers, [Pa]. It was demonstrated that stability and load-bearing characteristics of the slender pillar depend on several factors, including mechanical properties of the roof strata, cover depth, panel geometry, etc. Yan [3] further indicated that the bearing capacity of the slender pillar would decrease because both edges of the pillar are fractured. Due to the effect of moving the front abutment pressure with the advance of an active panel, the pillar along the gob-side would be more deformed and fractured. With the fracturing of the main roof, the pressure and deformation of host rock will increase rapidly and it would be more difficult to keep the gob-side entry stable. Therefore, the slender coal pillar is still subjected to the loading resulting from the fractured main roof, even though the main roof and solid coal (where they are touching) is more or less compressed, deformed, or crushed, as shown in Figure 2a.
In order to increase the stability and load-bearing capability of the gob-side entry, several different types of artificial walls were proposed [3,[8][9][10][11]. Many destressing technologies for the gobside entry, along with supplementary supports were suggested [12]. Zhang [8] proposed a roadside backfill body (RBB) using high-water quick-setting materials for highly-gassy multi-entry for LTCC. RBB thicknesses of 1.0 m, 1.2 m, 1.5 m, 1.8 m, and 2.0 m were analyzed using numerical models and it was eventually determined that the best choice was 1.5 m. Yang [9] proposed a similar approach. Zhang [10,11] presented an artificial wall built along the gob using a cementing material, fly ash, detrital refuse, and concrete. Tan [3] proposed a new type of artificial composite wall that consists of an upper flexible layer and a lower hard layer which can sustain large compression and provide strong supporting intensity.
Overall, Zhang's [8] and Zhang's [10,11] methods are positive approaches which intend to provide high abutment-load carrying capacity of the artificial pillar, while Tan's [3] artificial pillar is a "semi-positive" and "semi-passive" approach, the pillar can deform to some extent with the deformation and rotation of the roof beam structure. No matter what pillar it is, natural or artificial, It was demonstrated that stability and load-bearing characteristics of the slender pillar depend on several factors, including mechanical properties of the roof strata, cover depth, panel geometry, etc. Yan [3] further indicated that the bearing capacity of the slender pillar would decrease because both edges of the pillar are fractured. Due to the effect of moving the front abutment pressure with the advance of an active panel, the pillar along the gob-side would be more deformed and fractured. With the fracturing of the main roof, the pressure and deformation of host rock will increase rapidly and it would be more difficult to keep the gob-side entry stable. Therefore, the slender coal pillar is still subjected to the loading resulting from the fractured main roof, even though the main roof and solid coal (where they are touching) is more or less compressed, deformed, or crushed, as shown in Figure 2a.
In order to increase the stability and load-bearing capability of the gob-side entry, several different types of artificial walls were proposed [3,[8][9][10][11]. Many destressing technologies for the gob-side entry, along with supplementary supports were suggested [12]. Zhang [8] proposed a roadside backfill body (RBB) using high-water quick-setting materials for highly-gassy multi-entry for LTCC. RBB thicknesses of 1.0 m, 1.2 m, 1.5 m, 1.8 m, and 2.0 m were analyzed using numerical models and it was eventually determined that the best choice was 1.5 m. Yang [9] proposed a similar approach. Zhang [10,11] presented an artificial wall built along the gob using a cementing material, fly ash, detrital refuse, and concrete. Tan [3] proposed a new type of artificial composite wall that consists of an upper flexible layer and a lower hard layer which can sustain large compression and provide strong supporting intensity.
Overall, Zhang's [8] and Zhang's [10,11] methods are positive approaches which intend to provide high abutment-load carrying capacity of the artificial pillar, while Tan's [3] artificial pillar is a "semi-positive" and "semi-passive" approach, the pillar can deform to some extent with the deformation and rotation of the roof beam structure. No matter what pillar it is, natural or artificial, the optimum width of the slender pillar or artificial wall should be designed based on the expected abutment pressure [2].
The gob-side entry designed above will still cause gateroads and slender pillars to be subjected to a complex dynamic loading environment during their service life [13]. Conventional roof supports may not be adequate to meet the requirements for mining operations and supplementary supports, such as ladder beam, counter-pulled bolts, and variations of pipe-shed support systems must be used for ground control, which can result in high production costs. Designing a gob-side entry still remains relatively complicated and challenging at present. This paper proposes a new and relatively simple approach for gob-side entry design for LTCC. In this approach to the gob-side entry, mining-induced stress distribution should be made most of, the complex dynamic loading environment must be avoided, and support design should be simpler.
The arrangement is listed as follows: In Section 2, the geological, geotechnical settings, and mining conditions and circumstances are presented. In Section 3, the new approach is introduced. In Section 4, numerical modelling, incorporating the angle of break and a double-yield model of the gob model is then carried out. In Section 5, field testing and observation is conducted. In Section 6, pros and cons of the approach are discussed. Conclusions are listed according to the study of the paper in Section 7.

Engineering Background of the Case Study Mine
The analysis is based on the mining conditions of the 18111 and 18111-1 panels of the #8 coal seam in Zhenchengdi Colliery, Shanxi Province, China. The panels are 600 m long along the strike and 120 m wide along the dip. The average coal seam thickness is 5 m thick. The average cover depth and seam dip angle are 200 m and 8 • , respectively. The generalized stratigraphy for panel 18111-1 is shown in Figure 5 in descending order. The mining height is 2.3 m, and the caving height is 2.7 m. All entry cross-sections are trapezoidal, 2.8 m wide on the top, 3.6 m wide on the bottom, and 2.5 m high (Figure 6a,b). The major tectonic structure in the area is a simple mono-synclinal structure.
Energies 2018, 11, x FOR PEER REVIEW 4 of 24 the optimum width of the slender pillar or artificial wall should be designed based on the expected abutment pressure [2]. The gob-side entry designed above will still cause gateroads and slender pillars to be subjected to a complex dynamic loading environment during their service life [13]. Conventional roof supports may not be adequate to meet the requirements for mining operations and supplementary supports, such as ladder beam, counter-pulled bolts, and variations of pipe-shed support systems must be used for ground control, which can result in high production costs. Designing a gob-side entry still remains relatively complicated and challenging at present. This paper proposes a new and relatively simple approach for gob-side entry design for LTCC. In this approach to the gob-side entry, mining-induced stress distribution should be made most of, the complex dynamic loading environment must be avoided, and support design should be simpler.
The arrangement is listed as follows: In Section 2, the geological, geotechnical settings, and mining conditions and circumstances are presented. In Section 3, the new approach is introduced. In Section 4, numerical modelling, incorporating the angle of break and a double-yield model of the gob model is then carried out. In Section 5, field testing and observation is conducted. In Section 6, pros and cons of the approach are discussed. Conclusions are listed according to the study of the paper in Section 7.

Engineering Background of the Case Study Mine
The analysis is based on the mining conditions of the 18111 and 18111-1 panels of the #8 coal seam in Zhenchengdi Colliery, Shanxi Province, China. The panels are 600 m long along the strike and 120 m wide along the dip. The average coal seam thickness is 5 m thick. The average cover depth and seam dip angle are 200 m and 8°, respectively. The generalized stratigraphy for panel 18111-1 is shown in Figure 5 in descending order. The mining height is 2.3 m, and the caving height is 2.7 m. All entry cross-sections are trapezoidal, 2.8 m wide on the top, 3.6 m wide on the bottom, and 2.5 m high (Figure 6a,b). The major tectonic structure in the area is a simple mono-synclinal structure.

A New Gob-Side Entry (NGE)
The proposed new gob-side entry (NGE) is illustrated in Figure 7. The development entries on either end of the panel are located at different elevations [14]. For the 18111 panel, the tailgate is driven along the floor while the maingate is driven along the roof. Thus, the longwall face has a gradually elevating section on one end of the 18111 panel. The elevating section slope is incrementally increased by adjusting the inclination of each section of the armored face conveyor (AFC), shields, and other production machines. Figure 7c shows that the length and inclination of each AFC section are typically 1.5 m and 3° [15]. For the adjacent 18111-1 panel, the tailgate, i.e., the gob-side entry, is located at the floor level and is offset horizontally with respect to the maingate of the 18111 panel. Similarly, the maingate of the 18111-1 panel is driven along the roof, thus forming the same geometry with the 18111 panel. The two adjacent panels are extracted sequentially. Therefore, there are several locations for the gob-side entry, as shown in Figure 8. The case study mine decided to select Entry 1 to extract more coal. The stress distribution for this new configuration is illustrated in Figure 9.

A New Gob-Side Entry (NGE)
The proposed new gob-side entry (NGE) is illustrated in Figure 7. The development entries on either end of the panel are located at different elevations [14]. For the 18111 panel, the tailgate is driven along the floor while the maingate is driven along the roof. Thus, the longwall face has a gradually elevating section on one end of the 18111 panel. The elevating section slope is incrementally increased by adjusting the inclination of each section of the armored face conveyor (AFC), shields, and other production machines. Figure 7c shows that the length and inclination of each AFC section are typically 1.5 m and 3 • [15]. For the adjacent 18111-1 panel, the tailgate, i.e., the gob-side entry, is located at the floor level and is offset horizontally with respect to the maingate of the 18111 panel. Similarly, the maingate of the 18111-1 panel is driven along the roof, thus forming the same geometry with the 18111 panel. The two adjacent panels are extracted sequentially. Therefore, there are several locations for the gob-side entry, as shown in Figure 8. The case study mine decided to select Entry 1 to extract more coal. The stress distribution for this new configuration is illustrated in Figure 9.
Compared with the conventional side abutment pressure distribution in LTCC (Figure 3), two additional zones are observed for this new geometry: zone V, stress-free zone; and zone VI, de-stressed gob edge zone. Typically, 3-6 m away from the gob edge to the gob, caved rocks are hardly seen. That is the reason for naming zone V the stress-free zone, and the concept is similar to the wisdom that the safest place in a room during an earthquake is a corner where a bridge protective structure is likely to form.   Compared with the conventional side abutment pressure distribution in LTCC (Figure 3), two additional zones are observed for this new geometry: zone V, stress-free zone; and zone VI, destressed gob edge zone. Typically, 3-6 m away from the gob edge to the gob, caved rocks are hardly seen. That is the reason for naming zone V the stress-free zone, and the concept is similar to the wisdom that the safest place in a room during an earthquake is a corner where a bridge protective structure is likely to form.     Compared with the conventional side abutment pressure distribution in LTCC (Figure 3), two additional zones are observed for this new geometry: zone V, stress-free zone; and zone VI, destressed gob edge zone. Typically, 3-6 m away from the gob edge to the gob, caved rocks are hardly seen. That is the reason for naming zone V the stress-free zone, and the concept is similar to the wisdom that the safest place in a room during an earthquake is a corner where a bridge protective structure is likely to form.     Compared with the conventional side abutment pressure distribution in LTCC (Figure 3), two additional zones are observed for this new geometry: zone V, stress-free zone; and zone VI, destressed gob edge zone. Typically, 3-6 m away from the gob edge to the gob, caved rocks are hardly seen. That is the reason for naming zone V the stress-free zone, and the concept is similar to the wisdom that the safest place in a room during an earthquake is a corner where a bridge protective structure is likely to form.  The proposed NGE approach has many advantages: (i) increased recovery of over 10% [16]; (ii) reduced risk of spontaneous combustion [17]; (iii) the 18111-1 tailgate (the gob-side entry) is located in Energies 2018, 11, 1292 7 of 24 the de-stressed zone along the gob edge, i.e., zone V in Figure 9, and this zone is de-stressed and is only subjected to the load of caved immediate roof fragments below the main roof protective structure after the 18111 panel is mined out; and (iv) the complexity in determining the size of the slender pillar is avoided, and no artificial wall needs to be constructed, which reduces the mining cost.
Many researchers have focused on locating the gob-side development entry in zone I shown in Figure 3. However, another de-stressed zone was overlooked and that was the gob. Yavuz [18] derived the formula for the stress at a distance x in the gob below: where σ t is the compressive strength (MPa); b is the bulking factor of the caved roof; S x is the subsidence of surface at distance x (m); x is the horizontal distance from the origin (m), which is located at the center of the subsidence profile; and c 3 and c 4 are coefficients depending on lithology.
If the vertical expansion and sagging of the strata in the fractured zone are ignored, then Equation (3) is reduced to: Based on these equations, an in-depth study of a number of parameters, including the bulking factor of the caved material, the strength of immediate and main roof strata, the excavation height, etc., was performed [18]. The results suggest that for all variations of those parameters, the stress for the gob edge is always less than any other parts in the entire panel layout system and is far less than the pre-mining stress.
In the NGE, stress distribution for the 18111-1 panel after the excavation of the 18111 panel and is shown in Figure 10 (after Brady [19]) and Figure 11. As 18111-1 tailgate inner offset more than one entry with respect to the 18111 maingate (Figure 7b), the solid coal mass on the right of the 18111 maingate influenced by the side abutment pressure is located approximately 4 m away from the gob-side entry. Hence, the side abutment pressure "steps over" the gob-side entry and the stress concentration that occurs at the T-junction due to the superposition of the front abutment, and side abutment pressure is located within the solid coal of the 18111-1 panel. Thus, the 18111-1 tailgate is independent of both side abutment pressure and front abutment pressure. The proposed NGE approach has many advantages: (i) increased recovery of over 10% [16]; (ii) reduced risk of spontaneous combustion [17]; (iii) the 18111-1 tailgate (the gob-side entry) is located in the de-stressed zone along the gob edge, i.e., zone V in Figure 9, and this zone is de-stressed and is only subjected to the load of caved immediate roof fragments below the main roof protective structure after the 18111 panel is mined out; and (iv) the complexity in determining the size of the slender pillar is avoided, and no artificial wall needs to be constructed, which reduces the mining cost.
Many researchers have focused on locating the gob-side development entry in zone I shown in Figure 3. However, another de-stressed zone was overlooked and that was the gob. Yavuz [18] derived the formula for the stress at a distance x in the gob below: where σt is the compressive strength (MPa); b is the bulking factor of the caved roof; Sx is the subsidence of surface at distance x (m); x is the horizontal distance from the origin (m), which is located at the center of the subsidence profile; and c3 and c4 are coefficients depending on lithology.
If the vertical expansion and sagging of the strata in the fractured zone are ignored, then Equation (3) is reduced to: Based on these equations, an in-depth study of a number of parameters, including the bulking factor of the caved material, the strength of immediate and main roof strata, the excavation height, etc., was performed [18]. The results suggest that for all variations of those parameters, the stress for the gob edge is always less than any other parts in the entire panel layout system and is far less than the pre-mining stress.
In the NGE, stress distribution for the 18111-1 panel after the excavation of the 18111 panel and is shown in Figure 10 (after Brady [19]) and Figure 11. As 18111-1 tailgate inner offset more than one entry with respect to the 18111 maingate (Figure 7b), the solid coal mass on the right of the 18111 maingate influenced by the side abutment pressure is located approximately 4 m away from the gobside entry. Hence, the side abutment pressure "steps over" the gob-side entry and the stress concentration that occurs at the T-junction due to the superposition of the front abutment, and side abutment pressure is located within the solid coal of the 18111-1 panel. Thus, the 18111-1 tailgate is independent of both side abutment pressure and front abutment pressure.  The fall height for the caving rocks above the elevating section is relatively small compared with the other parts of the gob. Esterhuizen [20] indicated that, as the fall height reduces, the potential for fragments to rotate diminishes and the bulking factor of the caved rocks is reduced. Thus, the caved rock fragments have a low void ratio and are able to accept the loads more quickly [21]. Therefore, pressure transferred to the abutment close to the elevating section should be less, resulting in lower side abutment pressure for the adjacent panel. Since the geometry of the panel is asymmetrical, the stress distribution must be consequently asymmetrical. Side abutment pressure on the right side (red line) is gentler than that on the left (blue line).
On the other hand, compression arches commonly develop across longwall panels due to coal extraction, provided the panel is narrow enough or deep enough for both ends of the arch to span the panel width and bear on rock. These arches are zones of tangential compressive stress where some of the weight of the overburden overlying the arch can be transferred onto the abutments. Some, or all, of the downward deflected rock under the arch will bear on the collapsed rock under the arch. If the width of a longwall panel is too wide or too shallow for the stress arch to span the panel width, a smaller stress arch will form, with one end of the arch bearing on and compressing the collapsed gob [22]. However, since coal mines in China are getting deeper, the critical panel width for most longwall panels is not reached. Therefore, the full weight of the overburden strata is not transferred onto the gob due to bridging of strata between abutments [18]. That is not the case for NGE, fortunately. Since conventional rectangular pillars between adjacent panels are eliminated, the roof strata above the adjacent panels are connected and deform as a single unit similar to a super-critical LTCC panel, as shown in Figure 12. Caved rock fragments are more compacted and more overburden weight is picked up by the gob rather than the abutments leading to lower side abutment pressure for the adjacent active or future panel. Lastly, different locations for the gob-side entry generate different ground control issues. The effects of the different development entry locations on the stress distribution are shown in Figure 13.  The fall height for the caving rocks above the elevating section is relatively small compared with the other parts of the gob. Esterhuizen [20] indicated that, as the fall height reduces, the potential for fragments to rotate diminishes and the bulking factor of the caved rocks is reduced. Thus, the caved rock fragments have a low void ratio and are able to accept the loads more quickly [21]. Therefore, pressure transferred to the abutment close to the elevating section should be less, resulting in lower side abutment pressure for the adjacent panel. Since the geometry of the panel is asymmetrical, the stress distribution must be consequently asymmetrical. Side abutment pressure on the right side (red line) is gentler than that on the left (blue line).
On the other hand, compression arches commonly develop across longwall panels due to coal extraction, provided the panel is narrow enough or deep enough for both ends of the arch to span the panel width and bear on rock. These arches are zones of tangential compressive stress where some of the weight of the overburden overlying the arch can be transferred onto the abutments. Some, or all, of the downward deflected rock under the arch will bear on the collapsed rock under the arch. If the width of a longwall panel is too wide or too shallow for the stress arch to span the panel width, a smaller stress arch will form, with one end of the arch bearing on and compressing the collapsed gob [22]. However, since coal mines in China are getting deeper, the critical panel width for most longwall panels is not reached. Therefore, the full weight of the overburden strata is not transferred onto the gob due to bridging of strata between abutments [18]. That is not the case for NGE, fortunately. Since conventional rectangular pillars between adjacent panels are eliminated, the roof strata above the adjacent panels are connected and deform as a single unit similar to a super-critical LTCC panel, as shown in Figure 12. Caved rock fragments are more compacted and more overburden weight is picked up by the gob rather than the abutments leading to lower side abutment pressure for the adjacent active or future panel. Lastly, different locations for the gob-side entry generate different ground control issues. The effects of the different development entry locations on the stress distribution are shown in Figure 13. The fall height for the caving rocks above the elevating section is relatively small compared with the other parts of the gob. Esterhuizen [20] indicated that, as the fall height reduces, the potential for fragments to rotate diminishes and the bulking factor of the caved rocks is reduced. Thus, the caved rock fragments have a low void ratio and are able to accept the loads more quickly [21]. Therefore, pressure transferred to the abutment close to the elevating section should be less, resulting in lower side abutment pressure for the adjacent panel. Since the geometry of the panel is asymmetrical, the stress distribution must be consequently asymmetrical. Side abutment pressure on the right side (red line) is gentler than that on the left (blue line).
On the other hand, compression arches commonly develop across longwall panels due to coal extraction, provided the panel is narrow enough or deep enough for both ends of the arch to span the panel width and bear on rock. These arches are zones of tangential compressive stress where some of the weight of the overburden overlying the arch can be transferred onto the abutments. Some, or all, of the downward deflected rock under the arch will bear on the collapsed rock under the arch. If the width of a longwall panel is too wide or too shallow for the stress arch to span the panel width, a smaller stress arch will form, with one end of the arch bearing on and compressing the collapsed gob [22]. However, since coal mines in China are getting deeper, the critical panel width for most longwall panels is not reached. Therefore, the full weight of the overburden strata is not transferred onto the gob due to bridging of strata between abutments [18]. That is not the case for NGE, fortunately. Since conventional rectangular pillars between adjacent panels are eliminated, the roof strata above the adjacent panels are connected and deform as a single unit similar to a super-critical LTCC panel, as shown in Figure 12. Caved rock fragments are more compacted and more overburden weight is picked up by the gob rather than the abutments leading to lower side abutment pressure for the adjacent active or future panel. Lastly, different locations for the gob-side entry generate different ground control issues. The effects of the different development entry locations on the stress distribution are shown in Figure 13.    According to the slip line field theory [23], the floor in Zone I in Figure 14a is in the active Rankine state, which tends to move toward the excavation. The angle formed by the horizontal line and slip-line is 45° +/− φ/2, with the major principal stress depending upon the sign convention; Zone III is in the active Rankine state, the angle formed by the horizontal line and slip-line is 45° −/+ φ/2; Zone II is the transition zone and the curve "AB" may be taken as the logarithmic spiral with the equation: where, r is the radius of the spiral, r0 is the length of OA (i.e., starting radius at θ = 0) (m); θ is the angle formed by OA and OB (i.e., r and r0) (°); φ is angle of internal friction of roof strata (°); thus: From the geometric equivalence relation: By substituting Equation (7) into Equation (6): , h reaches the maximum, hence: (9) that is: Thus; As: According to the slip line field theory [23], the floor in Zone I in Figure 14a is in the active Rankine state, which tends to move toward the excavation. The angle formed by the horizontal line and slip-line is 45 • +/− φ/2, with the major principal stress depending upon the sign convention; Zone III is in the active Rankine state, the angle formed by the horizontal line and slip-line is 45 • −/+ φ/2; Zone II is the transition zone and the curve "AB" may be taken as the logarithmic spiral with the equation: where, r is the radius of the spiral, r 0 is the length of OA (i.e., starting radius at θ = 0) (m); θ is the angle formed by OA and OB (i.e., r and r 0 ) ( • ); φ is angle of internal friction of roof strata ( • ); thus: From the geometric equivalence relation: By substituting Equation (7) into Equation (6): When dh dθ = 0, h reaches the maximum, hence: that is: Thus; As: Hence, the maximum depth of the yield zone is: According to [23] and the above analysis, for a cover depth of 200 m (the average cover depth of the case study longwall panels), the stress concentration distribution within the floor strata is shown in Figure 14b. It shows that the floor strata on the gob side are in a destressed state, while the floor strata on the solid coal side are in a stressed state. The NGE is located right within the destressed zone of the floor strata. Therefore, the stress environment for the NGE is improved.
According to [23] and the above analysis, for a cover depth of 200 m (the average cover depth of the case study longwall panels), the stress concentration distribution within the floor strata is shown in Figure 14b. It shows that the floor strata on the gob side are in a destressed state, while the floor strata on the solid coal side are in a stressed state. The NGE is located right within the destressed zone of the floor strata. Therefore, the stress environment for the NGE is improved. A structural model developed for the roof strata movement using the NGE approach is shown in Figure 15. The top part of the conventional triangular slender pillar seems to be struck off creating a low-stress interspace for the stress and energy around NGE to dissipate, and the pressure of the masonry beam is transferred to the gob and abutment to the right of the NGE. Thus, the 18111-1 tailgate is only subjected to the weight of caved rock fragments near the gob edge, which is very low. An additional advantage is that the inner stress field discussed by Li [7] (Figure 4) that was caused by the deformed contact, is now located in the solid coal mass on the right of the 18111 maingate rather than next to the gob-side entry. A structural model developed for the roof strata movement using the NGE approach is shown in Figure 15. The top part of the conventional triangular slender pillar seems to be struck off creating a low-stress interspace for the stress and energy around NGE to dissipate, and the pressure of the masonry beam is transferred to the gob and abutment to the right of the NGE. Thus, the 18111-1 tailgate is only subjected to the weight of caved rock fragments near the gob edge, which is very low. An additional advantage is that the inner stress field discussed by Li [7] (Figure 4) that was caused by the deformed contact, is now located in the solid coal mass on the right of the 18111 maingate rather than next to the gob-side entry.
The horizontal span L b of the Block B [24] is related to the periodic roof weighting interval of the main roof L p and the panel width L, which is expressed by: Energies 2018, 11, 1292 11 of 24 In conventional LTCC, since a rectangular pillar is left between adjacent panels, the mining system for each panel is isolated. Therefore, before the first periodic weighting, the main roof strata acts like a plate clamped rigidly at four edges, as shown in Figure 16a. After the first periodic weighting, the main roof strata acts like a plate clamped rigidly at three edges with one free edge behind the working face, as shown in Figure 16b. However, for the 18111-1 panel using the NGE, before the first weighting, the edge of the 18111 panel side is free, as shown in Figure 16c. After the first weighting, the main roof strata acts like a plate clamped rigidly at two edges with the other two free edges (one on the 18111 panel side and one behind the working face), as shown in Figure 16d. The horizontal span Lb of the Block B [24] is related to the periodic roof weighting interval of the main roof Lp and the panel width L, which is expressed by: In conventional LTCC, since a rectangular pillar is left between adjacent panels, the mining system for each panel is isolated. Therefore, before the first periodic weighting, the main roof strata acts like a plate clamped rigidly at four edges, as shown in Figure 16a. After the first periodic weighting, the main roof strata acts like a plate clamped rigidly at three edges with one free edge behind the working face, as shown in Figure 16b. However, for the 18111-1 panel using the NGE, before the first weighting, the edge of the 18111 panel side is free, as shown in Figure 16c. After the first weighting, the main roof strata acts like a plate clamped rigidly at two edges with the other two free edges (one on the 18111 panel side and one behind the working face), as shown in Figure 16d.   [2] and Tan [4]).
The horizontal span Lb of the Block B [24] is related to the periodic roof weighting interval of the main roof Lp and the panel width L, which is expressed by: In conventional LTCC, since a rectangular pillar is left between adjacent panels, the mining system for each panel is isolated. Therefore, before the first periodic weighting, the main roof strata acts like a plate clamped rigidly at four edges, as shown in Figure 16a. After the first periodic weighting, the main roof strata acts like a plate clamped rigidly at three edges with one free edge behind the working face, as shown in Figure 16b. However, for the 18111-1 panel using the NGE, before the first weighting, the edge of the 18111 panel side is free, as shown in Figure 16c. After the first weighting, the main roof strata acts like a plate clamped rigidly at two edges with the other two free edges (one on the 18111 panel side and one behind the working face), as shown in Figure 16d.

Numerical Modeling
A few researchers [8,9] studied the gob-side entry using numerical modeling and gob behavior was involved. Some researchers specifically studied the gob behavior [18,20,25,26], but none considered the influence of the angle of break. Programs used by researchers included such software

Numerical Modeling
A few researchers [8,9] studied the gob-side entry using numerical modeling and gob behavior was involved. Some researchers specifically studied the gob behavior [18,20,25,26], but none considered the influence of the angle of break. Programs used by researchers included such software as PAFEC (Program for Automatic Finite Element Calculations), FLAC (Fast Lagrangian Analysis of Continua), and FLAC 3D (Fast Lagrangian Analysis of Continua in 3 Dimensions). FLAC [27] which has a higher solution speed and lower memory requirement is used for two-dimensional plane-strain modeling of the NGE. The built-in double-yield constitutive model in FLAC represents a material in which there may be significant irreversible compaction in addition to shear yielding, such as hydraulically-placed backfill or lightly-cemented granular material [27].
The model grid and assigned materials are shown in Figure 17. The dimensions of the model are 380 m (length) × 75.71 m (height). At the top of the model, a uniform vertical stress of 162.08 m × 0.027 MPa/m = 4.38 MPa was applied to simulate the overburden weight. An angle of break of 30 • obtained by physical modeling was incorporated. The horizontal sides were roller constrained and the bottom boundary was fixed both horizontally and vertically. Mohr-Coulomb constitutive properties were used for different lithologies and the double-yield model was used to define the gob. Rock mass engineering properties used for numerical modeling are given in Table 1.  [27] which has a higher solution speed and lower memory requirement is used for two-dimensional plane-strain modeling of the NGE. The built-in double-yield constitutive model in FLAC represents a material in which there may be significant irreversible compaction in addition to shear yielding, such as hydraulically-placed backfill or lightly-cemented granular material [27]. The model grid and assigned materials are shown in Figure 17. The dimensions of the model are 380 m (length) × 75.71 m (height). At the top of the model, a uniform vertical stress of 162.08 m × 0.027 MPa/m = 4.38 MPa was applied to simulate the overburden weight. An angle of break of 30° obtained by physical modeling was incorporated. The horizontal sides were roller constrained and the bottom boundary was fixed both horizontally and vertically. Mohr-Coulomb constitutive properties were used for different lithologies and the double-yield model was used to define the gob. Rock mass engineering properties used for numerical modeling are given in Table 1.  Table 1).   Table 1). A double-yield model was used to simulate the compaction of the gob material in FLAC. Pappas and Mark [21] studied the behavior of simulated longwall gob material using laboratory tests and concluded that Salamon's equation closely follows the laboratory results. Therefore, the gob model calibrations were done [28] using equivalent gob elements that follow Salamon's model (Equation (15)):  Table 2 and are expressed by: A test was carried out to obtain a set of reasonable parameters of the double-yield constitutive gob model. A simple model with dimensions 1 m × 2 m was constructed. Loading was simulated by applying a velocity on the top surface with fixed sides along the X axis, and fixed boundaries both in the X and Y directions. The input parameters were fitted by an iterative change in the bulk and shear modulus, the angle of dilation, the angle of friction, and the density of the gob material. The volumetric strain, vertical stress contours, and FLAC numerical stress-strain curves are shown in Figure 18 to obtain the fitted parameters. The final properties are given in Table 3. A comparison between the numerically-obtained strain-hardening curve and Salamon's theoretical curve is shown in Figure 19. The interface parameters used in modeling were: normal stiffness of 1.0 × 10 9 Pa; shear stiffness of 1.0 × 10 9 Pa; cohesion of 1 × 10 5 Pa; angle of internal friction of 20 • ; and the shear bond is 100 times the tensile bond [27]. Table 3. Parameters for gob material.        Figure 20 shows vertical stress distribution after excavation of the 18111 panel. The distance between the gob edge and the location where pre-mining stress is reached in the gob is about 51 m. This agrees with the rule of thumb proposed by Wilson [29] that the pre-mining stress is achieved 0.2-0.3 times the panel overburden depth. It is also observed that the elevating section is located in the destressed zone. Therefore, the gateroads located here would suffer lower ground pressure and the ground control problems would be, consequently, minimal. According to records of several panels in the Zhenchengdi coal mine, dynamic disasters, such as rockbursts or coal and gas outbursts, never occurred. In addition, Figure 20b indicates that the abutment pressure on the side of the elevating section is smaller and the gob pressure increases faster. This means that the load borne by the gob is larger, especially for the elevating section.  Figure 20 shows vertical stress distribution after excavation of the 18111 panel. The distance between the gob edge and the location where pre-mining stress is reached in the gob is about 51 m. This agrees with the rule of thumb proposed by Wilson [29] that the pre-mining stress is achieved 0.2-0.3 times the panel overburden depth. It is also observed that the elevating section is located in the destressed zone. Therefore, the gateroads located here would suffer lower ground pressure and the ground control problems would be, consequently, minimal. According to records of several panels in the Zhenchengdi coal mine, dynamic disasters, such as rockbursts or coal and gas outbursts, never occurred. In addition, Figure 20b indicates that the abutment pressure on the side of the elevating section is smaller and the gob pressure increases faster. This means that the load borne by the gob is larger, especially for the elevating section.    Figure 20 shows vertical stress distribution after excavation of the 18111 panel. The distance between the gob edge and the location where pre-mining stress is reached in the gob is about 51 m. This agrees with the rule of thumb proposed by Wilson [29] that the pre-mining stress is achieved 0.2-0.3 times the panel overburden depth. It is also observed that the elevating section is located in the destressed zone. Therefore, the gateroads located here would suffer lower ground pressure and the ground control problems would be, consequently, minimal. According to records of several panels in the Zhenchengdi coal mine, dynamic disasters, such as rockbursts or coal and gas outbursts, never occurred. In addition, Figure 20b indicates that the abutment pressure on the side of the elevating section is smaller and the gob pressure increases faster. This means that the load borne by the gob is larger, especially for the elevating section.   Readers may also note that the gob pressure distributions of physical modeling and numerical modeling are different. This is because that the roof subsidence does not reach the surface in physical modeling, and not all the overburden weight is sustained by the gob at the center of the panel, rather, the pressure is transferred to the abutments. This is also demonstrated by the fact that the abutment pressure in physical modeling is larger than that in numerical modeling.

Support Design for the New Gob-Side Entry
The layout of the gob-side entry determines its unique support system. For the NGE, experience from the multi-slicing longwall mining technology was used to decide on the support system.
During the extraction of the 18111 panel, #10 steel wire mesh was used along the elevating section of the 18111 panel (about 8-10 m wide). The wire mesh allows the development of the artificial roof for the gob-side entry of 18111-1 panel. Over six months after the extraction of the 18111 panel, the gob-side entry was excavated. An ELMB-75C machine (Hang Tian Chen Guang Co., Ltd., Nanjing, China) was employed for tunneling along with SGW-40T (Hua Yun An Te Co., Ltd., Beijing, China) type conveyor for transportation of coal. Each cut width was 0.8 m, and the cutting height was controlled to avoid damaging the roof steel wire mesh. Trapezoidal steel yieldable sets made of #11 I-shaped mine steel spaced at 800 mm were used for roadway support. The leg height was 2700 mm, and the top beam was 3100 mm long. Four laggings were used on the top pf the steel sets, and three laggings for the coal ribs. Figure 23 presents the detailed design. Readers may also note that the gob pressure distributions of physical modeling and numerical modeling are different. This is because that the roof subsidence does not reach the surface in physical modeling, and not all the overburden weight is sustained by the gob at the center of the panel, rather, the pressure is transferred to the abutments. This is also demonstrated by the fact that the abutment pressure in physical modeling is larger than that in numerical modeling.

Support Design for the New Gob-Side Entry
The layout of the gob-side entry determines its unique support system. For the NGE, experience from the multi-slicing longwall mining technology was used to decide on the support system.
During the extraction of the 18111 panel, #10 steel wire mesh was used along the elevating section of the 18111 panel (about 8-10 m wide). The wire mesh allows the development of the artificial roof for the gob-side entry of 18111-1 panel. Over six months after the extraction of the 18111 panel, the gob-side entry was excavated. An ELMB-75C machine (Hang Tian Chen Guang Co., Ltd., Nanjing, China) was employed for tunneling along with SGW-40T (Hua Yun An Te Co., Ltd., Beijing, China) type conveyor for transportation of coal. Each cut width was 0.8 m, and the cutting height was controlled to avoid damaging the roof steel wire mesh. Trapezoidal steel yieldable sets made of #11 I-shaped mine steel spaced at 800 mm were used for roadway support. The leg height was 2700 mm, and the top beam was 3100 mm long. Four laggings were used on the top pf the steel sets, and three laggings for the coal ribs. Figure 23 presents the detailed design. Readers may also note that the gob pressure distributions of physical modeling and numerical modeling are different. This is because that the roof subsidence does not reach the surface in physical modeling, and not all the overburden weight is sustained by the gob at the center of the panel, rather, the pressure is transferred to the abutments. This is also demonstrated by the fact that the abutment pressure in physical modeling is larger than that in numerical modeling.

Support Design for the New Gob-Side Entry
The layout of the gob-side entry determines its unique support system. For the NGE, experience from the multi-slicing longwall mining technology was used to decide on the support system.
During the extraction of the 18111 panel, #10 steel wire mesh was used along the elevating section of the 18111 panel (about 8-10 m wide). The wire mesh allows the development of the artificial roof for the gob-side entry of 18111-1 panel. Over six months after the extraction of the 18111 panel, the gob-side entry was excavated. An ELMB-75C machine (Hang Tian Chen Guang Co., Ltd., Nanjing, China) was employed for tunneling along with SGW-40T (Hua Yun An Te Co., Ltd., Beijing, China) type conveyor for transportation of coal. Each cut width was 0.8 m, and the cutting height was controlled to avoid damaging the roof steel wire mesh. Trapezoidal steel yieldable sets made of #11 I-shaped mine steel spaced at 800 mm were used for roadway support. The leg height was 2700 mm, and the top beam was 3100 mm long. Four laggings were used on the top pf the steel sets, and three laggings for the coal ribs. Figure 23 presents the detailed design. The support system should have capacity greater than the following equation to ensure stability: where HI is the height of the immediate roof, and γI is the unit weight of the immediate roof, 0.027 MN/m 3 (note that immediate roof here is a general definition which refers to the one or several roof strata that will become fragmented rocks within the caved zone).
From an engineering background we know that HI = 17.9 m, thus, PT = 0.48 MPa, which agrees very well with numerical modeling data that the vertical stress is only about 0.5 MPa and the stress concentration factor is less than 0.1. The numerical global model is, therefore, validated as well. A safety factor of 3 is used and, therefore, PT = 0.15 is required.

Support Resistance of Shields
Along the 18111-1 working face, the gob-side entry side to the maingate side, No. 5, 30, 59, and 76 shields located right over the elevating section were installed with two monitoring recorders for each shield. The number of support cycles varied from 6-8 each day with 0.6 m of face advance for each cycle. Due to large quantity of data, only 292 cycles of data were analyzed. Initial setting pressure P0, final pressure Pm, time-weighted average resistance (TWAR) Pt [30], and Pm' vs. the number of shield supporting cycles together with the amount of face advance are plotted in Figure 24, as well as an illustration of the locations of the monitored shields and: The support system should have capacity greater than the following equation to ensure stability: where H I is the height of the immediate roof, and γ I is the unit weight of the immediate roof, 0.027 MN/m 3 (note that immediate roof here is a general definition which refers to the one or several roof strata that will become fragmented rocks within the caved zone).
From an engineering background we know that H I = 17.9 m, thus, P T = 0.48 MPa, which agrees very well with numerical modeling data that the vertical stress is only about 0.5 MPa and the stress concentration factor is less than 0.1. The numerical global model is, therefore, validated as well. A safety factor of 3 is used and, therefore, P T = 0.15 is required.

Support Resistance of Shields
Along the 18111-1 working face, the gob-side entry side to the maingate side, No. 5, 30, 59, and 76 shields located right over the elevating section were installed with two monitoring recorders for each shield. The number of support cycles varied from 6-8 each day with 0.6 m of face advance for each cycle. Due to large quantity of data, only 292 cycles of data were analyzed. Initial setting pressure P 0 , final pressure P m , time-weighted average resistance (TWAR) P t [30], and P m ' vs. the number of shield supporting cycles together with the amount of face advance are plotted in Figure 24, as well as an illustration of the locations of the monitored shields and: p t = 1 2 (p 0 + p 1 )t 1 + 1 2 (p 1 + p 2 )t 2 + · · · + 1 2 (p n−1 + p n )t n t 1 + t 2 + · · · + t n (18) where p m is the average TWAR, σ pm is the mean square error, and P 1 , P 2 , . . . , P n is the support resistance during the corresponding period of t 1 , t 2 , . . . , t n , as shown in Figure 25.  (19) where m p is the average TWAR, pm  is the mean square error, and P1, P2, …, Pn is the support resistance during the corresponding period of t1, t2, …, tn, as shown in Figure 25.
where m p is the average TWAR, pm  is the mean square error, and P1, P2, …, Pn is the support resistance during the corresponding period of t1, t2, …, tn, as shown in Figure 25.   Taking pressures larger than p m as the signs of periodic roof weighting, there was a periodic roof weighting every 15-20 cycles. Therefore, the periodic roof weighting interval varies from 9 to 12 m, and according to Equation (14):

Pressure Distribution in the Elevating Section
The gradually elevating section may be analyzed as a triangular section, as shown in Figure 26. The length of the elevating section is about 15 m, height of the entry is 2.5 m, thus, the inclination of the elevating section is arctan (2.5/15) = 9.46 • . Therefore, the mining height M h of the elevating section is estimated as: Energies 2018, 11, x FOR PEER REVIEW 19 of 24 Taking pressures larger than m p  as the signs of periodic roof weighting, there was a periodic roof weighting every 15-20 cycles. Therefore, the periodic roof weighting interval varies from 9 to 12 m, and according to Equation (14): Thus, the width of Block B was calculated as 10 to 13 m.

Pressure Distribution in the Elevating Section
The gradually elevating section may be analyzed as a triangular section, as shown in Figure 26. where σt is the compressive strength (MPa); b is the bulking factor of the caved roof; Sx is the subsidence of surface at distance x, (m); x is the horizontal distance from the origin, (m), which is located at the center of the subsidence profile; and c3 and c4 are the coefficients depending on lithology and l is the distance between a certain point along the bottom of the elevating section to the vertex, (m). Since subsidence data are currently not available, it is proposed that this equation be validated in future papers.

Pressure and Convergence in Gateroads of the 18111-1 Panel
Pressure and convergence monitoring were carried out in both gateroads, that is, 18111-1 tailgate (the gob-side entry) and 18111-1 maingate. Immediately after the steel sets were installed on drivage, automated recording load cells were inserted on the top of the steel sets ahead of the face within the gateroads to collect pressure data. The convergences of the gateroads were monitored, as shown in Figure 27. The locations of the measurement points were selected in the middle of the two adjacent 9.3°l  (21) where σ t is the compressive strength (MPa); b is the bulking factor of the caved roof; S x is the subsidence of surface at distance x, (m); x is the horizontal distance from the origin, (m), which is located at the center of the subsidence profile; and c 3 and c 4 are the coefficients depending on lithology and l is the distance between a certain point along the bottom of the elevating section to the vertex, (m).
Since subsidence data are currently not available, it is proposed that this equation be validated in future papers.

Pressure and Convergence in Gateroads of the 18111-1 Panel
Pressure and convergence monitoring were carried out in both gateroads, that is, 18111-1 tailgate (the gob-side entry) and 18111-1 maingate. Immediately after the steel sets were installed on drivage, automated recording load cells were inserted on the top of the steel sets ahead of the face within the gateroads to collect pressure data. The convergences of the gateroads were monitored, as shown in Figure 27. The locations of the measurement points were selected in the middle of the two adjacent steel sets where the roof was clearly exposed and ready to be marked. Permanent pins were installed in the roof, floor, pillar rib, and longwall panel rib for each measurement station, respectively. Convergences were measured by a flexible tape. Differences of "AB" and "CD" with respect to their initial values at a certain location represent rib-to-rib convergence and roof-to-floor convergence, respectively. The data obtained are plotted in Figures 28 and 29.
Energies 2018, 11, x FOR PEER REVIEW 20 of 24 steel sets where the roof was clearly exposed and ready to be marked. Permanent pins were installed in the roof, floor, pillar rib, and longwall panel rib for each measurement station, respectively. Convergences were measured by a flexible tape. Differences of "AB" and "CD" with respect to their initial values at a certain location represent rib-to-rib convergence and roof-to-floor convergence, respectively. The data obtained are plotted in Figures 28 and 29.   steel sets where the roof was clearly exposed and ready to be marked. Permanent pins were installed in the roof, floor, pillar rib, and longwall panel rib for each measurement station, respectively. Convergences were measured by a flexible tape. Differences of "AB" and "CD" with respect to their initial values at a certain location represent rib-to-rib convergence and roof-to-floor convergence, respectively. The data obtained are plotted in Figures 28 and 29.   steel sets where the roof was clearly exposed and ready to be marked. Permanent pins were installed in the roof, floor, pillar rib, and longwall panel rib for each measurement station, respectively. Convergences were measured by a flexible tape. Differences of "AB" and "CD" with respect to their initial values at a certain location represent rib-to-rib convergence and roof-to-floor convergence, respectively. The data obtained are plotted in Figures 28 and 29.   The support load in the 18111-1 tailgate (gob-side entry) is much lower than that in the 18111-1 maingate and less than the pre-mining stress. This corroborates that the gob-side entry employing the NGE is located within the destressed zone. Steel sets were spaced at 800 mm along the longitudinal direction of the entry, hence, 400 mm of the roof on each side is carried by one set, and the roof is 3.1 m long, thus, the weight carried by one set is 3.1 × 0.8 × 1.79 × 0.027 = 0.12 MN. As the flange width of #11 I-shaped steel is 90 mm, it is assumed that all the weight of the caved rock fragments of the immediate roof acts on the top beam, thus, the average stress on the top of the beam is 0.12/(0.09 × 3.1) = 0.43 MPa, which also agrees very well with the numerical modeling results. The average value is lower than 0.43 MPa. It is believed that this is due to the influence of the angle of break. This is similar to multi-slice longwall mining where the panel gateroads in the lower seam are located directly underneath the gob of the immediate super-adjacent seam. Thus, the gateroads in the lower slice are located within the destressed zone of the super-adjacent panel [30]. In addition, the data show that the roof pressure is maintained at about 0.42 MPa without much fluctuation. This proves that the gob-side entry is independent of the side and front abutment pressure, which validates the theoretical analysis and numerical modeling. Therefore, the NGE has significant advantages from a ground control point of view.
Convergence for the two gateroads of 18111-1 panel shows that roof-to-floor convergence is similar and is small. Rib-to-rib convergence indicates that the gob-side entry converged less laterally than the non-gob-side entry which means the lateral pressure acting on the gob-side entry is not large either. In fact, according to engineering experience, the most roof to floor convergence result from the deformation of the pseudo-roof. The value of the deformation depends on the meshing quality and support design of the support system, such as the interval of adjacent steel sets, the density of lagging, etc. Nearly no convergences were observed in the gob side entry of the 22204 panel at this coal mine. Therefore, it demonstrates that the overall stress environment in the gob-side entry employing the NGE is improved. Thus, ground control problems, such as roof falls, bumps, and outbursts, are minimized and operation costs should be lower. Tangshan, Huafeng, and Dayuan coal mines, among others, are using the technology for deep inclined thick coal seams [31][32][33].

Pros and Cons of the NGE
The two adjacent panels must be extracted sequentially. As the gob-side entry of the panel is driven directly under the artificial roof or pseudo-roof of the gob edge of the panel, the development of the gob-side entry should be done at least six months (preferably one year) after extraction of the 18111 panel. This time is required to allow the caved and fractured roof strata to compact [30]. This point must be considered when a coal mine plans to use the NGE approach. Fortunately, this case study coal mine uses a two-wing longwall system as shown in Figure 30. The panels are mined in numbered order. For instance, #1 panel is first mined, then #2 panel is mined. Thus, it has been over six months after the #1 panel is mined when the #3 panel starts to be mined. If a coal mine uses one-wing system, then mining every other panel can be adopted for NGE. The support load in the 18111-1 tailgate (gob-side entry) is much lower than that in the 18111-1 maingate and less than the pre-mining stress. This corroborates that the gob-side entry employing the NGE is located within the destressed zone. Steel sets were spaced at 800 mm along the longitudinal direction of the entry, hence, 400 mm of the roof on each side is carried by one set, and the roof is 3.1 m long, thus, the weight carried by one set is 3.1 × 0.8 × 1.79 × 0.027 = 0.12 MN. As the flange width of #11 I-shaped steel is 90 mm, it is assumed that all the weight of the caved rock fragments of the immediate roof acts on the top beam, thus, the average stress on the top of the beam is 0.12/(0.09 × 3.1) = 0.43 MPa, which also agrees very well with the numerical modeling results. The average value is lower than 0.43 MPa. It is believed that this is due to the influence of the angle of break. This is similar to multi-slice longwall mining where the panel gateroads in the lower seam are located directly underneath the gob of the immediate super-adjacent seam. Thus, the gateroads in the lower slice are located within the destressed zone of the super-adjacent panel [30]. In addition, the data show that the roof pressure is maintained at about 0.42 MPa without much fluctuation. This proves that the gob-side entry is independent of the side and front abutment pressure,s which validates the theoretical analysis and numerical modeling. Therefore, the NGE has significant advantages from a ground control point of view.
Convergence for the two gateroads of 18111-1 panel shows that roof-to-floor convergence is similar and is small. Rib-to-rib convergence indicates that the gob-side entry converged less laterally than the non-gob-side entry which means the lateral pressure acting on the gob-side entry is not large either. In fact, according to engineering experience, the most roof to floor convergence result from the deformation of the pseudo-roof. The value of the deformation depends on the meshing quality and support design of the support system, such as the interval of adjacent steel sets, the density of lagging, etc. Nearly no convergences were observed in the gob side entry of the 22204 panel at this coal mine. Therefore, it demonstrates that the overall stress environment in the gob-side entry employing the NGE is improved. Thus, ground control problems, such as roof falls, bumps, and outbursts, are minimized and operation costs should be lower. Tangshan, Huafeng, and Dayuan coal mines, among others, are using the technology for deep inclined thick coal seams [31][32][33].

Pros and Cons of the NGE
The two adjacent panels must be extracted sequentially. As the gob-side entry of the panel is driven directly under the artificial roof or pseudo-roof of the gob edge of the panel, the development of the gob-side entry should be done at least six months (preferably one year) after extraction of the 18111 panel. This time is required to allow the caved and fractured roof strata to compact [30]. This point must be considered when a coal mine plans to use the NGE approach. Fortunately, this case study coal mine uses a two-wing longwall system as shown in Figure 30. The panels are mined in numbered order. For instance, #1 panel is first mined, then #2 panel is mined. Thus, it has been over six months after the #1 panel is mined when the #3 panel starts to be mined. If a coal mine uses onewing system, then mining every other panel can be adopted for NGE.   Coal mines currently using the NGE are all single-entry or two-entry systems. NGE is in use only in China or Russia. This practice may not meet the statutory requirements in many countries around the world. Therefore, additional research should be done to assess the performance of three-, four-, and more-entry panel systems.
Leakage of air from the NGE into the gob is another area of concern. Therefore, the NGE is preferably used for coal seams with low methane content and gob water amount. For the case study coal mine, as the coal seam averages about 5 m, the gob-side entry is totally connected to the gob. By grouting slime water into the gob from the 18111 maingate to the gob edge during extraction of the 18111 panel, the degree of consolidation is increased so that the air leakage meets the ventilation requirement [15]. In addition, grouting slime reduces the spontaneous heating of loose coal in the gob area which means that the prevention and extinguishing of a fire caused by residual loose coal are also rather good. Other proper measures, such as pre-drainage of gas through gob wells or cross-measure drill holes, can also be adopted to control the leakage. Employing an outer offset configuration of the NGE with a slender pillar in between or leaving a coal sheet between the two adjacent panels is shown in Figure 31, then it can be applied for coal seams with a higher gas content. Coal mines currently using the NGE are all single-entry or two-entry systems. NGE is in use only in China or Russia. This practice may not meet the statutory requirements in many countries around the world. Therefore, additional research should be done to assess the performance of three-, four-, and more-entry panel systems.
Leakage of air from the NGE into the gob is another area of concern. Therefore, the NGE is preferably used for coal seams with low methane content and gob water amount. For the case study coal mine, as the coal seam averages about 5 m, the gob-side entry is totally connected to the gob. By grouting slime water into the gob from the 18111 maingate to the gob edge during extraction of the 18111 panel, the degree of consolidation is increased so that the air leakage meets the ventilation requirement [15]. In addition, grouting slime reduces the spontaneous heating of loose coal in the gob area which means that the prevention and extinguishing of a fire caused by residual loose coal are also rather good. Other proper measures, such as pre-drainage of gas through gob wells or crossmeasure drill holes, can also be adopted to control the leakage. Employing an outer offset configuration of the NGE with a slender pillar in between or leaving a coal sheet between the two adjacent panels is shown in Figure 31, then it can be applied for coal seams with a higher gas content.

Concluding Remarks
This paper has presented a novel design approach to gob-side entry which is relatively easy. In this approach, determining the size of the slender pillar for the gob-side entry is avoided and no artificial wall needs to be constructed. Development entries on either end of the panel employing the NGE approach are located at different vertical levels. The gob-side entry is located at the gob edge and is offset horizontally with respect to the previous mined-out panel. The basic mining principle is similar to multi-slice longwall mining which has been improved using the NGE.
We carried out the theoretical analysis, which shows that the gob-side entry in this approach is located within the most destressed zone in the entire panel system. It is almost in a stress-free state, which avoids dynamic disasters such as coal bumps and coal outbursts. Therefore, the support and maintenance requirements for the gob-side entry are small. The unique geometry also creates a smaller side abutment pressure for the adjacent active panel.
Numerical modeling incorporating the angle of break was carried out. The double-yield model for the gob was calibrated using a simple model by fitting gob strain-hardening behavior with Salamon's theoretical model. Pre-mining stress in the gob is achieved about 51 m away from the gob edge. Side abutment pressure within the solid coal close to the gob-side entry is smaller than the other side of the gob. The stress concentration factor for the roof on the gob-side entry is less than 0.1. The stress around the gob edge near the gob-side entry is relatively low compared with the other side of the gob of the 18111 panel. The destressing effect of the NGE is evident. Numerical modeling results corroborate the theoretical analyses.
We carried out field testing and field observations. Support design for the NGE is also presented. In summary, the theoretical analysis, numerical modeling, and field observations are consistent and validate the scientific foundation of the new technology.

Concluding Remarks
This paper has presented a novel design approach to gob-side entry which is relatively easy. In this approach, determining the size of the slender pillar for the gob-side entry is avoided and no artificial wall needs to be constructed. Development entries on either end of the panel employing the NGE approach are located at different vertical levels. The gob-side entry is located at the gob edge and is offset horizontally with respect to the previous mined-out panel. The basic mining principle is similar to multi-slice longwall mining which has been improved using the NGE.
We carried out the theoretical analysis, which shows that the gob-side entry in this approach is located within the most destressed zone in the entire panel system. It is almost in a stress-free state, which avoids dynamic disasters such as coal bumps and coal outbursts. Therefore, the support and maintenance requirements for the gob-side entry are small. The unique geometry also creates a smaller side abutment pressure for the adjacent active panel.
Numerical modeling incorporating the angle of break was carried out. The double-yield model for the gob was calibrated using a simple model by fitting gob strain-hardening behavior with Salamon's theoretical model. Pre-mining stress in the gob is achieved about 51 m away from the gob edge. Side abutment pressure within the solid coal close to the gob-side entry is smaller than the other side of the gob. The stress concentration factor for the roof on the gob-side entry is less than 0.1. The stress around the gob edge near the gob-side entry is relatively low compared with the other side of the gob of the 18111 panel. The destressing effect of the NGE is evident. Numerical modeling results corroborate the theoretical analyses.
We carried out field testing and field observations. Support design for the NGE is also presented. In summary, the theoretical analysis, numerical modeling, and field observations are consistent and validate the scientific foundation of the new technology.
Funding: This research was funded by the National Natural Science Foundation of China, Excellent Young Scientists, grant number 51422404 and the NSFC-Shanxi Joint Fund for Coal-Based Low-Carbon Technology, grant number U1710258.