Experiment and Application of Coalcrete on Roadway Stability: A Comparative Analysis

Improving roadway stability in deep underground mines is quite challenging, as the conventional support structures easily fail. Roadway collapse and large deformation occur just several months after tunnel excavation. In this study, a relatively new prereinforcement technique, the jet grouting (JG), is introduced to improve roadway stability. A ﬁeld test was performed for evaluating the practicability and applicability of JG in soft coal mass. A series of laboratory tests were conducted to assess the properties of coalcrete (coal-grout after JG treatment). A two-dimensional numerical model was established for validating the input parameters. Based on the veriﬁed model, three JG support cases for roadway were modeled and compared with a conventional support case, namely, the currently used support in this mine “rock bolts+U-shaped steel set +shotcrete.” The results show that the proposed prereinforced JG support structures can considerably control the deformation and failure zone of the roadway and improve the bearing capacity of coal mass. The mechanism of maintaining roadway stability using JG techniques is further revealed. Some suggestions are further presented to improve the stability of the jet-grouted roadway.


Introduction
e stability of the roadway is tightly related to coal production and the safety of miners. To maintain roadway stability, some typical methods including rock bolts system [1][2][3] and U-shaped steel set [4] have been successfully applied in shallow coal mass. However, with the increase of mining depth, the aforementioned support methods may fail to keep roadway stable due to high geostress and tectonic stress and soft coal mass [5]. In these cases, large deformation of roadway, e.g., floor heave, roof subsidence, and sidewall convergences, appears commonly [6]. Many engineering practices have proved that using a single support method is extremely difficult or sometimes impossible to control the large deformation especially in the deep underground roadways with weak surrounding rock mass [7].
Consequently, compound support techniques such as "boltscables-shotcrete-square confined arch" [8] and "bolts-cables-U-shaped steel set-shotcrete" were designed and applied in poor geological conditions. ese support schemes can reduce roadway deformation to some extent, while they cannot prevent the long-duration deformation of roadway effectively due to the rheological properties of soft rock mass under the high-stress state [9].
To deal with this problem, grouting, as an additional reinforcement method to main support systems, was utilized for considerably improving the mechanical properties (cohesion and friction angle), integrity, and stress state of the surrounding rock mass [10,11]. Generally, the delayed grouting on the roadway was widely employed to improve the strength and deformation-resisting ability of soft and fractured surrounding coal mass, when other primary support structures have been installed [12]. Based on this grouting method, pregrouting before roadway excavation was proposed and applied in many underground projects such as reinforcement of headgates and tailgates [13], enhancement of heading face [14], and water inflow prevention [15]. In these projects, the pregrouting technique increased the overall quality of coal mass and, therefore, improved its self-bearing capacity. However, it should be noted that the currently used pregrouting approach is limited by the groutability of materials [16]. So, it could only be used as a supplement for the main support system. Also, it is noted that the coupling support system, e.g., "bolt-cable-U-shaped steel set-shotcrete-grouting" is too complex to operate conveniently. Constructing each part of the system is timeconsuming and strenuous, which would decrease tunneling speed and increase costs. Furthermore, according to the fundamental idea of the New Austrian Tunneling Method (NATM), the aim of support system for roadway is to improve the self-stabilization ability or load-bearing capacity of the soft coal seam. us, it is necessary to propose a support method to simplify the compound support system and draw on the advantages of pregrouting.
Jet grouting (JG) technique, a prereinforcement approach, is normally used for soil ground improvement [17]. Currently, the subhorizontal jet grouting has been used in tunnel construction including tunnel canopy [18], collapsed tunnel recovery [19], and jet-grouted umbrellas in a tunnel [20]. is method reinforces the surrounding loose soil or fragmented rock ahead of tunnel face with jet-grouted columns, which effectively controls the convergence and extrusion of the tunnel. Based on different geological environments, the product of JG is generally called soilcrete or sandcrete, i.e., a mix of cement grout with soil or sand, possessing some special advantages including low deformability and high strength.
e mechanical properties of soilcrete or sandcrete played a crucial role in the stability of a jet-grouted tunnel, which has been studied by scholars [21][22][23][24]. Inspired by this, the coalcrete generated by JG technique for improving underground soft coal mass was reported recently [25]. e JG processes can be simply formulated as follows [26]. e high-pressure grout with a high-velocity jet is injected into soft coal mass through a small nozzle, and then, with the rotation of drill stem, the grout is mixed with coal particles to create coal-cement composite, i.e., coalcrete; a coalcrete column is then formed as the drill withdraws continually; then along the roadway profile, the next coalcrete column overlaps the former one with a certain length; finally, a JG coalcrete prereinforced support system is formed. en, the shotcrete and U-shaped steel sets could be used for further strengthening the support [27,28]. e control effect of JG coalcrete columns on roadway stability should be analyzed systematically. In addition, due to the influence of JG pressure, coal mass strength, and grout materials, the diameter of coalcrete column varies. erefore, the size effect on the mechanical properties of coalcrete should be considered in this study, though the traditional mechanical parameters of coalcrete such as strength and Young's modulus have been studied based on small samples in previous studies.
is study investigates the support effect of JG coalcrete on roadway stability and its mechanism. Based on the field observation of the JG test, the applicability of JG in soft coal mass was evaluated. Laboratory tests of coalcretes considering size effects were conducted to determine mechanical parameters. A 2D numerical model was established to verify the input parameters by comparing the modeling results with filed measurements. Furthermore, by the verified numerical method, three typical JG support cases and one conventional support case were compared and analyzed deeply in terms of roadway deformation, failure zone, and stress state. Based on the results, some suggestions were further listed.
is pioneering work proposes a workflow for designing JG parameters in practice and facilitates the application of JG presupport system in the roadway with soft coal mass.

Field JG Test.
A field trial of JG was conducted in Guobei coal mine. e thickness of the coal seam in this mine is about 10 m and the roadway was normally driven along the floor. Due to the low-strength coal mass, the collapse and spalling in the driving face normally appeared (Figures 1(a) and 1(b)). Moreover, with the conventional compound support system (U-shaped steel set + rock bolts + shotcrete), large roadway deformation still occurred just 1 month after tunneling work started.
Considering these cases, the jet grouting technique was first tried in soft coal mass to improve its self-stabilization ability and form a reinforced support system for the roadway. e schematic diagram of JG treatment for coal mass is shown in Figure 2. In the field, when the drill stem withdrew from the hole in coal mass, the high-pressure cement grout (20∼26 MPa) rotated with the rod and therefore cut coal mass continuously. e fragment coal particles mixed with cement grout and then a coalcrete column was formed. Specifically, the rotating speed and the reverse speed of the rod were 18 r/min and 20 cm/min, respectively. When the processes of jet grouting finished, the mechanical properties of loose raw coal were entirely changed. e strength and deformation-resisting ability of the coal were improved. e diameter of formed coalcrete in the field is normally around 400 mm∼600 mm. According to the engineering practice of controlling a roadway in soft coal mass, the largest deformation normally occurs at the surface of the roadway because the broken coal mass surrounding the roadway is basically in the extent of 400 mm based on the borehole observation. Also, practice experience has proven that controlling the fractured coal mass around roadway profile is critical for maintaining long-term roadway stability. As described above, the prereinforced JG coalcrete column (diameter from 400 mm to 600 mm) has higher strength and rigidity than raw coal mass, which means that the fractured coal in the corresponding range from the surface is replaced and controlled by coalcrete before roadway excavation.
is meets the requirements of a roadway support system by such coalcrete columns and 2 Advances in Materials Science and Engineering indicates that the jet grouting technique in soft coal mass is practicable and promising.

Laboratory Tests of the Mechanical Properties of
Coalcrete. In practice, the coalcrete columns can improve the shear strength of soft coal mass before excavation and resist the surrounding stress after tunneling by constructing a stable arch structure [25]. Hence, it is necessary to evaluate the strength behaviors of coalcrete. A series of shear tests and uniaxial compressive tests considering the size effect were performed to assess the mechanical properties of coalcrete. To select reasonable mechanical parameters for the following numerical simulation, coalcrete specimens with different sizes were set as 50 mm, 100 mm, 150 mm, and 200 mm cube. A direct shear test machine (YZW500, Jinan, China) was used for obtaining the shear parameters of coalcrete according to the designed shear tests scheme (given in Table 1). As for the uniaxial tests of the specimens, a displacement control with the loading rate 0.5 mm/min was applied by a compression testing machine (MTS C64, Minnesota, USA).

Coal Measure Rock Mass
Properties. e surrounding rock mass in the roadway is in complex geological conditions such as groundwater conditions, fractures, and discontinuities. In this paper, the geological strength index (GSI) was applied to obtain the parameters of rock mass [29]. Furthermore, the rock mass rating (RMR) classification system was applied to determine the GSI; as RMR is calculated from a series of ratings of various parameters including UCS of rock, RQD, and joint spacing, it is more quantitative and practicable [30]. e empirical correlation between GSI and RMR is given in equation (1) [31]. Based on the RMR values obtained from the field observations, the GSI values of various rock masses are summarised in Table 2. Moreover, the intact rock properties such as compressive strength (σ ci ), m i constant, were tested in the laboratory. en, the rock mass properties (strength and deformation modulus, shown in Table 2) were calculated based on the equations in the literature [32,33]: (1) e Hoek-Brown criterion is given as follows: In equation (2), m b , s, and a are the rock mass constants.
In order to calculate the rock mass deformation modulus, GSI-based empirical equation is given as where E i means the deformation modulus of intact rock.

Modeling Framework and
Procedure. e numerical method has been successfully applied in many cases to analyze the roadway stability nowadays. e primary objective of this work is to analyze the control effect of JG columns on roadway stability by comparing the conventional support schemes with the JG support system. Although the displacement and stress fields during tunneling should be analyzed using three-dimensional numerical methods, two-dimensional methods were used successfully to analyze the relevant problems such as surface settlements and tunnel convergence [34], the effect of construction procedures [35], etc. Besides, compared with the three-dimensional model, the 2D model was more efficient and less time-consuming, which has been widely used for analyzing the effect of JG application in tunnels as well [18,36]. In this study, the software RS2 was introduced to reveal the effect of JG coalcrete on a representative roadway cross section.
Same as the size in the field, the width and height of the roadway in the numerical model are 5.6 m (symbolized D) and 4.8 m, respectively (shown in Figure 3). To reduce the influence of mesh size and boundary conditions, a series of trial analyses were performed. Finally, the model domain is 5D × 5D and the appropriate mesh discretization is 3-noded triangular elements. e geometry of the model is shown in Figure 3(a) and the enlarged rectangular region represents different support systems, i.e., conventional support system (rock bolts + U-shaped shed + shotcrete) (Figure 3(b)) and JG prereinforced support system (coalcrete + shotcrete) ( Figure 3(c)). e established model was fixed in the horizontal direction at the right and the left side. e bottom parts of the boundary were pinned. ese models contained approximately 22,000 elements with the optimum mesh size after testing. Due to the effect of last working face retreat on the stresses at the position of the studied roadway, the principal stresses would deflect. erefore, based on the previous study of 3D global-2D local modeling method [37], the stresses for the 2D model were transformed from the three-dimensional stresses. e applied vertical stress and horizontal stress are 19.23 MPa and 11.21 MPa, respectively.
Selecting a proper constitutive model for materials is essential. In this study, the Mohr-Coulomb (MC) model was utilized for modeling the rock mass and the coalcretes, which has been commonly used in coal measures rocks and jet-grouted materials (soilcrete, sandcrete, etc.) reported in the literature [38,39]. Furthermore, the mechanical behavior of coalcrete is studied in Section 4.1, in which the results show that coalcretes' failure suited the Mohr-Coulomb criterion well. A linear elastic constitutive model is used for modeling shotcrete and U-shaped shed by structural elements according to recommendations [40,41]. Several key steps were performed for establishing various models as follows: (1) e boundaries were set based on Figure 3(a), and stresses were initialized according to calculation. (2) e mechanical parameters were allocated to rock mass ( Table 2) and coalcretes ( Table 5 in Section 4.1).
It should be noted that the coalcrete columns surrounding the roadway profile were installed subhorizontally before tunnel face excavation.
(3) e arched roadway section was excavated and supported simultaneously. In this stage, models were divided into two types (see Figures 3(b) and 3(c)). e conventional support system, i.e., the rock bolts, U-shaped shed, and shotcrete, was installed at this stage. e support materials are consistent with the field and their properties are summarised in Table 3. By comparison, the proposed JG support system just needs to spray shotcrete since the prereinforced coalcrete columns have been already generated in stage 2.
(4) Different model simulations were run for various purposes. One model is used for corroborating input parameters (conventional support, verification in Section 3.2). Other models are employed for investigating the effect of the JG support system on roadway stability (comparison in Section 4.3).

Verification of the Established Model.
e 2D numerical analysis should consider the three-dimensional face effect after tunneling, as the stress redistribution and deformation occur during tunnel excavation. ere are some techniques used in the 2D plane strain model, such as stress relief method [42], stiffness reduction method [43], and disk calculation method. Among these approaches, the stress relief method, also called convergence-confinement method, was commonly utilized for modeling around tunnel face [34]. Referring to this method, a displacement release method was presented [44]. In this study, the displacement release ratio was calculated based on the previous recommendation [45].    Advances in Materials Science and Engineering e displacement release ratio for models increased from 0% to 100%, and each model has run for enough steps to reach the balance. When the ratio was set as 100%, the model could be regarded as that the tunnel face was far away from the studied cross section [46]. e displacement release ratio is related to real time, which means the stress acting on roadway can be expressed as the function of time, as the in situ displacement and its corresponding time were recorded simultaneously. By this method, the real stresses on roadway increase over time, and therefore the roadway deformation exhibits rheology property, which has been confirmed in deep tunnel in the literature [47]. In practice, convergence-measuring pins were adopted to measure the sidewall convergence with time (shown in Figure 4). en, the displacement/stress release rate was calculated, and eight representative results were further selected (as shown in Table 4).
Correspondingly, a 2D numerical model was constructed according to the modeling procedures described above. In the model, there are eight stages of applied stresses, which can simulate the roadway deformation over time. Figure 5 shows the evolution of sidewall convergence with time. Furthermore, the monitored deformation by the model was compared with the monitoring results in the field ( Figure 4). As can be seen, the trend and magnitude of these values are very similar, confirming that the model is reasonable and accurate.

Mechanical Test Results of Coalcrete.
e shear test results of coalcretes of various sizes under different normal stresses are shown in Figure 6. It can be seen that the shear stress increased linearly with the increase of normal stresses.
us, the Mohr-Coulomb failure criterion is used for fitting the data obtained from the tests, as given in the following equation:

Advances in Materials Science and Engineering
where τ, φ, c, and σ are the peak shear stress, friction angle, cohesion, and normal stress component, respectively.
It is clear that all coefficients exceeded 0.97, which means that there was a strong correlation between shear stress and normal stress of coalcrete specimens under different sizes, and the failure of coalcrete obeyed the Mohr-Coulomb criterion well. e calculated cohesion and friction angle are presented in Table 5. Moreover, deformation modulus and uniaxial compressive strength (UCS) based on test results are also given in Table 5. It can be due to the fact that the mechanical parameters of coalcrete such as cohesion, friction angle, deformation modulus, and UCS exhibit size effect, which are consistent with the experience and literature. In accordance with the recommendation in the studies [48], regression analysis with the fitting functions (shown in Table 5) was conducted. More specifically, the values of fitting coefficients R 2 were larger than 0.9, demonstrating that the function had a good agreement with experimental results.
In practice, the diameter of coalcrete columns depended on jet-grouting pressure, coal mass conditions, rod withdraw speed, etc. In general, it varied from 400 mm to 600 mm in soft coal mass. erefore, in this study, the typical diameter of the coalcrete columns, 300, 500, and 700 mm, was selected for further analysis. According to the regression equations and Table 5, it is reasonable to determine the   Advances in Materials Science and Engineering mechanical parameters of larger coalcrete. In this study, the tensile strength of coalcrete was regarded as 0.3 times that of UCS [49], which is different from the one-tenth of the compressive strength of rock materials [50]. e value of Poisson's ratio was set as 0.24 based on the test results. Normally, the values of Poisson's ratio with different sizes are similar regardless of the size.

Failure Mechanism of the Soft Coal
Roadway. e mechanism of deformation and failure of roadway in conventional support scheme was analyzed by displacement monitoring results and the evolution of the failure zone of surrounding coal measure rock mass. e results are depicted in Figure 7. As can be seen in Figure 7(a), after roadway excavation, deformation occurred in a short time with a relatively small extent of the damage zone, in which the shear failure dominated, whereas the tensile failure appeared at the vault and bottom corner of roadway, similarly with the studies presented in [51]. In this stage, the maximum extent of yielding zone just exceeded the length of rock bolts for the top two in support, which means rock bolts could be effective to some extent. As time increased, the shear failure zone expanded to the space in soft coal mass, especially along the left and right spandrel of the roadway. Correspondingly, the roadway displacement of the roof,   Advances in Materials Science and Engineering floor, and the sidewall increased significantly (Figure 7(b)). On the 30th day, the convergence of the roof and floor and the two sides reached 360 mm and 450 mm, respectively, accounting for more than half of the whole deformation in monitoring, which can be found in some studies having similar on-site monitoring tendency in deep coal mines [52]. e failure band developed sharply and propagated upwards and downwards, until it reached the top mudstone stratum and bottom sandstone horizon and stopped since the strength of these two strata was higher. Two shear zones along the spandrel of roadway formed in soft coal mass, and a fully damaged area was found in the immediate floor, i.e., mudstone lithology (Figure 7(c)). After 150 days, the failure zone further enlarged and its extent exceeded the control limit of all rock bolts, indicating that the support system was unstable and ineffective. Large deformation of the roadway such as roof subsidence, floor heave, and considerable sidewall convergence was observed at this stage (Figure 7(d)). Predictably, the convergence of roadway could develop further when subjected to stresses induced by longwall panel retreat, as the surrounding rock mass totally failed in a large range. Under such harsh geoengineering conditions, geohazards are very likely to occur. e tunnel needs to be maintained frequently.
Based on the above analysis, the roadway in such conditions cannot maintain stability, and the possible reasons are as follows. (1) e surrounding coal mass with poor quality cannot bear load effectively [53]. e dramatic stress redistribution and stress concentration due to roadway excavation resulted in a large failure zone in soft coal seam around the roadway section. is causes the deterioration of the bearing capacity of coal mass. en, the loading transferred to floor stratum (a relatively low-strength mudstone lithology) and broke the floor completely. (2) e currently used support method (i.e., rock bolts, shotcrete, and U-shaped steel set) cannot form a coordinated and stable control system in the roadway. e passive support component, U-shaped shed, cannot provide active support effect after excavation at the primary stage, and it usually bucked and failed because of uneven contact with surrounding coal mass in practice [54]. In terms of the active support component, boltshotcrete, their bearing capacity and confinement were too low to increase the residual strength of rock mass significantly, resulting in a large yielding extent. (3) e worse situation is that the failure zone extended rapidly with stress release, whereas the bearing capacity of surrounding rock mass became insufficient and the support structure became unstable, which further deteriorated the stability of roadway.
As discussed above, the coal mass properties play a vital role in the stability of the roadway. Improving the quality of soft coal mass indicates not only the reinforcement of coal but also the enhancement of integrity of support structure. Hence, the jet grouting support system is designed and compared with the conventional support schemes.

Comparative Studies of Various Support Cases.
Case studies of the jet grouting effect on roadway stability were carried out in three conditions (diameter of coalcrete column of 300 mm, 500 mm, and 700 mm), compared with the conventional support case mentioned above, i.e., without coalcrete. e detailed support materials in different cases were summarised in Table 6. e roadway section, rock mass properties, and the stresses acting on the model were identical in these numerical cases. It is noted that in practice the final status of the supported roadway is always a matter of concern to engineers. Accordingly, in this study, stresses acting on these models were the totally released stresses, i.e., the maximum stresses based on the calculation in Section 3.1. e strength parameters of coalcrete of different sizes were input to models according to the calculated results described above. In the following, the effect of various support systems on roadway stability was investigated, including deformation, failure zone, and stress distribution, and then recommendations of JG schemes were given based on the comparison results.

Deformation of the Roadway for Various Cases.
e recorded deformation results of the roadway are comparatively shown in Figure 8. It is clearly that case 1 had a relatively even and large convergence, i.e., the similar value of roof subsidence, floor heave, and sidewall deformation. By comparison, cases 2, 3, and 4 had a relatively larger displacement at the lower wall sides than any other position of the roadway. Also, case 1 had the largest roof subsidence and sidewall convergence, followed by case 2, case 3, and case 4. Moreover, case 2, case 3, and case 4 had similar values on floor heave, while case 1 had larger value on floor heave. e similar trend for all models was that the sidewall convergence was larger than the roof-to-floor convergence. It is clear that the roof-to-floor convergence and sidewall convergence reduced with the increase of diameter of coalcrete. More concretely, compared with case 1, the roof-to-floor convergence declined by 71%, 72%, and 75% for case 2, case 3, and case 4, respectively. In terms of sidewalls convergence, it decreased by 46%, 55%, and 63% for case 2, case 3, and case 4, compared with case 1. Figure 9 illustrates the failure zone distribution of roadway under four support conditions by an intuitively comparative way. Obviously, all cases had the irregular shape of failure zone and their extents were different, but all extended to three lithologies, i.e., coal mass layer, mudstone in the immediate floor, and sandstone in the main floor. Case 1 had the largest failure zone, followed by case 2, case 3, and case 4, which means that with the increase of the diameter of coalcrete, the failure zone of roadway reduced significantly. In order to evaluate the extent of failure zone quantitatively, some lines such as L w , L e , L f , and L r were defined to represent the distance of failure zone ( Figure 9). As can be seen, the characteristic of "L e > L w > L f > L r " indicates that the failure on roadway shoulder is larger than other areas. In comparison with case 1, the JG support cases (cases 2, 3, and 4) significantly reduced the failure zone on roadway shoulders, roof, floor, and sidewalls. More specifically, the length of the failure zone on roadway shoulder (L e ) decreased by 47%, 56%, and 61% for cases 2, 3, and 4, respectively. e failure distance at the roof decreased by nearly 100% for jet grouting support schemes.

Comparison of the Extent of Failure Zone.
It is widely accepted that the elements in shear or tension failure zone can represent the development of fractures or cracks in rock materials [55]. Figure 10 depicts the failure zone distribution around roadway under various cases. More tensile cracks in case 1 occurred at the roof in coal mass, floor, and corner in mudstone (see Figure 10(a)), whereas cases 2, 3, and 4 had tension cracks mainly concentrating on bottom sidewalls of coalcrete column and the corner of roadway (Figures 10(b)-10(d)). In general, combined with the number of shear and tensile cracks, the region with more failure led to coal mass and coalcrete more dilation and therefore caused deformation and roadway convergence (see the failure distribution and the deformed profile by grey lines in Figure 10). For example, case 1 had large deformation, as it had a lot of failure cracks around the  Advances in Materials Science and Engineering 9 roadway section. After the application of coalcrete columns (cases 2, 3, and 4), there was a significant reduction of shear and tension cracks on the vault and subsequently less deformation on the roof. Moreover, with the increase of the diameter of coalcrete column, the bottom and sidewall convergence declined considerably since the total number of failure cracks decreased. ere was similar floor heave for cases 2, 3, and 4 because the total number of cracks was close. After roadway excavation, though some failure cracks appeared on coalcrete because the redistributed stress exceeded the ultimate strength of coalcrete, the overall failure zone just extended into coal mass in a smaller extent. is is because though the failure occurred at the bottom sidewall of coalcrete columns (see Figures 10(b)-10(d)), the relatively higher residual strength of coalcrete still had a better self-supporting capacity, which indicates that better-confined stress still acted on surrounding coal mass and thus the self-stabilization ability of surrounding coal mass could be improved notably. Consequently, a smaller damage zone would extend into coal mass than conventional support (comparing case 1 with cases 2, 3, and 4 in Figure 9).

Stress Distribution Analysis.
It should be pointed out that the evolution of stress is closely related to the propagation of failure of the rock mass [31,56]. e relationship between each other can be explained as follows. For the excavation of underground constructions such as tunnels, roadway, and openings in rock mass, the failure normally first occurred in rock at shallow surface as the deviator stress reached its ultimate strength rapidly; the failing rock entered the postfailure stage but still had residual strength, causing a decline in the bearing capacity (self-supporting capacity) of shallow rock; without support, the stress further transferred and propagated into the deeper surrounding rock mass that had relatively larger confined stress meaning a higher ultimate strength of rock mass, until the concentrated stress was lower or equal to the ultimate strength of rock mass. e above processes of failure in rock mass lead to a large stress relaxation zone.
To show the effect of coalcrete on stress distribution at sidewall comparatively, observation lines were used for obtaining the maximum principal stress in the sidewall. All monitored results are shown in Figure 11. According to these, when the stress around the roadway is lower than the in situ stress (around 19 MPa), the areas are defined as stress release zones. us, there are stress release zones around the roadway, but the extent of them was different. e JG coalcrete columns can optimize the stress state to make major parts of support in better stress conditions, while some areas were still in poor force state (e.g., the stress concentration zone at the corner). When additional stress such as induced stress acted on the roof and sidewalls, stress may further transfer from the sides to the floor. e superposition effect of these stresses in roadway corners was negative for stabilizing floor and the whole JG structure, and, therefore, more attention should be paid to these areas in roadway.
It is clear that the stresses near the sidewall surface in supported cases (i.e., cases 1, 2, 3, and 4) were larger than the no-support case (approximately 0 MPa), which means that the bearing capacity of coal mass was increased by support systems. However, by further analysis, when the distance of monitoring points to the surface was larger than 3 m, the stress distribution of case 1 was similar to no-support case, which signifies that the conventional support cannot significantly improve the load-bearing capacity of deep coal mass. In addition, compared with case 1, the stress in sidewall surface in cases 2, 3, and 4 increased by 50%, 60%, and 68%, respectively, indicating that the bearing capacity of coalcrete was much larger than coal and it also increased with the diameter of coalcrete. ough the coalcrete columns have failed (shown in Figures 10(b)-10(d)), the stress in coalcrete (see the region marked by a blue dotted circle in Figure 11) still increased greatly with the distance from the roadway surface, which means the failed coalcrete still bore the load effectively. Interestingly, a stress dissipation phenomenon in cases 2, 3, and 4 appeared beyond the coalcrete column, and all values of dropped stress in coal mass were larger than the stress at the same position in case 1 or nosupport case. at means the residual strength of coal mass was enhanced after coalcrete application.

Selection Principles and Suggestions to the JG Support
System. Based on the results and analysis above, it is clear that the jet grouting support system ("coalcrete + shotcrete") is promising for stabilizing roadway. In comparison with the conventional support system, i.e., "bolt + shotcrete + Ushaped steel set," JG support system has many benefits on roadway stability such as less roadway deformation, the smaller extent of failure zone, and optimal stress states. Moreover, with the increase of the diameter of jet-grouted coalcrete column, the control effect gets better. However, the diameter of approximately 400 mm to 600 mm is more economic and practicable in the field, as bigger coalcrete column needs higher jet pressure and more cement grout, which is money consuming and uneconomic [57,58]. e selection of the coalcrete diameter depends on the  Advances in Materials Science and Engineering 11 requirements of control effect and geological conditions of roadway in practice. Also, it should be pointed out that there are still some defects of JG support system, e.g., some failure at bottom sidewalls and stress concentration at the corner of roadway. Subsequently, to keep the stability of JG roadway in the long term, some suggestions are additionally proposed based on previous successful experience in coal mines [59]: (1) e rock bolts or cables are recommended to install along roadway shoulder and sidewalls. By the jet grouting support method, though failure zone of the rock mass is in a controllable extent, some commonly used support components such as cable bolts and rock bolts would be used for further reinforcement. e bolts can effectively bond the yielding coal to intact coal mass to avoid the sudden collapse of the roadway surface. (2) e inverted arch is suggested for optimizing the stress conditions in the floor and bottom corner of the roadway. Due to the angle effect of the JG support system on the roadway, the high-stress concentration zone is in the corner, which is negative for the long-term stability of the whole support structure. e inverted arch can not only improve the stress state but also provide a certain bearing capacity, which could effectively transfer the stress into deep rock and reduce floor heave.

Concluding Remarks
is paper presented a relatively new prereinforcement support system of the roadway with soft coal mass, i.e., jet grouting (JG) method. Based on the field observation of JG in coal mass, this method met the requirements in techniques and could be used for stabilizing roadway in deep mines. According to the data collected from the field and laboratory, the GSI was used for estimating coal measure rock mass properties. A numerical model was further established and then verified, which confirms the reliability of input parameters. e mechanical properties of coalcrete produced from JG in coal mass were obtained by a series of tests and their size effect was considered. Furthermore, according to the verified model and the determined mechanical parameters of coalcrete, three JG numerical models were established and compared with the conventional support case (the currently used support approach, "rock bolts + U-shaped steel set +shotcrete").
e results show that the JG method could promote roadway stability by reducing the deformation, failure zone and optimizing the stress conditions.
In addition, the mechanism of this new support on stabilizing roadway was revealed. e generated coalcrete columns have higher strength, and larger deformationresisting ability can form a stable arch structure. Some advice such as installing rock/cable bolts and inverted arch on floor were further proposed for the stability of roadway by JG treatment. is study mainly focuses on the comparison of conventional support and JG support on roadway stability by numerical modeling methods, and in the future, the more detailed application information of this innovative technique for roadway would be studied and analyzed.

Data Availability
e Microsoft Excel Worksheet data used to support the findings of this study are available from the corresponding author (yuantiansun@cumt.edu.cn) upon request.