Abstract

The relationship between slope and prevention structures under the action of rainfall is complicated. In this paper, a series of model tests were carried out to study the coevolution mechanism of slope and prevention structures under rainfall. It was found that the prestressed grouting steel anchor tube frame beam composite structure has a good reinforcement and protective effect on the slope. Specifically, less damage degree is observed for a slope with prestressed grouting steel anchor tube frame beam, while the slope of without prestressed grouting steel anchor tube frame beam suffers from severe deformation. Under the effect of rainfall, the evolution process of slope and prevention structure cooperative deformation can be divided into three stages: creep stage, extrusion stage, and acceleration stage. The antisliding force provided by the prevention structure is positively correlated with the slope deformation. Moreover, as an initiating factor for changes in various monitoring parameters, rainfall plays a dominant role in the variation of earth pressure reinforcement by steel anchor tube.

1. Introduction

Rainfall is one of the most important factors that induce landslides, and there are countless landslides caused by rainfall all over the world [1, 2]. About 90% of landslide disasters in China alone are directly or indirectly related to rainfall by records [36]. The occurrence of landslides seriously endangers the safety of human life and property and causes the loss of national economic property [7]. Therefore, the protection and treatment of landslides have always been a concern of people. Advanced and effective treatment measures can effectively treat landslides and avoid disasters caused by landslides, which is the goal of continuous research and exploration. A large number of scholars have studied the reinforcement mechanism of different reinforcement methods and the coevolution mechanism of different prevention structures and slopes by using numerical simulation and model tests [811]. However, the synergistic relationship between the slope and the prevention structure is very complicated, and the infiltration of rainfall into the soil further aggravates this complicated relationship. Understanding the synergistic relationship between the slope and the prevention structure is conducive to further improving and perfecting the reinforcement technology. Therefore, an experimental study on the mechanism of synergistic evolution of slopes and protective structures under the action of rainfall was carried out. This research adopts a new type of prevention structure, namely, prestressed grouting steel anchor tube frame beam, and studied the coevolution mechanism between it and the slope by model tests. Additionally, setting up a control test to verify its reinforcement performance.

The prestressed grouting steel anchor tube is an improved structure based on prestressed anchor cable, and the prestressed anchor cables are usually used with antislide piles. The formation of antislide piles with prestressed anchor cables mainly starts with antislide piles, besides, as an excellent measure for slope reinforcement, antislide piles are widely used in slope reinforcement projects [46]. The synergy between antislide piles and slopes has been studied at home and abroad [10, 1216]. For example, Chen et al. [17] studied the failure mechanism of antisliding piles on a road cutting slope and believed that the direct removal of the slope in front of the pile without using anchor cables was the direct cause of the damage of antisliding piles. Moreover, continuous rainfall over several days aggravated the landslide damage because of the increase of the self-weight and degradation of the mechanical parameters of the slope materials. As a result, the damage of the antislide piles is aggravated. However, the local reinforcement of the slope and the protection of the slope surface by antislide piles are big defects. Therefore, researchers continue to improve on the basis of the antislide pile reinforcement method and gradually form a prestressed anchor cable antislide pile frame beam slope reinforcement structure, which gets a lot of attention and researches [1820]. Wu et al. [21] and Fan et al. [22] studied the dynamic response of slope reinforced by prestressed anchor cable and antislide pile by indoor shaking table test. Fan et al. [23], Branco et al. [24], and Bathurst et al. [25] studied the control effect of anchor cable antislide pile reinforcement on slope sliding, showing that the protective structure has good reinforcement and protection performance on slopes and slope surfaces. In addition, some scholars have studied the reinforcement characteristics of anchor cables in the roadway, surrounding rock, carrying out optimization research on the reinforcement parameters in these projects [46, 26].

In recent years, with the development of slope protection technology towards lighter weight and durability, the previous protection technology cannot meet the durability requirements. Therefore, the grouting steel anchor tube reinforcement technology emerged at the times required, becoming a new research boom in slope reinforcement technology in recent years. Wang et al. [27, 28] have studied the horizontal bearing capacity of vertical steel anchor tube micropile through a large-scale model groove test. Chen et al. [29] carried out a full-scale model test on the antisliding characteristics of the high-pressure grouted steel anchor tube micropile and discussed the antisliding characteristics of the high-pressure grouted steel anchor tube micropile under real stress conditions. Zhang et al. [30] used on-site large-scale model tests to study the new protection technology of grouting steel anchor tube pile group structure and its mechanism. Zhou et al. [31] used on-site prototype tests to study the antisliding performance of the grouted steel anchor tube single pile from the aspects of landslide thrust, pile side earth pressure, pile body bending moment, etc., analyzing the failure of the grouted steel anchor tube.

In conclusion, the grouting steel anchor tube reinforcement technology has gradually become a hot issue, while most of the current research is mainly limited to the active reinforcement performance of the grouted steel anchor tube structure itself. The research on the synergy between slope and protective structure is especially scarce. This article relies on the entity engineering of the diseased slope treatment of DaYangyun Expressway, combining the prestressed grouting steel anchor tube with the frame beam structure and carrying out the model test study of the coevolution of the slope and the prestressed steel anchor tube frame beam under the action of rainfall. The experiment results compared the deformation and failure of the side slope without the prestressed steel anchor tube frame beam to verify the reinforcement performance of the prestressed steel anchor tube frame beam. It is expected to provide a reference for follow-up theoretical research and engineering application and the research on the coevolution mechanism between the new technology of the prestressed grouting steel anchor tube and the slope.

2. Analysis of Engineering Geological Characteristics

The study is based on the diseased slope of section of DaYangyun Expressway, located in XiShan Township, ErYuan County, Yunnan Province. It belongs to the middle-low mountain geomorphological region of tectonic erosion, altitude between 2240~2290 m with sparse vegetation on the slope. Overall, the slope is relatively flat, between 15°~30°, the slope aspect is SW5°, and the route direction is NW84°. There is no fault structure in this area, the schistosity development of rock mass, and the dominant occurrence of strata is 123°∠ 32°. The slope, a red-bedded geological body, overlies gravel soil of 2-10 m and underlies highly weathered argillaceous siltstone and highly weathered silty mudstone, and the joint fissure has developed. Therefore, the engineering geological characteristics of rock mass itself are relatively weekly. Due to the site investigation time as rainy season, the foot of the slope accumulated water, and the surface drainage was slow, causing part of the surface water to seep through the surface overburden to the foundation interface, which softens the rock mass. The factors provide favorable conditions for the development of slope damage.

An overview of the diseased slope area is shown in Figure 1. The proposed road passes through the foot of the slope in the form of excavation, resulting in the free face being formed at the front edge of the slope. Because of the strong construction disturbance, the original stress balance of the soil is destroyed. Under the interaction of groundwater, rainfall, and other factors, it is easy to form shear cracks at the excavation site of slope construction. According to the field investigation, there are some collapses at the front edge of the slope shown as cracks at the rear edge and dislocation steps and appeared obvious shear cracks on the firstly-level slope. Combined with the phenomenon of drilling core exposure: during core drilling, which was very fast at 18.5 ~ 20 m, the core is disordered and accompanied by signs of rubbing, which presumably a sliding belt, it is the representation of slope deformation sliding.

Based on the above analysis, the slope is currently in an unstable state, it is speculated that further engineering disturbance, heavy rainfall, or continuous rainfall will induce landslides. The length of the diseased slope is about 97.5 m, the width is about 160 m, the sliding depth is about 20 m, and the volume of the landslide is , which belongs to a large middle-level landslide, the main sliding direction is 185°.

3. Reinforcement Mechanism of Prestressed Grouting Steel Anchor Tube Frame Beam Composite Structure

Prestressed grouting steel anchor frame beam composite structure is a new reinforcement technology that combined “splitting grouting” and “prestressed anchorage.” The steel anchor tube is not only used as a protective measure but also as a channel for grouting, and the rock and soil around the drilling are compacted by splitting grouting, forming a grouting inclusion of “tree root” [30, 32]. Increasing the bonding strength between the grouting body and the rock mass and improving the anchoring force between the grouting body and the rock mass, the steel anchor tube was integrated with the surrounding soil, improving the strength of the surrounding soil. By implanting steel strands inside the steel anchor tube and applying to prestress to the steel strands, the bond strength between splitting grouting and rock mass was fully exerted, and the anchoring force is further improved. At the same time, it can also increase the shear strength of the anchor cable and improve the strength parameters of the sliding zone rock mass. In addition, the frame beams are casted by reinforced concrete and fixed on the slope surface, which mainly plays the role of surface protection and transmission of landslide thrust to deep anchor solids. Compared with the traditional slope anchoring method, this new structural reinforcement technology greatly improves the safety and durability of anchorage. The structure of prestressed the grouting steel anchor tube and its reinforcement mechanism is shown in Figure 2.

4. Design of Model Test

4.1. Design of Similitude Relationship

The rational use of similar materials plays a decisive role in the success of model tests [33]. The parameters of similar materials in the landslide physical model mainly include geometric similarity ratio, density, cohesion, internal friction angle, bulk density, and permeability coefficient [27, 28]. Due to the limitations of existing test conditions, this study mainly considers the geometrical similarity parameters, cohesion, internal friction angle, density, and bulk density of materials. Since the longitudinal height of the selected section is 65 m, combined with the size of the existing model box (), the similarity coefficient of geometric conditions was which is the geometric similarity coefficient, is the model size, and denotes the prototype size of the side slope.

According to the similarity theory, and referring to previous studies of scholars [27, 28], the density similarity ratio and bulk density are both 1, and the geometric similarity ratio () is 50 in this landslide physical model. Based on the dimensional analysis method, the similarity ratios of the angle of internal friction and cohesion (c) are determined as ,

Based on the above derivation results of similar materials, combined with [34] the orthogonal design of the rock similar material ratio research, determination of the similar material proportion of model sliding bed, sliding zone, and sliding body by orthogonal matching ratio and indoor direct shear test, coarse sand: is used to simulate sliding bed; quartz sand: is used to simulate the sliding belt; quartz sand: is used to simulate sliding body. The specific parameters are shown in Table 1.

It is necessary to find appropriate structurally similar materials before carrying out model tests. Because of the difference in mechanical properties of different simulated objects, the selection of similar materials should comply with different similarity requirements, according to the main mechanical properties of the supporting structure in the reinforcement process. In this model test, the tensile strength of the steel anchor tube and anchor cable is used as the control parameter of similar material, and the bending strength of the frame beam and the antislide pile is used as the control parameter of similar material.

Because it is difficult to realize the absolute similarity selection of the mechanical parameters, the selection of similar materials for protection structures in this model test is mainly based on the geometric similarity constant of the material. To consider the mechanical parameters of materials, the protection structure similar to material selection was determined. The mechanical properties and geometrically similar parameters of the protection structure are shown in Table 2.

4.2. Design of Rainfall

According to China’s National Emergency Broadcasting and the Forecast information of Yunnan Meteorological Observatory, the daily rainfall can be as high as 80~ 120 mm in some areas. Refer to He et al. [35] for the conversion method of rainfall similarity relationship, model rainfall intensity : prototype rainfall intensity : , get the model rainfall intensity . According to the model rainfall time similarity relation : , the model rainfall time was obtained. Considering the weak response of slope to rainfall under strong support and combining with the rainfall intensity design value of [36] in the model test research (the rainfall intensity design value is 30 mm/h). This test enlarges rainfall intensity and shortens rainfall time appropriately. Therefore, the average rainfall intensity of the test design was 32 mm/h, the accumulated rainfall time was about 3 h, and the accumulated rainfall was 96 mm. We adopt the method of controlling the total rainfall constant and intermittent rainfall to better restore the real situation of rainfall.

4.3. The Layout of Measuring Instruments

The sliding displacement of the slope was recorded by using dial gauge displacement meters. On the side of the composite structure with prestressed grouting steel anchor tube frame beam, along the axis of the T4 antislide pile, deploy dial gauge displacement meters on the top of the T4 antislide pile, the fourth-level slope, and the fifth-level slope, respectively. On the side of without the prestressed grouting steel anchor frame beam, along the axis of the T7 antislide pile, deploy dial gauge displacement meters on the top of the T7 antislide pile, the fourth-level slope and the fifth-level slope, respectively. Four earth pressure boxes were installed at different depths behind T3, T5, T8, and T10 antislide piles to monitor the change of earth pressure behind the piles. And the earth pressure boxes were arranged under the steel anchor tubes of the third-level and fourth-level slopes to monitor the change of earth pressure under the steel anchor tubes. The layout of the sensor is shown in Figure 3.

4.4. Model Building

To better restore the prototype of the slope, the process of filling the model from bottom to top is followed, and the slope is divided into two sides along the central line and filled synchronously. One side is equipped with the prestressed grouting steel anchor tube frame beam. And the other side is not equipped with prestressed grouting steel anchor tube frame beam. In addition, the sample material was filled in layers and filled into the model box in the form of 10 cm thick soil layers, and each layer was compacted in the same way to ensure that the model was uniform. The specific process of model making is as follows: preparation of similar materials ⟶ fill the slide bed ⟶ install the prestressed anchor cable antislide pile and the earth pressure after the pile ⟶ fill the slide belt ⟶ fill the slide body ⟶ install the steel anchor tube and the earth pressure at the third-level and fourth-level slope ⟶ fill the slide body of the fifth-level and sixth-level slope ⟶ install the frame beam ⟶ grouting for steel anchor pipe ⟶ maintenance of grouting ⟶ prestressing for steel anchor pipe and anchor cable ⟶ install the dial indicator⟶ install the rainfall device ⟶ connect the sensor to the data collector ⟶ debugging test equipment. The process of model filling was shown in Figure 4. And the complete model was shown in Figure 5.

5. Results and Discussion

5.1. Analysis of the Experimental Phenomenon

The appearance of the progressive destruction process phenomenon of the slope under the action of rainfall was shown in Figure 6. As can be seen from Figure 6, the deformation and destruction of the slope gradually intensified with the continuous progress of rainfall, at 14 : 41, the accumulative rainfall was 24.5 mm, and microcracks first appeared on the slope of the fifth level. The rain lasted for five minutes until 14:26, and the accumulative rainfall reached 29.5 mm. A small collapse and grooves occurred on the side of the prestressed grouting steel anchor tube frame beam on the fourth-level slope. At 15 : 15, the accumulative rainfall has increased to 49.5 mm and then microtension cracks appeared on the platform of the six-level slope. At 15 : 49, the accumulated rainfall amounted to 74 mm on the left side of the frame beam of the third-level slope, local slump, and grooves developed. At the same time, the side grooves of the prestressed grouting steel anchor tube on the third-level slope were extremely developed, and the dissolution gullies 1 and 2 are gradually formed. On the fourth-level slope, the soil loss on the side without the prestressed grouting steel anchor tube gradually developed to form a gully 3. At 17 : 30, a large amount of soil collapsed on the side of the third-level slope without prestressed grouting steel anchor tubes, transporting and depositing on the slope foot under the action of rainfall. And the phenomenon of ground settlements appeared on the slope surface. Obviously, a significant difference in the damage degree of the slope on the side where the prestressed grouting steel anchor tube frame beam was installed and the side where was not installed. Mainly as follows, on the side of without the prestressed grouting steel anchor tube frame beam: rain washes and erodes slopes ⟶ crack appeared on the slope surface ⟶ the slope has developed dissolution grooves and ground settlements ⟶ grooves develop to form gullies ⟶ partial gullies development caused slumping ⟶ the pile is deformed by the extrusion of the soil; on the side of the prestressed grouting steel anchor tube frame beam: rain washes and erodes slopes ⟶ a few microcracks appear on the slope surface ⟶ the slope has developed dissolution grooves ⟶ local collapse and ground settlements occurred. Obviously, the deformation and failure degrees of the side slope were small, which indicates that the prestressed grouted steel anchor tube frame beam had a good effect on slope reinforcement and protection.

5.2. Analysis of Earth Pressure behind Prestressed Anchor Antislide Pile

To understand the response of prestressed anchor cable antislip piles under prestressed grouted steel anchor tube frame beams to slope compression during rainfall. We draw the curves of the peak earth pressure after T3, T5, T8, and T10 piles with the change of buried depth as shown in Figure 7. From Figure 7, the earth pressure behind the pile was negatively correlated with the buried depth and generally showed a trend of decreasing with the increase of the buried depth. And the peak earth pressure behind the pile near the surface was the largest, showing that under the action of rainfall, the soil bulk density and the slope weight and the diagonal component increase. And the near-surface section behind the prestressed anchor antislide pile responded most obviously to the sliding force of the slope, which was the main bearing part of the slope sliding force. It is suggested that the stiffness of the near-surface pile should be appropriately increased in practical engineering activities under similar conditions to improve its bearing capacity. In addition, from Figure 7, the peak earth pressure behind the pile near the sliding zone is significantly less than that near the surface, regardless of whether there is prestressed the grouting steel anchor tube frame beam on the slope. It is speculated that the reason for the above phenomena is as follows: in the process of rainfall, the soil bulk density near the surface increases first due to the influence of rainfall infiltration before the overall sliding of the slope occurs. As a result, the earth pressure box near the surface is first squeezed by the soil, so the response of the earth pressure near the surface behind the pile is the most positive.

5.3. Analysis of Earth Pressure under Prestressed Grouting Steel Anchor Tube

The change of the earth pressure curve under the steel anchor tube under the action of rainfall is shown in Figure 8. The earth pressure phases are clearly divided, and the earth pressure variation curve with time under the action of rainfall is mainly divided into four phases in Figure 8. Stage 1: the earth pressure is positively correlated with the accumulated rainfall, which is manifested as a continuous increase with the increase of the accumulated rainfall. The reason for this phenomenon is rainfall infiltration, causing the bulk density of the soil to increase, thereby increasing the oblique force component of the slope, the earth pressure box is increased by the extrusion effect, and the earth pressure curve shows a gradually increasing trend. Stage 2: the accumulated rainfall remains unchanged, meaning the rainfall is 0, and the earth pressure shows a trend of gradual decrease. Combined with the actual test phenomena observed, it can be seen that some rainfall did not penetrate into the soil during the test, but scoured the slope surface to form a gully and flowed away along the gully during the test. In addition, after the rainfall stopped, the rainwater on the slope surface and the local area within the slope continued to accumulate and run off in the gully. It is speculated that the reason for this phenomenon of stage 2 is that the rain stopped, which caused the soil bulk density to decrease, which in turn led to the decrease of the slope component force and the soil rebound. The squeezing effect decreases on the earth pressure box, so the earth pressure curve shows a decreasing trend. Stage 3: the rainfall continues and the earth pressure increases again, repeating the changing trend of stage 1. Stage 4: the accumulated rainfall remained unchanged, and the earth pressure showed a decreasing trend, which repeated the change characteristics of stage 2. Based on the above analysis and combined with the earth pressure curve can be seen, there is a positive correlation between earth pressure and accumulated rainfall. It shows that the magnitude of earth pressure is mainly controlled by rainfall, which is the initiating factor of earth pressure changes.

5.4. Analysis of Displacement Data

The synergistic deformation evolution process of slope and protective structure under rainfall was shown in Figure 9. The displacement monitoring point was shown in Figure 9(a). From Figure 9(b), the slope deformation process under rainfall can be divided into three stages, including the creeping stage: the slope deformation and displacement under the action of rainfall increased slowly. Extrusion stage: the accumulated rainfall continued to increase, the deformation displacement of the slope slowly increases, and the antislide body in the front of the landslide has been compressed. Fretting to accelerated deformation stage: when the accumulated rainfall reaches a certain value, the deformation degree intensifies and then enters the fretting to accelerated deformation stage. One of the special phenomena was that after a short period time from fretting to accelerated deformation, the phenomenon of measuring point 6 from fretting to accelerated deformation disappears. The reason for this phenomenon is that the measuring point is within the reinforcement range of prestressed grouting steel anchor tube frame beam, and the reinforcement of supporting structure plays a role. It is predicted that if there is no prestressed grouting steel anchor tube frame beam, measuring point 6 will show the development trend of the AB section. Showing that prestressed grouting steel anchor tube frame beam has a good inhibitory effect on slope deformation displacement. In addition, at the measuring point ⑨ (the top of the T4 pile), the pile tip was squeezed by soil due to the increase of soil weight from rainfall infiltration. And there was a little displacement change at the beginning; then, the soil in front of the pile was compacted to produce resistance. The reinforcement effect of the pile was exerted and the slope deformation was controlled. It is found that under the action of continuous rainfall compared ② and ④ with ⑥ and ⑧, respectively; there is a significant difference in the deformation and displacement of the side slope of the prestressed grouting steel anchor tube frame beam side and the nonprestressed grouting steel anchor tube frame beam side (Figures 8(c) and 8(d)). The maximum displacement difference of ⑥ and ⑧ is 0.693 mm, and the maximum displacement difference of ② and ④ is 0.826 mm. It is further proved that the superior reinforcement effect of prestressed grouting steel anchor tube frame beam composite structure is remarkable for controlling the deformation displacement of the slope.

5.5. Correlation Analysis of the Coevolution of Slope and Protection Structure Based on Multivariate Data Fusion

To further study the relationship between the slope and the prevention structure, we draw the typical deformation displacement curve of the slope, the earth pressure curve behind the pile, the earth pressure curve under the prestressed grouting steel anchor tube, and the cumulative rainfall curve (Figure 10). From Figure 10, the antisliding force provided by the prevention structure is different in the three-stage process of slope deformation under the action of rainfall. When the slope is in the creeping stage, the earth pressure on the side of the pile and under the steel anchor tube is relatively small, showing the antisliding force provided by the prevention structure is relatively small. The slope deformation entered the extrusion stage, the earth pressure behind the pile, and the earth pressure under the prestressed grouting steel anchor tube increased slowly, and there were obvious signs of increase relative to the creeping stage, indicating that under the action of rainfall the antisliding force increasing provided by the prevention structure is obviously when the slope deformation enters the compression stage. As the slope deformation entered the stage of fretting to accelerated deformation, the earth pressure behind the pile increases rapidly, from 0.29 kPa to 1.32 kPa, which was nearly 4.5 times that of the compression stage. Similarly, the increase of earth pressure under the steel anchor tube was relatively obvious, from 1.59 kPa to 2.57 kPa, meaning that under the action of rainfall after the slope deformation enters the stage of fretting to accelerated deformation, the prevention structure needed to provide greater antisliding force.

6. Conclusions and Recommendations

To study the coevolution mechanism between the slope and the prevention structure under the action of rainfall and verify the reinforcement effect of the prestressed grouting steel anchor tube frame beam. First, analyzing the engineering geological characteristics of the prototype diseased slope and concluding that the main inducing factor for the further development of the diseased slope to form a landslide was rainfall. Based on this, building the coevolution model of the slope and the prevention structure under the action of rainfall through indoor model tests, simulating the coevolution process of slope and control structure under the action of rainfall. The conclusion is as follows: (1)The prestressed grouting steel anchor tube frame beam had a good reinforcement and protective effect on the slope. In this case, the slope deformation displacement and damage degree are significantly different with or without prestressed grouting steel anchor tube frame beam. The deformation and failure degrees of the side slope of the frame beam structure with prestressed grouting steel anchor tube were weak, while the deformation and failure of the side slope without the structure were severe(2)The evolution process of the deformation displacement of the slope under the action of rainfall is divided into three stages: creeping stage ⟶ extrusion stage ⟶ fretting to accelerated deformation stage. When the slope deformation is in the creeping stage, the earth pressure behind the pile is relatively small, and the antisliding force provided by the prevention structure is relatively small. When the slope deformation enters the extrusion stage, the earth pressure behind the pile is increase slowly, and the antisliding force provided by the prevention structure also increases slowly. When the slope deformation enters the stage of fretting to accelerated deformation, the earth pressure behind the pile increases rapidly, which is about 4.5 times that of the extrusion stage, and it indicates that the prevention structure needs to provide greater antisliding force(3)The relationship between earth pressure and rainfall under prestressed grouting steel anchor tube in slope was revealed; rainfall was the initiating factor for changes in earth pressure, and the variation of earth pressure under prestressed steel anchor tube frame beam was mainly controlled by rainfall and was positively correlated with accumulated rainfall. In addition, under the action of rainfall, the earth pressure behind the prestressed anchor cable antisliding pile was negatively correlated with the buried depth, which decreased with the increase of the buried depth, and the peak value of the earth pressure behind the pile near the surface was the largest(4)Revealing the correlation between the deformation displacement of the slope and the protective structure under the action of rainfall. The antisliding force provided by the prevention structure was positively correlated with the slope deformation. The antisliding force provided by the prevention structure was increasing in varying degrees with the increasing deformation of the slope

Data Availability

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

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This research was supported by the National Key R&D Program of China (2018YFC1504904), the National Natural Science Foundation of China (42002280), the Science and Technology Program of Guizhou Province ([2019]1056), science and technology program of Gansu Province (no. 21JR7RA739), and Natural Science Foundation of Gansu Province (No.145RJZA068); and the authors sincerely appreciate the Broadvision Engineering Consultants for the geological survey information, Dr. Yuan Kun to provide modification advice, and Mr. Pang Jun to provide test help on this work.