Method for evaluating stability of covered slope of citrus park in hilly land
By acquiring and calculating multiple parameters in the citrus orchards in the hilly areas, the slope safety coefficient is evaluated, and the problem of inaccurate slope stability judgment in the prior art is solved, and a more accurate and reliable slope stability evaluation is achieved.
Patent Information
- Application Number
- CN202411889040.6
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-20
- Publication Date
- 2025-04-15
AI Technical Summary
The stability of the cover slopes of citrus orchards in the existing hilly areas is inaccurate, especially in areas with significant vegetation diversity and topographic changes, it is difficult to comprehensively analyze the comprehensive impact of vegetation on slope stability.
By obtaining the rainfall of the day, canopy water storage, soil and water parameters and vegetation parameters, the effective rainfall, topographic humidity index, additional cohesion provided by the plant root system and the vegetation’s own gravity load, combined with mechanical tests and survey results, the slope safety coefficient is determined, thereby evaluating slope stability.
This method can comprehensively analyze the comprehensive impact of vegetation on slope stability, improve the effectiveness and reliability of the stability assessment of citrus orchards covering slopes in hilly areas, and solve the problem of inaccurate judgments.
Smart Images

Figure CN119831333A_ABST
Abstract
Description
Technical Field
[0001] The invention relates to the technical field of slope stability judgment, and in particular to a method for evaluating the stability of slopes covered by citrus orchards in hilly areas. Background Art
[0002] Existing physics-based landslide susceptibility models effectively combine simulated rainfall data and demonstrate good prediction performance, but it is difficult to comprehensively analyze the comprehensive impact of vegetation on slope stability. In other words, they may face limitations when applied in areas with significant vegetation diversity and terrain changes, especially for hilly areas suitable for citrus cultivation, which leads to inaccurate judgments on the stability of slopes covered by citrus orchards in hilly areas. Summary of the invention
[0003] Technical problem solved by the present invention: The present invention provides a method for evaluating the stability of slopes covered by citrus orchards in hilly areas, which solves the problem of inaccurate stability judgment of slopes covered by citrus orchards in hilly areas.
[0004] The present invention solves the above technical problems by adopting a technical solution: a method for assessing the stability of a slope covered by a citrus orchard in a hilly area, comprising the following steps:
[0005] S1. Obtain the current day rainfall and the canopy water storage of the previous day in the citrus orchard in the hilly area, add the current day rainfall and the canopy water storage of the previous day to determine the current day canopy water storage, and determine the effective rainfall based on the current day canopy water storage and the maximum canopy water storage;
[0006] S2, quantify the saturated hydraulic conductivity of the root-soil complex based on mechanical tests and survey results of different vegetation types in citrus orchards in hilly areas;
[0007] S3, obtaining soil and water parameters and vegetation parameters of citrus orchards in hilly areas, wherein the soil and water parameters include sliding surface depth, sliding surface width, uphill drainage area, width of outflow boundary, slope area, soil unit width, soil shear distortion angle, slope angle, soil cohesion, soil internal friction angle, soil depth, water unit weight, soil unit weight above groundwater level, saturated water content and retained water content; and the vegetation parameters include vegetation diameter at breast height, tree height, forest density, forest density, root cross-sectional area, average tensile strength per unit root cross-sectional area and rooted soil area;
[0008] S4, determining a terrain moisture index according to the effective rainfall, the saturated hydraulic conductivity of the root-soil complex, the depth of the sliding surface, the upslope drainage area, the width of the outflow boundary and the shear distortion angle of the soil;
[0009] S5. Determine the additional cohesion provided by the plant root system according to the root cross-sectional area, the average tensile strength per unit root cross-sectional area, the internal friction angle of the soil body, the rooted soil area and the shear distortion angle of the soil body;
[0010] S6. Determine the vegetation's own gravity load based on the DBH, the tree height, the tree density, the forest density, the slope area, and the soil unit width;
[0011] S7, determining soil weight according to the water unit weight, the soil unit weight above the groundwater level, the saturated water content, the retained water content, the sliding surface depth, the sliding surface width, the slope angle and the terrain moisture index;
[0012] S8. determining the earthquake force according to the soil weight and the soil depth;
[0013] S9, determining the slope safety factor according to the soil cohesion, the additional cohesion provided by the plant roots, the sliding surface width, the vegetation gravity load, the soil weight, the water unit weight, the terrain moisture index, the sliding surface depth, the slope angle, the soil internal friction angle and the seismic force;
[0014] S10. Determine the slope stability of the citrus orchard in the hilly area based on the slope safety factor.
[0015] Furthermore, in S1, the calculation formula for effective rainfall is: P ef =max[0,Scan-Scan max ], where P ef represents effective rainfall, max represents the maximum value function, Scan represents the canopy water storage capacity on that day, Scan max Indicates the maximum water storage capacity of the canopy.
[0016] Furthermore, in S4, the calculation formula of terrain humidity index is: Where T = k s-rs h, m represents the terrain humidity index, P ef represents effective rainfall, T represents soil permeability, k s-rs represents the saturated hydraulic conductivity of the root-soil complex, h represents the depth of the sliding surface, and A up represents the upslope drainage area, B represents the width of the outflow boundary, θ represents the soil shear distortion angle, and min represents the minimum function.
[0017] Furthermore, in S5, the calculation formula for the additional cohesion provided by the plant roots is: C root =k'T rootRAR, where k'=cosθtanφ+sinθ, C root represents the additional cohesion provided by plant roots, k' represents the correction coefficient, θ represents the soil shear distortion angle, φ represents the soil internal friction angle, and A r represents the root cross-sectional area, A rs represents the rooted soil area, RAR represents the root area ratio, T root It represents the average tensile strength per unit root cross-sectional area.
[0018] Furthermore, in S6, the calculation formula for the vegetation's own gravity load is: Among them, W veg represents the vegetation's own gravity load, V represents the forest volume, ρ m represents the density of trees, g represents the acceleration of gravity, ρ l represents forest density, S represents slope area, and w cell represents the width of soil unit, and the calculation formula of forest volume is:
[0019] D represents the diameter at breast height, and L represents the tree height.
[0020] Furthermore, in S7, the soil weight is calculated as: W soil =hbcosβ[(1-m)γ m +mγ sat ], where W soil represents soil weight, h represents sliding surface depth, b represents sliding surface width, β represents slope angle, m represents terrain moisture index, γ m represents the unit weight of soil above the groundwater level, γ sat represents the unit weight of saturated soil, γ sat =γ w (θ s -θ i )+γ m , γ w represents the unit weight of water, θ s represents the saturated water content, θ i Indicates maintaining moisture content.
[0021] Furthermore, in S8, the calculation formula for earthquake force is: T E =k E W soil , where T E represents the earthquake force, W soil represents soil weight, k E Scaling factor representing the ratio of earthquake force to soil weight, PGA stands for peak seismic acceleration, g for gravitational acceleration, and H for soil depth.
[0022] Furthermore, in S9, the calculation formula for the slope stability of the citrus orchard in the hilly area is:
[0023] Among them, F s represents the slope safety factor, C soil represents soil cohesion, C root represents the additional cohesion provided by plant roots, b represents the width of the sliding surface, W veg represents the vegetation's own gravity load, W soil represents the soil weight, γ w represents the unit weight of water, m represents the terrain moisture index, h represents the depth of the sliding surface, T E represents the earthquake force, β represents the slope angle, represents the internal friction angle of soil.
[0024] Furthermore, in S10, the stability of the slope covering the citrus orchard in the hilly area is determined based on the slope safety coefficient, including: if the slope safety coefficient is greater than 1, the stability of the slope covering the citrus orchard in the hilly area is good; if the slope safety coefficient is not greater than 1, the stability of the slope covering the citrus orchard in the hilly area is poor.
[0025] Further, the slope stability of the citrus orchard in the hilly area is determined according to the slope safety factor, including: when F s >1.5, indicating that the landslide risk level is low and prone to occur. s >1.25, indicating that the landslide risk level is relatively low and prone to occur. s >1, indicating that the landslide risk level is moderately prone to occur. s >0.5, indicating that the landslide risk level is high and prone to occur. s When ≤0.5, it means that the landslide risk level is high and prone to occur. s Represents the slope safety factor.
[0026] Beneficial effects of the invention: The invention provides a method for assessing the stability of slopes covered in citrus orchards in hilly areas. Through the soil and water parameters and vegetation parameters of citrus orchards in hilly areas, the vegetation on the surface, as well as the slope moisture regulation in the underground part and the additional cohesion provided by the plant roots are quantitatively assessed. The comprehensive impact of vegetation on slope stability can be comprehensively analyzed, thereby improving the effectiveness and reliability of the assessment of the stability of slopes covered in citrus orchards in hilly areas, and solving the problem of inaccurate stability judgment of existing slopes covered in citrus orchards in hilly areas. BRIEF DESCRIPTION OF THE DRAWINGS
[0027] Figure 1It is a schematic flow chart of a method for assessing the stability of a slope covered by a citrus orchard in a hilly area provided by the present invention. DETAILED DESCRIPTION
[0028] Aiming at the problem of inaccurate judgment on the stability of slopes covered by citrus orchards in hilly areas, the present invention provides a method for evaluating the stability of slopes covered by citrus orchards in hilly areas. Figure 1 As shown, the following steps are included:
[0029] S1. Obtain the current day rainfall and the canopy water storage of the previous day in the citrus orchard in the hilly area, add the current day rainfall and the canopy water storage of the previous day to determine the current day canopy water storage, and determine the effective rainfall based on the current day canopy water storage and the maximum canopy water storage;
[0030] Specifically, the random simulation method of extreme rainfall scenarios, i.e., the Gumbel model and the bounded random cascade model, is used to reproduce rainfall events under the random rainfall mode, thereby obtaining the rainfall of the day in the citrus orchard in the hilly area. The canopy water storage on the previous day is the water retained in the plant canopy on the previous day, which is obtained through soil moisture monitoring, plant evaporation calculation or direct measurement. The maximum canopy water storage can be obtained through information such as plant characteristics, soil type and meteorological conditions. The calculation formula for effective rainfall is: P ef =max[0,Scan-Scan max ], where P ef represents effective rainfall, max represents the maximum value function, Scan represents the canopy water storage capacity on that day, Scan max Indicates the maximum water storage capacity of the canopy. Effective rainfall reflects the interception effect of the vegetation canopy on the surface. The canopy provides a buffer for the slope by intercepting rainfall and protects the stability of the slope.
[0031] S2, quantify the saturated hydraulic conductivity of the root-soil complex based on mechanical tests and survey results of different vegetation types in citrus orchards in hilly areas;
[0032] Specifically, the citrus orchards in hilly areas were taken as the research objects, and different vegetation types were found, such as citrus trees, shrubs and herbaceous plants, and representative soil and root samples were collected. The treated soil and root samples were tested using mechanical testing equipment to obtain relevant mechanical parameters. Combined with the vegetation type survey data and mechanical test results, quantitative analysis was performed to evaluate the water conductivity of the root-soil complex under different vegetation types. Based on the quantitative analysis results, the saturated water conductivity of the root-soil complex was calculated.
[0033] S3, obtaining soil and water parameters and vegetation parameters of citrus orchards in hilly areas, wherein the soil and water parameters include sliding surface depth, sliding surface width, uphill drainage area, width of outflow boundary, slope area, soil unit width, soil shear distortion angle, slope angle, soil cohesion, soil internal friction angle, soil depth, water unit weight, soil unit weight above groundwater level, saturated water content and retained water content; and the vegetation parameters include vegetation diameter at breast height, tree height, forest density, forest density, root cross-sectional area, average tensile strength per unit root cross-sectional area and rooted soil area;
[0034] Specifically, soil-water parameters mainly reflect the soil moisture status, topographic characteristics and physical and mechanical properties of the citrus orchard in the hilly area, while vegetation parameters mainly describe the growth status, distribution density and root characteristics of vegetation in the citrus orchard in the hilly area. Based on the width of the soil unit, the shear test can obtain the soil shear distortion angle and internal friction angle of the soil in the citrus orchard in the hilly area. The soil cohesion is the bearing capacity of the soil in the absence of water, the saturated water content is the water content of the soil in the saturated state, and the retained water content is the amount of water that the soil can retain in the dry state. The diameter at breast height of the vegetation is the diameter of the trunk at breast height (usually 1.3 meters above the ground), the forest density is the number of trees per unit area, and the forest density is the vegetation coverage of the forest open space. Field surveys were conducted in the citrus orchard in the hilly area, and appropriate measurement tools and techniques (such as GPS, rangefinder, soil sampler, tape measure, caliper, etc.) were used to collect soil-water parameter and vegetation parameter data.
[0035] S4, determining a terrain moisture index according to the effective rainfall, the saturated hydraulic conductivity of the root-soil complex, the depth of the sliding surface, the upslope drainage area, the width of the outflow boundary and the shear distortion angle of the soil;
[0036] Specifically, the calculation formula of terrain humidity index is: Where T = k s-rs h, m represents the terrain humidity index, P ef represents effective rainfall, T represents soil permeability, k s-rs represents the saturated hydraulic conductivity of the root-soil complex, h represents the depth of the sliding surface, and A upRepresents the upslope drainage area, B represents the width of the outflow boundary, θ represents the soil shear distortion angle, and min represents the minimum function. The terrain wetness index is a key parameter for quantifying pore water pressure. The pore water pressure is the water pressure in the soil below the groundwater level, acting on the soil in a direction perpendicular to the slope. The pore water pressure affects the effective stress of the soil on the slope, that is, the actual bearing capacity, which can offset part of the contact pressure between soil particles. When the pore water pressure increases (for example, after rainfall), it can lead to a decrease in the effective stress of the soil, thereby reducing the stability of the soil and increasing the risk of landslides. Conversely, lower pore water pressure can increase the contact force between soil particles, which helps to improve the stability of the soil. The terrain wetness index uses Darcy's law and phenomena discovered in field experiments to quantify the proportion of saturated soil layers. is the effective water flow through the sliding surface area, if If it is greater than 1, it means that the water supply exceeds the carrying capacity of the soil, and the terrain moisture index is 1. Less than or equal to 1, indicating that the soil can effectively receive and retain moisture, the terrain moisture index is By calculating the terrain wetness index, it is possible to assess whether the moisture status of the soil under specific rainfall conditions will lead to landslides or other soil instability.
[0037] S5. Determine the additional cohesion provided by the plant root system according to the root cross-sectional area, the average tensile strength per unit root cross-sectional area, the internal friction angle of the soil body, the rooted soil area and the shear distortion angle of the soil body;
[0038] Specifically, the calculation formula for the additional cohesion provided by the plant roots is: C root =k'T root RAR, where k'=cosθtanφ+sinθ, C root represents the additional cohesion provided by plant roots, k' represents the correction coefficient, θ represents the soil shear distortion angle, φ represents the soil internal friction angle, and A r represents the root cross-sectional area, A rs represents the rooted soil area, RAR represents the root area ratio, T root It indicates the average tensile strength per unit root cross-sectional area. The correction coefficient can adjust the shear strength of the soil through the shear distortion angle of the soil and the internal friction angle of the soil, reflecting the bearing capacity of the soil under different shear conditions. The root area ratio indicates the influence of the root system on the soil structure. The additional cohesion provided by the plant root system indicates the strengthening effect of the root system on the soil, that is, the shear strength increased by the root system in the soil. The higher the additional cohesion provided by the plant root system, the higher the shear strength of the soil.
[0039] S6. Determine the vegetation's own gravity load based on the DBH, the tree height, the tree density, the forest density, the slope area, and the soil unit width;
[0040] Specifically, the calculation formula for the vegetation's own gravity load is: Among them, W veg represents the vegetation's own gravity load, V represents the forest volume, ρ m represents the density of trees, g represents the acceleration of gravity, ρ l represents forest density, S represents slope area, and w cell represents the width of the soil unit, D represents the diameter at breast height, and L represents the tree height. When D ≥ 10, the growth characteristics of larger trees apply. When D < 10, it applies to smaller trees, and its volume is estimated by different mathematical formulas. This segmented calculation method takes into account the growth patterns of trees with different diameters. The self-weight of meadows has almost no effect on slope stability, so the self-weight of meadows is not considered. During the field survey, it was found that the self-weight of shrubs with a breast diameter of less than 10 cm had no significant effect on slope stability, so the self-weight of shrubs was not considered.
[0041] S7, determining soil weight according to the water unit weight, the soil unit weight above the groundwater level, the saturated water content, the retained water content, the sliding surface depth, the sliding surface width, the slope angle and the terrain moisture index;
[0042] Specifically, the soil weight is calculated as: W soil =hbcosβ[(1-m)γ m +mγ sat ], where W soil represents soil weight, h represents sliding surface depth, b represents sliding surface width, β represents slope angle, m represents terrain moisture index, γ m represents the unit weight of soil above the groundwater level, γ sat represents the unit weight of saturated soil, γ sat =γ w (θ s -θ i )+γ m , γ w represents the unit weight of water, θ s represents the saturated water content, θ i Indicates the moisture content. When m = 0, it means the soil is completely dry, and when m = 1, it means the soil is completely saturated.
[0043] S8. determining the earthquake force according to the soil weight and the soil depth;
[0044] Specifically, the calculation formula for earthquake force is: TE =k E W soil , where T E represents the earthquake force, W soil represents soil weight, k E Scaling factor representing the ratio of earthquake force to soil weight, PGA stands for peak seismic acceleration, g stands for gravitational acceleration, and H stands for soil depth. The deeper the soil depth, the smaller the scaling factor of the ratio of seismic force to soil weight, indicating that the response of deeper soil layers to seismic forces is attenuated.
[0045] S9, determining the slope safety factor according to the soil cohesion, the additional cohesion provided by the plant roots, the sliding surface width, the vegetation gravity load, the soil weight, the water unit weight, the terrain moisture index, the sliding surface depth, the slope angle, the soil internal friction angle and the seismic force;
[0046] Specifically, the calculation formula for the slope stability of citrus orchards in hilly areas is:
[0047] Among them, F s represents the slope safety factor, C soil represents soil cohesion, C root represents the additional cohesion provided by plant roots, b represents the width of the sliding surface, W veg represents the vegetation's own gravity load, W soil represents the soil weight, γ w represents the unit weight of water, m represents the terrain moisture index, h represents the depth of the sliding surface, T E represents the earthquake force, β represents the slope angle, It indicates the internal friction angle of the soil. The numerator indicates the anti-sliding force on the slope, and the denominator indicates the sliding force on the slope. The slope safety factor is the ratio of the anti-sliding force to the sliding force, indicating the stability of the slope. The larger the ratio, the more stable the slope.
[0048] S10. Determine the slope stability of the citrus orchard in the hilly area based on the slope safety factor.
[0049] Specifically, if the slope safety coefficient is greater than 1, the stability of the slope covered by the citrus orchard in the hilly area is better; if the slope safety coefficient is not greater than 1, the stability of the slope covered by the citrus orchard in the hilly area is poor.
[0050] It can also be subdivided into: when F s >1.5, indicating that the landslide risk level is low and prone to occur. s >1.25, indicating that the landslide risk level is relatively low and prone to occur. s>1, indicating that the landslide risk level is moderately prone to occur. s >0.5, indicating that the landslide risk level is high and prone to occur. s When ≤0.5, it indicates that the landslide risk level is high. Based on the observation data, susceptibility zoning mapping was carried out under extreme rainfall scenarios with different rainfall durations and different recurrence periods, and the corresponding visualization map of the safety and stability of the slopes covered by citrus orchards in hilly areas was obtained.
Claims
1. A method for assessing the stability of slopes covered by citrus orchards in hilly areas, characterized in that: The following steps are involved: S1. Obtain the current day rainfall and the canopy water storage of the previous day in the citrus orchard in the hilly area, add the current day rainfall and the canopy water storage of the previous day to determine the current day canopy water storage, and determine the effective rainfall based on the current day canopy water storage and the maximum canopy water storage; S2, quantify the saturated hydraulic conductivity of the root-soil complex based on mechanical tests and survey results of different vegetation types in citrus orchards in hilly areas; S3, obtaining soil and water parameters and vegetation parameters of citrus orchards in hilly areas, wherein the soil and water parameters include sliding surface depth, sliding surface width, uphill drainage area, width of outflow boundary, slope area, soil unit width, soil shear distortion angle, slope angle, soil cohesion, soil internal friction angle, soil depth, water unit weight, soil unit weight above groundwater level, saturated water content and retained water content; and the vegetation parameters include vegetation diameter at breast height, tree height, forest density, forest density, root cross-sectional area, average tensile strength per unit root cross-sectional area and rooted soil area; S4, determining a terrain moisture index according to the effective rainfall, the saturated hydraulic conductivity of the root-soil complex, the depth of the sliding surface, the upslope drainage area, the width of the outflow boundary and the shear distortion angle of the soil; S5. Determine the additional cohesion provided by the plant root system according to the root cross-sectional area, the average tensile strength per unit root cross-sectional area, the internal friction angle of the soil body, the rooted soil area and the shear distortion angle of the soil body; S6. Determine the vegetation's own gravity load based on the DBH, the tree height, the tree density, the forest density, the slope area, and the soil unit width; S7, determining soil weight according to the water unit weight, the soil unit weight above the groundwater level, the saturated water content, the retained water content, the sliding surface depth, the sliding surface width, the slope angle and the terrain moisture index; S8. determining the earthquake force according to the soil weight and the soil depth; S9, determining the slope safety factor according to the soil cohesion, the additional cohesion provided by the plant roots, the sliding surface width, the vegetation gravity load, the soil weight, the water unit weight, the terrain moisture index, the sliding surface depth, the slope angle, the soil internal friction angle and the seismic force; S10. Determine the slope stability of the citrus orchard in the hilly area based on the slope safety factor.
2. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S1, the calculation formula for effective rainfall is: P ef =max[0,Scan-Scan max ], where P ef represents effective rainfall, max represents the maximum value function, Scan represents the canopy water storage capacity on that day, Scan max Indicates the maximum water storage capacity of the canopy.
3. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S4, the calculation formula of terrain wetness index is: Where T = k s-rs h, m represents the terrain humidity index, P ef represents effective rainfall, T represents soil permeability, k s-rs represents the saturated hydraulic conductivity of the root-soil complex, h represents the depth of the sliding surface, and A up represents the upslope drainage area, B represents the width of the outflow boundary, θ represents the soil shear distortion angle, and min represents the minimum function.
4. The method for assessing slope stability of a citrus orchard in a hilly area according to claim 1, characterized in that: In S5, the calculation formula for the additional cohesion provided by the plant roots is: C root =k'T root RAR, among which, C root represents the additional cohesion provided by plant roots, k' represents the correction coefficient, θ represents the soil shear distortion angle, represents the internal friction angle of soil, A r represents the root cross-sectional area, A rs represents the rooted soil area, RAR represents the root area ratio, T root It represents the average tensile strength per unit root cross-sectional area.
5. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S6, the calculation formula for the vegetation's own gravity load is: Among them, W veg represents the vegetation's own gravity load, V represents the forest volume, ρ m represents the density of trees, g represents the acceleration of gravity, ρ l represents forest density, S represents slope area, and w cell represents the width of the soil unit, D represents the diameter at breast height, and L represents the tree height.
6. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S7, the soil weight is calculated as: W soil =hbcosβ[(1-m)γ m +mγ sat ], where W soil represents soil weight, h represents sliding surface depth, b represents sliding surface width, β represents slope angle, m represents terrain moisture index, γ m represents the unit weight of soil above the groundwater level, γ sat represents the unit weight of saturated soil, γ sat =γ w (θ s -θ i )+γ m , γ w represents the unit weight of water, θ s represents the saturated water content, θ i Indicates maintaining moisture content.
7. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S8, the formula for calculating the earthquake force is: T E =k E W soil , where T E represents the earthquake force, W soil represents soil weight, k E Scaling factor representing the ratio of earthquake force to soil weight, PGA stands for peak seismic acceleration, g for gravitational acceleration, and H for soil depth.
8. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S9, the calculation formula for the slope stability of citrus orchards in hilly areas is: Among them, F s represents the slope safety factor, C soil represents soil cohesion, C root represents the additional cohesion provided by plant roots, b represents the width of the sliding surface, W veg represents the vegetation's own gravity load, W soil represents the soil weight, γ w represents the unit weight of water, m represents the terrain moisture index, h represents the depth of the sliding surface, T E represents the earthquake force, β represents the slope angle, represents the internal friction angle of soil.
9. The method for assessing the stability of slopes covered by citrus orchards in hilly areas according to claim 1, characterized in that: In S10, the stability of the slope covering the citrus orchard in the hilly area is determined based on the slope safety coefficient, including: if the slope safety coefficient is greater than 1, the stability of the slope covering the citrus orchard in the hilly area is good; if the slope safety coefficient is not greater than 1, the stability of the slope covering the citrus orchard in the hilly area is poor.
10. The method for assessing slope stability of citrus orchards in hilly areas according to claim 1, characterized in that: In S10, the slope stability of the citrus orchard in the hilly area is determined according to the slope safety factor, including: when F s >1.5, indicating that the landslide risk level is low and prone to occur. s >1.25, indicating that the landslide risk level is relatively low and prone to occur. s >1, indicating that the landslide risk level is moderately prone to occur. s >0.5, indicating that the landslide risk level is high and prone to occur. s When ≤0.5, it means that the landslide risk level is high and prone to occur. s Represents the slope safety factor.