Abstract
Static fracturing technology uses chemical expansion agents to fracture roofs. With the aim of fracturing corner roofs on deep working faces, in this study, the static fracturing technology was investigated through theoretical analysis, laboratory experiments, numerical calculations, and field practice. The theoretical analysis and experiments demonstrated that the swelling force increased with a decrease in the fracturing hole spacing, and the optimal water-cement ratio was 0.33. Twelve groups of FLAC3D models were designed using SPSSAU. The results revealed that the optimal fracturing effect was achieved when the hole diameter was 60 mm, hole spacing was 40 cm, and hole depth was 6 m. The fracturing effect of hard corner roofs was monitored by peering into the borehole and evaluating the support resistance. Thus, it can be concluded that within the fracturing range, internal fissures in the rock stratum are developed and linked to each other. The support pressure was the highest, 7 h after grouting, with a value of approximately 26.1 MPa, and then decreased gradually to 17.58 MPa, indicating that the static fracturing technology attained the expected results.
1. Introduction
The shifting of working faces results in the formation of a “large space” structure inside the overburden, which can easily cause the roof to hang at the corners of the working face. Affected by the roof depth and plate hardness, a hanging roof may cause an increase in the supporting pressure of the working face, gas accumulation, and stress concentration of coal pillars, resulting in potential hazards for mining [1–3]. The conventional materials used in mining are explosives and water, which can also create hazards in mines. Therefore, the most pressing task is to identify a new type of fracturing material. In addition, research on hard roof fracturing technology is of immense significance for the improvement of the working face environment and for safe and green mining production.
China is rich in coal, and the conditions of each mine are considerably different. Different fracturing technologies, such as hydraulic fracturing [4, 5], blasting fracturing [6, 7], and CO2 presplitting [8], are used on a case-by-case basis and achieved highly desirable results. A static cracking agent, also known as a chemical expansive agent, is a nonexplosive, noncombustible, and nontoxic powdered chemical mixture with high expansion performance. Tanaka was the first to perform a volumetric expansion experiment using the oxidation processes of CaO and H2O in concrete fracturing [9]. After long-term development, Chinese experts have developed mixing materials for the expansive agent and have studied the factors that affect the fracturing effect [10–12]. Early Chinese research on expansive agent fracturing was mainly applied to the cutting of rock materials [13]. With the significant challenges caused by the deep mining of coal resources [14, 15], by exploring the high ground stress, high ground temperature, high water pressure, and strong production disturbance in the deep courses [16, 17], experts worldwide have analysed the large deformation of deep rock and features of the unloading failure zone of the arch form [18, 19] and development characteristics of fractures [20, 21] and have determined the energy characteristics and criteria for the impact of depth on lithomechanical behaviour and dynamic instability [22–24]. Despite long-term research and practice, few studies have been conducted on the application of static fracturing technology in mines and adjoining rock control. Furthermore, research has been conducted on the prevention and control of deep dynamic disasters [25–28]. Considering the relatively low reaction speed and environmentally friendly characteristics of static crushing agents, combined with deep adjoining rock control, static fracturing technology for hard roofs in deep working faces has been introduced.
To realise the safe, green, and efficient control of the hard roofs of deep working faces, in this study, an experiment was conducted by fracturing a roof at a corner. Static fracturing technology was used, and factors affecting the expansion effect were determined through theoretical analysis and laboratory tests. In addition, the fracturing parameters suitable for a mine were obtained, in combination with the affecting factors, and SPSSAU was used to design a numerical calculation scheme. Finally, the effect of applying the static fracturing technology was verified on a working site.
2. Project Background
2.1. Overview of Mine and Working Face Production
The roof structures of the Binchang, Huangling, and Weibei mining areas in Shaanxi Province are considerably different. However, their uniaxial compressive strengths are generally greater than 45 MPa. In addition, a mine owned by Shaanxi Coal and Chemical Industry Group Co., Ltd., with an approved production capacity of 2 Mt/a, is the main production unit of the group. It mainly mines at the No. 2 coal seam, which is within the Yan’an Formation of the middle Jurassic period, with an inclination angle of 1–5°. The coal seam ranges from stable to more stable, and the mine has a high gas content. The No. 21209 working face is 395–699 m deep and is adjacent to the No. 21207 gob in the southwest and the No. 21211 working face in the northeast. The rest of the coal seam is unmined. The surface elevation of the working face is +1175 m to +1424 m, and the underground elevation is +723 m to +784 m. The ground surface is covered with medium and low mountain forests, with no buildings and other facilities around the mine. Figure 1 depicts the location of the No. 21209 working face and the overburden characteristics.
For the No. 21209 working face, the designed strike length was approximately 1235 m, incline length was approximately 153 m, and average coal thickness was 1.1 m. The dimensions (width × height) of the air inlet and return roadways were 4.6 × 3.8 and 5.4 × 3.6 m, respectively. The coal mining method used was strike longwall retreat full-seam mining. The roof of the gob was treated using a fully caving method. Table 1 lists the characteristics of the adjoining rock in the No. 21209 working face. Lithomechanical experiments demonstrated that the natural compressive strength of the roof and coal seam floor of the No. 21209 working face was 38.6–57.5 MPa, predominantly in the range of 42.3–50.8 MPa, and the single shear strength was generally 4.22–6.4 MPa.
2.2. General Conditions of Mine
With the continuous alternation of the working face and the incomplete collapse of the adjacent working face, the roof of the working face is continuous and forms a “large space” structure inside the overburden with the collapse space. With the backstopping of the working face, the scale of the extraction area increases. The large space structure inside the overburden is easy to cause the stage stress failure of roof plate at the corner of working face, and the length of the overhanging roof at the upper corner is larger than that at the lower corner. With varying degrees of mining along the gob, it is easy to cause stage stress failure in the roof at the corners of the working face, and the suspended roof length of the upper corner is greater than that of the lower corner. As the working face continues to advance with the strike, the corner roof extends to the gob and flexes (ω). When the flexing of the roof reaches the bending limit, the roof at the corner collapses. The long-term suspended roof at the corner causes severe deformation of the adjoining rock of the two roadways, which can easily increase the support pressure of the working face, gas accumulation, and stress concentration of the coal pillar, gradually producing different strata behaviours (see Figure 2).
(a)
(b)
(c)
(d)
3. Mechanism of Static Fracturing Technology
3.1. Mechanism of Static Fracturing Rock Mass
Static fracturing is a chemical reaction that occurs after a chemical expansive agent is mixed with water, which can increase the reaction volume and heat. The expansive agent is a mixture mainly composed of CaO, Ca3SiO3, and MgO. The chemical equation is shown in
During the hydration reaction process of the chemical expansive agent, there are two volumetric model structures, i.e., free state and constrained state, as depicted in Figure 3. In the free state, the molecular structure is uniformly arranged and only expands in volume, without causing swelling pressure. If it is spatially constrained, the molecules are squeezed together, reducing the porosity. It expands instantly in a limited space and generates swelling pressure. If the swelling pressure is greater than the ultimate tensile strength of the constraint, the development of the internal fracture of the object accelerates until it bursts completely.
Constrained by the rock mass, the volumetric expansion produced by the hydration reaction process acts on the constraint. The expansive agent is also under an acting force with the same magnitude but opposite direction. The shape change process of the chemical expansive agent is based on rheology and elastic strain. Affected by the extremely small modulus of elasticity of the hardened solid formed after the hydration reaction of the chemical expansive agent, the swelling pressure generated by the chemical expansive agent under the restraint is caused by the elastic strain alone. In other words, as the proportion of elastic strain that accounts for the total strain increases, the rate of increase of the swelling pressure increases, and the effect of swelling and fracturing improves. In addition, the tensile stress inside the rock mass increases with the continuously increasing swelling pressure. When the ultimate tensile strength at the filling hole wall on the pivoting face is reached, several radial rupturing fissures appear. With the continuous increase in the swelling force, the fissures continue to develop until they are linked to each other, forming a continuous rupturing surface until the swelling and fracturing are completed.
3.2. Mechanical Mechanism of Static Fracturing
In static fracturing, a large swelling force generated by the volumetric expansion of the chemical expansive agent in the hydration reaction causes the constraint to reach the ultimate tensile strength and then rupture. When there exist multiple free faces of the constraint, the tensile stress generated by the volumetric expansion of the chemical expansive agent causes rock mass failure. In addition, rupturing fissures generally develop from the edge of the hole bottom to the free face. In general, porous swelling and fracturing technologies are adopted. Considering double-hole swelling and fracturing as an example, the effects of the physical and mechanical action are illustrated in Figure 4.
Assuming that there are two circular holes on an infinite elastomer under internal pressure, the acting effect of the two holes can be approximately regarded as the superposition of the respective effects on the elastomer.
Then,
From (3), by obtaining the reciprocal of r and equating it to zero, the following is obtained:
This point denotes the stress at the extreme point, which is
If the left boundary r = r1,
If the right boundary r = l − r2,
Then, the minimum stress on the line connecting the centres of the two holes O1 and O2 is
In engineering practices, it is generally considered that the blasting hole and explosives used are the same, and the difference in the charging time is ignored; i.e., when r = r1 = r2, q = q1 = q2.
According to the maximum principal shear stress theory, the crushing medium must meet the conditioni.e.,where [σ] is the ultimate tensile strength of the crushing medium.
Therefore, a higher stress can be generated when the hole spacing is as short as possible. If holes are equally spaced, then the compressive and tensile stresses are at their maximum.
4. Determination of Physical and Mechanical Parameters of Expansive Agent
The swelling force originates from the volumetric expansion caused by the hydration reaction. By comparing the volumetric expansion characteristics under different water-cement ratio conditions, the optimal range of the water-cement ratio is determined. Three indicators—swelling volume, swelling force, and temperature generated by the hydration reaction of the static fracturing of the rock stratum—are used as the basis for determining the parameters of the expansive agent.
4.1. Expansion Characteristics under Different Water-Cement Ratios
Uniformly graduated cylinders (500 ml) were prepared. The chemical expansive agent was mixed according to the four water-cement ratios, and they were stirred until no blocky particles were found. The expansion characteristics of the expansive agent with four different water-cement ratios are presented in Table 2. To ensure the reliability of the test, the test was performed in a room with a temperature change range of less than 2°C. The purified water required for the test was placed in the laboratory for 24 h in advance to ensure that the initial water temperature was consistent with the room temperature.
In the initial state, the slurry with a water-cement ratio of 1 : 1 to 1 : 3 exhibited good fluidity. At a water-cement ratio of 1 : 4, the fluidity was evidently insufficient, and the resistance during mixing was higher. After being rested for approximately 10 min, the slurry exhibited stratification. After continually being rested for 20 minutes, the surface of the sediment gradually dried and started to solidify. With regard to the horizontal swelling characteristics, the cracking became increasingly varied with the increase in the water-cement ratio. In addition, the surface of the 1 : 1 slurry dried and no longer changed after solidification; the 1 : 2 slurry began to bulge in the process and exhibited obvious cracking. The 1 : 3 slurry cracked gently, and the block became larger. The 1 : 4 slurry was too dense, and beaker expanded and cracked, making it impossible to observe the volumetric shape after expansion.
With regard to the vertical expansion characteristics, the expansive volume varied significantly with the increase in the water-cement ratio. The 1 : 3 slurry exhibited a maximum volumetric expansion of approximately 290%, of which the volume expanded by nearly three times. The 1 : 4 slurry demonstrated a volumetric expansion rate of 250%. Therefore, the optimal water-cement ratio was determined to be between 1 : 3 and 1 : 4.
4.2. Determination of Mechanical Properties of Expansive Agent
The optimal water-cement ratio of 1 : 3 to 1 : 4 was further categorised into five different cases, i.e., 0.25, 0.29, 0.33, 0.37, and 0.41, to monitor the swelling force mechanical behaviour to obtain the superior water-cement ratio, as depicted in Figure 5. Five different slurries were stirred and filled into uniform steel pipes for strain monitoring. The specifications of the steel pipes are listed in Table 3.
(a)
(b)
(c)
A temperature sensor and strain tester were used to monitor the strain of the steel pipes in real time during the reaction process until the strain and temperature of the steel pipes stabilised. The temperature and swelling force change curves of the expansive agent are illustrated in Figures 6 and 7 .
The entire temperature change process with different water-cement ratios basically tended to stabilise after 4 h, as depicted in Figure 6. The high temperature generated was approximately 180°C, which was considerably lesser than the ignition point of coal and, thus, could cause the spontaneous combustion of the coal bed. The maximum temperature decreased as the water-cement ratio increased. The maximum temperature was 184.8°C, and the minimum temperature was 144.6°C. Therefore, as the content of the expansive agent increased, the temperature generated during the hydration reaction also increased.
Figure 7 depicts the characteristics of the swelling force changes under the five water-cement ratio conditions with a hole diameter of 40 mm. For the water-cement ratios of 0.41, 0.37, 0.33, 0.29, and 0.25, the generated swelling forces were 25.50, 28.24, 39.14, 32.81, and 30.01 MPa, respectively. Therefore, the maximum swelling force was generated when the water-cement ratio was 0.33.
4.3. Mechanical Effect of Hole Diameter on Swelling Force
The generated swelling force varied with different explosive loads of various hole diameters. To further explore the influence of different hole diameters on the swelling force, the water-cement material that generated the largest swelling force was filled into pipes with hole diameters of 30, 36, 40, 45, and 60 mm to measure the swelling forces.
The swelling characteristics of five hole diameters after 700 min of continuous monitoring are depicted in Figure 8. The swelling force of each hole diameter reached its maximum and stabilised within 275 min. As the hole diameter increased, the time taken to reach the maximum swelling force decreased, and the corresponding swelling force increased. The maximum swelling forces for diameters of 30, 36, 40, 45, and 60 mm were 39.14, 52.24, 60.33, 68.44, and 74.12 MPa, respectively. The distribution characteristics of the maximum swelling forces of five hole diameters are depicted in Figure 9.
Thus, within the range of experimental testing, the hole diameter range is directly proportional to the hole diameter. According to the distribution characteristics of the hole diameter and swelling force, the fitting formula is shown in (12) and the degree of fitting (R2) is approximately 0.9987.
Therefore, if the hole diameter is 53.28 mm, the maximum swelling force is 72.57 MPa. After the hole diameter is increased to 53.28, the swelling force reduces with increasing hole diameter.
As a result, if the optimal water-cement ratio is 0.33, the swelling force is at its maximum when the hole diameter is 53.28 mm.
5. Determination of Static Fracturing Technical Parameters
5.1. Determination of Experimental Program
According to the actual production situation on-site, it was determined that the hole diameters were 42, 50, and 60 mm, the hole spacings were 40, 50, 60, 70, and 80 cm, and the hole depths were 2, 4, 6, and 8 m. The three factors (hole diameter, hole spacing, and hole depth) and levels were input into the orthogonal design software SPSSAU. Twelve sets of experimental programs were obtained, as presented in Table 4.
5.2. FLAC3D Modelling
According to the site conditions, combined with the orthogonal design experiments, the swelling and fracturing effects were studied with different hole diameters, depths, and spacings. The FLAC3D model of double-hole fracturing was established, as illustrated in Figure 10. The dimensions of the design model were 2 × 2 m. The depth (H) varied with the hole depth designed by the orthogonal experiment. Here, R1 and R2 denote the diameters of the two holes, and L denotes the hole spacing. Three indicators, the plastic zone width (ΔL), plastic zone volume (ΔV), and swelling stress (ΔU) among the holes, were determined, and the swelling and fracturing effects of each model were analysed.
(a)
(b)
A 0.5-MPa stress was applied in the x- and z-directions of the model to simulate the actual adjoining rock stress. Through comprehensive analysis of the stress, plastic zone, etc., perpendicular to the direction of the connecting line of the bores, the swelling and fracturing effect of each hole were determined, and the optimal parameters of the holes were determined. The roof parameters of the numerical computation are listed in Table 5.
At the beginning of the simulation, a radial stress with a water-cement ratio of 0.33 was applied to the inside of different hole diameters to simulate the swelling stress. According to fitting formula in (11), the obtained swelling forces corresponding to 42, 50, and 60 mm were 63.78, 71.39, and 74.12 MPa, respectively.
5.3. Distribution Characteristics of Swelling and Fracturing Stress
The fissures generated after the rock stratum was fractured, evolved along the two holes in a linear fashion. The stress disturbance characteristics of the five parameters obtained from the 42, 50, and 60 mm holes were similar, as depicted in Figure 11.
(a)
(b)
(c)
With the same hole diameter, the stress is greater at each hole and lowest on the outside of the hole, generally presenting a “double hump” distribution characteristic. The corresponding stress peak shifts toward both sides with increased spacing between the holes. In addition, the stress peak decreases gradually during the shift. A comparison of the stresses of different parameters with the three hole diameters reveals that the stress disturbance caused by the fracturing is the highest when the hole diameter is 60 mm and is relatively high when the hole spacing is 50 mm.
5.4. Plastic Failure Characteristics of Swelling and Fracturing
The swelling and fracturing of the three hole parameters caused plastic damage to the rock stratum, as presented in Table 6.
The fracturing effect is jointly affected by the hole diameter and hole spacing. The plastic failure range increases with increasing hole diameter and reduces with increasing hole spacing. When the hole diameter is 42 mm and the hole spacing is 40 cm, the plastic failure is linked by a width of approximately 1.09 m, and the failure characteristic presents a “∞” shape. As the hole spacing increases, the plastic failure range is abruptly reduced and not linked. With the increase in hole diameter, if the hole spacing is 40 cm, the plastic zone is linked and continuously expands to the periphery of the hole. For the other hole spacing, the plastic failure range gradually reduces with the increase in the hole spacing.
Therefore, under static fracturing, the fissures of the hole develop with a divergent shape. The swelling force increases with the increase in the hole diameter. The failure characteristic is the most obvious when the hole spacing is 40 cm.
5.5. Influencing Factors of SPSSAU on Fracturing
According to the numerical simulation computed by SPSSAU, the three evaluation indicators (ΔU, ΔL, and ΔV) are presented in Table 7. As theoretically analysed, the range results of the three indicators were obtained. The range expression is shown in (13). As the range increases, the influence of the changed level of a certain factor on the test indicator also increases.
5.5.1. Analysis of Stress Peak ΔU
The range of the tested stress peak ΔU is presented in Table 8. The inspection results of the three factors revealed the following order: factor A (hole diameter) > factor B (hole spacing) > factor C (hole depth) with the smallest range. This indicates that the bore diameter demonstrates the most significant influence on the model stress, followed by hole spacing and, finally, hole depth. Considering that as the stress applied on the rock mass increases, the fracturing effect improves, the optimal combination is determined as A3B1C3—that is, the optimal effect can be obtained with a hole diameter of 60 mm, hole spacing of 40 cm, and hole depth of 6 m.
5.5.2. Analysing Width of Plastic Zone (ΔL)
The ranges of each influencing factor on the width of the plastic zone (ΔL) are 0.52, 1.416667, and 0.4, respectively, as presented in Table 9. Evidently, factor B demonstrates the most significant influence on the experimental indicators. Therefore, factor B (hole spacing) is the main consideration, followed by factor A (hole diameter) and finally factor C (hole depth). The greater the width of the plastic zone is, the more easily the swelling and fracturing lines link to each other. Therefore, the optimal combination is A3B1C2—that is, the optimal effect can be obtained with a hole diameter of 60 mm, hole spacing of 40 cm, and hole depth of 4 m.
5.5.3. Analysing Volume of Plastic Zone ΔU
The influence of each influencing factor on the volume of the plastic zone ΔU is presented in Table 10. The range of the three factors can be arranged as follows: B > C > A. Therefore, hole spacing has the most significant influence. The influence of the hole diameter is slightly lesser than that of the hole depth. Therefore, under the conditions of these evaluation indicators, the optimal combination is A3B1C3 — that is, the optimal effect can be obtained with a hole diameter of 60 mm, hole spacing of 40 cm, and hole depth of 6 m.
5.6. Determination of Optimal Parameters
Twelve static fracturing cases performed using the SPSSAU orthogonal design software reveal that the swelling stress applied around the hole is significantly increased with an increase in the fracturing hole diameter. As the hole spacing decreases, the width and volume of the plastic zone increase. Tables 7–9 indicate that the three factors that influence the fracturing effect have their respective results, with hole spacing exhibiting the most significant influence. The hole depth also affects the volume range of the plastic zone. Therefore, the optimal hole parameters for swelling and fracturing were preliminarily obtained, a hole diameter of 60 mm, hole spacing of 40 cm, and hole depth of 6 m.
6. Engineering Practice
6.1. Static Fracturing Construction Scheme and Process
On No. 21209 working face, the average daily suspended roof area of the corner roof was approximately 10 × 3.5 m. Therefore, swelling and fracturing construction was performed in two grooves along the working face. The construction process is illustrated in Figure 12. The expansive agent and water were mixed at a mixing ratio of 0.33 and then injected into the fracturing hole with a grouting pump after stirring.
According to the width of the two roadways of the working face, 10 and 12 roof-cutting holes were designed, respectively. Two roadways of 50 m were placed in front of the working face for roof cutting. The expansive agent and water were mixed at a mixing ratio of 0.33. Considering the belt roadway as an example, the site construction scheme is illustrated in Figure 13.
The top of the roadway was designed with roof-cutting holes 0.5 m from both pit slopes. The remaining hole spacing was 0.4 m, hole length was 6 m, and inclination toward the gob was approximately 10°. The grouting liquid was directly injected into the hole using a grouting pump. To ensure the fracturing effect, a mine hole packer was used to seal the hole. The on-site grouting construction process is depicted in Figure 14.
(a)
(b)
6.2. Testing Fracturing Effect
Twenty-four hours after static fracturing, three monitoring methods, i.e., peering into the hole, evaluating the pressure distribution of the advanced support, and examining the bed separation of the roof, were adopted to analyse the effect on the physical and mechanical characteristics of the roof fracturing.
6.2.1. Physical Detection of Fracturing Effect
After static fracturing, detection was conducted within 0.4 m of the roof fracturing range. The fracturing effect is illustrated in Figure 15.
(a)
(b)
The fracturing inside the roof is depicted in Figure 15(a). The roof was blocky 0.5 and 1.5 m toward the inside. The fissure at 2.5 m was obviously expanded, and some blocks were present. The roof at 3.5–5.5 m exhibited bed separation and was partially blocky. In the entire fracturing range, the integrity of the roof was broken, and the fracture evolution became faster. The entire roof bottom was broken, as depicted in Figure 15(b). The roof was bending to the side of the gob and caving, and the rock mass was further broken into several larger rock blocks.
6.2.2. Distribution of Support Pressure
The support pressure conditions under the original state and static fracturing conditions were imposed at the same position of the belt roadway of the working face for two days, and they were continuously monitored for 12 h. The distribution characteristics of the support pressure are depicted in Figure 16.
In the original state, the working resistance of the support oscillated at approximately 24 MPa and remained stable. Six consecutive hours after static fracturing and grouting, the continuous hydration reaction started to occur. At the seventh hour, the support pressure increased to 26.1 MPa. After 8 h, the pressure decreased to 23 MPa. Over time, the support pressure continuously reduced. After 12 h, the support pressure decreased to 17.58 MPa. Thus, the corner roof could collapse in time, without causing a suspended roof.
Hence, during the entire process of chemical swelling and fracturing, there is no need to evacuate personnel, and other work processes are not adversely affected, thereby ensuring the continuity and efficiency of mine operations and improving the operation coefficient of safety.
7. Conclusions
(1)Static fracturing technology can increase the swelling force. The relationship between swelling force and fracturing hole spacing was investigated through theoretical analysis. It was concluded that if holes are equally spaced, then the compressive and tensile stresses are at their maximum.(2)After the expansive agent reacts chemically with water, the maximum swelling force is generated in the same hole if the water-cement ratio is 0.33. The swelling forces of different hole diameters demonstrate an inverted parabolic distribution characteristic. Fitting the relationship between hole diameter and swelling force demonstrated that the swelling force is the largest with a water-cement ratio of 0.33 and a hole diameter of 53.28 mm.(3)Twelve experimental cases were designed for three variants—hole diameter, hole depth, and hole spacing—with the SPSSAU orthogonal design software. The stress distribution and plastic failure characteristics under different fracturing conditions were computed using the corresponding FLAC3D numerical values. The optimal static fracturing parameters obtained include a hole diameter of 60 mm, hole spacing of 40 cm, and hole depth of 6 m.(4)The on-site monitoring results revealed that the hydration reaction occurs after 6 h. The support pressure increases up to 26.1 MPa at the seventh hour and then declines, and the support pressure continuously decreases to 17.58 MPa over time. The internal fissures in the rock stratum within the fracturing range are developed and linked to each other. Thus, the collapse of the corner roof can be achieved in time without causing a suspended roof.Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was sponsored by the National Natural Science Foundation of China (grant nos. 51874231 and 51504184), which is part of the Shaanxi Natural Science Fundamental Research Program Enterprise United Fund (grant no. 2019JLZ-04), and the Outstanding Youth Science Fund of the Xi’an University of Science and Technology (grant no. 2019YQ2-16).