In this study, the roof-cutting-type gob-side entry retaining is introduced, and its application in medium-thickness coal seams is studied. Based on the analysis of the construction procedure and principle, the mechanical model of the retained roadway structure and cantilever beam formed by roof cutting was established, and the support resistance and roof deformation were obtained. In addition, through technological design analysis and numerical simulation, the parameters of roof cutting were determined. The roof-cutting height and angle were designed to be 9 m and 15°, respectively. Flac3 D was used to analyze the stress evolution law under different mining conditions. The stress on the integrated coal side and roof subsidence was lower when the roof-cutting height was 8∼10 m and the cutting angle was 15°. Through field monitoring, the roof pressure, gob-side lateral gangue retaining pressure, anchor cable stress, and deformation of the surrounding rock eventually reached a stable state. This indicates that the roof cutting can effectively cut down the overlying strata over the gob and form a stable entry structure to meet the requirements of the next working face.

1. Introduction

The waste of coal energy has attracted attention due to the general increase in the degree of mining mechanization [13]. To ensure the safety of underground work, the traditional longwall mining method balances the roof pressure in the remaining coal pillars, but the waste of coal resources is more serious [4, 5]. Therefore, the traditional gob-side entry retaining (GSER) technology has been developed to solve the problem of wasting resources, that is, use certain technical means to retain the crossheading, and serve the next working face [6], as shown in Figure 1(b). This technology realizes nonpillar mining, increases the recovery rate, reduces the development ratio, and provides advantages in ventilation and gas accumulation control [79]. Therefore, it is widely used.

Although the traditional GSER solves the resource waste problem [10, 11], it also has limitations: it has complex construction technology, it has high costs [9], most of the backfills are rigid materials and thus do not largely deform [12], and overwhelming under concentrated stress occurs, eventually resulting in instability of the roadway [13, 14]. At present, there are many studies on the stability of retained roadway and gob-side support body. The study of stability mainly focuses on the deformation characteristics of the retained roadway [15] and the stress distribution characteristics under the influence of mining. To understand its deformation and failure mechanism, supporting measures in roadways were performed [16]. Research on GSER mainly focuses on the material and its mechanical properties or builds the corresponding mechanical model and proposes a method to determine the main parameters of the supporting body [1719]. Some researchers have conducted simulation experiments to study the influence of the fracture of overlying strata on the stability of the retained roadway under different geological conditions, and simulation analysis of the related roadside filling technology [20, 21]. The mentioned studies have made important contributions to the development of GSER technology, but the essential problem remains unresolved. It mainly includes two aspects: one is the deformation of the surrounding rock caused by the rotation and subsidence of the strata above the retained roadway; the other is the fracture of the overlying hard roof into a structure of transferring stress through mutual extrusion, as shown in Figure 2(b). Therefore, the surrounding rock of the retained roadway is still under the stress environment of the traditional coal-mining mode.

To fundamentally improve the stress distribution of surrounding rock, an in-depth study was conducted based on the field application of roof-cutting-type GSER. Using theoretical analysis and numerical simulation, we analyzed the surrounding rock movement law and performed field monitoring verification. Our study contributes to a better understanding of the development of GSER technology in a medium-thick coal seam.

2. Gob-Side Entry Retaining by Roof Cutting

2.1. Construction Procedure

Figure 2 shows the main construction steps of the roof-cutting-type GSER [22, 23]. The first step: the constant resistance and large deformation (CRLD) anchor cables are adopted to strengthen the original support in the prereserved roadway area. The second step: the bidirectional energy-cavity blasting (BECB) technology is carried out along the roof beside the gob, shown in Figures 2(a) and 2(b). The third step: the gangue retaining support is carried out along the cutting-line with the advance of the working face, so as to prevent the collapsed rock from entering the reserved roadway area. Generally, the measures of retaining gangue are “single hydraulic pillar + I-shaped steel (or U-shaped steel + reinforcing mesh).” The last step: withdraw the temporary support when the support section is far away from the working face and the surrounding rock movement is stable. At this time, the retained roadway is completed, it can serve the next working face mining, and thus it realizes the mining mode of a single roadway without coal pillars.

2.2. Principle

Figure 1 shows the surrounding rock structure of the retained roadway under different retaining process conditions. As shown in Figure 1(c), the roof cutting interrupts the stress transfer path between the roof strata of the retained roadway and overlying strata above the goaf, shown in Figures 1(a) and 1(b), caused by mutual compression, and Kγh is the peak lateral support pressure. Moreover, the collapsed strata can use its bulk increase character to support the upper strata, reduce the rotary subsidence of the upper strata, and optimize the stress environment of the retained roadway [24]. Compared with traditional mining technology, this technology changes the stress environment in the area near the retained roadway and also changes the law of strata movement in the vicinity, reduces the deformation of surrounding rock, and guarantees the safe use of the retained roadway.

2.3. Key Technologies

According to the analysis of the construction procedure and principle, the roof-cutting-type GSER has two key technologies: the CRLD anchor cable reinforced support and the roof cutting by directional blasting.

2.3.1. Constant Resistance and Large Deformation Cable Support

As shown in Figure 3, the retained roadway generally experiences the influence of advanced stress, mining dynamic pressure, and reuse stage, and the surrounding rock deformation and stress of the anchor bolt continue to increase. At present, the elongation of the common anchor bolt is low, and it cannot adapt to the large deformation of the surrounding rock. The anchor bolt support strength considerably decreases, its impact resistance is significantly reduced, and it is not conducive to underground safe construction. Therefore, in the process of retaining the roadway, the CRLD anchor cable is used to reinforce the original support. Figure 1(b) shows the deformation curve of the surrounding rock and the stress of the ordinary anchor bolt in the process of the traditional retaining roadway. The deformation exceeds 500 mm, and the anchor bolt breaks after the stress reaches a maximum value of 180 kN.

Figures 4(a) and 4(b) show the main structure of the anchor cable and the constitutive model and working principle, respectively [25, 26]. There are two main stages of deformation: the elastic deformation and in sequence the constant resistance tensile deformation. When the surrounding rock deformation is small, the axial force acting on the bolt is less than the constant resistance, and thus, the rod body is in the elastic state and produces elastic deformation. The deformation of the surrounding rock is resisted by the elastic deformation of the rod body. When the deformation is large and the axial force is larger than the constant resistance, the cone slides through the sleeve and generates constant resistance; i.e., the CRLD cable can produce constant support resistance under the condition of large deformation. Therefore, it is used for retained roadway support. Figure 4(b) shows the stress value of the cable at different deformation stages, which can be calculated as follows [25]:

The deformation energy absorbed at different deformation stages is

2.3.2. Directional Blasting Technology

After strengthening the roadway support with the CRLD anchor cable, the next step is the roof cutting, which is also called bilateral cumulative tensile blasting technology [22]. Its main steps are to regularly drill holes along the roof and charge for blasting, cutting the overlying roof strata. However, when the conventional blasting method is employed (Figure 5(a)), the surrounding rock is damaged in large areas, and there is a serious risk of accidents, such as roof caving. Therefore, before charging, a guide pipe with two rows of guide holes should be installed in the borehole, such that the stress and energy generated by blasting can be transmitted from the guide holes to the surrounding rock surface, making the overlying strata crack along the set direction, and minimizing the damage of the rock mass, as shown in Figure 5(b).

From the above blasting analysis, it can be seen that reasonable control of the spacing of blasting holes is significantly important to maintain the stability of the surrounding rock. By studying cumulative blasting, [22] the peak stress generated can be deduced as follows:where is the radius of the cartridge, is the hole diameter of the blast hole, is the initial detonation pressure, is the hole depth, is the air interval length in the blast hole, is the plugging length of the blast hole, is the explosive density, is the rock mass density, and is the cumulative energy coefficient, .

According to the attenuation law of the stress wave in the rock mass, the calculation formula for the blasting stress damage range is as follows [24]:where is the lateral pressure coefficient, is the initial damage of the rock mass, is the tensile strength of the rock mass, is the attenuation index of the explosion stress wave in the rock mass, and .

The above analysis is generally used to design the spacing of blasting holes. In addition to controlling the proper spacing, the height and angle are also important design parameters. The inclination angle is determined by the site condition and numerical simulation and is commonly less than 30°. The determination of the roof-cutting height is related to the lithology and mining height of the working face. The calculation principle is shown in Figure 6. Assuming that there are m rock layers within the roof cutting, the rock layers are marked as “1, 2, 3,…,m” from bottom to top. Theoretically, the caving rocks should fill the goaf below; that is, the filling height of the collapsed rock should be equal to the sum of the thicknesses of the coal seam and the rock strata within the range of roof cutting.

The following can be obtained:

The roof-cutting height iswhere is the roof-cutting height; is the thickness of the coal seam; are the coefficients of broken expand of the first, second, …, m-th layers, respectively; are the rock thicknesses of the first, second, …, m-th layers, respectively; and is the average coefficient of broken expand.

2.3.3. Gangue Retaining Support

After mining, to prevent the caved rock from rushing into the retained roadway, gangue retaining support is needed. The main form of gangue retaining is usually “support bracket + retractable gangue retaining 29# U-steel + mesh reinforcement.” The length of 29# U-steel is 2 m and 2.5 m, and the spacing is 500 mm. It is connected by two pairs of flanges. The upper and lower ends of the flange are 50 mm away from the lap end of the U-steel, respectively. It is not less than 300 mm below the floor. Besides, within the influence range of dynamic pressure, steel mesh is used to fix the duct cloth for preliminary air leakage prevention, and spray treatment is used after stabilization.

2.4. Mechanical Model

Compared with the traditional filling type GSER, the surrounding rock structure fundamentally changes due to the retaining roadway process; therefore, it is necessary to study the deformation mechanism and then combine it with the field application to develop the corresponding supporting measures. As shown in Figure 1(c), owing to the roof cutting, a wedge-shaped confined zone is formed above the retained roadway. It can be regarded as a cantilever structure; therefore, the mechanical model of the retained roadway was established, as shown in Figure 7. The self-weight of the strata above the retained tunnel can be regarded as uniform load, and the support in the roadway can be regarded as uniform stress. The support of the caved rock mass in the goaf to the upper stratum can be regarded as uniform stress, and the concentrated force on the right side of the retained roadway is the support of the single pillar. It can be seen that the main parameters of the model are the support resistance of a single hydraulic prop and the deformation of the cantilever beam; thus, the solution of the model was analyzed.

The support resistance at the roof-cutting side of the reserved roadway is as follows [23]:where is the support resistance of the single pillar, is the ultimate bending moment of the upper strata, is the shear force, is the support force of the gangue in the goaf, q is the self-weight, is the self-weight of the direct roof, is the width of the limit equilibrium area, is the roadway width, is the length of block BC, is the sinking value of end C when block AC is cut off, is the subsidence of the B end, and is the thickness of the old roof rock block.

The displacement of the cantilever beam can be obtained from Figure 7 as follows:

Based on the above analysis of support resistance and surrounding rock deformation, the required support resistance and deformation tendency of the retained roadway can be preliminarily obtained and provide a certain reference for the support and gangue retaining support in the roadway.

3. Field Application

3.1. Engineering Background

The Dianping coal mine is located in Lvliang City, Shanxi province, as shown in Figure 8. The depth of 200 working faces is 360 m, the length of open-off cut is 220 m, the length of the crossheadings is 1088 m, the length of the open-off cut to the stop line is 942 m, the average coal seam dip angle is 4°, and the average thickness of the coal seam is 3.1 m. The designed retained roadway was used as the return air roadway. The layout of the working face is shown in Figure 8(b). Figure 8(c) shows the lithology of the roadway roof and floor. The direct floor of the working face is deep grey sandy mudstone with an average thickness of 3.3 m; the old floor is medium-grained sandstone with an average thickness of 2.8 m. From the bottom to the top of the coal seam roof, there are 2.1 m thick sandy mudstone, 6.1 m thick medium sandstone, 0.7 m thick sandy mudstone, 2.4 m thick coal interlayer, and 3.0 m thick limestone.

3.1.1. Support Design

Figure 9(a) shows the support section. The roof adopted a 20 × 2000 mm steel bolt, metal mesh, and steel bar ladder for combined support. The row and line space between the bolts was 940 × 1000 mm, with six bolts in each row. The two sides of the roadway were supported by a 16 × 1600 mm steel bolt and plastic mesh, and the row and line space between the bolts was 1200 × 1000 mm. An anchor cable with a size of Ф 18.9 × 6200 mm and channel steel was also used to strengthen the roof support, two rows of anchor cables were arranged along the roadway, and the row and line space were 2350 × 1000 mm. Based on the traditional support scheme, two rows of CRLD anchor cables were arranged: the first 500 mm away from the gob, with row spacing of 1000 mm; and the second in the middle of the roadway, with row spacing of 2000 mm. In the first row, cables were connected by a steel belt.

The self-circulation portal hydraulic support was also used for support in the retained roadway. Furthermore, the roof broken area increased the single pillar to strengthen the support. The length of the portal hydraulic support was 2.8 m, the working resistance was 2040 kN, the designed row distance was 2000 mm, and the working resistance of the single prop was 280 kN.

Figure 9(b) shows the charge structure of the blasting hole, and the depth of the blasting hole is 8 m. Based on the cumulative blasting analysis and in situ tests, the distance between blasting holes was 500 mm. The binding energy tube is made from a special PVC material to control blasting energy. The device surface has two distributed rows of energy-accumulated grooves at 180° from each other with a spacing of 8 mm. The length of the binding energy tube is 1.5 m, the optimal sealing clay length was 2 m, the cartridge size was Ф27 × 300 mm, and the quality was 200 g/volume.

25# U-shaped steel shrinkable support is used to support the gangue after the support. The retractable 25# U-shaped steel is overlapped by upper and lower sections, the overlapping length is not less than 1 m, the steel shed length is 2.0 m and 2.5 m (it can be adjusted and matched according to the roadway height), two pairs of Kalan are used for connection, and the distance between the upper and lower edges of Kalan and the overlapping end of U-shaped steel is 50 mm respectively. 25# U-shaped steel shed shall be buried no less than 300 mm below the bottom plate.

4. Numerical Simulation

4.1. Model Establishment

To further study the evolution law of stress and displacement of roadway surrounding rock after roof cutting, the FLAC3 D numerical software was used to simulate three types of GSER, that is, the remaining coal pillar, wall building, and roof-cutting type based on the original size of the site [2730]. As shown in Figure 10(a), a numerical calculation model was established. The model dimensions were as follows: x × y × z = 300 × 150 × 80 m, the excavation size of the simulated roadway was x × y × z = 5 × 150 × 3 m, and the working face excavation size was x × y × z = 220 × 120 × 3 m. The numerical models under the three processes were basically the same, except for the marked dotted line in the figure, which was different at the reserved roadway. The first type is coal pillar mining, the second type is traditional GSER, and the third type is roof-cutting-type GSER. Figure 10(b) shows the model boundary conditions. The left and right boundaries of the model limit the displacement in the x-direction. The anterior and posterior boundaries limit the displacement in the y-direction and set the horizontal stress gradient, the lower boundary limits the displacement in the z-direction, and the upper boundary imposes uniform self-weight stress. The model failure complies with the Mohr–Coulomb strength criterion.

Table 1 shows the physical and mechanical parameters of each rock layer in the numerical model. As the roof and floor rock layers are mainly dominated by mudstone and sandy mudstone, the relevant indoor rock mechanics tests have also been performed in the second section. Thus, the experimental data can be directly applied. The other mechanical parameters such as medium grain sandstone and limestone are provided by the mine side.

5. Stress Distribution Law under Different Mining Technologies

5.1. Surrounding Rock Stress

Figures 11(a)11(c) show the stress distribution of the roadway surrounding rock under the conditions of three retaining roadway technologies, and Figure 11(d) shows the solid coal side stress curve. The distance between the stress concentration area and the coal wall is relatively close under the three types of roadway retention processes, and the side abutment pressures are 15.9, 14.1, and 12.1 MPa under the conditions of retaining coal pillars, retaining roadways with the traditional filling body, and roof cutting, respectively. In addition, compared with the traditional stress nephogram, the low-stress area of the roof is larger, which indirectly indicates that the roof cutting improves the stress distribution of the surrounding rock.

5.2. Stope Stress

Figure 12(b) shows the stress nephogram and curve at the position of the 5 m monitoring line in front of the working face. The advanced stress distribution rule and stress value of the working face are close under the conditions of the remaining coal pillar and traditional roadway reservation. However, under the condition of roof cutting, the stress clearly decreases, and its transverse influence range is approximately 60 m; that is, the roof-cutting-type GSER improves the stress distribution and also reduces the advanced stress concentration to a certain extent, which is conducive to the safe operation of the working face.

5.3. Research on Roof Cutting
5.3.1. Different Roof-Cutting Heights

From the above analysis, the stress distribution of the surrounding rock under the condition of roof cutting has been optimized to a certain extent. To obtain the best optimization effect, different roof-cutting heights and angles were simulated and analyzed to obtain the best roof-cutting parameter. Figure 13 shows the distribution law of the retained roadway surrounding rock stress under different roof-cutting heights. When the roof-cutting height is 6, 8, and 10 m, the peak values of stress concentration are 13.0, 12.1, and 11.4 MPa, respectively. From 6 to 10 m, the stress decreases by approximately 11%, and the stress concentration area transfers to the deep position of the entity coal. From Figure 13(d), it can be seen that the larger the roof-cutting height, the smaller the peak stress, and the farther the stress concentration area is from the coal wall, the more beneficial it is to maintain the roadway. However, when the cutting height reaches a certain degree, it has little influence on stress reduction; therefore, the height of the roof cutting should be higher than or equal to 8 m but lower than 10 m.

5.3.2. Different Roof-Cutting Angles

Figure 14 shows the distribution law of the retained roadway surrounding the rock stress under different roof-cutting angles. It can be seen that when the roof-cutting angles are 0°, 15°, and 25°, the peak values of stress concentration are 12.1, 10.9, and 11.7 MPa, respectively. From Figure 13(d), when the angle is 0° and 25°, the stress value is relatively large but similar for 15°, and the peak value is the least. Thus, 15° provides the best results.

The above analysis indicates that the roof-cutting-type GBER can effectively interrupt the stress transfer path between the roadway and the gob roof strata, considerably reducing the peak value of stress concentration, distancing the concentration area from the coal wall, and transferring to the deep position of the entity coal. Owing to the roof cutting, the continuity of the roof over the roadway and the gob is interrupted; therefore, the roof deformation of the roadway is less affected by the movement of overlying strata over the gob.

5.4. In Situ Monitoring

After adopting the technology of GSER by roof cutting, it is necessary to monitor the deformation and stress of the roadway and working face to further formulate protective measures. The monitoring includes two aspects: the monitoring of the working face pressure and the appearance of the mine pressure monitoring in the retained roadway, including the stress of a single pillar and anchor cable, the pressure of retaining gangue, and the deformation of the roadway. Figure 15 shows the layout of the monitoring points. The working face hydraulic support is numbered from the cutting side to the rightmost side, starting from bracket 3#, and for every 10 brackets, a monitoring point is set up. Because the single prop is installed after mining, it lags behind the working face and sets a monitoring point every 20 m. The gangue retaining pressure monitoring points lag the working face and are arranged every 20 m. The layout of the CRLD cable stress monitoring points and roadway deformation monitoring points are the same. Starting from the position of the open-off cut, one monitoring point is set every 50 m.

5.5. Working Face Pressure Monitoring

Figure 16 shows the stress curves of the 5# hydraulic support. From the stress monitoring curve, when the roof cutting was performed, the stress value clearly decreased, the peak stress decreased from 56 MPa to 44 MPa, the average stress decreased from 43 MPa to 37.5 MPa, and its average periodic weighting length increased from 12 m to 21 m. The increase in the periodic weighting step indicates that under the influence of roof cutting, the overlying strata collapse height beside the cutting surface was large, the rock fragmentation was small, and the goaf filling effect was good. Therefore, there was less space for the main roof to rotate, the rotation angle was smaller, and the rotation deformation was also relatively small, resulting in the main roof difficult to fracture; that is, the fracture weighting step increased.

5.6. Mine Pressure Monitoring in Roadway
5.6.1. Stress Analysis of Single Hydraulic Prop and CRLD Cable

Figure 17 shows the stress curve of the CRLD anchor cable and single prop. The CRLD anchor cable measuring point is marked as M4 and is 230 m away from the open-off cut. The single prop measuring points are 690, 710, 730, 750, and 770 m away from the open-off cut. The cable curve can be divided into four stages: advanced impact, intense impact (the support resistance sharply increases), steadily increasing resistance, and stabilization stages. The intense impact zone is within 30 m before and after the working face, and the lead abutment pressure and roof movement caused by mining at this stage have the highest impact. When the measuring point lags the working face by approximately 100 m, the anchor cable enters the constant resistance stage. The stress variation of a single prop is consistent with that of the CRLD anchor cable. Within the range of 20 m behind the working face, the prop pressure sharply increases. After lagging the working face at 25 m, the prop pressure basically reaches a stable state, which indicates that the strata movement is initially stable.

5.6.2. Stress Analysis of Retaining Gangue

The caved gangue is used as part of the retained roadway, and the study of the gangue pressure on the lateral support system and the gangue collapse degree have important significance for the stability. Therefore, it is necessary to conduct dynamic monitoring and research on the gangue pressure on the lateral support and the coefficient of bulk increase in the field. Figure 18 shows the stress of retaining gangue and the coefficient of the bulk increase curve. When the measuring point lags behind the working face by approximately 100 m, the gangue retaining pressure and coefficient gradually reach a stable state and exhibit a certain synergy. The maximum pressure is 2.9 MPa, and the average value is 1.93 MPa; the final coefficient of bulk increase is approximately 1.32.

5.6.3. Deformation Analysis of Retained Roadway

Figure 19 shows the displacement curve of the 8# measuring point. The measuring point is 160 m distant from the open-off cut. From the monitoring curve, it can be seen that the amount of roof subsidence is approximately 200 mm, the amount of floor heave is approximately 313 mm, and the convergence between roof and floor is 513 mm. When a single pillar is withdrawn, the deformation rate of the roadway increases, but the increase is small. When the measuring point lags behind the working face by approximately 300 m, the deformation of the roadway gradually stabilizes. Figure 20 shows the effect of the retained roadways in this coal mine.

6. Conclusions

This study mainly analyzes the application of roof-cutting-type GSER in the medium-thick coal seam of the Dianping coal mine, in the Liliu mining area, and Shanxi province. First, the technological process and its technical principle were studied in detail, and mechanical analysis of the surrounding rock structure was performed. On this basis, the field parameter was designed, and the displacement and stress evolution law under the technical condition of the reserved roadway was studied using a numerical model. Finally, field monitoring was used to verify and evaluate the effect of the reserved roadway. The conclusions are as follows:(1)The principle of roof-cutting-type GSER is to cut off the voussoir beam structures which were used to transmit horizontal thrust above the retained roadway due to the coal mining to achieve the purpose of optimizing the stress distribution of surrounding rock. The roof-cutting-type GSER has three core technologies: CRLD anchor cable strengthening support, roof cutting by directional blasting, and gangue retaining support, and through the study of these technologies, the basic technical parameters were determined; i.e., two rows of CRLD anchor cables were arranged; with the spacing of 1000 mm, the roof-cutting height and angle were designed as 9 m and 15°, respectively(2)Based on the analysis of the surrounding rock structure formed by roof cutting, a mechanical model of the cantilever beam was established, and the calculation method of material mechanics was used to solve the support resistance and structural deformation. It is concluded that the closer to the roof-cutting side, the higher the displacement deformation(3)Using FLAC3 D numerical simulation, the stress distribution laws under the three processes of retaining coal pillar, traditional wall building type GSER, and roof-cutting-type GSER were compared and analyzed. Under the condition of the remaining coal pillar and traditional GSER, the advanced stress distribution and value of the working face are close to each other, while under the condition of roof cutting, the advanced stress at the roof-cutting side decreases obviously, and its transverse influence range is about 60 m. On the basis of the above analysis, a simulation study was carried out on the roof-cutting-type GSER with different heights and angles, it was verified that the optimal cutting height was greater than 8 m and less than 10 m, and the optimal cutting angle was 15°. Above the analysis shows that roof cutting can effectively cut off the stress transfer between the roadway roof and the goaf, reduce the stress concentration degree, make the stress concentration area transfer to the solid coal, that is, to the side of roof cutting, and enhance the stability of the roadway roof(4)On-site monitoring was performed, including the working face pressure, CRLD anchor cable, and single prop stress and roadway deformation. The field monitoring results showed that the preliminary weighting step of the roof-cutting side was larger than that of the uncut side, and the maximum mining face pressure was lower. The amount of roof subsidence was approximately 200 mm, the amount of floor heave was approximately 313 mm, and the overall deformation was small. This indicates that the effect of the retained roadway was good [31, 32]

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This research was funded by the National Natural Science Foundation of China, under grant no. 51574248, and the National Key Research and Development Program, under grant no. 2016YFC0600901.