Abstract

The water-blocking properties of the clay layer at the bottom of the Cenozoic overburden in China are an important factor influencing the safety of thin bedrock coal seam mining. Clay has remolding properties that are unlike the nonreversible characteristics of cracks in brittle rock, and failure cracks in clay can reclose or continue to expand under the influence of different external factors. In this work, the soil layer on top of thin bedrock is the research object, and the influences of the particle composition, water content, soil layer thickness, and crack width on the crack development-closure state of soil layer are analyzed by the orthogonal test method. Visual analysis shows that the order of influence of each factor on the stability of soil layer is the crack width, particle composition, soil layer thickness, and water content. The stability of soil layer decreases with increasing crack width and sand content and decreasing soil layer thickness; in addition, soil layer stability decreases first and then increases with increasing water content. Further variance analysis shows that the crack width and particle composition are key factors that impact the stability of soil layer and that the soil layer thickness has some influence, while the water content has little effect on the stability of soil layer. In addition, the crack will reclose when the sand content in soil is less than 50% and the crack width is less than or equal to 1.0 mm, and the soil layer is prone to further failure when the sand content in soil is more than 50% and the crack width is greater than or equal to 3.0 mm; when the soil layer thickness is 15.0 cm, its stability is better than when the soil layer thickness is 10.0 cm or 5.0 cm.

1. Introduction

Environmental problems, such as water resource shortages and surface degradation caused by coal mining, are becoming increasingly prominent as the intensity of coal mining in China increases. The aforementioned is especially true in midwestern areas, where underground latent water resources are widely distributed in the Cenozoic overburden. The damage to underground latent water resources is particularly evident in the coal seam mining process [14]. Therefore, ideas and methods, such as “green coal mining [57]” and the “coexistence of coal and water [811],” have been proposed. These ideas and methods have provided new ways of protecting underground latent water resources.

To find solutions to specific problems, such as maintaining the safety of coal seam mining and researching the water-blocking mechanisms of clay, scholars in China and abroad have carried out useful studies by establishing safety assessment methods, conducting field measurements, and running numerical simulations. For example, the concept of the risk coefficient of water inrush in loose porous aquifers and the clay at the bottom of the Quaternary system that has been put forward can be used as part of the mining protection layer [1215]. These studies provide guidance to improve the extent of coal seam mining under thin bedrock.

Simultaneously, much research has also been conducted on the protection of water resources in the process of mining coal resources. The influence of coal mining on underground latent water resources were divided into four categories, namely, serious, moderate, slight, and no water loss, and an aquifer-protection mining technique can be successfully applied by modifying a few mining parameters, such as the mining height or advance rate. Moreover, the key factor in protecting underground latent water resources is the thickness of the weathered bedrock located immediately below the aquifer [1620]. Mining activities can cause not only surface and subsurface water loss but also chemical, trace metal, and microbiological pollution of surface and subsurface water [2125].

The water flowing fractured zone usually penetrates through thin bedrock while mining coal seams under thin bedrock. At this time, the clay layer at the bottom of the Cenozoic overburden becomes an important water-blocking structure and plays an important role in preventing the water in porous aquifers and phreatic aquifers from flowing to the working face [26]. In addition, the properties, distribution characteristics, and mining failure characteristics of clay have important influences on its water-blocking performance [2730]. Unlike the development of dry crack or the development of crack in slopes, retaining walls, and other soil engineering structures [3135], the direction of mining-induced shear/tension crack in the overlying soil layer during thin bedrock coal seam mining is generally nearly vertical and the crack is subject to erosion and scouring by upper phreatic aquifer water [36]. In addition, clay has remolding properties, which are unlike the nonreversible crack in brittle rock [3740]. The crack in soil will reclose under the action of stress and water swelling property and other factors. For example, Huang et al. [41] found that the crack in the soil swelled with water and the crack reclosed when the amount of expansion exceeded the width of the crack. Zhang et al. [42, 43] studied the functional relationship between the soil suction and resilient modulus by repeated load triaxial test and considered that the resilient modulus increases with the increase in the minimum bulk stress, soil suction, and degree of compaction and decreases with increasing octahedral shear stress.

The characteristics of crack development-closure in soil layer are of great significance to the stability of the overlying soil and its water-blocking effect in thin bedrock coal seam mining. In this paper, four factors that affect the development-closure state of crack, namely, particle composition, soil layer thickness, crack width, and water content, are examined. The crack development-closure characteristics of soil layer are analyzed by the orthogonal test method. The aim of the work is to provide a reference for the stability analysis of overlying soil and the protection of underground latent water resources in thin bedrock geological regions during thin bedrock coal seam mining.

2. Hydrogeological Conditions of Thin Bedrock Coal Seam Mining

2.1. Thin Bedrock Definition

At present, there is no clear standard for the definition of thin bedrock, and researchers have different definitions of thin bedrock according to their own research purposes. For instance, thin bedrock can be defined according to the height of three zones based on the view of water prevention [44]. Similarly, researchers define thin bedrock according to whether there exists a key stratum based on the view of mine pressure control [15]. This article defines thin bedrock (as shown in Figure 1) in terms of the water prevention view [44], and (1) if the bedrock thickness Hb is smaller than the height of the caving zone Hc, it is called ultrathin bedrock; (2) if the bedrock thickness Hb is larger than the height of the caving zone Hc but smaller than the height of the water flowing fractured zone Hc + Hf, it is called thin bedrock; and moreover, (3) if the bedrock thickness Hb is larger than the height of the water flowing fractured zone Hc + Hf, it is called normal thickness bedrock. Ultrathin bedrock and thin bedrock are commonly referred to as thin bedrock.

2.2. Hydrogeological Conditions of Thin Bedrock

The Carboniferous and Permian (C3-P1) and Jurassic (J1-J2) are the main coal forming periods in China (as shown in Figure 2). Tectonism leads to interruption or denudation in the strata deposit after the deposition of coal measure strata has been completed. As a result, some of the strata between the coal measure strata and the Cenozoic overburden are missing, and the Cenozoic overburden is directly deposited on the coal measure strata (as shown in Figure 3).

The water flowing fractured zone generally penetrates through thin bedrock and into the Cenozoic overburden with the progression of coal mining under such geological conditions. Accordingly, the failure state and water-blocking properties of soil in the Cenozoic overburden are key factors in determining whether a hydraulic connection between the Cenozoic porous aquifers and the working face will occur. Typical hydrogeological conditions of thin bedrock are shown in Figure 3.

2.3. Failure Process of the Overlying Soil of Thin Bedrock

The overlying soil layer of thin bedrock will deform or even fail with the deformation and caving of the bedrock as the coal seam is continuously exploited, and the process of deformation and failure of the overlying soil layer can be divided into three stages, as shown in Figure 4. First, the coal seam is mined out. Second, the bedrock collapses or fractures. Third, the collapse or fracture of the thin bedrock leads to a reduction in the supporting space of the overlying soil layers, and the overlying soil layers will subside due to the self-weight load. Then, the overlying soil layers undergo deformation and failure. Consequently, the essential cause of deformation or failure of the soil is the reduction in the supporting space in its lower part.

The overlying soil will deform and produce cracks under the influence of coal seam mining. The cracks may reclose or continue to broaden by scouring and finally lead to further soil destruction under the action of water flow scouring in loose porous aquifer and other conditions. The key factors affecting the development-closure state of the cracks in soil layer are particle composition, dry-wet state (which can be expressed by the saturation or water content), crack dimension (crack length and width), and soil stress state.

The main external forces are the self-weight stress and phreatic aquifer hydrostatic pressure in the overlying soil of thin bedrock, and their influence mode is long-term slow action. Therefore, the influence of the dynamic change in external forces can be temporarily neglected. Then, four factors can be selected, namely, the soil particle composition, water content, crack length, and crack width, and influence analysis of these factors on the development-closure state of the cracks in soil layer can be performed through orthogonal experimental methods.

3. Experimental Device and Scheme

3.1. Experimental Device and Specimen Preparation

The crack development-closure experimental device for soil layer is shown in Figure 5.

The main functions of the experimental device are as follows: (1) the water tanks I and II, the water inlet pipe, and the overflow pipe keep the water level above the soil layer constant; that is, a constant hydrostatic pressure is maintained and (2) the volume of outflow water from the cracks is measured by using a measuring cylinder, and the flow volume can be calculated in terms of the volume of outflow water.

The change in soil particle composition is simulated by changing the ratio of bentonite clay to fine sand (0.075 mm < particle diameter ≤ 0.25 mm), which is shown in Figure 6. The higher the sand content is, the lower the content of hydrophilic minerals in the soil is, and the worse the hydrophilicity and water swelling properties of the soil are.

The water content is realized by adding water. First, the dry weight of the soil layer is calculated in terms of the soil density (1.884 g/cm3). Then, the quality of water needed to added can be calculated according to the water content. The crack length means the length of the crack in the vertical direction, and the cracks in the soil layer are actually not straight. The cracks made in experiments are straight cracks because it is difficult to make bending cracks in specimens; therefore, the thickness of the soil layer is equal to the crack length. A crack is formed when the iron wire penetrates the soil layer after the soil layer was compacted and formed according to the design thickness (as shown in Figure 7).

3.2. Experimental Scheme

The symbols of each factor in the experiments and their level values are shown in Table 1.

There are four factors in the experiments, and each factor has three levels. The purpose of the experiments is to analyze the influence of four factors on the development-closure state of the cracks in soil layer. Simultaneously, the effect of the interaction between the particle composition and the other factors on the crack development-closure state is emphatically investigated. The free degree of four experimental factors (each of which has three levels) and four sets of interactions is ; therefore, the L27(313) orthogonal table can be used to design the experimental scheme, which is shown in Table 2.

The dependent variable of the experiments is the development-closure state of the cracks in soil layer, that is, the change in crack dimension, and the water rate of flow Q can be used to indicate the change in soil crack dimension (the water rate of flow is defined as the volume of fluid flow through cracks per unit time, and the unit is m3/s). Then, the development-closure state of the cracks in soil layer can be determined according to the change trend of Q. (1) The water rate of flow Q decreases gradually to 0 m3/s, which indicates that the crack dimension gradually decreases to closure. (2) The water rate of flow Q gradually increases to  m3/s, which means that the crack gradually broadens due to water flow scouring and finally leads to soil further destruction. (3) The water rate of flow Q is between 0 m3/s and  m3/s, which indicates that the crack is ultimately neither reclosed nor broadened by water flow scouring.

4. Experimental Results and Discussion

4.1. Direct Analysis of the Experimental Results
4.1.1. Experimental Results

The water rate of flow results at 6 h for each experiment is shown in Table 3.

Table 3 shows that the water rate of flow for 6 of the 27 experiments (experiments 1, 4, 7, 10, 13, and 16) gradually decreases to 0 m3/s, which indicates that the cracks in the soil finally reclose under the experimental conditions of these 6 experiments. There are 19 experiments (experiments 2, 3, 5, 6, 8, 9, 11, 12, 14, 15, 17, 18, 19, 20, 22, 23, 25, 26, and 27) whose water rate of flow is between 0 m3/s and  m3/s, which indicates that the cracks are ultimately neither reclosed nor broadened by water flow scouring under these 19 experimental conditions. The water rate of flow for 2 of the 27 experiments (experiments 21 and 24) gradually increases to  m3/s, which indicates that the cracks gradually broaden due to water flow scouring and finally lead to soil destruction under the experimental conditions of these 2 experiments.

4.1.2. Analysis of Experiments That the Water Rate of Flow Decreases Gradually to 0 m3/s

One comparison of the initial cracks and reclosed cracks in the experiments in which the water rate of flow decreases gradually to 0 m3/s is shown in Figure 8.

The water rate of flow curves changes with time in the experiments in which the water rate of flow gradually decreases to 0 m3/s is shown in Figure 9.

As shown in Figure 9, the water rate of flow of the 6 groups gradually decreases to 0 m3/s (that is, the cracks reclosed in these experiments) under the same experimental conditions of or , which indicates that when the ratio of soil to sand is greater than or equal to 1 : 1 and the crack width is less than or equal to 1.0 mm, the cracks are reclosed regardless of the thickness and water content of the soil. At the same time, the period with a sharp decrease in the water rate of flow in the 6 experiments is less than 20 minutes, which shows that the duration of the main change in the soil swelling increment (or water swelling velocity rate ) and soil compaction amount (or compression velocity rate ) is less than 20 minutes, namely, the sum of the swelling increment and compaction amount is larger than the sum of the broadening amount and original crack width , which results in the cracks reclosed.

4.1.3. Analysis of Experiments That the Water Rate of Flow between 0 m3/s and  m3/s

A comparison of the initial and stable cracks in the experiments in which the water rate of flow is between 0 m3/s and  m3/s is shown in Figure 10.

The water rate of flow curves changes with time in the experiments in which the water rate of flow is between 0 m3/s and  m3/s is shown in Figure 11.

As shown in Figure 11, the water rate of flow of experiment 27 is the largest and is 12.167 (10−6 m3/s), while that of experiment 25 is the smallest and is 0.008 (10−6 m3/s). Intuitively, it cannot be determined which factor is the main factor and which factor is the secondary factor leading to the difference in the rate of flow results. Thus, further analysis of the experimental results is needed to determine the order of influence of each experimental factor on the water rate of flow.

4.1.4. Analysis of Experiments That the Soil Layer Destroyed

One comparison of the initial cracks and destroyed soil in the experiments in which the water rate of flow increases gradually to  m3/s is shown in Figure 12.

The water rate of flow curves changes with time in the experiments in which the water rate of flow increases to  m3/s is shown in Figure 13.

As shown in Figure 13, the condition for cracks to be broadened by scouring causing further soil destruction is , and experiment 21 shows that the time from broadening of the crack by scouring to soil destruction is 25 minutes, while that in experiment 24, it is 310 minutes, which indicates that increasing the soil layer thickness or increasing the water content or both together can effectively increase the antidestructive capacity of the soil layer. Further analysis is needed to determine which conditions change to improve the antidestructive capacity of soil layer.

4.1.5. Analysis of Experimental Results

The experimental results show that the crack will close when the sand content in soil is less than 50% and the crack width is less than or equal to 1.0 mm and the soil layer is prone to further failure when the sand content in soil is more than 50% and the crack width is greater than or equal to 3.0 mm. The experimental results preliminarily reveal that the greater the sand content and the wider the crack width and the smaller the thickness of the soil layer, the worse the stability of the soil layer, because the orthogonal test method is adopted instead of the comprehensive test method, and visual analysis and variance analysis are needed to get more accurate conclusions.

4.2. Visual Analysis of the Experimental Results

There are 27 groups of experimental results in Table 3, and no two experiments among the 27 groups have identical experimental conditions; thus, it is not feasible to compare any of the 27 groups of experimental results directly. However, if these 27 sets of experimental results are combined, comparability of the orthogonal test results will be achieved.

Taking Cc as an example and level 1 (Cc10) of experimental factor Cc appeared in experiment Nos. 1–9 in Table 3, the average of the rate of flow of these nine experiments islevel 2 (Cc11) of experimental factor Cc appeared in experiment Nos. 10–18 in Table 3, the average of the rate of flow of these nine experiments isand level 3 (Cc13) of experimental factor Cc appeared in experiment Nos. 19–27 in Table 3, the average of the rate of flow of these nine experiments is

The soil layer in experiment No. 21 and No. 24 is destroyed, and their experimental results are  m3/s. Because the value cannot participate in the calculation, the rate of flow result in experiment No. 21 and No. 24 is calculated to be a larger value of 100 (10−6 m3/s).

, namely, the rate of flow is the largest when the particle composition is Cc13, and the rate of flow is the smallest when the particle composition is Cc10. This shows that the soil layer is prone to further failure when the particle composition is Cc13, and the soil layer is prone to maintain stability when the particle composition is Cc10.

The above calculations can also be performed in the orthogonal table, as shown in Table 4. Ij represents the sum of the experimental rate of flow results corresponding to level 1 of column j, and is the average of Ij; IIj represents the sum of the experimental rate of flow results corresponding to level 2 of column j, and is the average of IIj; IIIj represents the sum of the experimental rate of flow results corresponding to level 3 of column j, and is the average of IIIj.

The experimental factor Cc is in the first column; therefore, is , is , and is . Similarly, Cl, Cw, and are in the second, fifth, and tenth columns, respectively. Hence, is , is , and is , while is , is , and is ; in addition, is , is , and is .

The , , and values of each experimental factor in Table 4 are plotted using the longitudinal coordinate system, as shown in Figure 14.

It is easy to intuitively observe that if one experimental factor has a great influence on the experimental result, then this experimental factor is the main factor. Namely, the difference between the average rate of flow corresponding to different levels of the experimental factors is large, which is reflected in Figure 14, and the points corresponding to the three levels of the experimental factors differ greatly in the longitudinal coordinates. In contrast, if one experimental factor has little influence on the experimental result, then this experimental factor is the secondary factor. That is, the difference between the average rate of flow corresponding to different levels of the experimental factors is small, which is reflect in Figure 14, and the points corresponding to the three levels of the experimental factors differ little in the longitudinal coordinates.

The absolute average values of the difference between , , and in Table 4 are

The maximum absolute average value of the difference between the three levels of each experimental factor is factor Cw, followed by factors Cc and Cl and factor . In addition, the influences of Cw and Cc on the experimental results are more notable than those of Cl and . The result indicates that the order of influence of each factor on the stability of soil layer is the crack width, particle composition, soil layer thickness, and water content and also indicates that the influence of crack width and particle composition on the failure of soil layer is greater than that of the soil layer thickness and water content.

As seen from Figure 14, the order of influence of the particle composition on the rate of flow is , which shows that the higher the sand content in the soil is, the worse the stability of the soil layer is. The order of influence of the soil layer thickness on the rate of flow is ; that is, the smaller the soil layer thickness is, the more likely it is that the soil layer will fail. The order of influence of the crack width on the rate of flow is , which indicates that the larger the crack width is, the worse the stability of the soil layer is. The order of influence of the water content on the rate of flow is , which indicates that the stability of the soil layer decreases first and then increases with increasing water content and that soil layer is prone to further failure when the water content is 20%.

In addition, Table 4 also shows that (III1) is significantly larger than (I1) and (II1), and (III5) is significantly larger than (I5) and (II5), while (III2) is significantly smaller than (I2) and (II2). These observations indicate that when the sand content in the soil is higher than 50% and the crack width is larger than 2.0 mm, the soil layer is prone to further failure. When the soil layer thickness is 15.0 cm, the stability of the soil layer is better than those of the 10.0 cm or 5.0 cm; that is, increasing the thickness of the soil layer can increase the ability of soil layer to resist further failure.

4.3. Variance Analysis of the Experimental Results
4.3.1. Sum of Squares of Deviations of Experimental Factors

Taking Cc as an example, the calculation method of the sum of squares of deviations of experimental factor is illustrated. Level 1 (Cc10) of experimental factor Cc appeared in experiment Nos. 1–9 in Table 3, and the sum of the rate of flow of the nine groups of experimental results islevel 2 (Cc11) of experimental factor Cc appeared in experiment Nos. 10–18 in Table 3, and the sum of the rate of flow of the nine groups of experimental results islevel 3 (Cc13) of experimental factor Cc appeared in experiment Nos. 19–27 in Table 3, and the sum of the rate of flow of the nine groups of experimental results is

The soil layer in experiment No. 21 and No. 24 is destroyed, and their experimental results are  m3/s. Because the value cannot participate in the calculation, the rate of flow result in experiment No. 21 and No. 24 is calculated to be a larger value of 100 (10−6 m3/s).

Hence, the sum of squares of deviations of experimental factor Cc iswhere is the average rate of flow of the 27 groups of experimental results and .

Expansion of formula (8) results in

Let be factor G and be factor CT; then, formula (9) can be simplified as follows:

The data in Table 3 are substituted into formula (10) to obtain .

In the same way, the sum of squares of deviations of experimental factors Cl, Cw, and ω can be calculated as

The sum of squares of deviations of the above experimental factors can also be calculated on the orthogonal table, which is shown in Table 5.

As shown in Table 5, Ij represents the sum of the experimental results corresponding to level 1 of column j, and IIj represents the sum of the experimental results corresponding to level 2 of column j, while IIIj represents the sum of the experimental results corresponding to level 3 of column j. Therefore, I1 is , II1 is , and III1 is . Sj represents the sum of squares of deviations of the experimental results in column j; hence, S1 is SCc, S2 is SCl, S5 is SCw, and S10 is Sω.

4.3.2. Sum of Squares of Deviations of the Interaction between the Experimental Factors

The interaction between two experimental factors accounts for two columns in the orthogonal table; therefore, the sum of squares of deviations of the interaction between two factors is the sum of the sum of squares of deviations of the two columns. Therefore, the sum of squares of deviations of SCc×Cl is

Similarly, , , and can be obtained:

The sum of squares of deviations of the interaction between experimental factors can also be calculated on the orthogonal table, which is shown in Table 6.

4.3.3. Sum of Squares of Deviations of the Experimental Error

The sum of squares of deviations of the experimental error can be calculated by the sum of squares of deviations of the blank columns without experimental factors arranged in the orthogonal table. In this experiment, the error can be estimated by the sum of squares of deviations in column 13 of Table 3 because column 13 does not arrange the experimental factor, and the sum of squares of deviations of column 13 does not include the difference between the levels of experimental factors but only reflects the magnitude of experimental errors.

In addition, if the sum of squares of deviations of the other columns is close to the sum of squares of deviations of the blank column (column 13); then, the sum of squares of deviations of this column can be combined with the sum of squares of deviations of the blank column as an error estimate, which can make the error estimate more accurate.

According to Tables 5 and 6, the sum of squares of deviations of columns 4, 8, and 9 in this experiment is close to the sum of squares of deviations of the blank column (column 13). Therefore, S4, S8, and S9 can also be classified as errors as follows:

The data in Tables 5 and 6 are substituted into formula (14) to obtain .

4.3.4. Saliency Detection of Experimental Factors and Their Interactions

According to formula (8) for the calculation the sum of squares of deviations, it can be seen that the sum of squares of deviations may be large if the amount of data is large. That is, the sum of squares of deviations is related not only to the change in data itself but also to the amount of data. Therefore, the influence of the amount of data on the sum of squares of deviations is eliminated by introducing the free degree in orthogonal experiment, namely, comparing the average sum of squares of deviations with , where fA is the free degree of the sum of squares of deviations of experimental factor A and fe is the free degree of the sum of squares of deviations of experimental error.

For orthogonal experiments, there arewhere is the total number of experiments and is the level amount of experimental factor A.

Therefore, for this experiment, , , , , , , , , and .

According to formula (14),

The ratio of the average sum of squares of deviations of the experimental factors () to the average sum of squares of deviations of the errors () is called the F ratio, that is,

The F ratio of the experimental factors and the interaction between the experimental factors is calculated in terms of formula (17) as follows

Look up the F ratio table, and it shows the following: , , , , , and . Thus, the following conclusions can be drawn:, which indicates that the change in the experimental factor of soil particle composition Cc has a highly salient influence on the experimental results, which is recorded as ., which indicates that the change in the experimental factor of soil layer thickness Cl has some influence on the experimental results, which is recorded as ., which indicates that the change in the experimental factor of crack width Cw has a highly salient influence on the experimental results, which is recorded as ., which indicates that the change in the experimental factor of soil water content has a nonsalient influence on the experimental results, which is recorded as ⊗., which indicates that the interaction between soil particle composition Cc and soil layer thickness Cl has some influence on the experimental results, which is recorded as ., which indicates that the interaction between soil particle composition Cc and crack width Cw has a highly salient influence on the experimental results, which is recorded as ., which indicates that the interaction between soil particle composition Cc and soil water content has a nonsalient influence on the experimental results, which is recorded as ⊗., which indicates that the interaction between soil layer thickness Cl and crack width Cw has a nonsalient influence on the experimental results, which is recorded as ⊗.

The saliency detection of the above experimental factors and their interactions can be summarized as an analysis of variance table, as shown in Table 7.

5. Conclusions

In this work, based on the crack development-closure experiments, the following conclusions are obtained.

Firstly, soil particle composition and crack width have a highly salient influence on the stability of soil layer, and soil layer thickness has some influence on its stability, while water content has a nonsalient influence on the stability of soil layer. As the crack width and sand content increase and the soil layer thickness decreases, the stability of the soil layer decreases.

Secondly, the crack will reclose when the sand content in soil is less than 50% and the crack width is less than or equal to 1.0 mm, and the soil layer is prone to further failure when the sand content in soil is more than 50% and the crack width is greater than or equal to 3.0 mm.

Thirdly, the interaction of experimental factors Cc and Cw has a highly salient influence on the stability of soil layer, while interaction of Cc and Cl has some influence on the stability of soil layer, and interaction of Cc and ω and interaction of Cl and Cw have nonsalient influence on the stability of soil layer.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (41977238 and 51804339) and the China Postdoctoral Science Foundation (2019T120715 and 2018M640760).